Full genome sequencing of archived wild type and vaccine rinderpest virus isolates prior to their destruction

Nature

Full genome sequencing of archived wild type and vaccine rinderpest virus isolates prior to their destruction"


Play all audios:

Loading...

ABSTRACT When rinderpest virus (RPV) was declared eradicated in 2011, the only remaining samples of this once much-feared livestock virus were those held in various laboratories. In order to


allow the destruction of our institute’s stocks of RPV while maintaining the ability to recover the various viruses if ever required, we have determined the full genome sequence of all our


distinct samples of RPV, including 51 wild type viruses and examples of three different types of vaccine strain. Examination of the sequences of these virus isolates has shown that the


African isolates form a single disparate clade, rather than two separate clades, which is more in accord with the known history of the virus in Africa. We have also identified two groups of


goat-passaged viruses which have acquired an extra 6 bases in the long untranslated region between the M and F protein coding sequences, and shown that, for more than half the genomes


sequenced, translation of the F protein requires translational frameshift or non-standard translation initiation. Curiously, the clade containing the lapinised vaccine viruses that were


developed originally in Korea appears to be more similar to the known African viruses than to any other Asian viruses. SIMILAR CONTENT BEING VIEWED BY OTHERS GLOBAL GENOMIC SURVEILLANCE OF


MONKEYPOX VIRUS Article Open access 23 October 2024 PHYLOGENETIC MOLECULAR EVOLUTION AND RECOMBINATION ANALYSIS OF COMPLETE GENOME OF HUMAN PARECHOVIRUS IN THAILAND Article Open access 21


April 2021 IDENTIFICATION OF CRF89_BF, A NEW MEMBER OF AN HIV-1 CIRCULATING BF INTERSUBTYPE RECOMBINANT FORM FAMILY WIDELY SPREAD IN SOUTH AMERICA Article Open access 01 June 2021


INTRODUCTION Rinderpest (RP) was one of the most severe diseases of cattle ever recorded, with high morbidity rates, and mortality rates of 80% to 90% in naïve populations. The disease was


declared eradicated in 20111, thus becoming the second viral disease, after smallpox, to be eradicated, with global benefits estimated to be in the billions of dollars2. The RP virus (RPV)


itself has not been entirely eliminated, with a number of laboratories known to have samples of wild type RPV. Accidental release of RPV from such a laboratory is thought to be the most


likely pathway by which the virus might re-enter the environment3,4, although it might also be deliberately released as an act of sabotage or bioterrorism. The member states of the World


Organisation for Animal Health (OIE), and the Food and Agricultural Organisation of the United Nations (FAO) agreed to restrict all work with the virus and to allow the storage of the virus


only in highly secure Rinderpest Holding Facilities (RHFs) that have been inspected and approved jointly by OIE and FAO. The FAO-OIE RHF in the UK is the Pirbright Institute which, as the


Institute for Animal Health, and before that the Animal Virus Research Institute, was a centre for research on RPV since the 1960s. The institute developed the monoclonal antibody-based


competition ELISA (cELISA) used extensively for surveillance5, as well as the most widely used RT-PCR assay for RPV6 and the concept of different geographic lineages of the virus based on


the sequence of the product from the RT-PCR7. The first RPV genome sequences were determined at Pirbright8,9 and the system for recovering RPV from a cDNA copy of the genome was developed


there10. Because of this history, and its links to many of the countries where RPV was still endemic in the latter half of the 20th century, the institute had accumulated a significant


number of RPV isolates from a range of countries. Some of these were tissue samples while other isolates had been grown in cell culture for various research purposes. Most of these isolates


had not been subjected to extensive characterisation, either as pathogens or at the molecular level. Full genome sequences have only been determined for the “Plowright” tissue


culture-attenuated vaccine strain8 and the virulent virus from which it was derived9. Simply destroying all stocks of wild type virus would pose the risk that information would be lost that


might one day be useful, since other uncharacterised morbilliviruses are known to exist11 and may have the potential to move into the environmental niche presented by a global population of


cattle lacking immunity to these viruses. One way to mitigate this risk would be to sequence the genomes of these viruses prior to their destruction. The system for recovery of live RPV from


a copy of its genome10 is well-established, and has been used to create a large number of recombinant RPVs over the years e.g.12,13,14,15,16. Current DNA synthesis techniques are such that


a complete RPV genome could be built into the appropriate plasmid, as was recently done to recover live peste des petits ruminants virus (PPRV) from a cDNA copy of the genome built entirely


_de novo_17. Determining the full genome sequence of viruses in our existing archive would enable any of those viruses to be recreated should they ever be required in the future, meaning it


was no longer necessary to keep that actual virus. It would also provide a database of sequences that might be useful in tracing the origin of any potential future outbreak of RP, as well as


information about the evolution of the virus over time. The project to determine these sequences, followed by destruction of the virus samples was therefore proposed to, and approved by,


OIE and FAO. We present here the results of that project, which have identified several new and unexpected features of the genome of some RPV isolates, as well as improving our understanding


of the spread of the virus in Africa. RESULTS AND DISCUSSION SEQUENCING LIBRARIES All samples were first screened by reverse transcription real time PCR (RT-qPCR) specific for RPV18.


RPV-positive RNA samples were processed to create sequencing libraries which were sequenced using an Illumina MiSeq. RPV is a morbillivirus, an enveloped RNA virus with a negative sense


genome of 15882 bases8. As with all the morbilliviruses, it is pleiomorphic and cannot be easily purified free of host cell material19,20. Both tissue samples and cultured virus were


expected to contain a high percentage of host cell RNA. Libraries prepared using standard Nextera kits from cell-cultured RPV contained a highly variable fraction of RPV RNA (median: 5%;


range: 0.01–70%), while for tissue samples, even those giving similarly low Ct in the RT-qPCR assay (indicating high RPV content), the fraction of total RNA derived from RPV was lower


(median: 0.16%; range: 0–84%). Because of the uneven distribution of reads along the genome (Fig. 1), at least 3000 read pairs were required to give good coverage of the majority of the


genome, and most tissue samples gave many fewer RPV-specific reads than this. Several techniques were used in attempts to improve the fraction of total tissue RNA that was derived from RPV.


The most effective technique was amplification of total RNA using a single primer isothermal amplification system (SPIA), followed by partial depletion of host rRNA using a human


sequence-optimised system (see Methods), giving an approximately 35-fold improvement (mean 34.72, s.d. = 18.7) in the specific RPV content in the MiSeq libraries for the samples that were


analysed by both methods. Supplementary Table S1 gives the sample preparation method and sequencing results (total number of reads and number matching RPV) for each sample sequenced.


ASSEMBLING RPV GENOME SEQUENCES RPV genome assembly was performed by mapping the sequence data to an existing RPV genome, the wild type RPV Kabete ‘O’ sequence9 (RPV-KO). The program


_bowtie2_21 was the most effective at identifying reads mapping to heterologous RPV isolates, while _bwa-mem_22 was more effective at identifying and incorporating data from reads that were


derived from viral copy-back RNAs. Since all reads covering the ends of the genome were derived from such copy-back RNAs, each program provided information not available with the other. Both


mappers were therefore used, combining the information to give the final consensus sequence. Assembling the RPV genome by _de novo_ assembly using any of several existing programs was not


as effective. Average sequence coverage from the NGS data is shown in Fig. 1. For almost all isolates, one or two regions of the long GC-rich M-F UTR were not determined from the NGS data,


