BACKGROUND
apple stem grooving virus  is a member of the genus capillovirus in the family betaflexiviridae  <cit> .
asgv has been most commonly identified from apple, european pear, japanese pear, and citrus trees  <cit> .
in addition, asgv has been identified in lily  <cit>  and kiwi  <cit> , and it infects several virus indicator plants, including chenopodium, cucumber, nicotiana, phaseolus, and vigna species  <cit> .
asgv infection in fruit trees is usually latent without disease symptoms  <cit> ; however, asgv sometimes causes serious viral diseases  <cit> .
in many cases, fruit trees are co-infected by different viruses and viroids.
for instance, apple trees showing fruit deformation, leaf deformation, and mosaic, chlorosis, and rusting symptoms in india were co-infected by apple chlorotic leaf spot virus , apple mosaic virus , asgv, apple stem pitting virus , and apple scar skin viroid   <cit> .
the viral particles of asgv are flexuous filaments 620–680 nm long and 12 nm wide  <cit> .
asgv has a single-stranded  positive-sense monopartite rna genome containing 5′ capping and a poly tail at the 3′ region  <cit> .
the genome size of asgv is about  <dig>  ~  <dig>  nucleotides , and it encodes two overlapping open reading frames , orf <dig>  and orf <dig>   <cit> .
orf <dig> encodes a polyprotein containing a replicase and coat protein , while orf <dig> encodes a movement protein  that is overlapped with the replicase and cp regions  <cit> .
a previous study demonstrated that asgv mutants with a stop codon between the replicase and cp coding regions were capable of systemic infection with decrease of pathogenicity  <cit> .
this result revealed that expression of asgv cp via a subgenomic rna  was sufficient for viability of asgv.
furthermore, mutational analysis revealed core promoter sequences required for the sgrna transcription of asgv and potato virus t, which were conserved among viruses in the families alphaflexiviridae and betaflexiviridae  <cit> .
next-generation sequencing  produces huge amounts of sequencing data, which facilitate the identification of known and novel viruses and viroids in a wide range of plant species  <cit> .
in addition, ngs can be applied in plant virus diagnostics  <cit>  and virus ecology  <cit> .
several types of ngs platforms—including hiseq systems by illumina,  <dig> flx systems by roche, and solid systems by ab—have been developed  <cit> .
each ngs system has advantages and disadvantages  <cit> .
the selection of proper ngs platforms is dependent on the purposes of the study.
hiseq systems produce high throughput with a relatively shorter read length, whereas  <dig> flx systems generate low throughput with a longer read length.
for example, the identification and diagnostics of known and novel viruses can be conducted by hiseq systems  <cit> , and viral genome sequencing using extracted viral rnas can be performed by  <dig> flx systems  <cit> .
moreover, ngs systems are useful for virus–host interaction studies.
for instance, srna-seq has been used for virus-derived sirnas  of asgv from asgv-infected samples  <cit> .
this study showed an increase in sirna production towards the 3′ end of asgv and several trna-derived srnas were differentially regulated by asgv infection.
a previous study identified  <dig> conserved and  <dig> novel mirnas of pear associated with asgv infection and found several mirnas in response to high temperature, which was used to reduce asgv titers in the shoot meristem tip  <cit> .
pear transcriptome analysis between asgv-infected and asgv-free apple samples has been conducted and identified  <dig> up-regulated and  <dig> down-regulated genes in asgv infected shoot culture as compared to asgv-free shoot culture  <cit> .
several approaches to detect asgv have been developed, such as long-distance pcr  to amplify the complete genome of asgv  <cit> , multiplex reverse transcriptase -pcr for major apple viruses  and pear viruses  <cit> , and immunochromatographic assays by monoclonal antibodies specific for cp  <cit> .
moreover, using available genome sequences for asgv, two phylogenetic groups and four recombinants of  <dig> asgv isolates have been identified  <cit> , and the molecular evolution of subgenomic rna of asgv has been studied  <cit> .
several recent studies have demonstrated that many plant transcriptomes contain viral sequences that could be applied to studies associated with virus identification and viral genome assembly  <cit> .
in this study, we conducted in silico analyses using publicly available transcriptome data for viral genome assembly and identification using asgv as an exemplar virus.
we showed the application of transcriptome data for the analysis of single nucleotide variations  on the asgv genome.
moreover, the two viral genomes obtained were successfully applied in the phylogenetic and recombination analyses of known asgv genomes.
RESULTS
identification and de novo assembly of asgv genome from asgv-infected apple mrna transcriptome
of the known viruses, we selected asgv, which mostly infects fruit trees, including apple  and pear .
due to the clonal propagation of fruit trees, the possibility of virus infection is very high.
we screened several apple and pear transcriptomes and selected the transcriptomes infected by asgv for further study .
of the several previously reported apple transcriptomes in response to asgv infection, we selected two studies: one that performed mrna sequencing   <cit>  and one that performed srna-seq  <cit> .
both sets of samples included asgv-infected and asgv-free apple plants .
the first study was conducted to examine the expression profiles of apple trees infected with asgv without any disease symptoms using mrna-seq  <cit> .
two libraries from asgv-infected and asgv-free shoot cultures were constructed.
we de novo assembled transcriptomes of the two libraries using the trinity program .
the  <dig>  and  <dig>  contigs obtained from asgv-infected and asgv-free shoot cultures, respectively, were blasted against a viral reference genome database.
as a result, we identified  <dig> and  <dig> asgv-associated contigs from asgv-infected and asgv-free samples, respectively .
interestingly, the  <dig> identified asgv-associated contigs mostly covered the complete genome of asgv.
thus, we assembled a nearly complete asgv genome with  <dig>  nt.
the newly assembled asgv genome was referred to as asgv isolate fuji with the accession number ku <dig>  asgv isolate fuji contained two genes encoding the 241 kda polyprotein and 36 kda protein, respectively .
the 241 kda polyprotein is associated with viral methyltransferase, duf <dig>  helicase, rna-dependent rna polymerase , and the cp at the 3′ end, while the 36 kda protein is known as a viral mp .
the previous study sequenced the complete genome of asgv  from the asgv-infected sample by the sanger sequencing method.
to compare the asgv genome sequences obtained by the sanger method and de novo assembly, we aligned the two genome sequences by clustalw.
the two genome sequences were almost identical except the 3′ region, which showed many polymorphisms, indicating the presence of asgv variants in the transcriptomes .table  <dig> blast results to identify asgv-associated contigs from asgv-infected apple mrna-seq data
the megablast results with the best hits were listed.
subject ids indicates the identity number of an individual assembled contig.
subject ids indicate the best matched viral genome
fig.
<dig> identification of de novo viral genome assembly and snvs for asgv from rna-seq data.
a genome structure of asgv isolate fuji.
the conserved domains were identified by the smart program .
abbreviations: mt , hel , rna-dependent rna polymerase , movement protein , and coat protein .
b alignment of raw data against genome of asgv isolate fuji by bwa was visualized by tablet program.
c positions of identified snvs in asgv-infected apple transcriptome were visualized by tablet program.
d identified sequence reads from asgv-free apple sample, which were associated with asgv by bwa alignment.
e blast results showing sequence reads from asgv-free apple sample matched to asgv genome.
f alignment of raw data using asgv-infected srna data from cultivar gd against reference asgv genome by bwa was visualized by tablet program.
g positions of identified snvs in asgv-infected apple srna transcriptome were visualized by tablet program.
h alignment of raw data from pear sample against genome asgv isolate cuiguan by bwa was visualized by tablet program.
i positions of identified snvs of asgv in pear mrna transcriptome were visualized by tablet program
it is well known that rna viruses have a quasispecies nature with a high mutation rate within infected hosts.
thus, we analyzed the snvs of asgv in the asgv-infected sample.
we mapped raw data on the genome of asgv isolate fuji, and interestingly, reads were highly mapped on the regions for cp and mp .
using the samtools program, we identified  <dig> snvs.
in particular, many snvs were identified in the 5′ and 3′ regions of the asgv genome .
in many previous studies, the assembled contigs or transcripts were frequently used to identify viruses or viroids in the host transcriptome  <cit> .
although the assembled contigs did not contain any viral sequences in the asgv-free sample, it is possible that the raw sequence data contained viral sequences.
the single- or paired-end mrna sequencing by hiseq <dig> produces raw sequence data up to 101 bp in size.
therefore, the raw data can also be successfully applied to identify viral sequences in the host transcriptome data.
we aligned a raw fastq file from the asgv-free sample on the genome of asgv isolate fuji using the bwa program.
as shown in fig.
1d,  <dig> sequenced reads were mapped on the genome of asgv isolate fuji.
to confirm the alignment results, we blasted the fasta converted sequences against the asgv genome.
we found that  <dig> sequenced reads were aligned along the asgv genome .
the mapping and blast results using sequenced raw data clearly demonstrated the presence of asgv viral sequences in the asgv-free sample.
identification and de novo genome assembly of asgv from asgv-infected srna transcriptomes
previous studies have demonstrated that both mrna-seq and srna-seq are useful for virus identification  <cit> .
to validate the utility of srna-seq data for the de novo assembly of the asgv genome, we used srna data from a previous study that conducted apple leaf srna sequencing using samples from the apple cultivar golden delicious   <cit> .
the data were composed of  <dig> libraries from asgv-infected and asgv-free samples .
moreover, two different types of libraries were generated according to size fraction  <cit> .
the six libraries from asgv-infected samples were subjected to de novo transcriptome assembly using the trinity program followed by a blast search to identify viral contigs.
however, we obtained only  <dig> contigs with 425 bp of n <dig> value, and no asgv-associated contigs were identified by the blast search.
it seems that the trinity program was not optimal for de novo transcriptome assembly using srna data.
thus, we used the velvet program, which is well known for srna transcriptome assembly  <cit> .
the velvet assembler assembled a total of  <dig>  contigs, which were blasted against a plant viral database identifying  <dig> contigs associated with asgv .
we mapped the identified asgv-associated contigs on the reference genome of asgv .
the  <dig> contigs covered about 30 % of the asgv genome and displayed many gaps along the genome.
in order to confirm that srna reads cannot cover the complete genome of asgv, we mapped srna raw data on the asgv reference genome .
we found that several regions of asgv were not mapped by srna sequences.
based on the mapping results, we also identified  <dig> snvs from srna data by the sam toolkit .
identification and de novo genome assembly of asgv from pear mrna transcriptome
we used pear transcriptome data from a previous study that did not include any information on the virus infection.
the transcriptome data  was derived from a mixture of nine different fruit developmental stages of the pyrus pyrifolia cultivar cuiguan.
the transcriptome was initially assembled by soapdenovo2; however, we performed de novo transcriptome assembly again using the trinity program.
a total of  <dig>  transcripts were assembled .
assembled sequences were subjected to a blast search against a viral reference database.
we found nine contigs associated with asgv ranging from 222 bp to  <dig>  bp .
of the nine contigs associated with asgv, a single contig with  <dig>  bp was a nearly complete genome sequence of asgv.
after removing poly tails from the contig, we obtained a sequence with  <dig>  nt referred to as asgv isolate cuiguan .
in order to identify additional viruses infecting pears, all raw data converted to fasta format were blasted against the viral reference database.
interestingly, we found many additional viruses infecting pears .
of  <dig> viruses, six viruses including apricot latent virus, grapevine fleck virus, rupestris stem pitting associated virus- <dig>  aclsv, grapevine pinot gris virus, and zucchini yellow mosaic virus with very small numbers of reads were identified.
based on our knowledge, it seemed that the six identified viruses were not likely viruses infecting pears.
they might have been sequences that were partially homologous to host genes or other viral genomes.
in addition, associations of the six viruses with pears have not been reported.
the sequence reads associated with potato leafroll virus were identified as sequences from the host.
of four identified viruses infecting pears, asgv was dominant followed by prunus virus t , apple green crinkle associated virus , and aspv.table  <dig> identification of viruses from raw mrna-seq data of pear transcriptome by blast search
apple stem grooving virus
potato leafroll virus
prunus virus t isolate aze239
apple green crinkle associated virus
apple stem pitting virus
apricot latent virus
grapevine fleck virus
rupestris stem pitting associated virus-1
apple chlorotic leaf spot virus
grapevine pinot gris virus
zucchini yellow mosaic virus
we examined snvs for asgv isolate cuiguan within the pear transcriptome after alignment of the raw data on the asgv isolate cuiguan .
we found  <dig> snvs in the whole asgv genome .
interestingly, snvs were only identified in the replicase region containing helicase, rdrp .
however, snvs were not detected in the region of mp or cp.
of the identified nucleotide changes, c to t  was dominant followed by t to c , g to a , and a to g .
comparison of de novo sequence assemblers for viral genome assembly
in this study, we used two different programs for de novo transcriptome assembly, trinity and velvet.
to find the advantages and disadvantages of the two programs, we compared the number of total contigs and sizes of viral contigs.
we first compared the number of contigs of two different mrna libraries assembled by the two programs .
the number of contigs assembled by velvet was more than  <dig>  to  <dig>  times that assembled by trinity .
in addition, the number of identified viral contigs by velvet was more than four times that identified by trinity.
however, the portion of viral contigs in the transcriptome assembled by trinity was higher than that assembled by velvet .
moreover, the viral contigs assembled by trinity were much bigger than those assembled by velvet .
as a result, the velvet assembler assembled large numbers of contigs with relatively short lengths, while the trinity assembler assembled a few contigs with relatively long lengths.
for example, in the transcriptome from the srr <dig>  the longest contigs assembled by trinity were  <dig>  bp, while the longest contigs assembled by velvet were 646 bp .
furthermore, the velvet assembler assembled seven contigs associated with aspv and agcav that were not assembled by trinity .table  <dig> comparison of de novo transcriptome assemblers for assembly of viral contigs
srr <dig> and srr <dig> were data derived from pear and apple samples, respectively
fig.
<dig> comparison of two de novo assemblers based on number and sizes of assembled contigs and phylogenetic tree of  <dig> asgv isolates.
size distribution of identified viral contigs from asgv-infected apple sample assembled by trinity  and velvet .
size distribution of identified viral contigs from pear sample assembled by trinity  and velvet .
the green-, blue-, and red-colored bars indicate aspv, agcav, and asgv, respectively.
the sizes of only the longest and the shortest contigs in each transcriptome are indicated.
the phylogenetic tree was constructed based on the genome sequences of  <dig> asgv isolates and  <dig> ctlv isolates.
we followed the original annotations for ctlv and pbnlsv, which were highly homologous to asgv.
the accession number and the name of each isolate were indicated.
detailed information for each isolate can be found in additional file 10
phylogenetic analysis of asgv isolates
several previous studies have reported that asgv is closely related to citrus tatter leaf virus   <cit> .
to confirm previous results, we blast identified two asgv genomes in this study against the ncbi nucleotide database.
the blast results confirmed that ctlv is closely grouped with asgv isolates in the genus capillovirus.
from the genbank, we retrieved all asgv-associated sequences as well as ctlv-associated sequences.
after removing partial sequences, we collected a total of  <dig> genomes of asgv and ctlv isolates, including two asgv isolates in this study.
the host ranges of ctlv were mostly from citrus species as well as lilium species .
pear black necrotic leaf spot virus  isolated from pear was an isolate of asgv according to the annotation in genbank  <cit> .
most asgv isolates were isolated from apple and pear, and some isolates, such as asgv isolates matsuco and li- <dig>  were identified from citrus tamurana and lily, respectively.
to reveal phylogenetic relationships, we aligned genome sequences displaying high sequence similarity of asgv and ctlv.
sequence alignment and a phylogenetic tree using genome sequences of asgv and ctlv identified two largely divided clades .
the first clade contained  <dig> genomes, while the second clade included only pbnlsv and asgv isolate kfp.
the first clade could be further divided into three groups.
group a consisted of six ctlv isolates and a single asgv isolate, while group b contained only asgv isolates.
group c was the largest, including seven asgv isolates and two ctlv isolates.
recombination analysis for  <dig> asgv isolates
we analyzed recombination events among  <dig> asgv isolates.
the aligned genome sequences were subjected to the rdp <dig> program, which includes nine different algorithms for recombination detection.
the rdp <dig> program detected a total of  <dig> recombination events.
of them, we selected six recombination events that were supported by at least five recombination algorithms .
for example, pbnlsv contains two recombination sequences from asgv isolate li- <dig>  three isolates—asgv isolate chn, asgv isolate hh, and ctlv isolate mth—include recombination sequences from asgv isolate fuji in the 5′ region .
recombination events  <dig> and  <dig> were supported by seven algorithms.
the major parent of the recombinant sequence for asgv isolate ytg was asgv isolate fuji .
the major parent of the recombinant sequences for the three isolates—asgv isolate hh, asgv isolate chn, and ctlv isolate mth—were asgv isolate li- <dig> and ctlv isolate lily .fig.
<dig> identification of recombination events by rdp <dig> program.
a the positions of identified recombinants were depicted with the respective names of parental sequences.
the individual genome of the asgv isolate was indicated by a different colored bar.
the identified recombination events, including number  <dig>  and number  <dig> , were rechecked by the rdp <dig> program and visualized by plot data with pairwise identity information.
detailed information on the identified recombination events is provided in table 4
event
minor parent parent contributing smaller fraction of sequence, major parent parent contributing larger fraction of sequence, unknown only one parent and a recombinant need be in the alignment for a recombination event to be detectable.
the sequence listed as unknown was used to infer the existence of a missing parental sequence, ns no significant p-value was recorded for this recombination event using this method
athe actual breakpoint position is undetermined
discussion
the rapid development of ngs is enabling virologists to find viruses from numerous species  <cit> .
ngs-based approaches have identified not only known viruses but also novel viruses  <cit> .
in fact, many horticultural plants are frequently infected by viruses and viroids  <cit> .
in particular, fruit trees usually propagated by grafting and cuttage are reservoirs of various plant viruses and viroids  <cit> .
in addition, the big data produced by ngs techniques has prompted virus identification in silico  <cit> .
here, we discussed the library types, sequencing methods, and de novo assembler for virus identification and viral genome assembly.
the majority of plant viruses are composed of rna genomes, and dna viruses also replicate via an rna intermediate  <cit> .
thus, rna-based transcriptome libraries are preferable to dna-based genome libraries for virus identification.
in the current study, we used published plant transcriptome data.
to enrich viral rnas, ribosome-deleted libraries are usually prepared using extracted total rnas from virus-infected samples  <cit> .
however, we demonstrated that the mrna libraries using oligo d were successfully applied for virus identification.
of course, rna viruses with poly tails such as asgv are also easily identified by mrna libraries.
similarly, several polyadenylated rna viruses have been identified from sweet potato transcriptomes  <cit> .
several recent studies have also demonstrated that ribosome-deleted rna libraries as well as plant mrna libraries are suitable for the identification of viruses without poly tails or viroids  <cit> .
therefore, it might be ideal to use ribosome-deleted libraries for studies only focused on viruses.
in the case of studies of both viruses and host plants, mrna libraries can be usefully applied  <cit> .
in this study, we used data from two different library types, including mrna and srna libraries that were single-end sequenced by the hiseq <dig> system.
according to many recent studies, viral genomes have been de novo assembled from mrna as well as srna data  <cit> .
in our study, we assembled nearly complete genomes of two asgv isolates from the mrna data; however, the srna data could cover only 30 % of the asgv genome.
we compared the numbers of sequencing reads between the mrna and srna data.
however, the numbers of sequence reads between mrna and srna were very similar, indicating that the sequencing amount is not an important factor for viral genome assembly.
in fact, when the number of sequencing reads is increased, the number of viral-associated reads is increased.
therefore, the quantity of the sequenced data might play an important role in de novo genome assembly.
the number of pear  transcriptomes was about ten times that of apple transcriptomes .
the sequence reads associated with asgv were  <dig>  viral reads out of  <dig> , <dig> reads for the apple sample and  <dig>  viral reads out of  <dig> , <dig> reads for the pear sample.
although the number of total sequence reads in the apple sample was much smaller than that in the pear sample, the number of sequence reads associated with asgv was about  <dig>  times higher.
this result suggests the amount of viral replication in the host might be also an important factor in de novo viral genome assembly.
the portion of viral nucleic acids in the sample infected by virus is often low suggesting enrichment of virions prior to ngs  <cit> .
for example, purification of double-stranded  rnas from the prunus species followed by  <dig> pyrosequencing enabled to assemble four complete genomes of asian prunus virus  <dig> , apv <dig>  and apv <dig>  <cit> .
this study demonstrated successful application of dsrna purification for virus genome assembly using ngs technique.
in the case of srna, two different types were prepared based on size fraction  <cit> .
the libraries without size fraction contain a large number of asgv-associated reads, but the libraries with size fraction contain very few reads associated with asgv.
of course, the srna libraries were targeted for the identification of viral srnas.
we suppose that the small number of srnas might be related to the ability of the rna silencing machinery in the host.
in any case, a sufficient number of viral-associated reads is necessary for viral de novo genome assembly.
in addition, sequencing methods are important for virus identification and viral genome assembly.
in this study, all transcriptome data were single-end sequenced by hiseq <dig>  as compared to single-end sequencing, paired-end sequencing provides sequences from both ends of a fragment and generates high-quality and alignable sequence data.
the advantages of paired-end sequencing have been previously reported  <cit> .
thus, paired-end sequencing was far superior for the identification and genome assembly of the target virus.
for virus identification, assembled contigs are frequently used.
therefore, the choice of de novo assembler affects the quality and quantity of virus identification.
for instance, mrna data were very efficiently assembled by trinity; however, few and low-quality contigs were assembled from the srna data by trinity.
our comparative studies between the two de novo assemblers suggest trinity and velvet for de novo assembly of mrna data and srna data, respectively.
the obtained viral contigs assembled by trinity from mrna data were low in number but long in length, while the viral contigs assembled by velvet were high in number but short in length.
for the de novo assembly of a target virus with high-quality mrna data, trinity is ideal.
velvet cannot assemble a nearly complete viral genome, but it assembled many contigs, which enabled us to identify additional viruses, for example, viruses in the pear transcriptomes.
recently, several programs iva, price, and vicuna for de novo assembly of rna virus genome have been developed .
the choice of optimal de novo assembler might be dependent on researchers and purposes.
it is well known that rna viruses have a quasispecies nature within the host  <cit> .
however, to date, most studies have shown the variants and mutation rates of target viruses using cloning-based sanger sequencing methods  <cit> .
in this study, we successfully demonstrated the usefulness of plant transcriptome data for revealing the snvs of asgv.
in fact, it is quite difficult to find virus variants using transcriptome data, while cloning-based sequencing methods might reveal variants.
however, the cloning-based approaches require a rt-pcr amplification procedure to amplify full-length viral genomes.
practically, the amplification of full-length viral genomes is not easy even though plant viruses are relatively small.
we showed the presence of asgv variants in the transcriptome by comparing the asgv genome from the cultivar fuji derived from the sanger-sequencing method and de novo assembly.
we did not judge which asgv genome was the dominant asgv genome; however, it is highly likely that the de novo-assembled asgv was a consensus genome sequence of asgv.
the mutation rates of identified asgv genomes were varied:  <dig>  %  in the fuji, 1 %  in the gd, and  <dig>  %  in the cuiguan.
we suppose that several factors—including hosts, viral replication, and environmental cues—might affect the mutation rates.
the association of viral mutation rate with other factors will be an interesting subject for further study  <cit> .
CONCLUSIONS
taken together, our study showed the successful application of plant transcriptome data for virus identification, viral genome assembly, and viral mutation rates.
in addition, we discussed several factors, including library preparation, ngs systems, de novo assemblers, and sample conditions for virus identification and genome assembly.

 apple stem grooving virus is a member of the genus capillovirus in the family betaflexiviridae. asgv has been most commonly identified from apple, european pear, japanese pear, and citrus trees. in addition, asgv has been identified in lily and kiwi, and it infects several virus indicator plants, including chenopodium, cucumber, nicotiana, phaseolus, and vigna species. asgv infection in fruit trees is usually latent without disease symptoms; however, asgv sometimes causes serious viral diseases. in many cases, fruit trees are co-infected by different viruses and viroids. for instance, apple trees showing fruit deformation, leaf deformation, and mosaic, chlorosis, and rusting symptoms in india were co-infected by apple chlorotic leaf spot virus, apple mosaic virus, asgv, apple stem pitting virus, and apple scar skin viroid. the viral particles of asgv are flexuous filaments 620–680 nm long and 12 nm wide. asgv has a single-stranded positive-sense monopartite rna genome containing 5′ capping and a poly tail at the 3′ region. the genome size of asgv is about  ~  nucleotides, and it encodes two overlapping open reading frames, orf and orf. orf encodes a polyprotein containing a replicase and coat protein, while orf encodes a movement protein that is overlapped with the replicase and cp regions. a previous study demonstrated that asgv mutants with a stop codon between the replicase and cp coding regions were capable of systemic infection with decrease of pathogenicity. this result revealed that expression of asgv cp via a subgenomic rna was sufficient for viability of asgv. furthermore, mutational analysis revealed core promoter sequences required for the sgrna transcription of asgv and potato virus t, which were conserved among viruses in the families alphaflexiviridae and betaflexiviridae. next-generation sequencing produces huge amounts of sequencing data, which facilitate the identification of known and novel viruses and viroids in a wide range of plant species. in addition, ngs can be applied in plant virus diagnostics and virus ecology. several types of ngs platforms—including hiseq systems by illumina, flx systems by roche, and solid systems by ab—have been developed. each ngs system has advantages and disadvantages. the selection of proper ngs platforms is dependent on the purposes of the study. hiseq systems produce high throughput with a relatively shorter read length, whereas flx systems generate low throughput with a longer read length. for example, the identification and diagnostics of known and novel viruses can be conducted by hiseq systems, and viral genome sequencing using extracted viral rnas can be performed by flx systems. moreover, ngs systems are useful for virus–host interaction studies. for instance, srna-seq has been used for virus-derived sirnas of asgv from asgv-infected samples. this study showed an increase in sirna production towards the 3′ end of asgv and several trna-derived srnas were differentially regulated by asgv infection. a previous study identified conserved and novel mirnas of pear associated with asgv infection and found several mirnas in response to high temperature, which was used to reduce asgv titers in the shoot meristem tip. pear transcriptome analysis between asgv-infected and asgv-free apple samples has been conducted and identified up-regulated and down-regulated genes in asgv infected shoot culture as compared to asgv-free shoot culture. several approaches to detect asgv have been developed, such as long-distance pcr to amplify the complete genome of asgv, multiplex reverse transcriptase -pcr for major apple viruses and pear viruses, and immunochromatographic assays by monoclonal antibodies specific for cp. moreover, using available genome sequences for asgv, two phylogenetic groups and four recombinants of asgv isolates have been identified, and the molecular evolution of subgenomic rna of asgv has been studied. several recent studies have demonstrated that many plant transcriptomes contain viral sequences that could be applied to studies associated with virus identification and viral genome assembly. in this study, we conducted in silico analyses using publicly available transcriptome data for viral genome assembly and identification using asgv as an exemplar virus. we showed the application of transcriptome data for the analysis of single nucleotide variations on the asgv genome. moreover, the two viral genomes obtained were successfully applied in the phylogenetic and recombination analyses of known asgv genomes.  identification and de novo assembly of asgv genome from asgv-infected apple mrna transcriptome of the known viruses, we selected asgv, which mostly infects fruit trees, including apple and pear. due to the clonal propagation of fruit trees, the possibility of virus infection is very high. we screened several apple and pear transcriptomes and selected the transcriptomes infected by asgv for further study. of the several previously reported apple transcriptomes in response to asgv infection, we selected two studies: one that performed mrna sequencing and one that performed srna-seq. both sets of samples included asgv-infected and asgv-free apple plants. the first study was conducted to examine the expression profiles of apple trees infected with asgv without any disease symptoms using mrna-seq. two libraries from asgv-infected and asgv-free shoot cultures were constructed. we de novo assembled transcriptomes of the two libraries using the trinity program. the and contigs obtained from asgv-infected and asgv-free shoot cultures, respectively, were blasted against a viral reference genome database. as a result, we identified and asgv-associated contigs from asgv-infected and asgv-free samples, respectively. interestingly, the identified asgv-associated contigs mostly covered the complete genome of asgv. thus, we assembled a nearly complete asgv genome with  nt. the newly assembled asgv genome was referred to as asgv isolate fuji with the accession number ku asgv isolate fuji contained two genes encoding the 241 kda polyprotein and 36 kda protein, respectively. the 241 kda polyprotein is associated with viral methyltransferase, duf helicase, rna-dependent rna polymerase, and the cp at the 3′ end, while the 36 kda protein is known as a viral mp. the previous study sequenced the complete genome of asgv from the asgv-infected sample by the sanger sequencing method. to compare the asgv genome sequences obtained by the sanger method and de novo assembly, we aligned the two genome sequences by clustalw. the two genome sequences were almost identical except the 3′ region, which showed many polymorphisms, indicating the presence of asgv variants in the transcriptomes.table blast results to identify asgv-associated contigs from asgv-infected apple mrna-seq data the megablast results with the best hits were listed. subject ids indicates the identity number of an individual assembled contig. subject ids indicate the best matched viral genome fig. identification of de novo viral genome assembly and snvs for asgv from rna-seq data. a genome structure of asgv isolate fuji. the conserved domains were identified by the smart program. abbreviations: mt, hel, rna-dependent rna polymerase, movement protein, and coat protein. b alignment of raw data against genome of asgv isolate fuji by bwa was visualized by tablet program. c positions of identified snvs in asgv-infected apple transcriptome were visualized by tablet program. d identified sequence reads from asgv-free apple sample, which were associated with asgv by bwa alignment. e blast results showing sequence reads from asgv-free apple sample matched to asgv genome. f alignment of raw data using asgv-infected srna data from cultivar gd against reference asgv genome by bwa was visualized by tablet program. g positions of identified snvs in asgv-infected apple srna transcriptome were visualized by tablet program. h alignment of raw data from pear sample against genome asgv isolate cuiguan by bwa was visualized by tablet program. i positions of identified snvs of asgv in pear mrna transcriptome were visualized by tablet program it is well known that rna viruses have a quasispecies nature with a high mutation rate within infected hosts. thus, we analyzed the snvs of asgv in the asgv-infected sample. we mapped raw data on the genome of asgv isolate fuji, and interestingly, reads were highly mapped on the regions for cp and mp. using the samtools program, we identified snvs. in particular, many snvs were identified in the 5′ and 3′ regions of the asgv genome. in many previous studies, the assembled contigs or transcripts were frequently used to identify viruses or viroids in the host transcriptome. although the assembled contigs did not contain any viral sequences in the asgv-free sample, it is possible that the raw sequence data contained viral sequences. the single- or paired-end mrna sequencing by hiseq produces raw sequence data up to 101 bp in size. therefore, the raw data can also be successfully applied to identify viral sequences in the host transcriptome data. we aligned a raw fastq file from the asgv-free sample on the genome of asgv isolate fuji using the bwa program. as shown in fig. 1d, sequenced reads were mapped on the genome of asgv isolate fuji. to confirm the alignment results, we blasted the fasta converted sequences against the asgv genome. we found that sequenced reads were aligned along the asgv genome. the mapping and blast results using sequenced raw data clearly demonstrated the presence of asgv viral sequences in the asgv-free sample. identification and de novo genome assembly of asgv from asgv-infected srna transcriptomes previous studies have demonstrated that both mrna-seq and srna-seq are useful for virus identification. to validate the utility of srna-seq data for the de novo assembly of the asgv genome, we used srna data from a previous study that conducted apple leaf srna sequencing using samples from the apple cultivar golden delicious. the data were composed of libraries from asgv-infected and asgv-free samples. moreover, two different types of libraries were generated according to size fraction. the six libraries from asgv-infected samples were subjected to de novo transcriptome assembly using the trinity program followed by a blast search to identify viral contigs. however, we obtained only contigs with 425 bp of n value, and no asgv-associated contigs were identified by the blast search. it seems that the trinity program was not optimal for de novo transcriptome assembly using srna data. thus, we used the velvet program, which is well known for srna transcriptome assembly. the velvet assembler assembled a total of contigs, which were blasted against a plant viral database identifying contigs associated with asgv. we mapped the identified asgv-associated contigs on the reference genome of asgv. the contigs covered about 30 % of the asgv genome and displayed many gaps along the genome. in order to confirm that srna reads cannot cover the complete genome of asgv, we mapped srna raw data on the asgv reference genome. we found that several regions of asgv were not mapped by srna sequences. based on the mapping results, we also identified snvs from srna data by the sam toolkit. identification and de novo genome assembly of asgv from pear mrna transcriptome we used pear transcriptome data from a previous study that did not include any information on the virus infection. the transcriptome data was derived from a mixture of nine different fruit developmental stages of the pyrus pyrifolia cultivar cuiguan. the transcriptome was initially assembled by soapdenovo2; however, we performed de novo transcriptome assembly again using the trinity program. a total of transcripts were assembled. assembled sequences were subjected to a blast search against a viral reference database. we found nine contigs associated with asgv ranging from 222 bp to  bp. of the nine contigs associated with asgv, a single contig with  bp was a nearly complete genome sequence of asgv. after removing poly tails from the contig, we obtained a sequence with  nt referred to as asgv isolate cuiguan. in order to identify additional viruses infecting pears, all raw data converted to fasta format were blasted against the viral reference database. interestingly, we found many additional viruses infecting pears. of viruses, six viruses including apricot latent virus, grapevine fleck virus, rupestris stem pitting associated virus- aclsv, grapevine pinot gris virus, and zucchini yellow mosaic virus with very small numbers of reads were identified. based on our knowledge, it seemed that the six identified viruses were not likely viruses infecting pears. they might have been sequences that were partially homologous to host genes or other viral genomes. in addition, associations of the six viruses with pears have not been reported. the sequence reads associated with potato leafroll virus were identified as sequences from the host. of four identified viruses infecting pears, asgv was dominant followed by prunus virus t, apple green crinkle associated virus, and aspv.table identification of viruses from raw mrna-seq data of pear transcriptome by blast search apple stem grooving virus potato leafroll virus prunus virus t isolate aze239 apple green crinkle associated virus apple stem pitting virus apricot latent virus grapevine fleck virus rupestris stem pitting associated virus-1 apple chlorotic leaf spot virus grapevine pinot gris virus zucchini yellow mosaic virus we examined snvs for asgv isolate cuiguan within the pear transcriptome after alignment of the raw data on the asgv isolate cuiguan. we found snvs in the whole asgv genome. interestingly, snvs were only identified in the replicase region containing helicase, rdrp. however, snvs were not detected in the region of mp or cp. of the identified nucleotide changes, c to t was dominant followed by t to c, g to a, and a to g. comparison of de novo sequence assemblers for viral genome assembly in this study, we used two different programs for de novo transcriptome assembly, trinity and velvet. to find the advantages and disadvantages of the two programs, we compared the number of total contigs and sizes of viral contigs. we first compared the number of contigs of two different mrna libraries assembled by the two programs. the number of contigs assembled by velvet was more than to times that assembled by trinity. in addition, the number of identified viral contigs by velvet was more than four times that identified by trinity. however, the portion of viral contigs in the transcriptome assembled by trinity was higher than that assembled by velvet. moreover, the viral contigs assembled by trinity were much bigger than those assembled by velvet. as a result, the velvet assembler assembled large numbers of contigs with relatively short lengths, while the trinity assembler assembled a few contigs with relatively long lengths. for example, in the transcriptome from the srr the longest contigs assembled by trinity were  bp, while the longest contigs assembled by velvet were 646 bp. furthermore, the velvet assembler assembled seven contigs associated with aspv and agcav that were not assembled by trinity.table comparison of de novo transcriptome assemblers for assembly of viral contigs srr and srr were data derived from pear and apple samples, respectively fig. comparison of two de novo assemblers based on number and sizes of assembled contigs and phylogenetic tree of asgv isolates. size distribution of identified viral contigs from asgv-infected apple sample assembled by trinity and velvet. size distribution of identified viral contigs from pear sample assembled by trinity and velvet. the green-, blue-, and red-colored bars indicate aspv, agcav, and asgv, respectively. the sizes of only the longest and the shortest contigs in each transcriptome are indicated. the phylogenetic tree was constructed based on the genome sequences of asgv isolates and ctlv isolates. we followed the original annotations for ctlv and pbnlsv, which were highly homologous to asgv. the accession number and the name of each isolate were indicated. detailed information for each isolate can be found in additional file 10 phylogenetic analysis of asgv isolates several previous studies have reported that asgv is closely related to citrus tatter leaf virus. to confirm previous results, we blast identified two asgv genomes in this study against the ncbi nucleotide database. the blast results confirmed that ctlv is closely grouped with asgv isolates in the genus capillovirus. from the genbank, we retrieved all asgv-associated sequences as well as ctlv-associated sequences. after removing partial sequences, we collected a total of genomes of asgv and ctlv isolates, including two asgv isolates in this study. the host ranges of ctlv were mostly from citrus species as well as lilium species. pear black necrotic leaf spot virus isolated from pear was an isolate of asgv according to the annotation in genbank. most asgv isolates were isolated from apple and pear, and some isolates, such as asgv isolates matsuco and li- were identified from citrus tamurana and lily, respectively. to reveal phylogenetic relationships, we aligned genome sequences displaying high sequence similarity of asgv and ctlv. sequence alignment and a phylogenetic tree using genome sequences of asgv and ctlv identified two largely divided clades. the first clade contained genomes, while the second clade included only pbnlsv and asgv isolate kfp. the first clade could be further divided into three groups. group a consisted of six ctlv isolates and a single asgv isolate, while group b contained only asgv isolates. group c was the largest, including seven asgv isolates and two ctlv isolates. recombination analysis for asgv isolates we analyzed recombination events among asgv isolates. the aligned genome sequences were subjected to the rdp program, which includes nine different algorithms for recombination detection. the rdp program detected a total of recombination events. of them, we selected six recombination events that were supported by at least five recombination algorithms. for example, pbnlsv contains two recombination sequences from asgv isolate li- three isolates—asgv isolate chn, asgv isolate hh, and ctlv isolate mth—include recombination sequences from asgv isolate fuji in the 5′ region. recombination events and were supported by seven algorithms. the major parent of the recombinant sequence for asgv isolate ytg was asgv isolate fuji. the major parent of the recombinant sequences for the three isolates—asgv isolate hh, asgv isolate chn, and ctlv isolate mth—were asgv isolate li- and ctlv isolate lily.fig. identification of recombination events by rdp program. a the positions of identified recombinants were depicted with the respective names of parental sequences. the individual genome of the asgv isolate was indicated by a different colored bar. the identified recombination events, including number and number, were rechecked by the rdp program and visualized by plot data with pairwise identity information. detailed information on the identified recombination events is provided in table 4 event minor parent parent contributing smaller fraction of sequence, major parent parent contributing larger fraction of sequence, unknown only one parent and a recombinant need be in the alignment for a recombination event to be detectable. the sequence listed as unknown was used to infer the existence of a missing parental sequence, ns no significant p-value was recorded for this recombination event using this method athe actual breakpoint position is undetermined discussion the rapid development of ngs is enabling virologists to find viruses from numerous species. ngs-based approaches have identified not only known viruses but also novel viruses. in fact, many horticultural plants are frequently infected by viruses and viroids. in particular, fruit trees usually propagated by grafting and cuttage are reservoirs of various plant viruses and viroids. in addition, the big data produced by ngs techniques has prompted virus identification in silico. here, we discussed the library types, sequencing methods, and de novo assembler for virus identification and viral genome assembly. the majority of plant viruses are composed of rna genomes, and dna viruses also replicate via an rna intermediate. thus, rna-based transcriptome libraries are preferable to dna-based genome libraries for virus identification. in the current study, we used published plant transcriptome data. to enrich viral rnas, ribosome-deleted libraries are usually prepared using extracted total rnas from virus-infected samples. however, we demonstrated that the mrna libraries using oligo d were successfully applied for virus identification. of course, rna viruses with poly tails such as asgv are also easily identified by mrna libraries. similarly, several polyadenylated rna viruses have been identified from sweet potato transcriptomes. several recent studies have also demonstrated that ribosome-deleted rna libraries as well as plant mrna libraries are suitable for the identification of viruses without poly tails or viroids. therefore, it might be ideal to use ribosome-deleted libraries for studies only focused on viruses. in the case of studies of both viruses and host plants, mrna libraries can be usefully applied. in this study, we used data from two different library types, including mrna and srna libraries that were single-end sequenced by the hiseq system. according to many recent studies, viral genomes have been de novo assembled from mrna as well as srna data. in our study, we assembled nearly complete genomes of two asgv isolates from the mrna data; however, the srna data could cover only 30 % of the asgv genome. we compared the numbers of sequencing reads between the mrna and srna data. however, the numbers of sequence reads between mrna and srna were very similar, indicating that the sequencing amount is not an important factor for viral genome assembly. in fact, when the number of sequencing reads is increased, the number of viral-associated reads is increased. therefore, the quantity of the sequenced data might play an important role in de novo genome assembly. the number of pear transcriptomes was about ten times that of apple transcriptomes. the sequence reads associated with asgv were viral reads out of, reads for the apple sample and viral reads out of, reads for the pear sample. although the number of total sequence reads in the apple sample was much smaller than that in the pear sample, the number of sequence reads associated with asgv was about times higher. this result suggests the amount of viral replication in the host might be also an important factor in de novo viral genome assembly. the portion of viral nucleic acids in the sample infected by virus is often low suggesting enrichment of virions prior to ngs. for example, purification of double-stranded rnas from the prunus species followed by pyrosequencing enabled to assemble four complete genomes of asian prunus virus, apv and apv. this study demonstrated successful application of dsrna purification for virus genome assembly using ngs technique. in the case of srna, two different types were prepared based on size fraction. the libraries without size fraction contain a large number of asgv-associated reads, but the libraries with size fraction contain very few reads associated with asgv. of course, the srna libraries were targeted for the identification of viral srnas. we suppose that the small number of srnas might be related to the ability of the rna silencing machinery in the host. in any case, a sufficient number of viral-associated reads is necessary for viral de novo genome assembly. in addition, sequencing methods are important for virus identification and viral genome assembly. in this study, all transcriptome data were single-end sequenced by hiseq as compared to single-end sequencing, paired-end sequencing provides sequences from both ends of a fragment and generates high-quality and alignable sequence data. the advantages of paired-end sequencing have been previously reported. thus, paired-end sequencing was far superior for the identification and genome assembly of the target virus. for virus identification, assembled contigs are frequently used. therefore, the choice of de novo assembler affects the quality and quantity of virus identification. for instance, mrna data were very efficiently assembled by trinity; however, few and low-quality contigs were assembled from the srna data by trinity. our comparative studies between the two de novo assemblers suggest trinity and velvet for de novo assembly of mrna data and srna data, respectively. the obtained viral contigs assembled by trinity from mrna data were low in number but long in length, while the viral contigs assembled by velvet were high in number but short in length. for the de novo assembly of a target virus with high-quality mrna data, trinity is ideal. velvet cannot assemble a nearly complete viral genome, but it assembled many contigs, which enabled us to identify additional viruses, for example, viruses in the pear transcriptomes. recently, several programs iva, price, and vicuna for de novo assembly of rna virus genome have been developed. the choice of optimal de novo assembler might be dependent on researchers and purposes. it is well known that rna viruses have a quasispecies nature within the host. however, to date, most studies have shown the variants and mutation rates of target viruses using cloning-based sanger sequencing methods. in this study, we successfully demonstrated the usefulness of plant transcriptome data for revealing the snvs of asgv. in fact, it is quite difficult to find virus variants using transcriptome data, while cloning-based sequencing methods might reveal variants. however, the cloning-based approaches require a rt-pcr amplification procedure to amplify full-length viral genomes. practically, the amplification of full-length viral genomes is not easy even though plant viruses are relatively small. we showed the presence of asgv variants in the transcriptome by comparing the asgv genome from the cultivar fuji derived from the sanger-sequencing method and de novo assembly. we did not judge which asgv genome was the dominant asgv genome; however, it is highly likely that the de novo-assembled asgv was a consensus genome sequence of asgv. the mutation rates of identified asgv genomes were varied:  % in the fuji, 1 % in the gd, and  % in the cuiguan. we suppose that several factors—including hosts, viral replication, and environmental cues—might affect the mutation rates. the association of viral mutation rate with other factors will be an interesting subject for further study.  taken together, our study showed the successful application of plant transcriptome data for virus identification, viral genome assembly, and viral mutation rates. in addition, we discussed several factors, including library preparation, ngs systems, de novo assemblers, and sample conditions for virus identification and genome assembly. 
