- Open Access
A reverse genetics system for avian coronavirus infectious bronchitis virus based on targeted RNA recombination
Virology Journal volume 14, Article number: 109 (2017)
Avian coronavirus infectious bronchitis virus (IBV) is a respiratory pathogen of chickens that causes severe economic losses in the poultry industry worldwide. Major advances in the study of the molecular biology of IBV have resulted from the development of reverse genetics systems for the highly attenuated, cell culture-adapted, IBV strain Beaudette. However, most IBV strains, amongst them virulent field isolates, can only be propagated in embryonated chicken eggs, and not in continuous cell lines.
We established a reverse genetics system for the IBV strain H52, based on targeted RNA recombination in a two-step process. First, a genomic and a chimeric synthetic, modified IBV RNA were co-transfected into non-susceptible cells to generate a recombinant chimeric murinized (m) IBV intermediate (mIBV). Herein, the genomic part coding for the spike glycoprotein ectodomain was replaced by that of the coronavirus mouse hepatitis virus (MHV), allowing for the selection and propagation of recombinant mIBV in murine cells. In the second step, mIBV was used as the recipient. To this end a recombination with synthetic RNA comprising the 3′-end of the IBV genome was performed by introducing the complete IBV spike gene, allowing for the rescue and selection of candidate recombinants in embryonated chicken eggs.
Targeted RNA recombination allowed for the modification of the 3′-end of the IBV genome, encoding all structural and accessory genes. A wild-type recombinant IBV was constructed, containing several synonymous marker mutations. The in ovo growth kinetics and in vivo characteristics of the recombinant virus were similar to those of the parental IBV strain H52.
Targeted RNA recombination allows for the generation of recombinant IBV strains that are not able to infect and propagate in continuous cell lines. The ability to introduce specific mutations holds promise for the development of rationally designed live-attenuated IBV vaccines and for studies into the biology of IBV in general.
Avian coronavirus infectious bronchitis virus (IBV) primarily infects the upper respiratory epithelium of chickens, causing a respiratory disease that is frequently complicated by secondary bacterial pathogens . In addition, some IBV strains affect the renal tubuli, oviduct and parts of the gastrointestinal tract, leading to pathological lesions in these organ systems, with subsequent reduced weight gain and a drop in egg production. The virus has a worldwide presence in both commercial and backyard chickens, appearing in a wide variety of geno-, sero- and protectotypes . IBV is currently regarded as one of the economically most relevant viral pathogens in the poultry industry.
Infectious bronchitis virus is the prototype gammacoronavirus in the family Coronaviridae, order Nidovirales . The enveloped virus particles have a positive-sense RNA genome of 27.6 kb (Fig. 1a) . The 5′ two-third of the viral genome comprises gene 1, divided into two large open reading frames 1a and 1b, which code for 15 nonstructural proteins primarily involved in RNA replication and transcription. The 3′ one-third of the viral genome codes for structural proteins: spike protein (S, encoded by gene 2), envelope protein (E, encoded by gene 3c), and membrane protein (M, encoded by gene 4), each located in the viral envelope. The nucleocapsid protein (N, encoded by gene 6) occurs in the ribonucleoprotein core . Interspersed between the structural genes, coronaviruses carry a variable number of genus specific accessory genes . Most of their gene products are nonstructural, and their expression is not essential for virus replication in vitro [7,8,9,10,11,12]. The IBV genome contains the accessory genes 3 and 5, encoding the proteins 3a and 3b, and proteins 5a and 5b, respectively . In addition, an open reading frame located in the intergenic region was identified between genes 4 and 5 .
The typical coronavirus spikes are formed by trimers of the type 1 membrane protein S, which is often proteolytically cleaved into two subunits, S1 and S2 [4, 14]. The glycosylated S1 domain forms the ‘head’ of the spike and contains the receptor-binding domain . Avian gammacoronaviruses typically interact with glycans on the host cell surface. IBV in particular requires α2,3-linked sialic acids for attachment and entry [16,17,18]. The S2 domain builds the remaining part of the ectodomain (the ‘stalk’), the transmembrane domain and the internally located endodomain. The S protein is the main determinant of the coronaviral host species tropism .
Many mammalian coronaviruses of the genera Alphacoronavirus and Betacoronaviruses can be propagated in cultured cells, unlike most avian coronaviruses of the genus Gammacoronavirus. IBV can, however, readily be propagated in, and isolated from, embryonated fowl eggs. During passaging in embryonated eggs adaptation occurs, often leading to attenuation. For example, IBV strain H52 represents the 52nd serial passage of a Massachusetts-like IBV strain isolated in The Netherlands , which causes embryonic death within 48 h post infection (hpi), and still has a residual virulence in young chickens. Another 68 passages resulted in the IBV strain H120, which is more attenuated and has a lower pathogenicity in young chicks. Similar serial passaging of another IBV strain of the Massachusetts serotype isolated in the USA resulted in the generation of the non-pathogenic and cell-culture adapted IBV strain Beaudette .
In order to study its characteristics, several research groups have independently developed a reverse genetics system (RGS) for IBV which allow the manipulation of its genome [22,23,24,25,26]. All these systems are based on the non-pathogenic cell-culture adapted IBV strain Beaudette, or the highly attenuated IBV vaccine strain H120. A major drawback of the use of the non-virulent IBV strains Beaudette and H120 [20, 21] is the inability to provide insights in the infection process in chickens, as these strains no longer cause a clinically relevant phenotype in vivo. Yet, the RGS has provided significant insight in the fundamentals of avian gammacoronavirus replication. Key findings include that IBV cell tropism is determined by the spike gene [27, 28], that the low virulence of IBV Beaudette is caused by changes in the replicase gene , and that one or more of the IBV accessory gene products interfere with the hosts’ interferon response [30,31,32].
Targeted RNA recombination is another reverse genetics approach, so far only developed for mammalian coronaviruses from the genera Alphacoronavirus and Betacoronavirus [11, 19, 33]. This system is based on the exchange of the spike gene by that of a coronavirus with a different host tropism, which enables subsequent selection on cells susceptible to the heterologous species . As a consequence, manipulation is limited to the last third of the coronavirus genome, covering all genes encoded that are located 3′ of gene 1, starting with the spike gene. Targeted RNA recombination has been shown to be easy in use and to allow the rescue of highly defective mutants [11, 19, 33]. However, the system is based on the ability to propagate both the donor and the recipient coronavirus in cell culture, and is hence not implementable for pathogenic IBV. This problem was solved by transfecting IBV genomic RNA into otherwise non-susceptible cells, exchanging the IBV spike gene by that of the mouse hepatitis virus (MHV) provided as part of a synthetic RNA, and by subsequently rescuing recombinant IBV from infected/transfected cells in embryonated eggs (Fig. 1). This system was successfully established to introduce marker mutations in the last one-third of the genome of IBV. The resulting recombinant viruses demonstrated growth kinetics in ovo and the in vivo phenotypic characteristics in one-day-old chickens similar to IBV wild type. The results presented here demonstrate for the first time a host species switch for an avian gammacoronavirus by exchanging the spike gene with that of the highly divergent betacoronavirus MHV. This RGS enables the manipulation of the structural and accessory protein genes from the genome of virulent IBV.
Cells, eggs, viruses & antibodies
Baby hamster kidney (BHK-21) cells (ATCC CCL-10) and murine LR7 cells KUO2000 were cultured in Dulbecco’s Modified Eagle Medium (DMEM) (BioWhittaker), supplemented with 4 mM L-glutamine (Lonza, Basel, Switzerland), 10% fetal bovine serum (FBS) (BioWhittaker) and 0.05 mg/ml gentamicin (Gibco Invitrogen), at 37 °C and 5% CO2. Virus titers in cells were obtained by determining the 50% tissue culture infective dose (TCID50) per ml at 2 days post inoculation (p.i.) according to the Spearman–Kärber method .
Fertilized specific pathogen free (SPF) white leghorn eggs (Animal Health Service, Deventer, The Netherlands) were incubated at 37.5 °C and 45–65% relative humidity. Eight-day- embryonated chicken eggs (ECE) were inoculated via the allantoic cavity unless stated otherwise, and candled twice daily. Upon embryonic death or no later than 7 days p.i., eggs were transferred to 4 °C for 16–24 h prior to allantoic fluid (AF) and chorio-allantoic membrane (CAM) collection. Virus titration in ovo was based on the determination of the 50% embryonic infectious dose (EID50) per ml, as determined at day 7 p.i. according to Reed and Muench . For the production of a virus stock, ten 8-day-old ECE were inoculated with 100 EID50, incubated for 24 h, and subsequently cooled for 16–24 h before the AF was harvested and pooled.
IBV strain H52 (Boehringer Ingelheim (BI), Ingelheim, Germany) was propagated in embryonated SPF eggs and titrated. IBV strain Beaudette (Animal Health Service, Deventer, The Netherlands) was propagated and titrated in BHK-21 cells. Mouse hepatitis virus (MHV) strain A59 was propagated and titrated in LR7 cells.
Monoclonal antibody (MAb) Ch/IBV 26.1 against the IBV S2 protein was obtained from Prionics (Thermo Fisher Scientific, Waltham, MA, USA) [37, 38]. The production of rabbit polyclonal antiserum k134 against MHV was described previously . Chicken polyclonal antiserum was derived from a SPF chicken vaccinated with IBV strain H120 (BI, Ingelheim, Germany). Secondary fluorescently-labeled antibodies Alexa Fluor 488 goat anti-chicken IgY, Alexa Fluor 568 goat anti-rabbit IgG, and Alexa Fluor 488 goat anti-mouse IgG (Invitrogen by Thermo Fisher Scientific) were stored in 50% glycerol at −20 °C.
Anti-IBV immuno histochemistry
CAMs were collected from ECEs, washed in PBS, fixed in neutral buffered 10% formalin in PBS for 24 h, stored in 70% ethanol and finally paraffin-embedded. Four micrometer sections of CAM were mounted on glass slides and subsequently deparaffinized and rehydrated in alcohol series. Next, the sections were subjected to endogenous peroxidase inactivation and antigen retrieval as described before . Sections were washed in phosphate buffered Normal Antibody Diluent (NAD, ScyTek Laboratories, Logan, USA) containing 0.1% Tween-20, and after primary antibody incubation with PBS 0.1% Tween-20. Sections were incubated for 60 min at room temperature with MAb Ch/IBV 26.1 diluted 1:100 in NAD. Antibody binding was detected by Dako Envision HRPO labeled polymer anti-mouse (Dako, by Agilent Technologies, Santa Clara, USA) diluted 1:1 in NAD, and visualized by 3-Amino-9-ethylcarbazole (AEC, Dako). Slides were counterstained with hematoxylin, mounted with Aquatex (Merck, Darmstadt, Germany), and viral antigen presence was assessed by light microscopy (BX60, Olympus, Tokyo, Japan).
RNA isolation, reverse transcription, PCR and RT-qPCR
RNA was isolated from harvested AF using the QIAamp viral RNA Mini Kit (Qiagen, Hilden, Germany) according to manufacturer’s protocol. Reverse transcription (RT) was performed using the Transcriptor First Strand cDNA Synthesis Kit (Roche, Basel, Switzerland) according to manufacturer’s protocol, with random hexamers for standard PCR, or with specific primers for sequencing and cloning purposes. PCR was performed with recombinant Taq DNA polymerase (Thermo Fisher Scientific) for plasmid characterization or with Phusion Hot Start II High-Fidelity DNA Polymerase (Thermo Fisher Scientific) for sequencing and cloning purposes.
One-step RT-qPCR was used to semi-quantitatively assess virus load in AF. Forward primer IBV.RdRp.F41 (3′-CATGCAGTTTGTTGGAGATCCT-5′) and reverse primer IBV.RdRp.R41 (3′-GTGACCTGGTTTTACCGTTTGA-5′) targeting the conserved region of gene 1b (nucleotide position 13,412 to 13,580 in IBV Beaudette GenBank accession number M95169.1) coding for the RNA-dependent RNA polymerase protein. Primers were obtained from Biolegio (Nijmegen, The Netherlands) and used at a final concentration of 300 nM each with the iTaq universal SYBR Green one-step kit (Bio-Rad Laboratories, Hercules, California, USA). The RT-qPCR reaction was carried out in a Bio-Rad CFX Connect real-time PCR system, starting with 10 min at 50 °C and 1 min at 95 °C, followed by 40 cycles of 10 s at 95 °C and 30 s at 60 °C, and ending with a dissociation step for the determination of the melting point of the obtained PCR fragment.
Construction of IBV donor plasmid
The complete genome sequence of IBV H52 BI was determined by Sanger sequencing using primers as described by Zhou et al. . The 5′- and 3′-UTR sequences were identified using the 2nd generation 5′/3′ RACE kit (Roche, Basel, Switzerland). The IBV H52 BI genome sequence was 27,640 nucleotides (nt) in length, including an annotated 10 nt polyA tail.
The design of the donor plasmids principally followed the strategy previously described by Kuo et al. . The final donor plasmid p-IBV was constructed from the stepwise ligation of fragments derived from five plasmids (Fig. 1b and Table 1, and described below in detail): Plasmid (p)IBV-5 comprises a T7 RNA polymerase promotor, 2 G nucleotides, and the near full-length 5′-untranslated region (UTR), with an unintended A to C substitution at position 54. Plasmid IBV-1b comprises the last 754 nt of the pol 1b gene, including the 50 nt overlap with the spike gene, and the first 66 nt of the spike gene, including the signal sequence. Plasmid IBV-S contains the near full length ectodomain of the spike gene, 3211 nt in length. Plasmid IBV-SIR comprises the last 212 nt of the spike gene (the transmembrane and the endodomain), the accessory gene 3, the envelope gene, the membrane gene and half of the intergenic region. Plasmid IBV-3 T comprises the 3′-terminal region of the IBV genome, including the second half of the IR, the accessory gene 5, the nucleocapsid gene, the 3′-UTR and a 100 nt poly-A sequence. All plasmids were generated by GenScript (Piscataway, NJ, USA) and provided in the plasmid pUC57-simple, a standard cloning plasmid with the polylinker removed.
In order to allow cloning of the fragments in a stepwise approach, naturally occurring restriction enzyme sites (RES) located in the viral cDNA were used, except for the BstBI site between p-IBV-5 and p-IBV-1b, which is only partly present in the 5′-UTR, and the XhoI site between p-IBV-1b and p-IBV-S, which was introduced without changing the amino acid sequence (silent mutation). Restriction enzyme sites were made unique by silently removing these RES from other parts of the genome included in the donor plasmid (Additional file 1: Table S1). In addition, semi-unique RES were introduced by silent mutations within 200 nt up- and downstream of the accessory genes 3 and 5. Finally, unique RES MssI and PacI were included after the poly-A sequence, allowing linearization of the plasmid by a single restriction enzyme digest. All genome fragments were ligated step-by-step into p-IBV-5 using the restriction enzymes specified in Table 1. Each ligation mixture was subsequently transfected into HB101 competent cells and plasmid DNA was isolated by performing midiprep DNA isolation (Qiagen, Hilden, Germany). The final plasmid consisted of p-IBV-5-1b-S-SIR-3 T, now called p-IBV (Fig. 1c). The composition of each plasmid was confirmed after each cloning step by PCR, restriction enzyme digestion and sequencing of each of the inserts (Macrogen, Amsterdam, The Netherlands).
Construction of mIBV donor plasmid
The ectodomain of the MHV A59 spike gene was amplified from pTUG  by PCR using primers with an XhoI overhang (Table 2) and ligated into pJet1.2 resulting in p-MHV-S. Site directed mutagenesis (SDM) with the Q5 SDM kit (New England Biolabs, Ipswich, USA) was used to silently remove an EcoRI and an XhoI RES interfering with subsequent cloning steps (Table 2). The ectodomain of MHV spike was ligated into p-IBV-5-1b, followed by subsequent cloning steps using the IBV fragments SIR and 3 T. This resulted in the plasmid p-IBV-5-1b-mhvS-SIR-3 T, now called p-mIBV (Fig. 1b).
Construction of N transcript plasmid
A plasmid comprising the nucleocapsid gene and 3′-UTR sequence of IBV H52 BI was generated by PCR amplifying the respective region using primers IBV-H52.N.ATG.FW and IBV-M41#2.IR.RV (Table 2). The amplicon was ligated into pJet1.2 downstream of the T7 promotor sequence, resulting in p-IBV-N, and the correctness of the insert was verified by sequencing.
In vitro transcription
Capped, run-off donor transcripts were synthesized from p-IBV, p-mIBV and p-IBV-N using the mMessage mMachine T7 kit (Ambion by Thermo Fisher Scientific). In brief, p-mIBV was PacI-linearized, and p-IBV and p-IBV-N were MssI-linearized. Linearized plasmid DNA was ethanol precipitated. Transcription reactions were prepared according to the manufacturer’s instructions, using 1.5 and 0.5 μg linearized DNA per 10 ul reaction for p-(m)IBV and p-IBV-N, respectively. After 1 h of incubation at 37 °C, production of RNA was verified by analyzing 1 μl of the reaction volume by gel electrophoreses. After an incubation of 2 h the reaction was stopped by transferring the reaction tubes to ice.
Targeted RNA recombination and rescue of mIBV
The IBV spike gene was replaced by a chimeric MHV-IBV spike gene in the IBV genome by targeted RNA recombination between p-mIBV generated donor RNA and recipient virus (IBV) RNA, as described before . IBV H52 viral RNA was transfected into BHK-21 cells, a cell line known to support replication of IBV , but not infection with IBV H52 (data not shown). Thus, 20 μl of IBV H52 BI RNA obtained from allantoic fluid, mixed with 10 μl transcript reaction mixtures of p-mIBV and p-IBV-N each were transfected into BHK-21 cells by electroporation using two pulses at 850 V and 25 μF in a Gene Pulser electroporation apparatus (Bio-Rad). Transfected BHK-21 cells were seeded onto monolayers of LR7 cells having an approximate confluence of 70–80% and incubated at 37 °C.
Two days after seeding, when syncytia in the LR7 monolayer were observed, the cell culture supernatant was harvested and rescued viruses were purified by two rounds of plaque purification on LR7 cells. Characterization of the last one-third of the genome of candidate recombinants was performed by RT-PCR and subsequent Sanger sequencing of the obtained cDNA-fragments, using the primer sets specified in Table 3. Murinized IBV (mIBV) strain #1B3-IIA was selected based on sequence analysis, and virus stocks were propagated and stored at −80 °C. Aliquots were titrated on LR7 cells.
Targeted RNA recombination and rescue of recombinant IBV
Recombinant IBV (rIBV) was generated by substituting the IBV spike ectodomain back into the mIBV genome by targeted RNA recombination between p-IBV-generated donor RNA and recipient virus mIBV. LR7 cells were infected with mIBV at a multiplicity of infection (MOI) of 0.4 for 4 h. Capped, run-off donor transcripts from p-IBV were transfected into the mIBV-infected LR7 cells by electroporation with two pulses at 850 V and 50 μF. Electroporated LR7 cells were resuspended in 2 ml DMEM (at 37 °C) and tenfold dilutions (up to 10−3) were prepared. Two hundred microliters of LR7 cell suspensions were inoculated into the allantoic cavity of 10-day-old ECEs, using 5 eggs per dilution. The eggs were candled twice daily and scored for embryonic death. Upon death, or at 7 days p.i., the eggs were transferred to 4 °C. Sixteen to twenty-four hrs later the AF was collected aseptically for RT-qPCR, and the CAMs were fixed in 10% formalin for immunohistochemistry (IHC). The AF from eggs inoculated with the highest dilution of electroporated LR7 cells, in which virus was detected by RT-qPCR and IHC, was subjected to two additional rounds of end-point dilution in 8-day-old ECE. Genetic characterization of candidate recombinants was performed by RT-PCR and subsequent Sanger sequencing of the region encoding the structural and accessory genes using the primer sets specified in Table 3.
Biological characterization of the chimeric nature of mIBV was performed by immunofluorescence (IF) double staining for IBV and MHV. BHK-21 and LR7 cells were grown on coverslips and inoculated with IBV Beaudette, and MHV A59 and mIBV #1B3-IIA, respectively, at an MOI of 1.0. Cells were fixed with PBS 4% paraformaldehyde (Aurion, Wageningen, The Netherlands) for 20 min at room temperature after 5¼, 8, and 10 hpi for MHV, IBV, and mIBV, respectively. Subsequently, cells were permeabilized with PBS containing 0.1% Triton X-100, blocked with goat serum (Gibco by Life Technologies), and incubated for 45–60 min with a combination of two primary antibodies in NGS; rabbit anti-MHV polyserum k134 diluted 1:400 and chicken anti-IBV-H120 serum diluted 1:400, or rabbit anti-MHV polyserum k134 diluted 1:400 and mouse MAb Ch/IBV 26.1 diluted 1:200. Cells were washed three times with PBS/0.05% Tween-20 and incubated in the dark for 45 min with a combination of two fluorescently labelled secondary antibodies diluted 1:200 in NGS: Alexa Fluor 488 goat anti-chicken IgY and Alexa Fluor 568 goat anti-rabbit IgG, or Alexa Fluor 488 goat anti-mouse IgG and Alexa Fluor 568 goat anti-rabbit IgG. Cells were washed three times and nuclei were stained with 300 nM DAPI in PBS for 5–10 min in the dark. Cells were washed once with milliQ and mounted with Fluorsave (Calbiochem by Merck Millipore, Billerica, MA, USA). Slides were viewed using an Olympus BX60 microscope with filters I3, A and N2.1 with a Leica DFC425C color CCD and Leica LAS-AF software (Leica Microsystems, Wetzlar, Germany).
In ovo kinetics of IBV and rIBV-wt
Eight-day-old ECEs were inoculated with 102 EID50 of IBV H52 BI or recombinant IBV wild-type (rIBV-wt). Eggs were candled twice daily and 6, 12, 24, 36, and 48 hpi. Five previously selected eggs per virus strain were transferred to 4 °C for 16–24 h, and AF was aseptically harvested and stored at −80 °C. For analysis, AF samples were thawed and tenfold diluted in PBS without Ca and Mg, and nucleic acids were extracted with the QIAamp DNA Blood Mini kit (Qiagen, Hilden, Germany) and the addition of carrier RNA, using the Hamilton Starlet pipet robot (Reno, Nevada, USA). Extracted nucleic acids were analyzed by RT-qPCR for the amount of IBV RNA following the protocol from Callison et al.  with small adaptations. Briefly, the same primers and probe were applied, while the thermoprofile was adapted for use of the ABI TaqMan Fast Virus 1-Step Master Mix (Applied Biosystems by Thermo Fisher Scientific) and the Roche 480 LightCycler. All nucleic acid samples were run and analyzed in triplicates using a tenfold dilution series of IBV H52 BI as reference for quantification. The embryonic survival at each of the time points was calculated according to the number of embryos alive at each time point compared to the total number of animals still in the experiment. A paired t-test was performed to analyse the differences in embryonic death between IBV H52 BI and rIBV-wt.
In vivo characteristics of IBV H52 BI and rIBV-wt
SPF layer-type chickens were used in the experimental infection study at Boehringer Ingelheim Animal Health Operations (BIAHO, Weesp, The Netherlands). The chickens were hatched from SPF eggs at BIAHO from eggs obtained from Charles River (BIOVO Kft, Mohács, Magyarország, Hungary). One-day-old chickens were kept in separate isolators under controlled housing conditions, including filtered supply and exhaust air.
Animals were housed in separate groups and inoculated via eye-drop with 103 EID50 in 0.1 ml of IBV H52 BI (n = 5), rIBV-wt (n = 5), or not inoculated (n = 5, negative control). Clinical symptoms monitored included ruffled feathers, decreased consciousness, depression, gasping, coughing, tracheal rales, and nasal discharge. Seven days p.i. animals were euthanized, and evaluated for their tracheal ciliary activity. To this end, the trachea was sliced into 10 transversal sections: 3 from the upper part, 4 from the middle part, and 3 from the lower part. Ciliary activity was examined by low-magnification microscopy within 2 h after sampling. Ciliostasis of each tracheal section was scored on a scale from 0 (100% ciliary activity) to 4 (no ciliary activity, i.e. complete ciliostasis), with the maximum score for each trachea being 40. Finally, the mean ciliostasis score for each group of animals was calculated. A unpaired t-test was performed to analyse whether there were differences in ciliostasis scores between IBV H52 BI and rIBV-wt.
Generation and antigenic characterization of mIBV and rIBV-wt
Viral RNA of IBV H52 BI and RNAs transcribed from plasmids p-mIBV and p-IBV-N were co-transfected into BHK-21 cells and seeded onto monolayers of LR7 cells. At 2 days post transfection, syncytia were observed in the LR7 monolayers, suggesting the successful generation of recombinant mIBV. After two rounds of plaque purification on LR7 cells, candidate recombinants were characterized antigenically and genetically.
IF staining of LR7 cells infected with mIBV showed positive staining with both anti-IBV and anti-MHV sera, indicating the chimeric nature of mIBV (Fig. 2a). IF staining with an anti-IBV-S2 MAb was positive for IBV Beaudette-infected BHK-21 cells (taken along as positive control for IF), but not for LR7 cells infected with mIBV, indicating the absence of IBV S protein in mIBV (Fig. 2b).
LR7 cells infected with mIBV and subsequently transfected with RNA transcribed from plasmid p-IBV were inoculated in tenfold dilution series into the allantoic cavity of 10-day-old ECEs. No embryonic death was observed up to 7 days p.i., but embryos in the lowest dilutions showed signs of stunting and curling typical for embryos infected with IBV. The presence of replicating recombinant IBV was demonstrated by RT-qPCR on viral RNA extracted from the AF (data not shown), and by IHC on CAM tissue (Fig. 2c). In contrast, CAMs of eggs inoculated with mIBV-infected LR7 cells (not transfected) did not show any positive signal for IBV in IHC. During the first and second passage of rIBV-wt in eight-day-old ECEs for end-point dilution purposes, infected embryos died between 2 and 3 days p.i.
Genetic characterization of mIBV and rIBV-wt
Using specifically located primers (Fig. 1 and Table 3), the intended genome structure and sequence of mIBV and recombinant IBV wild-type (rIBV-wt) were verified by RT-PCR (Fig. 3). The most important findings were that (1) murinized IBV contained the correct 5′ S gene sequence (Primer set [F] - with a forward primer located in the IBV 1b gene at a position upstream of the 1b sequence present in p-mIBV and a reverse primer located in the MHV spike gene - resulted in a detectable PCR product for mIBV, but not for IBV or p-mIBV); (2) murizined IBV contained the MHV spike gene at the location of the IBV spike gene (as both primer set [F] and [H] - each with a primer in the MHV spike gene and one in IBV gene 1b [F] or in the M gene [H] - resulted in a detectable PCR product for mIBV, but not for IBV or MHV); (3) the IBV spike gene was absent from the mIBV genome (as primer set [D] targeting the IBV spike gene resulted in a detectable PCR product for IBV, but not for mIBV or MHV(4); recombinant IBV was the result of recombination between genomic RNA from mIBV and RNA transcribed from p-IBV (as primer set [C] - with a forward primer in IBV gene 1b located upstream of the 1b sequence present in p-mIBV and a reverse primer located in the IBV spike gene - resulted in a detectable PCR product for rIBV-wt, but not for mIBV or p-IBV); (5) recombinant IBV contained the IBV spike gene at the location of the MHV-derived spike in mIBV (as both primer set [C] and primer set [E] - each with one primer in the IBV spike gene and the other in IBV ORF 1b [C] or the M gene [E] - resulted in a detectable PCR product for rIBV-wt, but not for mIBV or MHV); (6) the MHV spike gene was absent from the rIBV-wt genome (as primer set [G] targeting the MHV spike gene resulted in a detectable PCR product for MHV and mIBV, but not for rIBV-wt).
Sequence analysis of the 3′ 9 kb of the mIBV genome (starting 1 kb upstream of the BstBI RES in gene 1b, which marks the start of p-mIBV) confirmed the expected genetic identity of mIBV as observed after RT-PCR analysis (Additional file 2: Figure S1). The 3′ 9 kb of mIBV and rIBV-wt were exactly as designed, including the deliberate synonymous mutations listed in Additional file 1: Table S1. A single spontaneous silent mutation (T to C) was observed in the spike of rIBV-wt at position 22,644.
Recombinant IBV growth kinetics in embryonated eggs
In ovo growth kinetics of IBV and rIBV-wt were assessed after inoculating ECEs with 102 EID50 per egg by determination of the relative viral load in the AF of five eggs per virus at 6, 12, 24, 36 and 48 hpi by RT-qPCR. At 12 hpi, the parental wild-type IBV showed somewhat higher viral loads as compared to rIBV-wt, while from 24 hpi onwards, viral loads were comparable for both viruses (Fig. 4a). The virus titers remained at the same level until embryos started to die between 36 and 48 hpi. No differences in embryonic death between rIBV-wt and IBV H52 BI groups was observed (Fig. 4b; p > 0.05).
In vivo characteristics of rIBV-wt
The pathogenicity of rIBV-wt was compared to that of the parental IBV H52 BI strain by inoculating one-day-old SPF chickens. During the course of the infection, no clinical symptoms were observed in any of the groups (data not shown). At 7 days p.i. the animals were euthanized and the mean ciliostasis scores were determined as a readout for the ability of the respective viruses to infect and cause lesions in the primary target organ, the trachea. As expected, the negative control animals scored very low (Fig. 5), while IBV H52 BI infection resulted in a score of 26. rIBV-wt had a mean score of 19 (P > 0.05), indicating that rIBV maintained the ability to infect one-day-old chickens and induced lesions to a similar extent as the parental virus strain.
Here we developed a novel RGS based on targeted RNA recombination, which allows manipulation of the genome of virulent IBV. The resulting recombinant virus has the same characteristics as the wild type IBV H52 both in embryonated eggs and in one-day-old chickens.
By adapting the classical targeted RNA recombination approach [11, 19, 33] to IBV, the inability to culture IBV strains like H52 on continuous cell lines has been overcome. For this, a cell-line known to support the replication, but not the entry of, IBV was used. This observation was used to create, by co-transfection of IBV viral RNA and a transcript of chimeric IBV carrying the MHV spike gene, the recombinant mIBV virus. Upon transfection of synthetic rIBV donor RNA into mIBV-infected murine LR7 cells, subsequent infectious IBV virus particles could be rescued in ECE. The feasibility of the approach was demonstrated by the generation of recombinant rIBV virus carrying silent marker mutations.
The inability to select for individual IBV recombinants by plaque purification was circumvented by a combination of three approaches. First, the mIBV-infected and rIBV donor RNA-transfected LR7 cells were inoculated into ECEs by end-point dilution. Second, early RT-PCR and sequencing based screening of the genetic make-up of recombinants helped to identify and discard erroneous recombinants. Third, two subsequent end-point dilution series were executed in ECEs, each leading to the selection of rIBV-wt. Finally, the genetic identity of the 5′ 9 kb of rIBV-wt, i.e. the part of the IBV genome had been subject to manipulation, was confirmed by sequence analysis.
The replication and pathogenicity of rIBV-wt in ECE was comparable to that of the parental IBV H52 BI. Viral loads in the AF were similar with respect to maximum virus titers (Fig. 4) and embryonic death induced by both viruses did not differ (not shown). The pathogenicity of rIBV-wt in one-day-old chickens was also comparable to that of parental IBV H52 BI, as demonstrated by comparable mean ciliostasis scores at 7 days p.i. (Fig. 5). Taken together, rIBV-wt has the same properties as IBV H52 both in ovo and in vivo and can thus be used to provide insights in the infection process in chickens. Previously described RGS based on non-pathogenic IBV [21,22,23,24,25] can be used for in vivo studies only upon introduction of virulence factors including spikes from other IBV serotypes [26,27,28]. Our newly developed RGS for IBV-H52 directly allows the elucidation of factors that determine the pathogenicity of IBV, as well as studying its protective immunity in vivo.
Here, for the first time host species switching of an avian gammacoronavirus to mammalian cells was demonstrated, by exchanging the spike gene with that of the betacoronavirus MHV. Manipulation of coronavirus genomes by targeted RNA recombination using an interspecies chimeric coronavirus has been demonstrated for alpha- and betacoronaviruses, thereby switching species tropism between mammalian hosts. Our observation confirms the spike gene as the principal determinant of host species tropism of both avian and mammalian coronaviruses.
In summary, a novel reverse genetics system (RGS) based on targeted RNA recombination that allowed manipulation of the genome of virulent IBV was developed. This system makes use of an interspecies chimeric coronavirus, which is created by replacing the ectodomain of the IBV spike protein by that of MHV. The spike ectodomain exchange results in a host species tropism switch, which enables replication of mIBV in cell culture. Upon recombination of mIBV with synthetic donor RNA carrying the IBV spike gene, rIBV could be rescued in ECEs. In ovo growth kinetics and in vivo characteristics were comparable for rIBV-wt and parental IBV H52 BI, suggesting no attenuating effect of the recombination process or of the introduced synonymous marker mutations. This system will allow the introduction of mutations in the 3′ one-third of the IBV genome, allowing the manipulation of the structural and accessory genes. The use of this system for both fundamental and applied research is promising, and potentially enables the development of a new generation of rationally designed live-attenuated IBV vaccines.
Embryonated chicken eggs
Embryonic infectious dose
Hours post infection
Infectious bronchitis virus
Mouse hepatitis virus
Murinized infectious bronchitis virus
Restriction enzyme sites
Reverse genetics system
Recombinant infectious bronchitis virus
Specific pathogen free
Cook JK, Jackwood M, Jones RC. The long view: 40 years of infectious bronchitis research. Avian Pathol. 2012;41(3):239–50.
Sjaak de Wit JJ, Cook JK, van der Heijden HM. Infectious bronchitis virus variants: a review of the history, current situation and control measures. Avian Pathol. 2011;40(3):223–35.
de Groot RJ, Baker SC, Baric R, Enjuanes L, Gorbalenya AE, Holmes KV, et al. Family Coronaviridae in: Ninth Report of the International Committee on Taxonomy of Viruses. Oxford: Elsevier; 2012.
Cavanagh D. Coronavirus avian infectious bronchitis virus. Vet Res. 2007;38(2):281–97.
Masters P, Perlman S. Ch.28 Coronaviridae. In: Howley P, Knipe DM, editors. Fields Virology. 6th ed. Wolters Kluwer; 2013.
Liu DX, Fung TS, Chong KK, Shukla A, Hilgenfeld R. Accessory proteins of SARS-CoV and other coronaviruses. Antivir Res. 2014;109:97–109.
Casais R, Davies M, Cavanagh D, Britton P. Gene 5 of the avian coronavirus infectious bronchitis virus is not essential for replication. J Virol. 2005;79(13):8065–78.
Haijema BJ, Volders H, Rottier PJ. Live, attenuated coronavirus vaccines through the directed deletion of group-specific genes provide protection against feline infectious peritonitis. J Virol. 2004;78(8):3863–71.
Hodgson T, Britton P, Cavanagh D. Neither the RNA nor the proteins of open reading frames 3a and 3b of the coronavirus infectious bronchitis virus are essential for replication. J Virol. 2006;80(1):296–305.
Li C, Li Z, Zou Y, Wicht O, van Kuppeveld FJ, Rottier PJ, et al. Manipulation of the porcine epidemic diarrhea virus genome using targeted RNA recombination. PLoS One. 2013;8(8):e69997.
Yount B, Roberts RS, Sims AC, Deming D, Frieman MB, Sparks J, et al. Severe acute respiratory syndrome coronavirus group-specific open reading frames encode nonessential functions for replication in cell cultures and mice. J Virol. 2005;79(23):14909–22.
Bentley K, Keep SM, Armesto M, Britton P. Identification of a noncanonically transcribed subgenomic mRNA of infectious bronchitis virus and other gammacoronaviruses. J Virol. 2013;87(4):2128–36.
Wickramasinghe IN, van Beurden SJ, Weerts EA, Verheije MH. The avian coronavirus spike protein. Virus Res. 2014;194:37–48.
Promkuntod N, van Eijndhoven RE, de Vrieze G, Grone A, Verheije MH. Mapping of the receptor-binding domain and amino acids critical for attachment in the spike protein of avian coronavirus infectious bronchitis virus. Virology. 2014;448:26–32.
Winter C, Schwegmann-Wessels C, Cavanagh D, Neumann U, Herrler G. Sialic acid is a receptor determinant for infection of cells by avian Infectious bronchitis virus. J Gen Virol. 2006;87(Pt 5):1209–16.
Wickramasinghe IN, de Vries RP, Grone A, de Haan CA, Verheije MH. Binding of avian coronavirus spike proteins to host factors reflects virus tropism and pathogenicity. J Virol. 2011;85(17):8903–12.
Ambepitiya Wickramasinghe IN, de Vries RP, Weerts EA, van Beurden SJ, Peng W, McBride R, et al. Novel Receptor Specificity of Avian Gammacoronaviruses That Cause Enteritis. J Virol. 2015;89(17):8783–92.
Kuo L, Godeke GJ, Raamsman MJ, Masters PS, Rottier PJ. Retargeting of coronavirus by substitution of the spike glycoprotein ectodomain: crossing the host cell species barrier. J Virol. 2000;74(3):1393–406.
Bijlenga G, Cook JK, Gelb J Jr, de Wit JJ. Development and use of the H strain of avian infectious bronchitis virus from the Netherlands as a vaccine: a review. Avian Pathol. 2004;33(6):550–7.
Geilhausen HE, Ligon FB, Lukert PD. The pathogenesis of virulent and avirulent avian infectious bronchitis virus. Arch Gesamte Virusforsch. 1973;40(3):285–90.
Britton P, Evans S, Dove B, Davies M, Casais R, Cavanagh D. Generation of a recombinant avian coronavirus infectious bronchitis virus using transient dominant selection. J Virol Methods. 2005;123(2):203–11.
Casais R, Thiel V, Siddell SG, Cavanagh D, Britton P. Reverse genetics system for the avian coronavirus infectious bronchitis virus. J Virol. 2001;75(24):12359–69.
Fang S, Chen B, Tay FP, Ng BS, Liu DX. An arginine-to-proline mutation in a domain with undefined functions within the helicase protein (Nsp13) is lethal to the coronavirus infectious bronchitis virus in cultured cells. Virology. 2007;358(1):136–47.
Youn S, Leibowitz JL, Collisson EW. In vitro assembled, recombinant infectious bronchitis viruses demonstrate that the 5a open reading frame is not essential for replication. Virology. 2005;332(1):206–15.
Zhou YS, Zhang Y, Wang HN, Fan WQ, Yang X, Zhang AY, et al. Establishment of reverse genetics system for infectious bronchitis virus attenuated vaccine strain H120. Vet Microbiol. 2013;162(1):53–61.
Casais R, Dove B, Cavanagh D, Britton P. Recombinant avian infectious bronchitis virus expressing a heterologous spike gene demonstrates that the spike protein is a determinant of cell tropism. J Virol. 2003;77(16):9084–9.
Armesto M, Evans S, Cavanagh D, Abu-Median AB, Keep S, Britton P. A recombinant avian infectious bronchitis virus expressing a heterologous spike gene belonging to the 4/91 serotype. PLoS One. 2011;6(8):e24352.
Armesto M, Cavanagh D, Britton P. The replicase gene of avian coronavirus infectious bronchitis virus is a determinant of pathogenicity. PLoS One. 2009;4(10):e7384.
Kint J, Fernandez-Gutierrez M, Maier HJ, Britton P, Langereis MA, Koumans J, et al. Activation of the chicken type I interferon response by infectious bronchitis coronavirus. J Virol. 2015;89(2):1156–67.
Kint J, Dickhout A, Kutter J, Maier HJ, Britton P, Koumans J, et al. Infectious Bronchitis Coronavirus Inhibits STAT1 Signaling and Requires Accessory Proteins for Resistance to Type I Interferon Activity. J Virol. 2015;89(23):12047–57.
Kint J, Langereis MA, Maier HJ, Britton P, van Kuppeveld FJ, Koumans J, et al. Infectious Bronchitis Coronavirus Limits Interferon Production by Inducing a Host Shutoff That Requires Accessory Protein 5b. J Virol. 2016;90(16):7519–28.
Haijema BJ, Volders H, Rottier PJ. Switching species tropism: an effective way to manipulate the feline coronavirus genome. J Virol. 2003;77(8):4528–38.
Masters PS, Rottier PJ. Coronavirus reverse genetics by targeted RNA recombination. Curr Top Microbiol Immunol. 2005;287:133–59.
Spearman C. The Method of “Right and Wrong Cases” (Constant Stimuli) without Gauss’s Formula. Br J Psychol. 1908;2:227-227-242.
Reed LJ, Muench H. A simple method of estimating fifty per cent endpoints. Am J Hyg. 1938;27:493-493-497.
Koch G, Hartog L, Kant A, van Roozelaar DJ. Antigenic domains on the peplomer protein of avian infectious bronchitis virus: correlation with biological functions. J Gen Virol. 1990;71(Pt 9):1929–35.
De Wit JJ, Koch G, Kant A, Van Roozelaar DJ. Detection by immunofluorescent assay of serotype-specific and group-specific antigens of infectious bronchitis virus in tracheas of broilers with respiratory problems. Avian Pathol. 1995;24(3):465–74.
Rottier PJ, Horzinek MC, van der Zeijst BA. Viral protein synthesis in mouse hepatitis virus strain A59-infected cells: effect of tunicamycin. J Virol. 1981;40(2):350–7.
Promkuntod N, Wickramasinghe IN, de Vrieze G, Grone A, Verheije MH. Contributions of the S2 spike ectodomain to attachment and host range of infectious bronchitis virus. Virus Res. 2013;177(2):127–37.
Vennema H, Godeke GJ, Rossen JW, Voorhout WF, Horzinek MC, Opstelten DJ, et al. Nucleocapsid-independent assembly of coronavirus-like particles by co-expression of viral envelope protein genes. EMBO J. 1996;15(8):2020–8.
Callison SA, Hilt DA, Boynton TO, Sample BF, Robison R, Swayne DE, et al. Development and evaluation of a real-time Taqman RT-PCR assay for the detection of infectious bronchitis virus from infected chickens. J Virol Methods. 2006;138(1–2):60–5.
Woo PC, Lau SK, Lam CS, Lau CC, Tsang AK, Lau JH, et al. Discovery of seven novel Mammalian and avian coronaviruses in the genus deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene source of gammacoronavirus and deltacoronavirus. J Virol. 2012;86(7):3995–4008.
Geert de Vrieze, Maartje Woelders, Maloeke de Jong, and Alexandra Negatsch are acknowledged for excellent technical support.
This research was financially supported by Boehringer Ingelheim, Ingelheim, Germany.
Availability of data and materials
Data generated or analyzed during this study and described in this manuscript are included in this manuscript and its supplementary information files.
SJvB, EM, PJMR and MHV designed the studies. SJvB, AJB, AKK, GC, DS, and HCP carried out the experiments. SJvB wrote the manuscript. MHV, EM, and PR edited the manuscript. All authors read and approved the final manuscript.
SJvB and AJB perform contract research for Boehringer Ingelheim Animal Health at Utrecht University. AKK, HCP and EM are employed by Boehringer Ingelheim Veterinary Research Center, Hannover, Germany. GC and DS are employed by Boehringer Ingelheim Animal Health Operations BV, Weesp, The Netherlands.
Consent for publication
Ethics approval and consent to participate
Incubation of all embryonated chicken eggs used in this study was terminated before or at day 17 by cooling to 4 °C for 16-24 h, and embryos were destroyed before or at day 18 post incubation. No chicken embryos were grown to a viable age, and therefore no approval of the Utrecht University’s ethics committee was required.
The in vivo residual pathogenicity and efficacy study was carried out at Boehringer Ingelheim Animal Health Operations by GC and DS, The Netherlands, as permitted by the Dutch authority for animal experiments (Centrale Commissie Dierproeven) under project license number AVD224002015158 and by the Dutch authority for genetically modified organisms (Bureau GGO) under permit number IG 15–012.
Silent mutations introduced in rIBV. Nucleotide positions and sequences refer to the IBV H52 BI genome (see Additional file 2: Figure S1). Modified nucleotides in recombinant IBV wt are depicted in lower case; n.a. = not applicable. Purpose of introduction of restriction enzyme site are indicated for each site; in case of enzyme site removal the purpose was to create unique restriction enzyme sites for cloning; n.a. = not applicable in this study. (DOCX 19 kb)
Alignment of 3′ 9 kb of mIBV and rIBV-wt with IBV H52 BI. Alignment of the 3′ 9 kb of mIBV 1B3IIA P6 (excluding the MHV derived spike ectodomain sequence) and recombinant (r)IBV wild-type (wt) P4 with IBV H52 BI. Numbers refer to nucleotide positions in the IBV H52 BI genome. Restriction enzyme sites are highlighted in yellow, with the corresponding enzyme indicated above the sequences. An additional thymidine residue to keep the MHV spike gene ectodomain sequence in frame with the IBV spike gene signal sequence at position 20,385 is highlighted in green and marked with a # above the sequence. A spontaneous T to C silent substitution in the spike of rIBV-wt at position 22,644 is highlighted in red. (DOCX 41 kb)
About this article
- Avian coronavirus
- Infectious bronchitis virus
- Mouse hepatitis virus
- Targeted RNA recombination
- Reverse genetics system
- Vaccine development
- Embryonated eggs