Deep Sequencing of Foot-and-Mouth Disease Virus Reveals RNA Sequences Involved in Genome Packaging

Nonenveloped viruses protect their genomes by packaging them into an outer shell or capsid of virus-encoded proteins. Packaging and capsid assembly in RNA viruses can involve interactions between capsid proteins and secondary structures in the viral genome, as exemplified by the RNA bacteriophage MS2 and as proposed for other RNA viruses of plants, animals, and human. In the picornavirus family of nonenveloped RNA viruses, the requirements for genome packaging remain poorly understood. Here, we show a novel and simple approach to identify predicted RNA secondary structures involved in genome packaging in the picornavirus foot-and-mouth disease virus (FMDV). By interrogating deep sequencing data generated from both packaged and unpackaged populations of RNA, we have determined multiple regions of the genome with constrained variation in the packaged population. Predicted secondary structures of these regions revealed stem-loops with conservation of structure and a common motif at the loop. Disruption of these features resulted in attenuation of virus growth in cell culture due to a reduction in assembly of mature virions. This study provides evidence for the involvement of predicted RNA structures in picornavirus packaging and offers a readily transferable methodology for identifying packaging requirements in many other viruses.IMPORTANCE In order to transmit their genetic material to a new host, nonenveloped viruses must protect their genomes by packaging them into an outer shell or capsid of virus-encoded proteins. For many nonenveloped RNA viruses the requirements for this critical part of the viral life cycle remains poorly understood. We have identified RNA sequences involved in genome packaging of the picornavirus foot-and-mouth disease virus. This virus causes an economically devastating disease of livestock affecting both the developed and developing world. The experimental methods developed to carry out this work are novel, simple, and transferable to the study of packaging signals in other RNA viruses. Improved understanding of RNA packaging may lead to novel vaccine approaches or targets for antiviral drugs with broad-spectrum activity.

A universal protocol to generate consensus level genome sequences for foot-and-mouth disease virus and other positive-sense polyadenylated RNA viruses using the Illumina MiSeq

BACKGROUND

Next-Generation Sequencing (NGS) is revolutionizing molecular epidemiology by providing new approaches to undertake whole genome sequencing (WGS) in diagnostic settings for a variety of human and veterinary pathogens. Previous sequencing protocols have been subject to biases such as those encountered during PCR amplification and cell culture, or are restricted by the need for large quantities of starting material. We describe here a simple and robust methodology for the generation of whole genome sequences on the Illumina MiSeq. This protocol is specific for foot-and-mouth disease virus (FMDV) or other polyadenylated RNA viruses and circumvents both the use of PCR and the requirement for large amounts of initial template.

RESULTS

The protocol was successfully validated using five FMDV positive clinical samples from the 2001 epidemic in the United Kingdom, as well as a panel of representative viruses from all seven serotypes. In addition, this protocol was successfully used to recover 94% of an FMDV genome that had previously been identified as cell culture negative. Genome sequences from three other non-FMDV polyadenylated RNA viruses (EMCV, ERAV, VESV) were also obtained with minor protocol amendments. We calculated that a minimum coverage depth of 22 reads was required to produce an accurate consensus sequence for FMDV O. This was achieved in 5 FMDV/O/UKG isolates and the type O FMDV from the serotype panel with the exception of the 5' genomic termini and area immediately flanking the poly(C) region.

CONCLUSIONS

We have developed a universal WGS method for FMDV and other polyadenylated RNA viruses. This method works successfully from a limited quantity of starting material and eliminates the requirement for genome-specific PCR amplification. This protocol has the potential to generate consensus-level sequences within a routine high-throughput diagnostic environment.

Coronavirus NSP6 restricts autophagosome expansion

Autophagy is a cellular response to starvation that generates autophagosomes to carry long-lived proteins and cellular organelles to lysosomes for degradation. Activation of autophagy by viruses can provide an innate defense against infection, and for (+) strand RNA viruses autophagosomes can facilitate assembly of replicase proteins. We demonstrated that nonstructural protein (NSP) 6 of the avian coronavirus, infectious bronchitis virus (IBV), generates autophagosomes from the ER. A statistical analysis of MAP1LC3B puncta showed that NSP6 induced greater numbers of autophagosomes per cell compared with starvation, but the autophagosomes induced by NSP6 had smaller diameters compared with starvation controls. Small diameter autophagosomes were also induced by infection of cells with IBV, and by NSP6 proteins of MHV and SARS and NSP5, NSP6, and NSP7 of arterivirus PRRSV. Analysis of WIPI2 puncta induced by NSP6 suggests that NSP6 limits autophagosome diameter at the point of omegasome formation. IBV NSP6 also limited autophagosome and omegasome expansion in response to starvation and Torin1 and could therefore limit the size of autophagosomes induced following inhibition of MTOR signaling, as well as those induced independently by the NSP6 protein itself. MAP1LC3B-puncta induced by NSP6 contained SQSTM1, which suggests they can incorporate autophagy cargos. However, NSP6 inhibited the autophagosome/lysosome expansion normally seen following starvation. Taken together the results show that coronavirus NSP6 proteins limit autophagosome expansion, whether they are induced directly by the NSP6 protein, or indirectly by starvation or chemical inhibition of MTOR signaling. This may favor coronavirus infection by compromising the ability of autophagosomes to deliver viral components to lysosomes for degradation.