even when sequencing the RPV-KO isolate itself, so this was not a problem caused by mapping to a heterologous RPV isolate. Given that the sections of genome not found in the sequencing


library were among the most GC-rich sections of the genome sequence (Fig. 1), it is likely that this problem was due to a failure of cDNAs containing these sequences to be effectively


amplified during the PCRs used to attach barcodes and adapters during library preparation, as has previously been reported23. The genome sequence in these regions was therefore determined by


RT-PCR and sequencing the products by Sanger sequencing. Where the 5′ end or the 3′ end of the genome were not recovered from the Miseq dataset, the missing information was obtained by RACE


(see Methods) and Sanger sequencing. SEQUENCE FEATURES IN THE COMPLETED GENOMES A total of 121 full genomes were determined plus 2 more genomes lacking only ~40 bases at the 5′ end. Of


these 123, 10 were preparations of the RBOK vaccine strain from different sources and a further 11 genomes were preparations of lapinised RPV in use at different RP research laboratories,


specifically the East African Veterinary Research Organisation at Muguga, Kenya (EAVRO), the Plum Island Animal Disease Centre in the USA (PIADC) and the predecessor of the Pirbright


Institute, the Animal Virus Research Institute (AVRI), while 11 were preparations of goat-adapted vaccine from different sources. The remaining 91 genomes represented 51 discrete isolates of


wild type RPV. All the genomes were 15882 bases long, as previously reported8,9,24,25,26,27, except for 5 samples of goat-adapted vaccine virus (GtVacc), each of which had an extra 6 bases


in the long GC-rich 5′ UTR of the F gene (Fig. 2). These samples could be divided into two groups, GtVacc from Bangladesh (vaccine seed and production vaccine) and a sample of GtVacc from


India; for the latter, we sequenced some of the original material and also freeze-dried tissue from a UK goat that had been inoculated with this material. The Indian and Bangladeshi vaccines


had slightly different sequence modifications (Fig. 2), suggesting either two independent insertion events at the same point or an unstable insertion event during goat passage of the


vaccine virus which resolved in two different ways. Other samples of GtVacc prepared in Kenya, or grown at Pirbright from samples sent from Kenya, did not have this insertion, showing that


the insertion event occurred in India after the GtVacc was transferred to Kenya. Isolates of other morbilliviruses have been found with additional bases, always in multiples of six and


usually in the M-F UTR. A variant of peste des petits ruminants virus (PPRV) with an insertion of six bases in the F gene 5′ UTR was recently identified in China in 201328. Similarly,


several variants of measles virus (MV) have been found with a net insertion of six bases29,30,31. The mechanism of how these insertions and deletions occur is not yet clear, though it has


been suggested from studies of MV genome variants that it is the result of errors of the viral polymerase when transcribing regions with extended homopolymers30. Most of the regulatory


elements in the virus genome sequences (promoter sequences, gene start and stop sequences) were highly conserved. A notable variation was in a group of Middle Eastern isolates from the 1980s


and 1990s (RPV/Oman/79, Saudi/81, Yemen/81, Lebanon/82, Kuwait/83, Iraq/85, Turkey/92, Iran/94, Iran/95), where the otherwise conserved H-L intergenic region (CGT, 9197–9) had changed to


CTT, the same sequence as all the other intergenic regions in the virus. This may have had an effect on transcription of mRNA encoding the viral RNA-dependent RNA polymerase (L protein), and


it has been recorded that several members of this group showed particularly high virulence32. Another notable sequence variation was the absence of a classical start codon for the F protein


of a group of viruses found in sub-Saharan Africa in the period 1983–93. The start codon for the F protein is normally assumed to be that at 590–2 of the F gene transcript: this is


immediately followed by sequence encoding a classic hydrophobic signal peptide and cleavage site33, followed by the highly conserved sequence at the start of the F2 peptide of the F


protein34 (Supplementary Fig. S1). However, the genome sequences found in the viruses circulating in Nigeria in 1983, Egypt in 1984, Kenya in 1988–91 and Sudan in 1992–3 (Egypt/84,


Kenya/Kajiado/88, Kenya/Kiambu/88, Kenya/Ngong/88, Kenya/Olentoko/89, Kenya/Suswa/88, Kenya/WPokot/86, Kenya/WPokot/89, Kenya/WPokot/91, Nigeria/Tambo/83, Nigeria/Yankari/buffalo/83,


Sudan/93/RBsS, Sudan/Wakobu/92/RBS) have AUA (normally isoleucine) at this position instead of AUG. Interestingly, in these genomes there is no upstream AUG in the correct reading frame to


give rise to the F protein, although all of these isolates (and only these isolates) have an AUG at 514–6, which is in the wrong reading frame (−1 relative to the F protein ORF). This


contrasts with the F gene of RPV/Nigeria/Sokoto/1964, which has the even less efficiently used ACA codon at 590–2, but has an in-frame AUG codon just upstream at 545–7 (Supplementary Fig. 


S1). These data suggest that this group of viruses had to rely on abnormal translation or translation initiation in order to generate F protein. The initiation of translation of the F


protein from F gene mRNA is unusual in several morbilliviruses. The first RPV genomes sequenced8,9 had an additional in-frame AUG codon well upstream of the putative leader peptide sequence,


at 320–322, and it was not clear which AUG was used for translation. Some, but not all, RPV isolates also have upstream AUGs which might give rise to F proteins with extended peptides


before the signal peptide proper, e.g. India/Bison/89, SriLanka/87 and Korea/Fusan-B at 269–71, and the entire set of lapinised viruses, which have the first AUG at 152–4 of the F gene


transcript, translation from which would give rise to a very extended F signal peptide, but would avoid initiation at the AUG at 561–3 (reading frame +1 relative to the F protein), also


found in all the lapinised virus sequences. In the sequences presented here, about half the genomes have upstream AUGs that are in the wrong reading frame to give rise to the F protein; in


addition to the African isolates mentioned above, RPVs Afghanistan/95, India/Ajmer/54, India/Bangalore/72, India/Bison/89, India/Bombay/54, India/HillBull, India/Hissar/53, India/Ranipet/73,


India/Ranipet/80, Oman/79, Pakistan/85, Russia/89, Russia/Tuva/92, SriLanka/87, and Turkey/Pendik/49 all have an AUG codon at 91–3 (reading frame −1 relative to the F protein). In all these


cases, as with the African viruses from late 1980s/1990s, F protein expression would require either a translational frameshift35 ocurring between the out-of-frame AUG and the coding


sequence for the signal peptide, or leaky scanning to get past the incorrect AUG with, in some cases, translation initiation from an AUA codon in order to generate the F protein, despite the


low efficiency with which this codon is used36. A third possibility is that the long 5’UTR sequence of the F gene transcript has the ability to direct the ribosome to start translation from


a particular codon: work on MV37 showed that the long UTR directs translation initiation to the second available AUG in frame with the F protein ORF, ignoring the first. The dependence on


abnormal translation would be expected to lead to reduced expression of the F protein, which may be important for viral fitness. Although the RPV/Egypt/84 virus was recorded as being


particularly mild32, there is no indication from the literature that all of the isolates with the variant start codon were particularly attenuated. Studies on MV38 and on canine distemper


virus (CDV, another morbillivirus)39 showed that removing the long UTR in this region increased F protein expression. In the case of CDV it was shown that this also severely attenuated the


