key: cord-0971425-n9q99hdn authors: Laconi, Andrea; Berends, Alinda J.; de Laat, Esther C.H.; Urselmann, Tara A.P.M.P.; Verheije, Hélène M. title: Infectious bronchitis virus Mass-type (GI-1) and QX-like (GI-19) genotyping and vaccine differentiation using SYBR green RT-qPCR paired with melting curve analysis date: 2019-11-08 journal: J Virol Methods DOI: 10.1016/j.jviromet.2019.113771 sha: 44c73e155f1c27e1d6fcbff9751faf8eb76a9755 doc_id: 971425 cord_uid: n9q99hdn Infectious Bronchitis Virus (IBV) is a highly contagious virus of chicken, causing huge economic losses in the poultry industry. Many genotypes circulate in a given area, and optimal protection relies on vaccination with live attenuated vaccines of the same genotype. As these live vaccines are derived from field viruses and circulate, understanding the prevalence of different IBV genotypes in any area is complex. In a recent study, the genome comparison of an IBV QX vaccine and its progenitor field strain led to the identification of vaccine markers. Here we developed a simplex SYBRgreen RT-qPCR assay for differentiation between QX-like field and vaccine strains and a multiplex SYBRgreen RT-qPCR assay for IBV genotyping with melting curve analysis, as each virus produced distinct and reliable melting peaks. Both the simplex and the multiplex assays showed excellent efficiency, sensitivity and specificity representing a low cost diagnostic tool for IBV genotyping and vaccine differentiation. Infectious Bronchitis Virus (IBV) is a highly contagious avian coronavirus belonging to the family Coronaviridae, genus Gammacoronavirus. IBV infects primarily domestic fowl causing disease associated with respiratory signs, and in some cases with kidney lesions and proventriculitis, resulting in increased mortality in unprotected birds and in serious economic losses (Cook et al., 2012) . There is a wide variety of genetically distinct IBV genotypes; some have been continuously detectable since suitable typing techniques were developed, some have appeared, briefly caused a problem and then disappeared, some proved restrained to a particular geographical area, some widespread (Cook et al., 2012; Jackwood et al., 2012) . The control of the disease occurs mainly by vaccination using live attenuated homologous vaccines (Geerligs et al., 2011; Jackwood et al., 2003) , or a combination of two antigenically distinct live attenuated vaccines Jordan, 2017; Terregino et al., 2008) . Effective knowledge of the prevalence and circulation of different IBV genotypes is needed to decide upon the most effective vaccination protocol. However, this decision is complicated by the circulation and thus presence of the very live vaccines applied for disease control. For example, withdrawal of IBV 793/B (lineage GI-13) live vaccine in Northern Italy led to the disappearance of this previously thought ubiquitous genotype, suggesting that IBV 793/B virulent strains were absent in the area (Franzo et al., 2014a) . IBV detection is commonly performed using TaqMan® RT-qPCR assays (Callison et al., 2006) . IBVs' differentiation and genotyping rely on serological analyses and/or on RT-PCR followed by sequencing of full or partial S1 gene, which, due to its variability and biological function, is the region commonly targeted for this purpose (Cavanagh et al., 1999; Worthington et al., 2008) . However, they both represent expensive and/or time-consuming methodologies. SYBRgreen RT-qPCR paired with melting curve analysis represents a more economic and time-effective strategy for IBV genotyping and for the differentiation between vaccine and field strains in a single analysis. Using this approach genotyping and/or vaccine differentiation relies on sequence variants identification based on their individual melting temperatures, a strategy successfully applied for typing of other avian viruses (Franzo et al., 2014b; Pham et al., 2005; Steer et al., 2011) . However, due to the high number of IBV genotypes (or lineages) (Valastro et al., 2016) and the lack of vaccine markers for most of IBV vaccines (Listorti et al., 2017) the development of effective SYBRgreen RT-qPCR assays represents a challenge for any research group embarking in this task. Here analysis to a) distinguish between IBV field strains and homologous attenuated live vaccine and b) support multiplexing towards the differentiation between IBV genotypes. Comparison of the full genome sequences of vaccine Poulvac IB QX (Zoetis) and its progenitor 1148-A identified several vaccine markers, including three consecutive nucleotide changes in the S1 gene (Listorti et al., 2017) . Those unique vaccine markers affecting the GC content of the target region were targeted to develop a SYBRgreen RT-qPCR assay enabling the differentiation of vaccine Poulvac IB QX from QX-like field strains via melting curve analysis. Next we developed a multiplex assay for the simultaneous genotyping of IBV Mass-type (lineage GI-1) and IBV QXlike (lineage GI-19) strains, as those two genotypes represent a major threat for the poultry industry worldwide (Cook et al., 2012; de Wit et al., 2018; Jackwood et al., 2012) . To achieve this, specific primer sets were designed and RT-qPCR and SYBRgreen analysis settings were optimized in order to detect and differentiate the sequence variants of interest. Both the simplex and the multiplex approaches proved reliable and effective diagnostic tools, showing excellent sensitivity, specificity and repeatability within and between runs. The application of these protocols will aid in the understanding of the molecular epidemiology of IBV and they represent proves of principle for future developments and applications of melting curve analyses towards IBV diagnosis. IBV-D388 was isolated by GD Animal Health (Deventer, The Netherlands) in March 2004 from 19-day-old broiler breeders with respiratory signs and increased mortality due to renal failure and typed as IBV QX-like (Lineage GI-19) (de Wit et al., 2011) . Poulvac IB QX (Zoetis) is a live attenuated avian infectious bronchitis virus derived from QX-like field strain IB 1148-A after serial passages in embryonated chicken eggs (ECEs) (Geerligs et al., 2011) . IBV M41 strain was kindly provide by GD Animal Health (Deventer, The Netherlands); the virus was isolated from birds showing respiratory sings and drop in egg production. Nobilis® IB H120 (MSD/Animal Health, Boxmeer, The Netherlands) is a live attenuated avian infectious bronchitis Mass-type vaccine. For each virus a stock was produced. Briefly, ten 8-day-old specific pathogens free (SPF) ECEs were inoculated with a 50 % embryonic infectious dose (EID 50 ) of 100, incubated for 48 h, and cooled at 4°C for 24 h prior to harvesting and pooling of the allantoic fluid (AF). Viruses titration was performed in 8-day-old ECEs by determining the EID 50 per ml at 7 days post infection (dpi) according to the Reed and Muench method (Reed and Muench, 1938) . RNA was extracted from the allantoic fluid using QIAamp viral RNA minikit (QIAGEN, Hilden, Germany) following the manufacturer's protocol. All available IBV QX-like S1 sequences (96) deposited on ViPR (https://www.viprbrc.org) were downloaded on 30 November 2018 and aligned using MAFFT version 7 (https://mafft.cbrc.jp/alignment/ software/). QX-specific oligo pairs were designed to amplify the portion of S1 gene containing the two transversion (G→T) affecting the GC content of the gene, hence affecting also the melting temperature of the amplicon. To confirm the higher GC content in QX-like field strains detected worldwide, the multiple sequence alignment was trimmed to include solely the target region. For easiness of analysis, the dataset was processed with cd-hit-est test of the CD-HIT Suite (http://weizhongcluster.ucsd.edu/cdhit_suite/cgi-bin/index.cgi?cmd = cd-hit-est) (Li and Godzik, 2006) to cluster sequences that shared 100 % identity, such as each cluster represents a unique prototype sequence. The resulting unique prototype sequences were processed using uMELT Melting Curve Predictions Software (https://www.dna.utah.edu/umelt/umelt. html) to compare the predicted melting curve of vaccine Poulvac IB QX to those of the field strains. Eventually, the specificity of the oligo pairs was assessed using Primer-BLAST (https://www.ncbi.nlm.nih.gov/ tools/primer-blast/) (Jian et al., 2012) . Similarly, all available Masstype S1 sequences (166) deposited on ViPR (https://www.viprbrc.org) were downloaded on 30 November 2018. Using the same approach Mass-type specific oligo sets were designed, with the exception that after clustering a web-based Shannon Entropy calculator (https://www. hiv.lanl.gov/content/sequence/ENTROPY) was used to identify conserved region among the dataset. All RT-qPCRs were carried out in a Bio-Rad CFX Connect real-time PCR system using the iTaq universal SYBRgreen one-step kit (Bio-Rad Laboratories, Hercules, California, USA). Optimization of the assays was performed using ten-fold viral dilutions of Poulvac IB QX, IBV-D388, IBV M41 GD and Nobilis® IB H120, different primer concentrations, annealing temperatures, extension step lengths and temperature increment during the melting curve analysis. The efficiency of the simplex RT-qPCR was determined using serial dilution of Poulvac IB QX ranging from 10 7,83 EID 50 to 10 0,83 EID 50 and of IBV-D388 ranging from 10 7,31 EID 50 to 10 0,31 EID 50 . Similarly the efficiency of the multiplex assay was determined using the same dilution panels for the QX-like viruses, plus dilution panels for IBV M41 GD (10 7,63 EID 50 -10 0,63 EID 50 ) and Nobilis® IB H120 (10 7,75 EID 50 -10 0,75 EID 50 ). Two operators tested the lowest dilution detected as positive by the assays, using nine independent replicates; the lowest viral amount that can be detected in at least 50 % of replicates was considered the assays' limit of detection (LoD). Two operators using three viral dilutions for each virus, corresponding to high, medium and low dilution (one log higher than LoD) during three different experiments performed at weekly intervals, evaluated the repeatability. Inter and intra-run coefficient of variations (CV) were calculated to assess the robustness of the assay. The significance of Tm difference between vaccine and field strains and between genotypes among different runs and operators were evaluated by T-test using GraphPad Prism 7 (https://www.graphpad.com/scientificsoftware/prism/). Specificity of the assays was assessed against other IBV's genotype (IBV M41 GD, Nobilis® IB H120, H52 BI, IBV B1648, IBV 4/91 GD, Nobilis® IB IBV 4/91) and a panel of avian viruses, including Avian Influenza H9N2 (A/Chicken/Saudi Arabia/SP02525/3AAV/2000), Avian Metapneumovirus subtype A (Nobilis® TRT), Newcastle Disease Virus (Nobilis® ND LaSota), and Infectious Laryngotracheitis Virus (Nobilis® ILT). Differentiation of Poulvac IB QX and IBV-D388 by melting curve analysis using the simplex assay was performed testing three panels of 20 samples of known infectious status (allantoic fluid obtained from IBV, non-IBV and mock infected ECEs) prepared by one operator and blindly tested by two different operators at weekly intervals. The diagnostic sensitivity (DSe), the diagnostic specificity (DSp), the overall accuracy (Acc) and the vaccine differentiation were calculated. The ability of the multiplex RT-qPCR to differentiate between the two genotypes and identified vaccines Poulvac IB QX by melting curve A. Laconi, et al. Journal of Virological Methods 275 (2020) 113771 analysis was carried out following the experiment set up described above. In silico analysis confirmed the distinctiveness of the two transversion (G→T) occurred during the attenuation process of vaccine Poulvac IB QX (Fig. 1) ; as a consequence the target region of the vaccine strain has a GC content of 34,48 %, while for the field strains it ranged from 36.62 % to 41.38 % (supplementary materials 1). The predicted differences in Tm between the vaccine and the field strains based on uMelt analysis ranged from 1°C to 3°C. The GC content of the QX field strain used in this study was 37.93 %, and the predicted Tm difference between Poulvac IB QX and IBV-D388 was 1.5°C. The Mass-type specific oligo were designed to amplify a region longer than the QX-like one, characterized by a GC content ranging from 40.00 % of IBV M41 GD to 45.71 % of Nobilis® IB H120 (supplementary materials 1). uMelt analysis predicted a Tm difference between the two genotypes ranging from 3°C to 5°C, while the predicted Tm difference between the two Mass-type viruses used in this study was 1.5°C. The oligo pairs, their concentrations, and cycling protocol maximizing RT-qPCR performances, avoiding nonspecific amplification, was defined. Sequences of the selected oligo pairs, their concentration and the thermal profile are reported in Table 1 . The detection range of the RT-qPCR assay was from 10 7,83 EID 50 to 10 1,83 EID 50 for Poulvac IB QX and 10 7,31 EID 50 to 10 1,31 EID 50 for IBV-D388, showing excellent efficiency and linearity for both the viruses up to the dilution above the limit of detection (Fig. 2) , paired with excellent robustness, as CVs were consistently low both within and between runs ( Table 2) . The assay produced two distinguishable melting temperatures between the vaccine and the field strain, 72.2°C ( ± 0,2°C) and 73.2°C ( ± 0,2°C) respectively (Fig. 3A) . The melting temperatures observed were constant and reproducible among different replicates, operators and experiments, as depicted in Fig. 3B . T-tests confirmed the presence of a significant (P < 0.001) difference in the Tm between the two viruses, with a mean Tm of 72,28°C and of 73,25°C for Poulvac IB QX and IBV-D388, respectively. DSe was equal to 98.41 %, Dsp to 100 %, Acc to 98,86 % and vaccine differentiation accuracy to 100 %. To conclude, no increase in fluorescence associated with sigmoidal amplification curve was observed (data not shown) for no template controls, non QX-like IBVs nor other avian pathogens as listed in Materials and Methods. The detection range of the multiplex RT-qPCR was from 10 7,83 EID 50 to 10 2,83 EID50 for Poulvac IB QX, from 10 7,31 EID 50 to 10 2,31 EID 50 for IBV-D388, from 10 7,63 EID 50 to 10 2,83 EID 50 for IBV M41 GD and 10 7,75 EID 50 to 10 1,75 EID 50 for Nobilis® IB H120. Furthermore, the assay proved linear when tested against the four dilution panels up to the lowest dilution detected (Fig. 4) and robust, as CVs were consistently low both within and between runs ( Table 3) . The melting curve analysis confirmed the Tm previously observed for QX-like viruses, while the Tm detected for IBV M41 GD and Nobilis® IB H120 were 75.4°C ( ± 0,2°C) and 76.4°C ( ± 0,2°C) respectively (Fig. 5A) . The Tm proved to be constant and reproducible among different replicates, operators and experiment as depicted in Fig. 5B . T-test confirmed the presence of significant difference (P < 0.001) between the Tm of the two genotypes, with a mean Tm of 72,54°C and of 75,84°C for QX-like and Mass-type IBVs respectively. Furthermore, Ttest showed a significant difference (P < 0.001) between the Tm of vaccine Poulvac IB QX and its progenitor, and between vaccine Nobilis® IB H120 and Mass-type field strain IBV M41 GD. No increase in fluorescence associated with a sigmoidal amplification curve was observed for no template controls, non-QX-like and non-Mass-type IBVs nor other avian pathogens (data not shown). A slight increase in fluorescence was detected when processing an IBV B1648 (GI-14 lineage) isolate; however, it was associated with a non-sigmoidal amplification curve (linear shape) at Ct values above 38 and paired with no melting curves. In the diagnostic setting used, the assay showed DSe of 98,11 %, DSp of 100 %, Acc of 98,31 % and genotyping accuracy of 100 %. In the present study, we developed two SYBRgreen RT-qPCR assays paired with melting curve analysis for the differentiation between IBV vaccine and field strains and for genotyping. To this purpose we first developed and validated a protocol to differentiate vaccine Poulvac IB QX from QX-like field strains using three distinct vaccine makers arose during the attenuation process (Listorti et al., 2017) . Then we developed a multiplex protocol for IBV genotyping, focusing on QX-like and Mass-type (GI-1) strains, being two of the most widespread IBV genotypes (Cook et al., 2012; de Wit et al., 2018; Jackwood et al., 2012) . Both the assays showed high analytical sensitivity, however the limit of detection estimated for Poulvac IB QX and IBV-D388 differs between the two assays, with the simplex proving to be more sensitive than the multiplex approach. It is possible that the lower sensitivity of the multiplex assay is due to heterodimers assembling between QX-like and Mass-type oligo, reducing the availability of the QX-specific oligo during the annealing phase of the thermal cycle. The good efficiency and linearity at least up to one dilution above the limit of detection, together with the high repeatability observed for both the assays, suggest their suitability for viral quantification not only for diagnostic, but also for research purposes. To support, the simplex RT-qPCR has been successfully used in a study aiming to assess the tissue distribution of IBV-D388 at early time points after infection (manuscript in preparation). The same study showed that the efficiency of the protocol is not affected by the sample's matrix, as viral RNA was detected and quantified in tracheal and cloacal swabs, lungs and kidneys. The optimization of the designed oligo sets, reaction chemistries, thermal protocol and melting curve analyses avoided the presence of nonspecific amplifications, which, coupled with the in silico evaluation performed, ensured the high specificity of the assays, as demonstrated when tested against others IBV genotypes and different avian pathogens. Implementation of the melting curve analysis enables differentiation between vaccine Poulvac IB QX and field strains IBV-D388, as the melting temperatures of each virus proved constant within not only the same run, but also among different experiments and operators, and no overlaps between the melting temperatures of the two viruses were observed throughout the study. Noteworthy, the melting curve analysis proved extremely accurate, being able to identify sequence variants that differ only by three nucleotides, of which just two affect the GC content of the amplicon. QX-like field viruses will continuously evolve and the persistence in the field of Poulvac IB QX might lead to the introduction of mutations in the latter; however, the chance of the emergence of field variants showing a sequence identical to Poulvac IB QX by random mutations, or back mutations in the applied vaccine lessening the strength of the protocol, are extremely low (Laconi et al., 2019; Oade et al., 2019) . A constant difference in melting temperature was also observed between vaccine Nobilis® IB H120 and field strain IBV M41 GD; however, the unavailability of vaccine markers for Nobilis® IB H120 hinders to confidently identify it. In our study, the melting curve analysis towards vaccine differentiation proved an efficient and reliable tool, remarking the urgency of the availability of vaccine and progenitor sequences to determine markers for all IBV vaccines. This will be needed to obtain a precise knowledge of field viruses circulating in a certain area. Melting curve analysis was also successfully used to genotype QX- A. Laconi, et al. Journal of Virological Methods 275 (2020) 113771 like and Mass-type strains in a single run using a multiplex RT-qPCR approach. These two genotypes were chosen as target representing a major threat worldwide; however, this multiplex approach might be easily improved and adapted, via designing of specific oligo, to identify other IBV genotypes depending on the geographical area of interest and the epidemiological data available. Taken as a whole, the data gathered in the present study clearly demonstrate that melting curve analysis enables IBV genotyping and vaccine differentiation, representing a cost and time effective method in comparison to those in use in diagnostic laboratories worldwide. This research was financially supported by the research institute. The authors have no competing interests to declare. Fig. 4 . Standard curves obtained for the multiplex assay using serial ten-fold dilutions of vaccine Poulvac IB QX (A), field strain IBV D388 (B), Nobilis® IB H120 (C) and IBV M41 GD (D). Each point of the curves represents the mean of three replicates; branches represent standard deviation. Linear and R 2 are reported for each standard curve. Repeatability of the multiplex RT-qPCR assay for IBV genotyping. Mean Ct (three replicates), standard deviation (SD) and CV are reported for each virus according to virus, dilution, operator and week of experiment. Week 1 Week 2 Fig. 5. Melting curves observed for Poulvac IB QX, IBV-D388, Nobilis IB H120 and IBV M41 GD when tested using the multiplex approach for IBV genotyping (A). Melting temperatures among six independent runs (B). Each point represent the mean of nine replicates, branches represent standard deviation (SD) (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article). Development and evaluation of a real-time Taqman RT-PCR assay for the detection of infectious bronchitis virus from infected chickens Longitudinal field studies of infectious bronchitis virus and avian pneumovirus in broilers using type-specific polymerase chain reactions The long view: 40 years of infectious bronchitis research Detection of different genotypes of infectious bronchitis virus and of infectious bursal disease virus in European broilers during an epidemiological study in 2013 and the consequences for the diagnostic approach Induction of cystic oviducts and protection against early challenge with infectious bronchitis virus serotype D388 (genotype QX) by maternally derived antibodies and by early vaccination A sensitive, reproducible, and economic realtime reverse transcription PCR detecting avian metapneumovirus subtypes A and B Continued use of IBV 793B vaccine needs reassessment after its withdrawal led to the genotype's disappearance Efficacy and safety of an attenuated live QX-like infectious bronchitis virus strain as a vaccine for chickens Molecular evolution and emergence of avian gammacoronaviruses Attenuation, safety, and efficacy of an infectious bronchitis virus GA98 serotype vaccine author (s) Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction Vaccination against infectious bronchitis virus: a continuous challenge Two similar commercial live attenuated AMPV vaccines prepared by random passage of the identical field isolate, have unrelated sequences Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences Identification of IBV QX vaccine markers: should vaccine acceptance by authorities require similar identifications for all live IBV vaccines? Vaccine 35 Crossm attenuation of infectious bronchitis virus in eggs results in different patterns of genomic variation across multiple replicates Rapid detection and differentiation of Newcastle disease virus by real-time PCR with melting-curve analysis A simple method of estimating fifty per cent endpoints Application of high-resolution melting curve analysis for typing of fowl adenoviruses in field cases of inclusion body hepatitis Pathogenicity of a QX strain of infectious bronchitis virus in specific pathogen free and commercial broiler chickens, and evaluation of protection induced by a vaccination programme based on the Ma5 and 4/91 serotypes S1 gene-based phylogeny of infectious bronchitis virus: an attempt to harmonize virus classification A reverse transcriptase-polymerase chain reaction survey of infectious bronchitis virus genotypes in Western Europe from Supplementary material related to this article can be found, in the online version, at doi:https://doi.org/10.1016/j.jviromet.2019. 113771.