Visualizing the autophagy pathway in avian cells and its application to studying infectious bronchitis virus

Autophagy is a highly conserved cellular response to starvation that leads to the degradation of organelles and long-lived proteins in lysosomes and is important for cellular homeostasis, tissue development and as a defense against aggregated proteins, damaged organelles and infectious agents. Although autophagy has been studied in many animal species, reagents to study autophagy in avian systems are lacking. Microtubule-associated protein 1 light chain 3 (MAP1LC3/LC3) is an important marker for autophagy and is used to follow autophagosome formation. Here we report the cloning of avian LC3 paralogs A, B and C from the domestic chicken, Gallus gallus domesticus, and the production of replication-deficient, recombinant adenovirus vectors expressing these avian LC3s tagged with EGFP and FLAG-mCherry. An additional recombinant adenovirus expressing EGFP-tagged LC3B containing a G120A mutation was also generated. These vectors can be used as tools to visualize autophagosome formation and fusion with endosomes/lysosomes in avian cells and provide a valuable resource for studying autophagy in avian cells. We have used them to study autophagy during replication of infectious bronchitis virus (IBV). IBV induced autophagic signaling in mammalian Vero cells but not primary avian chick kidney cells or the avian DF1 cell line. Furthermore, induction or inhibition of autophagy did not affect IBV replication, suggesting that classical autophagy may not be important for virus replication. However, expression of IBV nonstructural protein 6 alone did induce autophagic signaling in avian cells, as seen previously in mammalian cells. This may suggest that IBV can inhibit or control autophagy in avian cells, although IBV did not appear to inhibit autophagy induced by starvation or rapamycin treatment.

Coronavirus nsp6 proteins generate autophagosomes from the endoplasmic reticulum via an omegasome intermediate

Autophagy is a cellular response to starvation which generates autophagosomes to carry cellular organelles and long-lived proteins to lysosomes for degradation. Degradation through autophagy can provide an innate defence against virus infection, or conversely autophagosomes can promote infection by facilitating assembly of replicase proteins. We demonstrate that the avian coronavirus, Infectious Bronchitis Virus (IBV) activates autophagy. A screen of individual IBV non-structural proteins (nsps) showed that autophagy was activated by IBV nsp6. This property was shared with nsp6 of mammalian coronaviruses Mouse Hepatitis Virus, and Severe Acute Respiratory Syndrome Virus, and the equivalent nsp5-7 of the arterivirus Porcine Reproductive and Respiratory Syndrome Virus. These multiple-spanning transmembrane proteins located to the endoplasmic reticulum (ER) where they generated Atg5 and LC3II-positive vesicles, and vesicle formation was dependent on Atg5 and class III PI3 kinase. The vesicles recruited double FYVE-domain containing protein (DFCP) indicating localised concentration of phosphatidylinositol 3 phosphate, and therefore shared many features with omegasomes formed from the ER in response to starvation. Omegasomes induced by viral nsp6 matured into autophagosomes that delivered LC3 to lysosomes and therefore recruited and recycled the proteins needed for autophagosome nucleation, expansion, cellular trafficking and delivery of cargo to lysosomes. The coronavirus nsp6 proteins activated omegasome and autophagosome formation independently of starvation, but activation did not involve direct inhibition of mTOR signalling, activation of sirtuin1 or induction of ER stress.

Analysis of Foot-and-mouth disease virus nucleotide sequence variation within naturally infected epithelium

An important epidemiological tool in the control of epidemics of Foot-and-mouth disease (FMD) is genetic tracing using complete virus genome sequence data. However to interpret these genetic data, it is important to quantify underlying variation present in FMDV populations from individual tissue samples. Cloned complete capsid sequences from two virus populations from epithelium from a cow (n=26) and from a sheep (n=15) infected during the UK 2001 outbreaks were generated. Genetic diversity of the two virus populations differed significantly, with sequences representing virus from the cow having a mutation frequency of 2.79 x 10(-4) mutations per nucleotide sequenced (mpns) and those from the sheep having 3.94 x 10(-4) mpns (chi(2)=8.24, P=0.004). The dN/dS ratio of sequences from the cow was higher (1.228) than that from the sheep (0.187) although not significantly so. The sequences from the cow epithelium exhibited significantly higher than expected number of changes within neutralising antigenic sites (P=0.0007). The performance of two different reverse transcriptase enzymes was found not to differ with respect to the frequency (P=0.559, chi(2)=0.341) or dN/dS ratio (P=0.863, chi(2)=0.03) of the mutations observed. These data provide insight into the population diversity that exists within a single lesion and help toward understanding the mechanisms that underpin sequence evolution of FMDV.

