mteb/CovidQA huggingface dataset:
Text: Bacterial artificial chromosomes (BACs) are ideally suited for the stable maintenance of large DNA sequences derived from viral genomes [1] . A considerable number of BAC systems have been established for large DNA viruses; in particular many different herpesvirus genomes have been cloned into BACs (for review see [2] ). The first BAC systems using RNA virus cDNAs were described for coronaviruses [3] [4] [5] [6] and recently the first BAC containing a full-length cDNA for a negative-stranded RNA virus was described [7] . Similarly, cDNAs corresponding to the full-length genomes of members of the Flaviviridae family (Japanese encephalitis virus [8] and Dengue virus [9] ) have been inserted into BACs.
BACs containing full-length cDNAs of pestiviruses (also within the Flaviviridae), including bovine viral diarrhea virus (BVDV) and classical swine fever virus (CSFV) have recently been established [10, 11] . Infectious pestiviruses can be rescued using RNA transcripts derived from these BACs. The pestiviruses have single stranded positive sense RNA genomes, about 12.3 kb in length, which includes a single long open reading frame, encoding a large polyprotein, flanked by 5′ and 3′ untranslated regions (UTRs) that are critical for autonomous replication of the genome [12, 13] . The polyprotein is cleaved by cellular and viral proteases into four structural proteins (nucleocapsid protein C, envelope glycoproteins E rns , E1 and E2) and eight nonstructural proteins (N pro , p7, NS2, NS3, NS4A, NS4B, NS5A and NS5B). The availability of genetically defined and stable pestivirus BACs facilitates the functional study of viral proteins or RNA structures and also the development of new marker
vaccine candidates. Several CSFV vaccines with marker properties based on chimeric pestiviruses have been developed over the years [14] . In particular, chimeric pestiviruses with substitution of the entire E2 protein have been described [15] [16] [17] but also mutants with more subtle modifications, such as the modification of the important TAV-epitope [18] within the CSFV-E2 protein [19, 20] are promising marker vaccine candidates.
Manipulation of BACs using traditional cloning procedures can be difficult (e.g. because of a lack of convenient restriction enzyme sites) and thus a range of methodologies that apply bacterial genetics, including homologous recombination (e.g. Red/ET homologous recombineering) within the E. coli host, have been developed (for review, see [21] ). The use of homologous recombination allows site-directed mutagenesis of BACs [22] and, by employing a counterselection scheme, specific modifications can be obtained without leaving residual "foreign" sequences [23] . The main advantage of this method is that there are no target limitations (e.g. based on size or location) and no need for suitable restriction sites. The integration of the modified sequence is performed in vivo (within E. coli) thereby potentially being more accurate than in vitro approaches like PCR-based methods. Although in vitro cloning approaches based on the use of high-fidelity polymerases for PCR amplification have
significantly improved in recent years, the use of in vivo approaches should allow a more accurate method of mutagenesis due to the use of the cells own high-fidelity replication system which includes proof reading. Whereas BAC recombination has been commonly used for modifying DNA viruses, there are only very few reports about the use of this technology for RNA viruses [7, 24, 25] .
Here, a generally applicable strategy for the manipulation and rescue of chimeric pestiviruses from BACs is described as a model, and the flexibility of this approach is demonstrated by generating different modifications in the viral cDNA of the new CSFV-BAC, pBeloR26, derived from the modified live vaccine strain "C-strain Riems". The targeted recombination-mediated mutagenesis described here includes the substitution of the 9 amino acid (aa) linear TAV-epitope (TAVSPTTLR) present in the E2 protein with the corresponding region (TTVSTSTLA) of a heterologous pestivirus (border disease virus, BDV, strain "Gifhorn") and also the replacement of the entire CSFV E2 protein coding region with the whole E2 coding region from the same BDV, to generate marked vaccine viruses that can be discriminated using specific anti-E2 monoclonal antibodies. The genetic stabilities of both the BAC constructs (within E. coli) and the rescued viruses have also been assessed.
Porcine kidney (PK15) and sheep fetal thymoid (SFT-R) cells were grown at 37°C (with 5% (v/v) CO 2 ) in Dulbecco's minimal essential medium (DMEM) supplemented with 5% (v/v) pestivirus-free fetal calf serum. Virus from a bait containing the modified live vaccine CSFV "C-strain Riems" (Riemser Arzneimittel AG, Germany) was propagated once in PK15 cells and termed vRiemser. RNA obtained from BDV strain "Gifhorn" [26] was used for amplification of the Gifhorn E2-coding sequence.

