Beyond the consensus: dissecting within-host viral population diversity of foot-and-mouth disease virus by using next-generation genome sequencing

The diverse sequences of viral populations within individual hosts are the starting material for selection and subsequent evolution of RNA viruses such as foot-and-mouth disease virus (FMDV). Using next-generation sequencing (NGS) performed on a Genome Analyzer platform (Illumina), this study compared the viral populations within two bovine epithelial samples (foot lesions) from a single animal with the inoculum used to initiate experimental infection. Genomic sequences were determined in duplicate sequencing runs, and the consensus sequence of the inoculum determined by NGS was identical to that previously determined using the Sanger method. However, NGS revealed the fine polymorphic substructure of the viral population, from nucleotide variants present at just below 50% frequency to those present at fractions of 1%. Some of the higher-frequency polymorphisms identified encoded changes within codons associated with heparan sulfate binding and were present in both foot lesions, revealing intermediate stages in the evolution of a tissue culture-adapted virus replicating within a mammalian host. We identified 2,622, 1,434, and 1,703 polymorphisms in the inoculum and in the two foot lesions, respectively: most of the substitutions occurred in only a small fraction of the population and represented the progeny from recent cellular replication prior to onset of any selective pressures. We estimated the upper limit for the genome-wide mutation rate of the virus within a cell to be 7.8 × 10(-4) per nucleotide. The greater depth of detection achieved by NGS demonstrates that this method is a powerful and valuable tool for the dissection of FMDV populations within hosts.     

Distinguishing low frequency mutations from RT-PCR and sequence errors in viral deep sequencing data

Background

RNA viruses have high mutation rates and exist within their hosts as large, complex and heterogeneous populations, comprising a spectrum of related but non-identical genome sequences. Next generation sequencing is revolutionising the study of viral populations by enabling the ultra deep sequencing of their genomes, and the subsequent identification of the full spectrum of variants within the population. Identification of low frequency variants is important for our understanding of mutational dynamics, disease progression, immune pressure, and for the detection of drug resistant or pathogenic mutations. However, the current challenge is to accurately model the errors in the sequence data and distinguish real viral variants, particularly those that exist at low frequency, from errors introduced during sequencing and sample processing, which can both be substantial.

Results

We have created a novel set of laboratory control samples that are derived from a plasmid containing a full-length viral genome with extremely limited diversity in the starting population. One sample was sequenced without PCR amplification whilst the other samples were subjected to increasing amounts of RT and PCR amplification prior to ultra-deep sequencing. This enabled the level of error introduced by the RT and PCR processes to be assessed and minimum frequency thresholds to be set for true viral variant identification. We developed a genome-scale computational model of the sample processing and NGS calling process to gain a detailed understanding of the errors at each step, which predicted that RT and PCR errors are more likely to occur at some genomic sites than others. The model can also be used to investigate whether the number of observed mutations at a given site of interest is greater than would be expected from processing errors alone in any NGS data set. After providing basic sample processing information and the site’s coverage and quality scores, the model utilises the fitted RT-PCR error distributions to simulate the number of mutations that would be observed from processing errors alone.

Conclusions

These data sets and models provide an effective means of separating true viral mutations from those erroneously introduced during sample processing and sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1456-x) contains supplementary material, which is available to authorized users.     

Variation in the ribosomal internal transcribed spacers (ITS1 and ITS2) among eight taxa of the Mimulus guttatus species complex

The internal transcribed spacer region (ITS1 and ITS2) of the 18S-25S nuclear ribosomal DNA sequence and the intervening 5.8S region were sequenced from three individuals in each of eight taxa of the Mimulus guttatus species complex. Three discrete variants, or "types," of ITS sequences were found, among which 30%-40% of sites differed, compared with 1%-2% within types. Dot plots indicate that these types were not related by conspicuous rearrangements or inversions. More than one ITS type was often found in the same taxon, and two of three ITS types span species boundaries, indicating their presence prior to speciation. These ITS sequences showed essentially no positional homology with the nearest sequenced relative, tomato. In contrast, the 5.8S region was relatively unvaried, with 8 of 162 sites varied in the sample among all eight taxa. The phylogeny inferred by the most common ITS sequence type, rooted by the two other ITS types, agreed with isozymes in showing the distinctness of M. nudatus, M. laciniatus, and M. tilingii from the other five taxa.      