virus, suggesting that limiting F protein expression is important for pathogenesis. The 3′ ends of the genome and antigenome act as the promoters (RNA polymerase binding sites and sites of


transcription initiation) for transcription of the antigenome and genome respectively and are referred to as the Genome Promoter (GP) and Antigenome Promoter (AGP). The GP also acts as the


promoter for the transcription of viral mRNAs. The AGP was highly conserved across all RPV strains sequenced, with the first 18 bases of the genome completely conserved in all the isolates


sequenced, as were 41 out of the first 50 bases. The GP is less conserved, with a variant base at position 5 and another at position 12, and only 34 conserved in the first 50 bases. The


vaccine strains of both RPV and PPRV have a G at position 5 of the GP, while their virulent parents and most virulent strains have A, leading to the suggestion that this is an attenuating


mutation40,41. However, we found a G in this position in the virulent viruses isolated in Russia in 1989 and 1992, and in the cattle-passaged Fusan parent of the Nakamura III lapinised


virus24, suggesting that this mutation is not attenuating by itself, or that it can be compensated for by other mutations elsewhere in the genome. We consistently observed small numbers of


changes in the viral sequence during cell culture passage; it was possible to observe a change in sequence from one base, though a mixture, to a different base. These changes were, however,


few in number, and usually silent: over 12 passages in cell culture, the RPV/Lap/AVRI strain showed only 11 consistent changes, and 8 of these were silent. The set of RBOK vaccine strain


sequences also showed few differences, perhaps because they were already adapted to cell culture. Out of 12 positions showing variations in at least 3 samples, 7 were silent, or led to


homologous changes in amino acid. All the other variations were in the F or H glycoproteins, and may reflect adaptation from the original bovine kidney cells42 to Vero cells (as used to grow


the vaccine in later stages of the eradication programme). COPY-BACK RNAS IN VIRUS SAMPLES The NGS datasets revealed that chimeric RNAs, that is RNA transcripts that mapped to more than one


part of the genome, were found in all preparations of RPV, whether cell cultured virus or infected tissue. These chimeric RNAs appeared to arise from copy-back events, as the supplementary


alignment was always on the opposite strand to the primary alignment (for examples, see Fig. 3). The fraction of RPV-derived cDNAs in the sequencing libraries that had supplementary


alignments varied from 0.9% to 16.4% (mean = 9.4%, s.d. = 2.7%), and this fraction did not differ significantly between cell cultured virus and infected tissues. The chimeric RNAs appeared


to be transcribed from all parts of the genome, the amount of chimeric RPV RNA closely following the overall pattern of RPV RNA in the samples (Fig. 4). However, some samples showed strong


peaks of chimeric RNAs at specific positions along the genome, notably at the AGP (compare Fig. 4a,c,d with Fig. 4b), suggesting that these samples contained a copy-back defective


interfering particle (DI) that included the AGP at each end and was thus replicated efficiently. Preparations containing a DI based on the GP could also be identified (Fig. 4e,f). Some


preparations also showed strong peaks at internal points along the genome (Fig. 4f); these regions do not contain promoters and so copy-back RNAs containing these sequences should not be


amplified, and further investigation will be required to identify the exact nature of the chimeric RNAs in these isolates, and whether this pattern is seen in related viruses such as measles


virus. Copy-back DIs have long been known to be produced during the replication of RNA viruses, and are normally seen in virus preparations that have been passaged in cell culture at too


high a multiplicity of infection (m.o.i.)43,44. Most of the samples found to contain the signature of a replicable copy-back DI were indeed from cell culture passaged virus; the exceptions


were the tissue samples taken from an animal infected with RPV/Saudi/81 (Fig. 4d), adding RPV to the list of viruses including influenza virus45, dengue virus46 and West Nile virus47 for


which DIs have been found in natural infection. The only published RT-qPCR assay for RPV18 was developed after the virus was declared eradicated, and so was not subject to the kind of global


testing undergone by the simple RT-PCR assays in use during the control and eradication programme6. It was therefore useful to assess this assay against the large set of RPV isolates in


this study. About a third of isolates had a G at position 10 of the forward primer instead of the published A/T (Fig. 5a). We also found mismatches in the reverse primer for


RPV/Kenya/kudu/95 and RPV/SriLanka/87 and in the probe for RPV/Kuwait/83, Russia/89 and Russia/92, and for RPV/Sudan/Nyala/Reedbuck/85 (Fig. 5a). This last mismatch, located close to the 3′


end of the probe, was the only one to have a major effect on the assay (Fig. 5b), possibly because it was coupled in this virus with a mismatch near the 3′ end of the forward primer; the


reaction was obviously much less efficient, giving a very high Ct despite the NGS sequencing library having a high content of RPV sequence; this deviation from the consensus did not


completely prevent the detection of this isolate of RPV, but did reduce the sensitivity of the assay when that isolate was the target. The viruses which had a mismatch closer to the 5′ end


of the probe (RPV/Kuwait/83, Russia/89 and Russia/92) also showed lower efficiency in the RT-qPCR (decreased slope of the amplification plots) (Fig. 5c), but the differences here were


relatively minor. All other isolates showed essentially normal amplification efficiency (not shown) despite a difference of one base in the forward or reverse primer. This assay has been


adopted at Pirbright in its role as OIE Reference Laboratory for Rinderpest and FAO World Reference Laboratory for Ruminant Morbilliviruses, albeit with modification to improve probe binding


(Methods). The assay clearly works for almost all virus isolates, but the possibility of a naturally occurring variant having decreased detection efficiency should be borne in mind.


PHYLOGENETIC ANALYSIS OF THE RPV GENOMES The genome sequence was partitioned into coding sequence (CDS) and untranslated regions (UTR), and the CDS further partitioned by codon position,


giving a total of 25 partitions which were then grouped using _PartitionFinder2_48 into 6 groups of partitions, the members of each group having essentially the same parameters for the best


fit evolutionary model (see Methods for details). The partitions, and the groups they were assigned to, are given in Tables 1 and 2. The UTRs were mostly grouped together, the consistent


exception being UTR4, encompassing the long GC-rich region between the M protein CDS and the F protein CDS; in addition to having a high C/G content (Fig. 1), this region was notable for a


very strong strand-specific A/T bias, with an A:T ratio of 4.8 for the antigenome strand, compared to ~1 for the other UTRs (Group 1) and for the groups containing CDS codon positions 2 and


3 (Groups 4 and 5); the CDS codon positions 1 showed the slight bias towards A over T (A:T = 1.8) that has previously been reported for a large number of eukaryotic CDS49. A notable


exception to the partitioning of codon positions 1 and 2 was found for CDS2. CDS2 encodes the P protein (a structural protein which links the nucleocapsid protein (N) to the viral polymerase


(L)) and also two non-structural proteins, C and V. The C protein is encoded in an alternate reading frame such that codon position 2 for the P protein open reading frame is codon position


1 for the C protein open reading frame; these overlapping reading frames are probably the reason why CDS2_pos1 and CDS2_pos2 group together, and separately from most of the other CDS


position 1s. For the phylogenetic analysis, duplicate sequences were removed and independent samples of the same isolate replaced with a consensus sequence; 8 Asian RPV genome sequences that


have been previously published24,26,27 or simply deposited in the public sequence databases were included for comparison. The final alignment included 70 genomes, which were analysed by