Oligonucleotide primers used are listed in Additional file 1: Table S1 .
The BAC construct, pBeloR26, was constructed using the long RT-PCR method as previously described [11] using RNA derived from the "C-strain Riems". Briefly, full-length viral cDNAs flanked by NotI sites were amplified by long RT-PCR using primers 5′Cstrain_T7_Not1 (which includes a T7 promotor for in vitro transcription, a NotI site and a region corresponding to the first 44 nt of the genome) and 3′CSFV_Not1 (that contains a NotI site and sequence complementary to the 3′-terminal 35 nt of the genome that are conserved among many CSFVs including the Cstrain). The product (ca. 12.3 kbp) was digested with NotI and inserted into similarly digested pBeloBAC11 (New England Biolabs, GenBank accession U51113). All BACs were modified and maintained in E. coli DH10B cells (Invitrogen) grown at 37°C in LB medium containing chloramphenicol (Cam, 15 μg/ml). The electroporation of bacteria was performed in 0.1 cm cuvettes using 1 pulse at 1800 V, 25 μF and 200 Ω in a Gene Pulser Xcell (Bio-Rad).
BACs to be used as templates for long PCR or for screening by restriction enzyme digestion were purified from 4 ml overnight cultures of E. coli DH10B using the ZR BAC DNA Miniprep Kit (Zymo Research). BACs required for direct genome sequencing were purified from 500 ml cultures using the Large-construct kit (Qiagen).
Modifications to the full-length CSFV cDNA were accomplished in E. coli DH10B (streptomycin resistant, Strep R ) using the Counter Selection BAC Modification Kit (Gene Bridges, Heidelberg, Germany).