Induction of haemolytic activity in Escherichia coli by the slyA gene product

The Salmonella typhimurium protein SlyA_ST, originally described as a cytolysin, shows sequence similarities to several known bacterial regulatory proteins. A homologue to the slyA_St gene has been localised to min 37 of the Eschericia coli K-12 chromosome and has been designated slyA_EC When introduced in trans on a plasmid, the slyA_EC gene conferred a haemolytic phenotype on wild-type but not clyA-knockout strains of E. coli K-12. The clyA gene encodes a novel haemolysin that is not expressed by wild-type E. coli under tested laboratory conditions. Western and Northern blot analyses, and DNA-band-shift assays support a model whereby the SlyA_EC protein activates clyA expression by binding to the clyA promoter region, thereby supporting the sequence similarity data in suggesting that SlyA_ST is a haemolysin activator rather than being a haemolysin per se.     

Acute toxicity effects of ibuprofen on behaviour and haematological parameters of African catfish Clarias gariepinus (Burchell, 1822)

Indiscriminate discharge of pharmaceutical waste into the aquatic ecosystem may pose serious health challenges to aquatic biota. The effect of acute exposure to ibuprofen was evaluated using changes in behaviour and haematological parameters under static bio-assay method in Clarias gariepinus. Test specimens were exposed to acute concentrations of ibuprofen (0.28, 0.33, 0.38, 0.43 and 0.48 mg l^?1) for 24, 48, 72 and 96 h durations respectively. Behavioural and phenotypic changes were observed in surviving fish. There were significant (p < 0.05) concentration and duration-dependent increases in erythrocyte (RBC), haemoglobin (Hb), pack cell volume (PCV) and leukocytes (WBC) in treated fish compared to the control. Insignificant decreases (p > 0.05) in mean corpuscular volume (MCV), mean corpuscular haemoglobin (MCH) and mean corpuscular haemoglobin concentration (MCHC) were observed in treated fish compared to the control. Ibuprofen elicited dose and duration- dependent decrease in neutrophil counts with the decreases being significant (p < 0.05) in the higher doses of 0.43 and 0.48 mg l^?1. Ibuprofen did not elicit any significant changes in monocytes, basophils and eosinophils. Changes observed in this study showed that ibuprofen negatively affected the health of the fish and we recommend that discharge of ibuprofen into the aquatic environment should be monitored and controlled.     

Isocitrate dehydrogenase 1 R132C mutation occurs exclusively in microsatellite stable colorectal cancers with the CpG island methylator phenotype

The CpG Island Methylator Phenotype (CIMP) is fundamental to an important subset of colorectal cancer; however, its cause is unknown. CIMP is associated with microsatellite instability but is also found in BRAF mutant microsatellite stable cancers that are associated with poor prognosis. The isocitrate dehydrogenase 1 (IDH1) gene causes CIMP in glioma due to an activating mutation that produces the 2-hydroxyglutarate oncometabolite. We therefore examined IDH1 alteration as a potential cause of CIMP in colorectal cancer. The IDH1 mutational hotspot was screened in 86 CIMP-positive and 80 CIMP-negative cancers. The entire coding sequence was examined in 81 CIMP-positive colorectal cancers. Forty-seven cancers varying by CIMP-status and IDH1 mutation status were examined using Illumina 450K DNA methylation microarrays. The R132C IDH1 mutation was detected in 4/166 cancers. All IDH1 mutations were in CIMP cancers that were BRAF mutant and microsatellite stable (4/45, 8.9%). Unsupervised hierarchical cluster analysis identified an IDH1 mutation-like methylation signature in approximately half of the CIMP-positive cancers. IDH1 mutation appears to cause CIMP in a small proportion of BRAF mutant, microsatellite stable colorectal cancers. This study provides a precedent that a single gene mutation may cause CIMP in colorectal cancer, and that this will be associated with a specific epigenetic signature and clinicopathological features.     

Specificity of siderophore receptors in membrane vesicles of Bacillus megaterium.