both maximum likelihood (ML) (Fig. 6) and Bayesian methods (Supplementary Fig. S2) in order to avoid any errors linked to the known weaknesses of either method50. In fact, the phylogenetic


trees produced by the two methods differ only in minor details of the relative placement of the set of very closely related Indian sequences. In both cases, the tree branches were very


strongly supported by the estimates of robustness, i.e. the results of the bootstrap (ML tree) or the posterior probability values (Bayesian tree). The phylogenetic trees produced by both


analyses are unrooted (i.e. no molecular clock is suggested), although for clarity in identifying the RPV isolate at each tip we have displayed them in the style of rooted trees. The root of


the RPV tree was determined by carrying out phylogenetic analyses on the same set of RPV genomes but with the inclusion of a MV genome sequence to act as an outgroup (Supplementary Fig. 


S3), and by assuming that the ancestral node closest to the MV sequence would be the root of the RPV tree. This placed the root on the section of the tree between the RPV/KabeteO clade and


that containing the goat-adapted vaccines, the two clades derived from the oldest isolates of RPV. As expected, there was a clear differentiation between isolates from countries in Asia and


the Middle East and those from countries in Africa (Fig. 6). Previous studies of the evolutionary relationships of different RPV isolates7,51 divided the viruses into three lineages, one


covering all isolates from Asia, the Middle and Near East, and two containing all the isolates from Africa (Africa 1 and Africa 2). These analyses were based on a relatively short (322


bases) stretch of the F CDS, and many of the sequences used in those analyses were from PCR products obtained from diagnostic samples; the associated viruses were not isolated and were


therefore not available to us for full genome sequencing. Africa 1 contained (among others) RPVs Egypt/84, Sudan/92 and the Kenyan isolates from 1988, 1989 and 1991, which are all still


closely linked when their full genome sequences are analysed. Africa 2 contained (among others) Tanzania/61/RBT1, Kenya/62/RGK1, Kenya/Kudu/93 and Nigeria/Sokoto/64; however, the full genome


sequence of these isolates shows that their most recent common ancestor is very distant from all of the other available sequences. It may not be reasonable to consider these viruses a


single clade or lineage, especially as Tanzania/61/RBT1 clearly came from a branch that separated off from the others before the split between “Africa 1” and the rest. Historically, RPV is


thought to have made a single incursion into sub-Saharan Africa in 188752, and spread rapidly over most of the continent. It was always a significant puzzle as to why there were two distinct


clades in Africa, or three if one includes the clade containing RPV/KabeteO and its derivatives, given that all African RPV isolates are thought to come from this single epizootic. The new


data suggest that, based on their evolutionary distance from each other, the African isolates fall into at least 5 clades (Tanzania/61/RBT1&2; Kenya/62/RGK1, Kenya/Kudu/93 and


Kenya/Tala/66; Sudan/Nyala/Reedbuck/85; Nigeria/Sokoto/64; the rest) which differ from each other by more than the distances separating any of the available Asian isolates. Alternatively,


they can be seen as a single large clade, with multiple branches that probably represent the isolation of the virus at very different times and in different geographical areas, a pattern


that is more in conformity with the available historical evidence. For comparison with previous publications, we generated a similar phylogenetic tree using all the available 322-base


sections of the RPV F gene (Supplementary Fig. S4), including those derived from diagnostic samples. It is clear from these data that, while there is a strong clade containing many mostly


East African samples, the Nigeria/58, Tanzania/61 and Nigeria/Sokoto/64 isolates each form a separate branch, as does the group of Kenyan isolates including Kenya/Tala/66 and Kenya/62/RGK1.


Because of the limited amount of data for each isolate, the support values for many of the branches are not as strong as those for the full genome sequences, but the overall pattern is very


similar. In the whole genome analyses, the clade containing the Kabete ‘O’ isolate and related viruses clustered with the other African isolates, although on a distinct branch, reflecting


that the parent of these viruses was originally isolated in Kenya in 1910. The Kabete ‘O’ challenge strain, from which Plowright derived his vaccine by repeated passage in BK cells53, was


originally maintained by cattle passage for use as a standard challenge virus54, and this virus was clearly shared with other countries, as the challenge virus in use in Egypt in the 1980s


was a variant of the Kenyan Kabete’O’ virus, and not a wild type virus from Egypt. An unexpected finding was that the entire group of lapinised viruses related to the original Nakamura III


vaccine55, including the Korea/Fusan-B which is thought to be the virus from which Nakamura derived his vaccine24, was clustered with the African viruses. This is not a reflection of how the


figure is drawn, but of the fact that the evolutionary distance from the Nakamura/Fusan clade to the African viruses, as estimated by both ML and Bayesian methods, is less than the distance


from the RPV/KabeteO clade to the rest of the African viruses. This finding was unexpected as there is no history suggesting a link between scientists working in Africa and the Japanese


scientists working at Fusan in Korea at the time (approximately 1934–8) when Nakamura was beginning his work on adapting RPV to rabbits56. Unfortunately, the origin of the virus used in


those studies was not given, nor was it referred to as anything other than “Laboratory Strain”56, although Nakamura himself refers to it in a later work as the “laboratory “O” strain”57,


reminiscent of the name (Kabete ‘O’) of the strain used in Kenya at that time. Further research may clarify the exact origin of the virus used in Korea for these studies. It is interesting


to note that the lapinised RPV held at the EAVRO in Kenya (and at some time transferred to AVRI in the UK for the preparation of anti-RPV sera for use in diagnosis) branches off the


lapinised virus line before the extra adaptation steps that lead to the current set of Japanese/Korean vaccines. This virus may thus represent something closer to the original Nakamura III


virus56 than the sample sequenced under that name, which had undergone many further passages in rabbits since the original vaccine was created. The goat-adapted viruses, originally developed


for use as a vaccine during the late 1920s in India58, form a distinct clade which, as is the case with the Kabete ‘O’/RBOK-vaccine viruses, probably reflects both the long time since their


isolation in the wild and the artificial way in which the virus was maintained. An interesting observation in this group of sequences is that the so-called Kabete-Adapted Goat (KAG) vaccine


virus, originally reported as having been developed from the Kabete ‘O’ wild type virus59, is clearly not related to RPV/KabeteO at all but is closely related to the Indian goat-adapted


viruses. The sequence of the RPV/KAG virus was the same in all the samples sequenced, including unopened vials prepared at EAVRO in Kenya. The attenuation of the Kabete ‘O’ virus over 250


passages in goats is on record60,61. These data, however, suggest that at some point the Edwards goat-adapted virus, sent to Kenya in 193662 was switched with the African virus being


passaged there in goats and, in the absence of any way of distinguishing strains at that time, this mistake was then propagated. The sequences of a large selection of Indian samples of RPV,


recorded as having been prepared as challenge virus from samples taken in different places at different times (See Supplementary Table S1), were essentially identical. The sample labelled


“Hill Bull” is a sample of one of the standard challenge virus used at Mukteswar in India for many decades; though it is unknown exactly which one is represented by RPV/India/HillBull. While


it was reasonable to expect the Ajmer, Hissar and Bombay isolates, all from around 1953–4, to be similar to each other, the Bangalore isolate of 1972 and those from Ranipet in 1973 and 1980


appear to be essentially identical to those viruses from 20–30 years previously. It is possible that RPV in India had become completely stable over a long period, having fully adapted to