The Red/ET recombination involved three steps (i-iii).
Step i) the temperature-sensitive pRedET expression plasmid (Gene Bridges) was introduced into electroporationcompetent E.coli DH10B cells containing the parental BAC (phenotype Cam R , Strep R ). The pRedET expresses the phage lambda proteins redα, redβ and redγ, under control of the arabinose-inducible pBAD promoter, allowing homologous recombination to occur. Immediately after electroporation, pre-warmed LB medium without antibiotics (1 ml) was added to the cells which were then incubated at 30°C for 1 hour, prior to spreading onto agar plates containing Cam (15 μg/ml) and tetracycline (Tet) (3 μg/ml) and then incubated at 30°C overnight to maintain the pRedET. The presence of the pRedET plasmid (conferring Tet R ) was verified by visual inspection of BAC-DNA preparations from the Cam R /Tet R colonies using agarose gel electrophoresis.
Step ii) counter-selection marker cassettes with an extra NotI site for screening purposes (rpsL-neo, 1325 bp) were amplified by PCR using primers with 30 nt or 50 nt extensions that were homologous to the target site in the BAC using the rpsL-neo plasmid (Gene Bridges) as template and the Phusion hot start II HF DNA polymerase (Thermo Scientific) with cycling conditions as follows: 98°C for 30s, followed by 35 cycles of 98°C for 10s, 60°C for 20s, 72°C for 60s, and 1 cycle at 72°C for 4 min. The PCR products (ca. 1400 bp) were isolated on 1% (w/v) TBE agarose gels and purified using a GeneJET gel extraction kit (Thermo Scientific). Samples (30 μl), from an E. coli culture containing pRedET and the parental BAC grown overnight at 30°C in LB media (Cam, Tet), were used to inoculate 1.4 ml of fresh LB media with the same antibiotics to obtain exponentially growing bacteria at 30°C. Red/ET recombination proteins were induced by adding 50 μl of 10% (w/v) L-arabinose (Sigma). The PCR
product (200 ng) containing the rpsL-neo cassette was introduced into these bacteria using electroporation (as above). Following electroporation, the cells were grown at 37°C for 70 min (to allow recombination) and then selected on plates containing Cam (15 μg/ml), Tet (3 μg/ml) and kanamycin (Kan, 15 μg/ml) overnight at 30°C to maintain the pRedET. Note, the rpsL cassette confers Streptomycin sensitivity (Strep S ) onto the resistant DH10B strain and the neo confers Kanamycin resistance (Kan R ). The correct phenotype (Cam R , Kan R , Tet R , Strep S ) of the resulting colonies was confirmed by streaking the colonies onto plates containing Cam (15 μg/ml), Tet (3 μg/ml) and Kan (15 μg/ml) and grown at 30°C. Importantly, for the third step, the replacement of the rpsL-neo cassette (using counter-selection), the selected colonies were also streaked onto plates containing Cam (15 μg/ml) plus Strep (50 μg/ml) and shown to be Strep S indicating incorporation of a functional rpsL gene. The
structures of the intermediate BACs were verified by restriction enzyme analysis and sequencing around the inserts.
Step iii) the replacement of the rpsL-neo selection cassettes from the intermediate constructs using linear DNA fragments was achieved through counter-selection and Red/ET recombination. Again, the homologous sequences at the ends of the DNA fragment were used for Red/ET mediated recombination events to replace the rpsL-neo cassette with the sequence of interest. Counterselection against the rpsL-neo cassette (phenotype Cam R , Kan R , Tet R , Strep S ) was employed using media containing Cam (15 μg/ml) and Strep (50 μg/ml) to isolate the required derivatives (phenotype Cam R and Strep R ).
Initially, the intermediate construct, pBeloR26_E2rpsLneo ( Figure 1 ), was generated using Red/ET recombination by insertion of the rpsL-neo cassette with an extra NotI site for screening purposes which was amplified using primers Criems-TAVfor and Criems-TAVrev (Additional file 1: Table S1 ) in place of the TAVSPTTLR coding sequence (27 nt) . Secondly, the rpsL-neo cassette in this intermediate construct was then replaced using counter-selection Red/ ET recombination using a single-stranded oligonucleotide, Riems_TAV_Gifhorn (Additional file 1: Table S1 ) with the same homology arms as used for the rpsL-neo cassette, to introduce the coding sequence for the BDV "Gifhorn" epitope sequence (TTVSTSTLA). The resulting construct was named pBeloR26_TAV (Figure 1 ). The initial intermediate construct (with rpsL-neo) was then used to produce the pBeloR26_E2gif construct ( Figure 1 ). For this, the E2 coding sequence was amplified from cDNA prepared from BDV "Gifhorn" RNA using two different
primer pairs, one set with 50 nt homology arms (Criems_E2_gifFlong/Criems_ E2_gifRlong) and another with 30 nt homologous sequences (Criems_E2_gifF/Criems_E2_gifR).
For generation of BACs with substitution of the entire E2 coding sequences, PCR products consisting of the sequence of interest flanked with homology arms identical to the target area were generated by PCR (as for the rpsLneo cassette). For making constructs with substitution of shorter sequences (e.g. the TAV-epitope), the recombination was achieved using synthetic single stranded oligonucleotides rather than PCR products. Pre-heating of single stranded oligonucleotides at 95°C for 2 min followed by snap-freezing, prior to electroporation, empirically showed the best results. In each case, the DNA molecules were introduced into E. coli containing the BAC derivatives including the rpsL-neo cassettes together with the pRedET plasmid by electroporation as described above. The structures of the modified BACs were verified by restriction enzyme analysis and subsequent full-genome sequencing (see below).
BAC DNA (1 μg) was linearized with NotI or 1 μl BAC DNA was used as template for long PCR amplification using primers 5′C-strain_T7_Not1 and 3′CSFV (Additional file 1: Table S1 ). Linearized BACs or PCR products were purified with the GeneJet PCR purification kit (Thermo Scientific) and transcribed in vitro using a Megascript T7 kit (Invitrogen). Viruses were rescued from RNA transcripts (1 to 5 μg) by electroporation of porcine (PK15) or ovine (SFT-R) cells essentially as described previously [24] . Cells were analysed using immunofluorescence microscopy (typically after 3 days) for the expression of NS3 and E2 proteins using specific monoclonal antibodies (mAbs), these were anti-NS3 (WB103/105, pan-pestivirus), anti-CSFV E2 (WH211, WH303, both CSFV specific) and anti-BDV E2 (WB166, BVDV/BDV specific) (AHVLA Scientific, United Kingdom) together with Alexa 488 conjugated goat antimouse IgG antibody (Molecular Probes, Invitrogen). The nuclei of cells were visualized using DAPI (Vector
Laboratories) and images were recorded using a BX63 fluorescence microscope (Olympus). For peroxidase staining, cells were fixed and stained for the presence of pestivirus antigens using biotinylated pig anti-CSFV/BVDV polyclonal IgG followed by avidin-conjugated horseradish peroxidase (eBioscience) as previously described [27] . The same staining procedure was also performed using the anti-E2 mAbs. Samples containing virus-positive cells were passaged onto new cells. Virus growth curves were generated as previously described [24] . Briefly, PK15 or SFT-R cells were infected at a multiplicity of infection (MOI) of 0.1 pfu/cell and grown for three days.
BAC DNAs (5 μg), purified using the Large-construct kit (Qiagen), or PCR products (1 μg) amplified from viral cDNA or from BACs using the long PCR method (as above) were consensus sequenced using a 454 FLX (Roche) or an Ion PGM (Life Technologies). Both Newbler (Roche) and the bwa.bwasw alignment algorithm [28] were used for mapping the reads to the expected sequence. A combination of Samtools [29] and LoFreq SNV-caller [30] was used for downstream single nucleotide variant (SNV) analysis. Finally, clone consensus sequences were aligned using MAFFT in the Geneious software platform (Biomatters).