a phylogenetic tree revealed that asgv and citrus tatter leaf virus are the same virus, which was divided into two clades. in addition, we compared the usability of two de novo assemblers, trinity and velvet, for virus identification and genome assembly. taken together, our in silico analyses using ngs data provide a successful application of plant transcriptomes to reveal extensive information associated with viral genome assembly, snvs, phylogenetic relationships, and genetic recombination. apple stem grooving virusde novo genome assemblyrecombinationrna-seqsingle nucleotide variationtranscriptomehttp://dx.doi.org/ /501100003627rural development administrationpj01186102cho won kyong issue-copyright-statement© the author 2016 we de novo assembled nearly complete genomes of asgv isolates fuji and cuiguan from apple and pear transcriptomes, respectively, and identified single nucleotide variations of asgv within the transcriptomes. we demonstrated the application of ngs raw data to confirm viral infections in the plant transcriptomes. recombination analyses identified six recombination events from viral genomes. in this study, we performed integrated analyses using plant transcriptome data for plant virus identification using apple stem grooving virus as an exemplar virus.
 next-generation sequencing provides many possibilities for plant virology research. in this study, we performed integrated analyses using plant transcriptome data for plant virus identification using apple stem grooving virus as an exemplar virus. we used publicly available transcriptome libraries from three different studies, two mrna-seq studies and a small rna-seq study.  we de novo assembled nearly complete genomes of asgv isolates fuji and cuiguan from apple and pear transcriptomes, respectively, and identified single nucleotide variations of asgv within the transcriptomes. we demonstrated the application of ngs raw data to confirm viral infections in the plant transcriptomes. in addition, we compared the usability of two de novo assemblers, trinity and velvet, for virus identification and genome assembly. a phylogenetic tree revealed that asgv and citrus tatter leaf virus are the same virus, which was divided into two clades. recombination analyses identified six recombination events from viral genomes.  taken together, our in silico analyses using ngs data provide a successful application of plant transcriptomes to reveal extensive information associated with viral genome assembly, snvs, phylogenetic relationships, and genetic recombination. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords apple stem grooving virusde novo genome assemblyrecombinationrna-seqsingle nucleotide variationtranscriptomehttp://dx.doi.org/ /501100003627rural development administrationpj01186102cho won kyong issue-copyright-statement© the author 2016 