the local hosts, and showing minimal sequence drift from 1950–80; however this does not accord with the continuous genetic drift seen in other RPV isolates from the same region and time,


such as the India/bison/89 and Sri Lanka/87 isolates. It is more likely that an established challenge strain was repeatedly re-isolated from cattle and relabelled as a new challenge strain.


Despite the intensive work on RPV carried out over many decades in India, primarily at Mukteswar, we have relatively few genome sequences from this part of the world. It will be useful if


any viruses still being held in India are similarly sequenced before being destroyed. In summary, we have sequenced, and now destroyed, almost all the RPV held at the Pirbright Institute.


The full genome sequences provide evidence supporting a single entry of the virus into sub-Saharan Africa and its expansion into multiple subclades. Further bioinformatic analyses may reveal


more detailed information about the growth and evolution of RPV. MATERIALS AND METHODS VIRUSES All the virus samples used in these studies were from the archive at the Pirbright Institute.


The complete list of RPV isolates sequenced, the methods of RNA extraction and purification, the method used to prepare the sequencing library and the fraction of RNA mapping to RPV in each


case is given in Supplementary Table S1, along with the sample history. RNA EXTRACTION Several methods were used to extract RNA from samples of cultured virus or from samples of tissue from


infected animals. Extraction with phenol-based reagents followed by ion-exchange spin-column purification was effective for cell culture supernatants containing cultured virus; for many


tissue samples, this method gave RNA that clearly contained an inhibitor of the RT or PCR step, as shown by an improved response in RT-qPCR after sample dilution. Extraction with the


Kingfisher automated magnetic bead system gave reproducibly cleaner results, although lower yield, and was adopted as the standard method by the end of the project. TRIzol LS was used to


dissolve cell culture samples, followed by RNA extraction using the Direct-zol RNA Miniprep kit or phase separation using the TRIzol protocol, when the RNA was either precipitated directly


using isopropanol, or extracted from the aqueous layer using the Zymo RNA Clean & Concentrator-5 kit. For extraction of RNA from freeze-dried tissue, the sample was resuspended directly


in 1 ml of TRIzol reagent and extracted using one of the methods above. Automated extraction of RNA was carried out using the LSI MagVet kit on a Kingfisher Flex Purification System.


PRELIMINARY SCREENING USING RT-QPCR The level of RPV-specific RNA in samples was estimated using a variation on the previously published RT-qPCR assay targeting the L gene18 in which the


probe contained a minor groove binder and non-fluorescent quencher to improve its effective melting temperature. Reactions were performed using 3 µl of sample RNA in a 20 µl reaction volume


on an AB7500 fast real-time PCR instrument: reverse transcription at 50 °C for 15 min, Taq activation at 95 °C for 20 s then 40 amplification cycles of 95 °C for 3 s and 60 °C for 30 s. NGS


LIBRARY PREPARATION Sequencing libraries were prepared using either transposon-based fragmentation of cDNA (Nextera XT DNA Library Prep kit, Illumina) or single primer isothermal


amplification (SPIA) (Trio RNA-Seq kit, NuGEN); in each case the reactions were carried out according to the manufacturer’s instructions. For library preparation using the Nextera kit, first


strand cDNA was generated from total RNA (0.4–4 µg depending on concentration) using random hexamer primers and SuperScript III reverse transcriptase according to the manufacturer’s


protocol. RNA was then digested at 37 °C for 20 min with 2U RNase H and double stranded cDNA synthesised using NEBNext Second Strand Synthesis enzyme mix and reaction buffer in a final


volume of 80 µl at 16 °C for 2.5 h. Double stranded cDNA was purified using the Illustra GFX DNA purification kit. The concentration of cDNA was assessed using the Qubit dsDNA HS assay kit


and adjusted to a final concentration of 0.2 ng/µl. Libraries were generated from 1 ng cDNA using the Nextera XT DNA Library Prep kit. For library preparation using SPIA, RNA was quantified


using the Qubit RNA HS quantification kit and 50 ng total RNA was used for library preparation with the Trio RNA-Seq kit. After amplification, enzymatic fragmentation and library


construction, host rRNA sequences were depleted as directed by the kit manufacturer (Trio RNA-Seq kit, NuGen). For all libraries, paired-end read sequencing was carried out using the


Illumina MiSeq platform and version 2 reagents. ANALYSIS OF NGS DATA All NGS datasets were analysed with a custom script in which the data was first quality trimmed using _Sickle_63 and then


mapped to the sequence of RPV/KabeteO (Accession numbers NC_006296/X98291) using _bowtie2_21 and _bwa-mem_22. In each case, duplicates were removed from the reads that mapped to the bait


sequence using the _SAMtools_ package64 and the consensus sequence was determined using the same package. The two consensus sequences were then compared, any disagreements resolved by manual


inspection of the read data, and the two merged into a single consensus. Regions at the genome ends or in the M-F GC-rich region that were not determined by the NGS data were filled by RACE


(genome ends) or a specific GC-rich PCR protocol. AMPLIFICATION AND SEQUENCING OF THE GC-RICH REGION The GC-rich region was amplified in two fragments using the primer pairs GC_F2/Frag3R


and Frag4F/GC_R2 (sequences of all named primers are given in Supplementary Table 2). cDNA template was prepared using Superscript III according to the manufacturer’s instructions. PCR was


carried out using KAPA HiFi reaction mix (Roche) with 10pmol each of forward and reverse primer and 2 µl of first strand cDNA template (50 µl final volume). The PCR cycling conditions were:


95 °C for 5 min, then 40 cycles of 98 °C for 20 s, 60 °C or 54 °C (for the GC_F2/Frag3R and Frag4F/R2 primer pairs, respectively) for 15 s, 72 °C for 1 min, and a final elongation step of 72


 °C for 1 min. In cases where the primers above did not result in amplification, isolate-specific primers with high melting temperatures were designed based on the next generation sequencing


data and PCR carried out as above but using a combined annealing/extension step of 68 °C for 45 s. In all cases, PCR products were purified using the Illustra GFX DNA purification kit,


according to manufacturer’s instructions. The purified PCR product was sequenced in both directions using the BigDye Terminator v3.1 reagents with the addition of dGTP to 1 µM. The cycle


sequencing conditions were: denaturation at 96 °C for 1 min followed by 30 cycles of 96 °C for 10 s, 50 °C for 5 s and 60 °C for 4 min. RAPID AMPLIFICATION OF SEQUENCE ENDS (RACE) To amplify


the 5′ end of the genome and/or antigenome in order to determine the terminal sequences of isolates, a variant of our previously published RACE protocol65 was used. RPV-specific cDNA was


prepared using SuperScript III and 2pmol each of RACE2, RACE3, RACE5 and RACE6 primers in a final volume of 20 µl. After incubation with 2U RNase H and 50U RNase 1 f (NEB) at 37 °C for 60 


min, single-stranded cDNA was purified using the Qiagen DyeEx 2.0 spin column purification kit. Purified cDNA (10 µl) was tailed in a 20 µl reaction containing 10U terminal deoxynucleotide


transferase (TdT), 0.1 mM dATP and 0.5X complete SuperScript III reaction buffer; the reaction was incubated at 37 °C for 5 min followed by inactivation of TdT at 70 °C for 10 min. The 5′