Generation of a BAC containing full-length cDNA corresponding to the modified live vaccine "C-strain Riems"
BACs containing the full-length cDNA corresponding to the parental vRiemser ("C-strain Riems") were constructed according to the method described previously for the "Paderborn" strain of CSFV [11] . BACs containing the complete CSFV cDNAs were identified by restriction Figure 1 Schematic representation of the CSFV genome organization and the BACs constructed and used in this study. Nucleotide (nt) and amino acid (aa) positions within R26 for the 5′ and 3′ termini together with the translational start and stop codons of the polyprotein coding region plus cleavage sites used to make the individual proteins (N pro , C, E rns , E1, E2, p7, NS2, NS3, NS4A, NS4B, NS5A and NS5B) are indicated. Insertion of the rpsL-neo in place of the TAV-epitope within CSFV E2 for the intermediate construct (R26_rpsLneo) and the subsequent replacement with the TTVSTSTLA sequence (R26_TAV) and the complete substitution of the E2 sequence (R26_E2gif) are shown. Names of BAC constructs begin with "pBelo" and
rescued viruses with "v" (e.g. pBeloR26 and vR26). Cell culture passage no. of virus is indicated with "/P" (e.g. vR26/P-4).
digest analysis and following linearization by NotI, RNA transcripts were produced and electroporated into PK15 cells. This screening resulted in the identification of a BAC containing a cDNA insert of 12316 nt, pBeloR26 (Figure 1) , which yielded infectious virus, termed vR26, that could be propagated in SFT-R cells (Figure 2 , upper panels) and in PK15 cells (Figure 3 ). The rescued vR26 displayed higher growth rate at the early stage (about 10fold difference in virus yield at 24 h) compared to the parental vaccine virus, but after 48 hours similar virus titres were obtained (Figure 3 ). Full-genome sequencing of the cloned BAC template, pBeloR26, revealed a number of differences throughout the genome when compared to the full-length consensus sequence of the cDNA used for the cloning procedure (see Table 1 ). These differences are non-representative variants within the cDNA. Overall, the BAC sequence differed from the cDNA sequence in 18 positions, 9 of these lead to predicted
amino acid substitutions within the polyprotein; one in each of N pro , E rns , E1, E2 and NS3 and four amino acid substitutions in NS5B (Table 1) . When compared to the published reference sequence (GenBank accession AY259122.1), the pBeloR26 BAC sequence differed at an additional 11 positions, 1 of these lead to a predicted amino acid substitution and there was one large insertion (27 nt) in the hypervariable region of the 3′-UTR (Additional file 2: Table S2 ).
To determine the utility of the targeted recombinationmediated mutagenesis system for pestiviruses, two different modifications of the E2 protein coding sequence within pBeloR26 were generated using the Red/ET recombination methodology. Initially, the sequence encoding the linear TAV-epitope (TAVSPTTLR) within the CSFV-E2 was substituted with the sequence encoding the corresponding region (encoding TTVSTSTLA) from the BDV strain "Gifhorn" as described in the Materials and Methods section. More than 90% of the colonies obtained using this procedure contained the required BAC
Anti-CSFV E2 (WH211) Figure 2 Antibody reaction patterns of pestivirus infected cells. SFT-R cells were infected with vR26 and its two derivatives vR26_E2gif and vR26_TAV plus vGifhorn [26] . After 72 h, the cells were fixed and stained with monoclonal antibodies against the NS3 protein (WB103/105, left column), the CSFV E2 protein (WH303 and WH211, middle columns) and the BDV E2 protein (WB166, right column) as indicated and viewed using a fluorescence microscope. structure as determined by NotI digestions. The complete genome sequences of the CSFV cDNA within two selected BACs, designated pBeloR26_TAV have been verified (data not shown). In addition, the complete coding sequence (1119 nt) for the CSFV-E2 protein was substituted by the corresponding sequence from BDV "Gifhorn". Again more than 90% of the colonies obtained contained the required BAC and the same proportion of correctly recombined BACs was obtained using either 30 nt or 50 nt homology arms. The chimeric BAC was
designated, pBeloR26_E2gif and the complete virus genome sequence (cDNA) was verified (data not shown).
After electroporation with RNA transcripts derived from either pBeloR26_TAV or pBeloR26_E2gif a large number of CSFV NS3-positive cells could be observed (data not shown) and chimeric virus stocks, termed vR26_TAV and vR26_E2gif, were generated after further passages in cells. Cells infected with these viruses and with the parental vR26 and vGifhorn strains were all stained with mAbs directed against the NS3 protein ( Figure 2 ). However, in contrast to the parental vR26 virus, the chimeric viruses rescued from the recombined BACs were not recognized by anti-E2 mAbs specific for the CSFV-E2 proteins ( Figure 2 ) and thus, consistent with their structure, displayed the same antibody reaction pattern as vGifhorn. Two different anti-CSFV E2 mAbs, WH211 and WH303, were used for the staining and the latter has been shown previously to target the TAV-epitope [18] . As anticipated, cells infected with either the vGifhorn or with the chimeric vR26_E2gif could be shown to express the "Gifhorn"
E2 protein using staining with an anti-BDV mAb ( Figure 2 ). The presence of the BDV epitope TTVSTSTLA in vR26_ TAV was insufficient to permit efficient recognition by this anti-BDV mab, although a weak signal was observed in some cells.
The BAC constructs pBeloR26 and pBeloR26_E2gif were analysed for the genetic stability of the cDNA to determine the suitability of the BAC vector for maintaining full-length pestivirus cDNAs. E. coli DH10B cells containing the BACs were passaged 15 times, by overnight growth, and the complete viral cDNAs within the BACs were sequenced after the 1st and the 15th passage. No mutations were observed within the 12316 nt virus cDNA sequences after this extensive propagation of the BACs in the bacterial host, indicating a highly stable system for the maintenance of complete pestivirus cDNA sequences.
The viruses, vR26 and vR26_E2gif, rescued from their respective BAC constructs, were also tested for their genetic stability within mammalian cells. Linearized BAC DNA was transcribed in vitro and the RNA was electroporated into PK15 cells. Three days after electroporation the cells were stained with the anti-NS3 antibody to detect the presence of replicating virus. Samples containing virus positive cells were passaged onto new cells, this process 