Membrane vesicles of Bacillus megaterium strains SK11 and Ard1 bound the ferrischizokinen and ferriferrioxamine B siderhores (iron transport cofactors). An approximately equimolar uptake of both labels of [3H, 59Fe]ferrischizokinen indicated binding of the intact chelate. Binding reached equilibrium in 2 to 5 min, was temperature independent, and was unaltered by the addition of several energy sources. A 91% dissociation of bound [Fe]ferrischizokinen was achieved in 60 s by the addition of excess ferrischizokinen. Ferriaerobactin, a siderophore which is structurally related to ferrischizokinen, caused no detectable release of bound [59Fe]ferrischizokinen. Of several other ferrigydroxamates tested, only ferriferrichrome A achieved the release (11%) of [Fe]ferrischizokinen. Rapid dissociation (92%) of bound [59Fe]ferriferrioxamine B by the addition of ferriferrioxamine B was observed, and a 67% release of [59Fe]ferriferrioxamine B was caused by ferriA2265, its structural relative. Ferrischizokinen, ferriferrichrome A, and ferrirhodotorulic acid produced a 6, 25, and 29% dissociation, respectively, of [59Fe]ferriferrioxamine B; ferriaerobactin caused no dissociation. [59Fe]ferriaerobactin was bound by the membranes, but its dissociation was not effected by unlabeled ferriaerobactin, suggesting no specific receptors for this chelate. The respective binding affinity constants and maximal binding capacities of membrane vesicles of strain SK11 were 2 x 10(7) M-1 and 280 pmol per mg of protein for ferrischizokinen and 7 x 10(7) M-1 and 37 pmol per mg of protein for ferriferrioxamine B. These values in strain Ard1 were, respectively, 1.4 x 10(7) M-1 and 186 pmol per mg of protein for ferrischizokinen and 11 x 10(7) M-1 and 23 pmol per mg of protein for ferriferrioxamine B. Separate, specific binding sites (receptors) for ferrischizokinen and ferriferrioxamine B exist on the vesicles. The ferrischizokinen receptors have a lower affinity but a higher binding capacity (eightfold) than that shown by the ferriferrioxamine B receptor. These receptors may be components of independent transport systems.     

Variation in heat shock proteins within tropical and desert species of poeciliid fishes

The 70-kilodalton heat shock protein (hsp70) family of molecular chaperones, which contains both stress-inducible and normally abundant constitutive members, is highly conserved across distantly related taxa. Analysis of this protein family in individuals from an outbred population of tropical topminnows, Poeciliopsis gracilis, showed that while constitutive hsp70 family members showed no variation in protein isoforms, inducibly synthesized hsp70 was polymorphic. Several species of Poeciliopsis adapted to desert environments exhibited lower levels of inducible hsp70 polymorphism than the tropical species, but constitutive forms were identical to those in P. gracilis, as they were in the confamilial species Gambusia affinis. These differences suggest that inducible and constitutive members of this family are under different evolutionary constraints and may indicate differences in their function within the cell. Also, northern desert species of Poeciliopsis synthesize a subset of the inducible hsp70 isoforms seen in tropical species. This distribution supports the theory that ancestral tropical fish migrated northward and colonized desert streams; the subsequent decrease in variation of inducible hsp70 may have been due to genetic drift or a consequence of adaptation to the desert environment. Higher levels of variability were found when the 30- kilodalton heat shock protein (hsp30) family was analyzed within different strains of two desert species of Poeciliopsis and also in wild-caught individuals of Gambusia affinis. In both cases the distribution of hsp30 isoform diversity was similar to that seen previously with allozyme polymorphisms.      

Neuronal synchrony mediated by astrocytic glutamate through activation of extrasynaptic NMDA receptors

Fast excitatory neurotransmission is mediated by activation of synaptic ionotropic glutamate receptors. In hippocampal slices, we report that stimulation of Schaffer collaterals evokes in CA1 neurons delayed inward currents with slow kinetics, in addition to fast excitatory postsynaptic currents. Similar slow events also occur spontaneously, can still be observed when neuronal activity and synaptic glutamate release are blocked, and are found to be mediated by glutamate released from astrocytes acting preferentially on extrasynaptic NMDA receptors. The slow currents can be triggered by stimuli that evoke Ca2+ oscillations in astrocytes, including photolysis of caged Ca2+ in single astrocytes. As revealed by paired recording and Ca2+ imaging, a striking feature of this NMDA receptor response is that it occurs synchronously in multiple CA1 neurons. Our results reveal a distinct mechanism for neuronal excitation and synchrony and highlight a functional link between astrocytic glutamate and extrasynaptic NMDA receptors.