end of the genome was amplified in a 50 µl reaction containing 1X KOD hot start master mix (Merck Millipore), 15pmol each of Q1 and RACE4a primers, 0.5pmol QT primer and 2 µl tailed cDNA.


Cycling conditions were: 95 °C for 2 min 20 s, 33 °C for 15 s, 70 °C for 15 s, followed by 40 cycles of 95 °C for 20 s, 53 °C for 10 s and 70 °C for 15 s. For the 5′ end of the antigenome, a


hemi-nested PCR protocol was used in which the first stage was as for the 5′ end of the genome, except for the use of RACE1 primer instead of RACE4a. The product of this reaction was


dominated by amplicons derived from the 5′ end of the N gene mRNA transcript, so the 5′ end of the antigenome was then amplified in a second PCR using 15pmol each Q1 and M13RACE7c primers


and 0.5 µl 1st round PCR product as template (50 µl final volume). Cycling conditions were: 95 °C for 2 min 20 s, followed by 40 cycles of 95 °C for 20 s, 63 °C for 10 s and 70 °C for 15 s.


PCR products were purified as above and sequenced using RACE4a for the genome 5′ end and the standard M13 forward primer for the antigenome 5′ end. All Sanger sequencing was carried out on


an Applied Biosystems 3730 DNA Analyser. PHYLOGENETIC DATA ANALYSIS Where an isolate had been sequenced from more than one sample, only the consensus sequence was used for phylogenetic


analysis. Sequences of the same isolate after sequential passage were also not included. In addition to 62 sequences from this study, 8 published full length RPV genomes were included:


cattle-passaged RPV/Korea/Fusan-B (AB547189) and the lapinised derivative RPV/Lap/NakamuraIII (AB547190)24; three further lapinised/avianised derivatives of Nakamura III namely RPV/Lap/L72


(JN234008), RPV/Lap/LA77 (JN234009) and RPV/Lap/LA96 (JN234010)26; an additional lapinised/avianised strain, RPV/Lap/LATC06 (GU168576)25; the current lapinised/avianised RPV vaccine kept in


Japan, RPV/Lap/LA-AKO (LC057619)27; an unpublished sequence of a lapinised strain from the University of Tokyo, RPV/Lap/Lv (LC168749). We did not include our previously published sequences


for the RBOK vaccine strain and the Kabete ‘O’ wildtype virus8,9 as these viruses were included in the set of genome sequences determined in this study and those sequences would be more


reliable than those originally determined from cDNA libraries. The optimal partitioning of the sequences was determined using PartitionFinder248. Optimisation was restricted to the most


general models, the General Time Reversible (GTR) model with or without a gamma-distributed set of variable rates (+G) and a fraction of the sequence that was completely invariant (+I). The


GTR + G + I model produced, overall, the better fit as judged by the Bayesian Information Criterion (BIC) values. This model was used for all subsequent analyses. The ML tree was determined


using RAxML66. The program was run 3 times, each time using 20 random trees as starting points, and the best fit tree taken from these results. Bootstrap support values were obtained from


1000 bootstraps using RAxML’s rapid bootstrapping algorithm. The best fit tree was also determined by Bayesian methods using MrBayes67. The default priors were used throughout, and the


search was run for 1,000,000 generations. Tree figures were prepared with FigTree v1.4.4. DATA AVAILABILITY All the genome sequences resulting from this study are available in the public


sequence databases; the relevant accession numbers are listed in Supplementary Table S1. REFERENCES * World Organisation for Animal Health (OIE). Final report of the 79th OIE General Session


p337-342 Resolution No. 18, Declaration of global eradication of rinderpest and implementation of follow-up measures to maintain world freedom from rinderpest.


http://www.oie.int/fileadmin/Home/eng/Media_Center/docs/pdf/RESO_18_EN.pdf (2011). * Roeder, P. L. & Rich, K. The global effort to eradicate rinderpest. IFPRI discussion paper 0923.


(International Food Policy Research Institute, 2009). * Beauvais, W. _et al_. Modelling the expected rate of laboratory biosafety breakdowns involving rinderpest virus in the


post-eradication era. _Prev. Vet. Med._ 112, 248–256 (2013). Article  CAS  PubMed  Google Scholar  * Fournie, G. _et al_. The risk of rinderpest re-introduction in post-eradication era.


_Prev. Vet. Med._ 113, 175–184 (2014). Article  PubMed  Google Scholar  * Anderson, J., McKay, J. A. & Butcher, R. N. In Seromonitoring of rinderpest throughout Africa: phase one.


Proceedings of the final research coordination meeting of the IAEA rinderpest control projects, Cote d’Ivoire 19–23 November 1990 IAEA-TECDOC-623 (International Atomic Energy Agency, Vienna,


1990). * Forsyth, M. A. & Barrett, T. Evaluation of polymerase chain reaction for the detection and characterisation of rinderpest and peste des petits ruminants viruses for


epidmiological studies. _Virus Res._ 39, 151–163 (1995). Article  CAS  PubMed  Google Scholar  * Chamberlain, R. W. _et al_. Evidence for different lineages of rinderpest virus reflecting


their geographic isolation. _J. Gen. Virol._ 74, 2775–2780 (1993). Article  CAS  PubMed  Google Scholar  * Baron, M. D. & Barrett, T. The sequence of the N and L genes of rinderpest


virus and the 5’ and 3’ extra-genic sequences: the completion of the genome sequence of the virus. _Vet. Microbiol._ 44, 175–186 (1995). Article  CAS  PubMed  Google Scholar  * Baron, M. D.,


Kamata, Y., Barras, V., Goatley, L. & Barrett, T. The genome sequence of the virulent Kabete ‘O’ strain of rinderpest virus: comparison with the derived vaccine. _J. Gen. Virol._ 77,


3041–3046 (1996). Article  CAS  PubMed  Google Scholar  * Baron, M. D. & Barrett, T. Rescue of rinderpest virus from cloned cDNA. _J. Virol_ 71, 1265–1271 (1997). Article  CAS  PubMed 


PubMed Central  Google Scholar  * Drexler, J. F. _et al_. Bats host major mammalian paramyxoviruses. _Nat. Commun._ 3, 796 (2012). Article  ADS  PubMed  CAS  Google Scholar  * Baron, M. D.,


Foster-Cuevas, M., Baron, J. & Barrett, T. Expression in cattle of epitopes of a heterologous virus using a recombinant rinderpest virus. _J. Gen. Virol_ 80, 2031–2039 (1999). Article 


CAS  PubMed  Google Scholar  * Baron, M. D. & Barrett, T. Rinderpest viruses lacking the C and V proteins show specific defects in growth and transcription of viral RNAs. _J. Virol._ 74,


2603–2611 (2000). Article  CAS  PubMed  PubMed Central  Google Scholar  * Das, S. C., Baron, M. D. & Barrett, T. Recovery and characterization of a chimeric rinderpest virus with the


glycoproteins of peste-des-petits-ruminants virus: homologous F and H proteins are required for virus viability. _J. Virol._ 74, 9039–9047 (2000). Article  CAS  PubMed  PubMed Central 


Google Scholar  * Walsh, E. P. _et al_. Recombinant rinderpest vaccines expressing membrane-anchored proteins as genetic markers: evidence of exclusion of marker protein from the virus


envelope. _J. Virol._ 74, 10165–10175 (2000). Article  CAS  PubMed  PubMed Central  Google Scholar  * Baron, M. D., Banyard, A. C., Parida, S. & Barrett, T. The Plowright vaccine strain