*Nt position 10665 in vR26/P-12 is reverted from A to G as in the parental cDNA.
was repeated for 12 separate passages (each of three days). The virus titre (as TCID 50 /ml) was determined for each passage. Passage of the rescued vR26_E2gif chimeric virus in PK15 cells resulted in rapidly decreasing virus titres and was discontinued after the 2nd passage ( Figure 4A ). Instead, further passage of this chimeric virus was performed in ovine SFT-R cells (the preferred cell type for BDV) and resulted in much higher titers of the chimeric virus. Virus titers reached more than 10 6 TCID 50 /ml after the 1st passage and remained stable for 12 passages ( Figure 4A ). The rescued vR26 was also efficiently propagated on the SFT-R cells but maintained a slightly lower titer than the vR26_E2gif chimeric virus ( Figure 4A ).
To check that the viruses retained their antibody reaction properties ( Figure 2 ) after these passages, cells were infected with viruses from the 12th SFT-R cell culture passage (termed vR26/P-12 and vR26_E2gif/P-12) and stained with a polyclonal anti-pestivirus serum and with specific mAbs directed against the CSFV-E2 and BDV-E2 proteins ( Figure 4B ). Cells infected with either the vR26/P-12 or the chimeric vR26_E2gif/P-12 were each detected by the polyclonal anti-pestivirus serum as expected. The anti-CSFV-E2 mAb specifically detected cells infected with vR26/P-12 but not cells infected by the chimeric virus containing the BDV-E2 protein (consistent with the results shown in Figure 2 ). In contrast, the anti-BDV-E2 mAb specifically detected infection by the vR26_E2gif/P-12 and did not recognize cells infected with vR26/P-12. Each result is in accord with the structure of the viruses. The 4th passage of vR26 (vR26/P-4) displayed a slower growth rate than the virus obtained after 12
passages (see Figure 5A ). It also had a reduced growth rate compared to both the vR26_E2gif/P-4 and vR26_E2gif/P-12. The fulllength sequence of pBeloR26 had revealed ten non-silent mutations compared to the reference sequence (AY25 9122.1) for this virus (Additional file 2: Table S2 ). Any of these mutations could be responsible for the impaired growth acting alone or in concert. For further investigation of this issue, full length cDNAs prepared from vR26/ P-4, vR26/P-12, vR26_E2gif/P-4 and vR26_E2gif/P-12 were deep-sequenced using both the 454 FLX and Ion PGM platforms for comparison and to determine the quasispecies distribution (Additional file 3: Figure S1 and Additional file 4: Figure S2 ). Sequencing data from both platforms revealed that both the vR26/P-12 and vR26_E2gif/P-12 were close to 100% changed at nt position A10665G compared to the BAC clones (resulting in the predicted amino acid substitution D3431G within the NS5B protein, the RNAdependent RNA polymerase, see
Figure 5B ). This adaptation is a reversion back to the consensus cDNA sequence of the parental vaccine virus, vRiemser (Additional file 2: Table S2 ). Additionally, vR26/P-4 and vR26_E2gif/P-4 already showed evidence for this reversion being present within the population. For vR26/P-4, the level of reversion was 57%, while for vR26_E2gif/P-4 the extent of change was 73% (see Figure 5B ).