Integrating genetic and epidemiological data to determine transmission pathways of foot-and-mouth disease virus

Estimating detailed transmission trees that reflect the relationships between infected individuals or populations during a disease outbreak often provides valuable insights into both the nature of disease transmission and the overall dynamics of the underlying epidemiological process. These trees may be based on epidemiological data that relate to the timing of infection and infectiousness, or genetic data that show the genetic relatedness of pathogens isolated from infected individuals. Genetic data are becoming increasingly important in the estimation of transmission trees of viral pathogens due to their inherently high mutation rate. Here, we propose a maximum-likelihood approach that allows epidemiological and genetic data to be combined within the same analysis to infer probable transmission trees. We apply this approach to data from 20 farms infected during the 2001 UK foot-and-mouth disease outbreak, using complete viral genome sequences from each infected farm and information on when farms were first estimated to have developed clinical disease and when livestock on these farms were culled. Incorporating known infection links due to animal movement prior to imposition of the national movement ban results in the reduction of the number of trees from 41472 that are consistent with the genetic data to 1728, of which just 4 represent more than 95% of the total likelihood calculated using a model that accounts for the epidemiological data. These trees differ in several ways from those constructed prior to the availability of genetic data.

Transmission pathways of foot-and-mouth disease virus in the United Kingdom in 2007

Foot-and-mouth disease (FMD) virus causes an acute vesicular disease of domesticated and wild ruminants and pigs. Identifying sources of FMD outbreaks is often confounded by incomplete epidemiological evidence and the numerous routes by which virus can spread (movements of infected animals or their products, contaminated persons, objects, and aerosols). Here, we show that the outbreaks of FMD in the United Kingdom in August 2007 were caused by a derivative of FMDV O₁ BFS 1860, a virus strain handled at two FMD laboratories located on a single site at Pirbright in Surrey. Genetic analysis of complete viral genomes generated in real-time reveals a probable chain of transmission events, predicting undisclosed infected premises, and connecting the second cluster of outbreaks in September to those in August. Complete genome sequence analysis of FMD viruses conducted in real-time have identified the initial and intermediate sources of these outbreaks and demonstrate the value of such techniques in providing information useful to contemporary disease control programmes.

Foot-and-mouth disease (FMD) outbreaks in the United Kingdom during August and September 2007 have caused severe disruption to the farming sector and cost hundreds of millions of pounds. Investigating and determining the source of these outbreaks is imperative for their effective management and future prevention. Foot-and-mouth disease virus (FMDV) has a high mutation rate, resulting in rapid evolution. We show how complete genome sequences (acquired within 24–48 h of sample receipt) can be used to track FMDV movement from farm to farm in real time. This helped to determine the most likely source of the outbreak, assisted ongoing epidemiological investigations as to whether these field cases were linked to single or multiple releases from the source, and predicted the existence of undetected intermediate infected premises.

Molecular epidemiology of the foot-and-mouth disease virus outbreak in the United Kingdom in 2001

The objective of this study was to quantify the extent to which the genetic diversity of foot-and-mouth disease virus (FMDV) arising over the course of infection of an individual animal becomes fixed, is transmitted to other animals, and thereby accumulates over the course of an outbreak. Complete consensus sequences of 23 genomes (each of 8,200 nucleotides) of FMDV were recovered directly from epithelium tissue acquired from 21 farms infected over a nearly 7-month period during the 2001 FMDV outbreak in the United Kingdom. An analysis of these consensus sequences revealed very few apparently ambiguous sites but clear evidence of 197 nucleotide substitutions at 191 different sites. We estimated the rate of nucleotide substitution to be 2.26 x 10(-5) per site per day (95% confidence interval [CI], 1.75 x 10(-5) to 2.80 x 10(-5)) and nucleotide substitutions to accrue in the consensus sequence at an average rate of 1.5 substitutions per farm infection. This is a sufficiently high rate showing that detailed histories of the transmission pathways can be reliably reconstructed. Coalescent methods indicated that the date at which FMDV first infected livestock in the United Kingdom was 7 February 2001 (95% CI, 20 January to 19 February 2001), which was identical to estimates obtained on the basis of purely clinical evidence. Nucleotide changes appeared to have occurred evenly across the genome, and within the open reading frame, the ratio of nonsynonymous-to-synonymous change was 0.09. The ability to recover particular transmission pathways of acutely acting RNA pathogens from genetic data will help resolve uncertainties about how virus is spread and could help in the control of future epidemics.