of rinderpest virus has attenuating mutations in most genes. _J. Gen. Virol_ 86, 1093–1101 (2005). Article  CAS  PubMed  Google Scholar  * Muniraju, M. _et al_. Rescue of a vaccine strain of


peste des petits ruminants virus: _In vivo_ evaluation and comparison with standard vaccine. _Vaccine_ 33, 465–471 (2015). Article  CAS  PubMed  PubMed Central  Google Scholar  * Carrillo,


C. _et al_. Specific detection of rinderpest virus by real-time reverse transcription-PCR in preclinical and clinical samples from experimentally infected cattle. _J. Clin. Microbiol._ 48,


4094–4101 (2010). Article  CAS  PubMed  PubMed Central  Google Scholar  * Plowright, W., Cruickshank, J. G. & Waterson, A. P. The morphology of rinderpest virus. _Virology_ 17, 118–122


(1962). Article  CAS  PubMed  Google Scholar  * Underwood, B. & Brown, F. Physico-chemical characterisation of rinderpest virus. _Med. Microbiol. Immunol. (Berl.)_ 160, 125–132 (1974).


Article  CAS  Google Scholar  * Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. _Nat. Methods_ 9, 357–359 (2012). Article  CAS  PubMed  PubMed Central  Google


Scholar  * Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 1303, 3997. Preprint at https://arxiv.org/abs/1303.3997 (2013). * Quail, M. A. _et al_.


Optimal enzymes for amplifying sequencing libraries. _Nat. Methods_ 9, 10–11 (2012). Article  CAS  Google Scholar  * Fukai, K., Morioka, K., Sakamoto, K. & Yoshida, K. Characterization


of the complete genomic sequence of the rinderpest virus Fusan strain cattle type, which is the most classical isolate in Asia and comparison with its lapinized strain. _Virus Genes_ 43,


249–253 (2011). Article  CAS  PubMed  Google Scholar  * Yeh, J. Y. _et al_. Genetic characterization of the Korean LATC06 rinderpest vaccine strain. _Virus Genes_ 42, 71–75 (2011). Article 


CAS  PubMed  Google Scholar  * Jeoung, H. Y. _et al_. Complete genome analysis of three live attenuated rinderpest virus vaccine strains derived through serial passages in different culture


systems. _J. Virol._ 86, 13115–13116 (2012). Article  CAS  PubMed  PubMed Central  Google Scholar  * Takamatsu, H., Terui, K. & Kokuho, T. Complete genome sequence of Japanese vaccine


strain LA-AKO of rinderpest virus. Genome Announc 3, (2015). * Bao, J. _et al_. Complete genome sequence of a novel variant strain of peste des petits ruminants virus, China/XJYL/2013.


Genome Announc 2, (2014). * Bankamp, B. _et al_. Wild-type measles viruses with non-standard genome lengths. _PLoS One_ 9, e95470 (2014). Article  ADS  PubMed  PubMed Central  CAS  Google


Scholar  * Ivancic-Jelecki, J., Slovic, A., Santak, M., Tesovic, G. & Forcic, D. Common position of indels that cause deviations from canonical genome organization in different measles


virus strains. _Virol. J._ 13, 134 (2016). Article  PubMed  PubMed Central  CAS  Google Scholar  * Gil, H. _et al_. Measles virus genotype D4 strains with non-standard length M-F non-coding


region circulated during the major outbreaks of 2011-2012 in Spain. _PLoS One_ 13, e0199975 (2018). Article  PubMed  PubMed Central  CAS  Google Scholar  * Taylor, W. P. Epidemiology and


control of rinderpest. _Revue Scientifique et Technique Office International des Epizooties_ 5, 407–410 (1986). Article  Google Scholar  * von Heijne, G. A new method for predicting signal


sequence cleavage sites. _Nucl. Acids Res_ 14, 4683–4690 (1986). Article  Google Scholar  * Evans, S. A., Baron, M. D., Chamberlain, R. W., Goatley, L. & Barrett, T. Nucleotide sequence


comparisons of the fusion protein gene from virulent and attenuated strains of rinderpest virus. _J. Gen. Virol._ 75, 3611–3617 (1994). Article  CAS  PubMed  Google Scholar  * Ketteler, R.


On programmed ribosomal frameshifting: the alternative proteomes. _Front Genet_ 3, 242 (2012). Article  PubMed  PubMed Central  Google Scholar  * Peabody, D. S. Translation initiation at


non-AUG triplets in mammalian cells. _J. Biol. Chem._ 264, 5031–5035 (1989). CAS  PubMed  Google Scholar  * Cathomen, T., Buchholz, C. J., Spielhofer, P. & Cattaneo, R. Preferential


initiation at the second AUG of the measles virus F mRNA: A role for the long untranslated region. _Virology_ 214, 628–632 (1995). Article  CAS  PubMed  Google Scholar  * Takeda, M. _et al_.


Long untranslated regions of the measles virus M and F genes control virus replication and cytopathogenicity. _J. Virol._ 79, 14346–14354 (2005). Article  CAS  PubMed  PubMed Central 


Google Scholar  * Anderson, D. E. & von Messling, V. Region between the canine distemper virus M and F genes modulates virulence by controlling fusion protein expression. _J. Virol._ 82,


10510–10518 (2008). Article  CAS  PubMed  PubMed Central  Google Scholar  * Banyard, A. C., Baron, M. D. & Barrett, T. A role for virus promoters in determining the pathogenesis of


rinderpest virus in cattle. _J. Gen. Virol_ 86, 1083–1092 (2005). Article  CAS  PubMed  Google Scholar  * Eloiflin, R. J. _et al_. Evolution of attenuation and risk of reversal in peste des


petits ruminants vaccine strain Nigeria 75/1. Viruses 11, (2019). * Plowright, W. & Ferris, R. D. Studies with rinderpest virus in tissue culture. II. Pathogenicity for cattle of


culture-passaged virus. _J. Comp. Pathol_ 69, 173–184 (1959). Article  CAS  PubMed  Google Scholar  * Johnston, M. D. The characteristics required for a Sendai virus preparation to induce


high levels of interferon in human lymphoblastoid cells. _J. Gen. Virol._ 56, 175–184 (1981). Article  CAS  PubMed  Google Scholar  * von Magnus, P. Propagation of the PR8 strain of


influenza A virus in chick embryos. III. Properties of the incomplete virus produced in serial passages of undiluted virus. _Acta Pathol. Microbiol. Scand_ 29, 157–181 (1951). Article 


Google Scholar  * Saira, K. _et al_. Sequence analysis of _in vivo_ defective interfering-like RNA of influenza A H1N1 pandemic virus. _J. Virol._ 87, 8064–8074 (2013). Article  CAS  PubMed


  PubMed Central  Google Scholar  * Li, D. _et al_. Defective interfering viral particles in acute dengue infections. _PLoS One_ 6, e19447 (2011). Article  ADS  CAS  PubMed  PubMed Central 


Google Scholar  * Pesko, K. N. _et al_. Internally deleted WNV genomes isolated from exotic birds in New Mexico: function in cells, mosquitoes, and mice. _Virology_ 427, 10–17 (2012).


Article  CAS  PubMed  Google Scholar  * Lanfear, R., Frandsen, P. B., Wright, A. M., Senfeld, T. & Calcott, B. PartitionFinder 2: new methods for selecting partitioned models of


evolution for molecular and morphological phylogenetic analyses. _Mol. Biol. Evol._ 34, 772–773 (2017). CAS  PubMed  Google Scholar  * Bofkin, L. & Goldman, N. Variation in evolutionary


processes at different codon positions. _Mol. Biol. Evol._ 24, 513–521 (2007). Article  CAS  PubMed  Google Scholar  * Yang, Z. & Rannala, B. Molecular phylogenetics: principles and


practice. _Nat. Rev. Genet._ 13, 303–314 (2012). Article  ADS  CAS  PubMed  Google Scholar  * Wamwayi, H. M., Fleming, M. & Barrett, T. Characterisation of African isolates of rinderpest


virus. _Vet. Microbiol._ 44, 151–163 (1995). Article  CAS  PubMed  Google Scholar  * Spinage, C. A. Cattle Plague: A History. (Kluwer Academic/Plenum, 2003). * Plowright, W. & Ferris,


R. D. Studies with rinderpest virus in tissue culture. The use of attenuated culture virus as a vaccine for cattle. _Res. Vet. Sci_ 3, 172–182 (1962). Article  Google Scholar  * MacOwen, K.


D. S. Department of Veterinary Services Annual Report 1955. 29 (The Government Printer, Colony and Protectorate of Kenya, Nairobi, 1956). * Nakamura, J., Kishi, S., Kiuchi, J. &


Reisinger, R. An investigation of antibody response in cattle vaccinated with the rabbit-passaged LA rinderpest virus in Korea. _Am. J. Vet. Res._ 16, 71–75 (1955). CAS  PubMed  Google


Scholar  * Nakamura, J., Wagatuma, S. & Fukusho, K. On the experimental infection with rinderpest virus in the rabbit. 1. Some fundamental experiments. _Journal of the Japanese Society


of Veterinary Science_ 17, 25–30 (1938). Google Scholar  * Nakamura, J. Two rinderpest live virus vaccines, “Lapinized” and “Lapinized-Avianized”. _Japan Agricultural Research Quarterly_ 1,


11–17 (1966). CAS  Google Scholar  * Edwards, J. T. In FEATM Seventh Congress Vol. III 699–706 (Thacker’s Press and Directories, Ltd., India, 1929). * Daubney, R. Rinderpest: a résumé of


recent progress in East Africa. _Journal of Comparative Pathology and Therapeutics_ 50, 405–409 (1937). Article  Google Scholar  * Daubney, R. In Rinderpest vaccines: their production and


use in the field (ed K. V. L. Kesteven) 6-18 (Food and Agriculture Organisation of the United Nations, 1949). * Daubney, R. In Proceedings of the 4th international congresses on tropical


medicine and malaria. 1358–1365 (U.S. Government Printing Office, 1948). * Taylor, W. P., Roeder, P. L. & Rweyemamu, M. M. In Rinderpest and peste des petits ruminants (eds T. Barrett,


P-P. Pastoret, & W.P. Taylor) Ch. 11, (Academic Press, 2006). * Sickle: A sliding-window, adaptive, quality-based trimming tool for FastQ files v. 1.33


(https://github.com/najoshi/sickle, 2011). * Li, H. _et al_. The sequence alignment/map format and SAMtools. _Bioinformatics_ 25, 2078–2079 (2009). Article  PubMed  PubMed Central  CAS 


Google Scholar  * Baron, M. D. & Barrett, T. Sequencing and analysis of the nucleocapsid (N) and polymerase (L) genes and the terminal extragenic domains of the vaccine strain of


rinderpest virus. _J. Gen. Virol._ 76, 593–602 (1995). Article  CAS  PubMed  Google Scholar  * Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large


phylogenies. _Bioinformatics_ 30, 1312–1313 (2014). Article  CAS  PubMed  PubMed Central  Google Scholar  * Ronquist, F. _et al_. MrBayes 3.2: efficient Bayesian phylogenetic inference and


model choice across a large model space. _Syst. Biol_ 61, 539–542 (2012). Article  PubMed  PubMed Central  Google Scholar  Download references ACKNOWLEDGEMENTS This work was carried out with


the financial support of the WMD Threat Reduction Program of Global Affairs Canada and the Defense Threat Reduction Agency of the United States Department of Defense. PRN and CB are


supported by the UK Department for Environment, Food and Rural Affairs and the UK Biotechnology and Biological Sciences Research Council (BBSRC). AUTHOR INFORMATION Author notes * Paolo


Ribeca Present address: Biomathematics and Statistics Scotland, JCMB, The King’s Buildings, Peter Guthrie Tait Road, Edinburgh, EH9 3FD, Scotland, UK AUTHORS AND AFFILIATIONS * The Pirbright


Institute, Ash Road, Pirbright, Surrey, GU24 0NF, UK Simon King, Paulina Rajko-Nenow, Honorata M. Ropiak, Paolo Ribeca, Carrie Batten & Michael D. Baron Authors * Simon King View author


publications You can also search for this author inPubMed Google Scholar * Paulina Rajko-Nenow View author publications You can also search for this author inPubMed Google Scholar *


Honorata M. Ropiak View author publications You can also search for this author inPubMed Google Scholar * Paolo Ribeca View author publications You can also search for this author inPubMed 


Google Scholar * Carrie Batten View author publications You can also search for this author inPubMed Google Scholar * Michael D. Baron View author publications You can also search for this


author inPubMed Google Scholar CONTRIBUTIONS S.K. carried out the major part of the RNA extraction, library construction and sequencing, both N.G.S. and Sanger; P.R.N. developed the initial


techniques and carried out the rest of the library construction and N.G.S. sequencing; H.M.R. assisted with all stages of the laboratory sequencing work; P.R. advised on the bioinformatics


analysis and data archiving; C.B. had overall responsibility for the project; M.D.B. initiated the project, developed the RACE protocol, set up the bioinformatics pipeline, carried out all


bioinformatics and phylogenetics analyses, researched the virus history, wrote the manuscript and prepared the figures. All authors reviewed the manuscript. CORRESPONDING AUTHOR


Correspondence to Michael D. Baron. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION PUBLISHER’S NOTE Springer Nature remains


neutral with regard to jurisdictional claims in published maps and institutional affiliations. SUPPLEMENTARY INFORMATION SUPPLEMENTARY TABLE. SUPPLEMENTARY INFORMATION. RIGHTS AND


PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any


medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The


images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not


included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly


from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE King, S.,


Rajko-Nenow, P., Ropiak, H.M. _et al._ Full genome sequencing of archived wild type and vaccine rinderpest virus isolates prior to their destruction. _Sci Rep_ 10, 6563 (2020).


https://doi.org/10.1038/s41598-020-63707-z Download citation * Received: 02 December 2019 * Accepted: 11 March 2020 * Published: 16 April 2020 * DOI:


https://doi.org/10.1038/s41598-020-63707-z SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not


currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative


Trending News

Full genome sequencing of archived wild type and vaccine rinderpest virus isolates prior to their destruction

ABSTRACT When rinderpest virus (RPV) was declared eradicated in 2011, the only remaining samples of this once much-feare...

Latests News

Full genome sequencing of archived wild type and vaccine rinderpest virus isolates prior to their destruction

ABSTRACT When rinderpest virus (RPV) was declared eradicated in 2011, the only remaining samples of this once much-feare...

Top