Transcriptome Analysis of the Small Brown Planthopper, Laodelphax striatellus Carrying Rice stripe virus.

Rice stripe virus (RSV), the type member of the genus Tenuivirus, transmits by the feeding behavior of small brown planthopper (SBPH), Laodelphax striatellus. To investigate the interactions between the virus and vector insect, total RNA was extracted from RSV-viruliferous SBPH (RVLS) and non-viruliferous SBPH (NVLS) adults to construct expressed sequence tag databases for comparative transcriptome analysis. Over 30 million bases were sequenced by 454 pyrosequencing to construct 1,538 and 953 of isotigs from the mRNA of RVLS and NVLS, respectively. The gene ontology (GO) analysis demonstrated that both libraries have similar GO structures, however, the gene expression pattern analysis revealed that 17.8% and 16.8% of isotigs were up- and down-regulated significantly in the RVLS, respectively. These RSV-dependently regulated genes possibly have important roles in the physiology of SBPH, transmission of RSV, and RSV and SBPH interaction.

Rice stripe virus (RSV), the prototype virus of genus tenuivirus, is a non-enveloped RNA virus which forms circular filamentous shape of ribonucleoprotein particles with four species of segmented single-stranded RNA genome of RNA1, 2, 3 and 4 according to the size of each RNA molecule (Hibino et al., 1985; Ishikawa et al., 1989). The largest segment RNA1 encodes a 336.8 kDa protein which is thought to be a RNA dependent RNA polymerase, RdRP (Toriyama et al., 1994). RNA2 encodes 94 kDa glycoprotein and NS2 (Takahashi et al., 1993), RNA3 encodes NS3 and CP (Zhu et al., 1991), and RNA4 encodes SP and NS4 (Zhu et al., 1992). RNA1 is a negative-sense single strand while RNA2, 3, and 4 are ambisense in their coding strategy (Hamamatsu et al., 1993). RSV is one of the most serious viral rice pathogens which occur in South Korea, China, and Japan. The typical symptoms of RSV infected rice plant are chlorotic stripes or necrotic streaks on leaves, and premature wilting. Some of the leaves emerge abnormally unfolded, and then elongate to become twisted and droop. The RSV often causes significant reduction in yield. Especially the rice plants inoculated with RSV in 11-leaf stage or earlier produce no grain (Hibino, 1996; Ling, 1972).

The mechanisms of the insect vector mediated plant virus transmission are diverse based on the differences of the duration of which the vector insects can harbor infectious viral particles. The transmission mechanisms can be classified into three categories; non-persistent (harboring viruses for a few seconds/minutes), semi-persistent (a few hours/days, but lost upon molting), and persistent (often throughout the lifespan of the vectors) (Ng and Falk, 2006). Most of the plant viruses are dependent on the vector insects for plant-to-plant spread, and one of the most important vector insects are hemipterans, which can transmit viruses during feeding with their piercing-sucking mouthpart in the phloem (Hogenhout et al., 2008). The small brown planthopper (SBPH), Laodelphax striatellus Fallén (Hemiptera: Delphacidae), is one of the vectors of RSV and classified as a vector of persistent transmission because RSV replicates persistently in the SBPH ovary to transmit virus to the offspring transovarially (Toriyama, 1986) as well as circulatively to proliferate and escape from the midgut, salivary gland and ovary barriers for host-to-host infection (Nault and Ammar, 1989). Many of the tenuiviruses rely on planthopper vectors for transmission, and the vector specificity is high (Falk and Tsai, 1998). The SBPH is one of the most important insect vectors of RSV, however, our knowledge of the SBPH biology is limited and most of the interactions between RSV and SBPH have received less attention so far.

To investigate the RSV-SBPH interaction, we employed next generation 454-pyrosequencing technology to construct expressed sequence tag (EST) libraries from SBPH because SBPH genome sequence is not available at this moment. The pyrosequencing technology (Ronaghi et al., 1996) on the Roche 454 GS-FLX platform was used to create EST databases of RSV-viruliferous SBPH (RVLS) and non-viruliferous SBPH (NVLS), and the gene expression pattern of RVLS and NVLS were compared based on the read counts of each isotigs. We also conducted quantitative real-time PCR (qPCR) experiments with selected genes from the EST database to validate the sequencing results.

RSV-viruliferous and non-viruliferous Laodelphax striatellus.

The RVLS and NVLS were collected from the RSV infected rice field and healthy rice seedlings, respectively. The RVLS and NVLS colonies were reared in the laboratory on 2 – 3 cm tall rice seedlings in glass vessels under a light:dark cycle of 16:8 h at 25°C and transferred to fresh seedlings every 10 – 14 days to assure sufficient nutrition. After a few generations, total RNA for nucleotide sequencing was extracted from RVLS and NVLS with TRI-Reagent (Invitrogen, USA) according to the manufacturer’s instruction, and RT-PCR was performed to detect RSV in the RVLS and NVLS. Two pairs of primers for specific amplification of RSV sequences were designed as follows; RSV-CP-F: 5′-GCT CAT AAA GCT GAT GTT G-3′ and RSV-CP-R: 5′-GTC TTC ACT TTC CCA TTG-3′ correspond to the nucleotide sequence of coat protein gene in RNA3 with 969 base pair (bp) of expected PCR product. RSV-NCP-F: 5′-GAA AGA TCC AGA GAG AGT CA-3′ and RSV-NCP-R: 5′-CAA AAT AAG TTG GCA TCT TC-3′ correspond to the nucleotide sequence of non-capsid protein gene in RNA4 with 537 bp of expected PCR product. Another two pairs of primers for specific amplification of SBPH sequences were designed as follows; CCOS-F: 5′-ATT GCT CAT AAA GCT GAT GTT G-3′ and CCOS-R primer: 5′-CCG TAC CAG GAC GAA TTA-3′ correspond to the nucleotide sequence of cytochrome c oxidase subunit II with 684 bp of expected PCR product. ATP-F: 5′-CCC TCA AGT CAC CTC TCT-3′ and ATP-R: 5′-ATG GAA TTT CTG AGG AA-3′ correspond to the nucleotide sequence of ATP synthase F0 subunit 6 with 655 bp of expected PCR product. cDNAs were synthesized with the reverse primer of each genes, and RT-PCR reactions were performed with AccuPower one-step RT-PCR premix (Bioneer, Korea). The RT-PCR result confirmed that only the RNA samples from the RVLS amplified viral sequences while both RVLS and NVLS samples amplified SBPH specific sequences (data not shown).

Transcriptome sequencing and EST library construction.

To conduct the transcriptome sequencing of RVLS and NVLS, total RNA samples were extracted with TRI-Reagent from approximately 500 SBPHs of RVLS and NVLS colonies, respectively. 200 ng of mRNA samples were eluted from the total RNA samples by using oligo-dT hybridization followed by heat fragmentation, and reverse transcription to synthesize cDNA with random hexamer. Second strand was synthesized with second strand enzyme followed by 2 hour incubation in buffer for recovery of double stranded cDNA, 5′-phosphorylation and addition of 3′-dA overhang. Finally, adaptor was ligated and the samples were amplified to be prepared for sequencing. Approximately 5 ug of each amplified cDNA samples from RVLS and NVLS were subjected to two of 1/8 plate sequencing runs with a Roche 454-GS-FLX Titanium sequencer. The sequence data were processed with GS-FLX software v2.5.3 to remove low quality sequences, and the resulting sequences were then screened against the NCBI UniVec database and Escherichia coli genome sequences to remove contaminant sequences, then the sequences shorter than 100 bp were discarded. The processed sequences were assembled into isotigs using Newbler v 2.5.3 software to construct two EST libraries of RVLS and NVLS. Both data sets of assembled contigs and singletons from RVLS and NVLS were pooled to construct a combined sequence dataset which was supposed to have more and longer SBPH contigs. The combined dataset was used to build a new reference EST library to compare with RVLS and NVLS, and to analyze differential gene expression pattern. The features of the three EST libraries are summarized in Table 1, and the EST libraries in FASTA format are available upon request. The sequence data of RVLS and NVLS were deposited in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) with accession numbers of SRR538840 and SRR538841, respectively.

Table 1

Summary of the features of Laodelphax striatellus EST libraries

	RSV-viruliferous L. strialellus	Non-viruliferous L. strialellus	Combined samples
Total bases (bp)	30,071,326	30,566,488	60,637,814
Number of reads	90,147	90,009	180,156
Average reads length (bp)	336	342	342
Number of isotigs	1,538	952	2,555
Average isotig length (bp)	720	699	745
Range of isotig length (bp)	15∼8,596	69∼9,265	82∼9,179
Number of singletons	25,855	14,733	33,808

Gene ontology analysis.

To obtain an overview of the transcriptome of SBPH, the isotigs of RVLS and NVLS were subjected to GenBank nr database with BLASTX and gene ontology (GO) analysis to annotate the isotigs into three main categories; biological process, cellular component, and molecular function (Ashburner et al., 2000). The GO distributions of the isotigs and singletons of both RVLS and NVLS were similar to each other (Table 2), and the detailed GO annotations of each subcategory also showed similar trends of distributions. This result implies that the RSV did not alter the gene expression attributes of SBPH significantly. The 2,555 isotigs of the combined dataset were analyzed by comparing with Drosophila proteins using BLASTX with a cut-off value of 1e-05 to identify their gene functions. Total of 1,467 isotigs of the BLASTX hits were subjected to Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology analysis (Kanehisa et al., 2004) to categorize their functional classes (Fig. 1A).

Table 2

BlastX search and Gene Ontology analysis of Laodelphax striatellus EST libraries

	RSV-viruliferous L. strialellus		Non-viruliferous L. strialellus
	Isotig	Singleton	Isotig	singleton
Total	1,538	25,855	952	14,733
NR BLASTX hit (e < 10−3)	927 (60.3%)	9,846 (61.9%)	595 (62.5%)	5,690 (38.6%)
GO hit	829 (53.9%)	8,232 (31.8%)	521 (54.7%)	4,684 (31.8%)
– Biological process	295 (19.2%)	2,800 (10.8%)	182 (19.1%)	1,629 (11.1%)
– Cellular component	203 (13.2%)	2,137 (8.3%)	137 (14.4%)	1,182 (8%)
– Molecular function	331 (21.5%)	3,295 (12.7%)	202 (21.2%)	1,873 (12.7%)
GO no-hit	709 (46.1%)	17,623 (68.2%)	431 (45.3%)	10,049 (68.2%)

Fig. 1

KEGG analysis of total (A), and UP and DOWN (B) isotigs. (A) Functional groups of the combined EST database defined by KEGG orthology were classified as carbohydrate metabolism, energy metabolism, lipid metabolism, nucleotide metabolism, amino acid metabolism, biosynthesis of other secondary metabolites, xenobiotics biodegradation and metabolism, transcription, translation, folding sorting and degradation, replication and repair, membrane transport, signal transduction, cellular processes, and organismal systems. (B) The isotigs which showed significantly different transcription level based on the read number (a binominal probability of < 0.1) were classified as UP or DOWN. UP and DOWN indicate the genes up- and down-regulated in RVLS, respectively.

Comparison of gene expression.

To compare the gene expression pattern of RVLS and NVLS, the total read counts of SBPH isotigs from RVLS and NVLS were normalized arithmetically, and then total of the 2,543 SBPH isotigs from the combined EST library were plotted based on the normalized read counts in RVLS and NVLS sequencing results. The significant differences of the transcription levels were determined by comparing the read counts of each isotig (binominal probability of < 0.1) to identify the isotigs which were up- or down-regulated by RSV. The distribution of isotigs showed that majority of the host genes were not significantly affected while 453 (17.8%) and 428 (16.8%) of isotigs were up- and down-regulated in RVLS, respectively (Fig. 2). The isotigs which showed significant differences of the transcription levels were also subjected to KEGG orthology anaylsis to show their functional classes. The isotigs related to folding, sorting and degradation, and cellular processes were up-regulated while the genes related to translation and xenobiotics metabolism were down-regulated compared with those of the NVLS. In particular, the genes related to translation were significantly down-regulated (Fig. 1B). Especially, the numbers of isotigs which showed more than ten normalized read count in either RVLS or NVLS with more than 10-fold of up-and down-regulation in RVLS were 234 and 308, respectively. These differentially expressed isotigs were subjected to GenBank insect nr database with BLASTX, and the total of 23 up-regulated and 22 down-regulated isotigs of BLASTX hits were summarized in Tables 4 and 5, respectively.

Fig. 2

Graph of the read counts of RVLS and NVLS for each isotig. The sequence read counts of RVLS and NVLS of each isotig were graphed on an x–y plot. For convenience, isotigs were ploted on two separated graphs. For the contigs which have less than 200 read counts are shown in the left panel, and the isotigs with more than 200 read counts were graphed in the write panel. Total of 2,543 SBPH isotigs were subjected to this plot, and 1,662 of the isotigs were not significantly affected (p > 0.1) while 453 and 428 of isotigs were up-and down-regulated, respectively. Up- and down-regulated isotigs are indicated as diamonds and squares, respectively.

Table 4

The genes down-regulated more than 10-fold in RSV viruliferous L. striatellus

Isotig No.	Ratio (NVLS/RVLS)	Accession No.	Gene description [species] (%identity)
isotig00992	RVLS=0	XP_001601896.1	Cleavage stimulation factor 64-kDa subunit [Nasonia vitripennis] (82%)
isotig01261	RVLS=0	XP_623653.3	Hypothetical protein LOC551259 [Apis mellifera] (68%)
isotig01263	RVLS=0	XP_001607475.1	Cop9 signalosome complex subunit [Nasonia vitripennis] (86%)
isotig01449	RVLS=0	XP_396542.4	Serrate RNA effector molecule homolog [Apis mellifera] (68%)
isotig01458	RVLS=0	EFZ09607.1	Hypothetical protein SINV_05505 [Solenopsis invicta] (52%)
isotig01601	RVLS=0	ADR73026.1	Putative acetylcholinesterase 1 [Laodelphax striatellus] (33%)
isotig01730	RVLS=0	NP_001155554.1	ACYPI003981 [Acyrthosiphon pisum] (89%)
isotig01842	RVLS=0	XP_002427287.1	Peroxisomal membrane protein 70 abcd3 [Pediculus humanus] (89%)
isotig01844	RVLS=0	EFY84384.1	Glyceraldehyde-3-phosphate dehydrogenase [Metarhizium acridum] (91%)
isotig01902	RVLS=0	ACY71296.1	Ribosomal protein L13A [Chrysomela tremulae] (76%)
isotig01972	RVLS=0	Q7Z1B8.1	Sec61p gamma subunit [Gryllotalpa orientalis] (94%)
isotig01982	RVLS=0	EFR24025.1	Hypothetical protein AND_11693 [Anopheles darlingi] (88%)
isotig01075	23.04	XP_967420.1	Hypothetical protein TcasGA2 [Tribolium castaneum] (92%)
isotig00365	15.02	EFN76208.1	Sodium/potassium-transporting ATPase subunit beta-2 [Harpegnathos saltator] (76%)
isotig01096	14.02	XP_001121309.1	Cullin-2 [Apis mellifera] (87%)
isotig01020	13.02	XP_002424439.1	26S proteasome non-ATPase regulatory subunit [Pediculus humanus] (94%)
isotig01563	12.02	XP_002430165.1	Eukaryotic translation initiation factor [Pediculus humanus] (87%)
isotig01265	11.02	XP_002434356.1	Retinal protein [Ixodes scapularis] (75%)
isotig01019	11.02	XP_002430627.1	D-beta-hydroxybutyrate dehydrogenase [Pediculus humanus] (64%)
isotig02418	10.02	XP_001983864.1	GH15337 [Drosophila grimshawi] (63%)
isotig01163	10.02	XP_974479.1	Hypothetical protein TcasGA2 [Tribolium castaneum] (61%)
isotig02288	10.02	AAP88325.1	Cytochrome c oxidase polypeptide VIII [Drosophila mauritiana] (45%)

Table 5

Primers used for quantitative real-time PCR expression profiling of putative transcripts generated from 454-pyrosequencing

Targeted contigs	Sequence of the primer pairs
Vitellogenin	Fw: 5′ GGACGAATACTGAAGCACACTTGG 3′Re: 5′ CAGCAGCAGCAGCAGCAAC 3′
Heat shock protein 1	Fw: 5′ GAAGAGAGCGTTTATTGAGTGTTGTG 3′Re: 5′ CTCCTAAATGGGCGATGACTATGATG 3′
Heat shock protein 2	Fw: 5′ GTCCTCCTCGTCAGCCTTGTC 3′Re: 5 GAAGAAGCACTCGCAGTTTATCGG 3′
Translation initiation factor 1	Fw: 5′ GCAAGGAAATCACGCCATACTGTC 3′Re: 5′ CAATGACTGGAGTGTGTTCTGTGG 3′
Translation initiation factor 2	Fw: 5′ TGTGTACTTGGTAGCCGATGACTG 3′Re: 5′ CGATGCCGAGAAGCGAAACG 3′
Translation initiation factor 3	Fw: 5′ TGCGTGAGATGTCGTTTGATGTG 3′Re: 5′ CGAGAGGTTTGGAAAGTGTGGTG 3′
Ribosomal protein 1	Fw: 5′ TGTCACAAAGCGTGCCAAGAAAC 3′Re: 5′ TGGTAAGGTCGGTCCTGTAATGC 3′
Ribosomal protein 2	Fw: 5′ CGGGCGTGATGAATGGATTTGG 3′Re: 5′ GCTCTCCTGCGTCGTCTCG 3′
Ribosomal protein 3	Fw: 5′ GTCACAAACGGAAGTATCGCTCAC 3′Re: 5′ ATCGGCTCGTCACAGTTCTCG 3′
Ribosomal protein 4	Fw: 5′ GTAGAAGCGAAACAGCCCAACTC 3′Re: 5′ GTGACCCTTCCTACCGAATCCG 3′

Verification of 454 pyrosequencing results using qPCR.

To validate the pyrosequencing results of the SBPH transcriptomes, the transcription levels measured based on the read counts were compared with those of qPCR. For qPCR, total of seven translation-related isotigs (translational initiation factor 1, 2, and 3, and ribosomal protein 1, 2, 3, and 4), which belong to the KEGG functional group of which showed the most significant difference of the transcription level between RVLS and NVLS, were selected. Also heat shock protein 1 and 2, and vitellogenin genes were selected for qPCR along with elongation factor 2 (EF2) gene as a reference for normalization (Tong et al., 2009). The cDNA was synthesized from the total RNA of RVLS and NVLS using QuantiTect Reverse Transcription Kit (QIAGEN, Germany) according to the manufacturer’s instructions, followed by qPCR conducted using EvaGreen qPCR Mastermix (Applied Biological Materials Inc, Canada) and CFX96™ Real-Time System (BIO-RAD, USA). The primers used for the qPCR are listed in Table 5. The qPCR results confirmed that the ten isotigs were indeed transcribed in both RVLS and NVLS, and the calculation of the relative transcription levels of the isotigs using 2−ΔCt method (Pfaffl, 2001) demonstrated that the seven genes (vitellogenin, heat shock protein 1 and 2, and ribosomal protein 1, 2, 3 and 4) and the other three genes (translation initiation factor 1, 2 and 3) were up- and down-regulated in RVLS, respectively. To calculate fold differences of the isotigs, the read counts of the isotigs in RVLS and NVLS which are the corresponding isotigs of the qPCR experiment were measured in reads per kilobase of exon per million mapped sequence reads (RPKM) (Mortazavi et al., 2008). The sequence reads of the five isotigs (heat shock protein 1 and 2, translation initiation factor 1, 2 and 3) were confirmed in both RVLS and NVLS while the other five isotigs (vitellogenin, ribosomal protein 1, 2, 3 and 4) were confirmed in RVLS exclusively. The expression levels measured by qPCR were in a good accordance with the ratios of read counts of the EST libraries (Fig. 3).

Fig. 3

Validation of the 454 pyrosequencing results by qPCR. The expression profiles of the ten isotigs were analyzed by qPCR to validate the 454 pyrosequencing results. The relative transcription levels are expressed in fold differences of RVLS to NVLS calculated by qPCR using 2−ΔCt method. Vitellogenin, heat shock protein 1 and 2, and ribosomal protein 1, 2, 3 and 4 were up-regulated while translation initiation factor 1, 2 and 3 were down-regulated in RVLS. The transcription level calculated with sequence read counts of the five of the ten genes, heat shock protein 1 and 2, translation initiation factor 1, 2 and 3, demonstrated that the transcription levels measured by genomics and experimental methods were in a good agreement. The other five genes, vitellogenin, ribosomal protein 1, 2, 3 and 4, which showed substantial read counts in RVLS, had no read counts in NVLS so that these genes were not able to calculate the fold differences (indicated with asterisk). NGS and qPCR indicate the fold differences calculated with the sequencing and qPCR results, respectively.

Rice stripe virus transcripts in viruliferous L. striatellus.

The RSV transcripts of three of the seven viral ORFs, NS3, CP and RdRP were confirmed in RVLS, and their sequence read counts were 293, 224 and 5, respectively. The total RSV read counts were 522, and this was 0.63% of the total of 83,041 valid sequence reads of RVLS contigs and singletons. The low read count of RdRP suggests that the replication of RSV is not highly active in SBPH while NS3 was the most abundant RSV transcript in RVLS which codes for viral suppressor of siRNA (Hemmes et al., 2009). CP was the second most abundant RSV transcript which codes for capsid protein (Hayano et al., 1990). Both NS3 and CP are the genes located on genomic and antigenomic strand of RSV RNA3, respectively. Therefore the similar levels of NS3 and CP transcripts possibly suggest that the genes in RNA3 are under the same gene expression regulation mechanism, however, qPCR experiments to quantify the RSV transcripts in RVLS need to be performed to prove the hypothesis.

The 454 pyrosequencing results demonstrated that the next generation sequencing (NGS) technique was able to reflect the characteristics of the transcriptome of RVLS. The genes in the RVLS which showed different transcription patterns indicated that the RSV may alter the gene expression of the SBPH which possibly affects the physiology, fitness and behavior of RVLS. The up-regulated heat shock protein genes might be one of the examples that the RSV may take the advantage of altering gene transcription to improve fitness of the vector because heat shock proteins play a significant role in generating tolerance to stresses (Krebs, 1999; Sorensen et al., 2003). It is also interesting to note that the transcription level of vitellogenin in RVLS was highly up-regulated, which means RSV might modulate the host vitellogenin pathway genes to promote the oogenesis or alter sex ratio of RVLS to facilitate its persistent transovarial transmission via the offsprings of RVLS. The highly up-regulated ribosomal protein genes were also the possible consequence of the up-regulation of vitellogenin because the ovarian cells need more amount of translational machinery, however, the down-regulated translation initiation factors are paradoxical because the recruitment of transcription initiation factors is a crucial step of protein synthesis. Many viruses can inhibit host translation initiation factors to take over host gene expression by which are usually viral protein mediated protein level interactions, rather than transcription regulation (Walsh and Mohr, 2011). It suggests that RSV possibly can alter the transcription of the genes for cellular translational machinery direct or indirect manner. Another example of viral modulation of cell physiology is the major viral transcript NS3, a RNAi suppressor which binds 21 nt siRNA, and it might be a viral strategy to evade defense system of SBPH and rice plant (Alvarado and Scholthof, 2009; Hemmes et al., 2007). Early studies on vector transmission of plant viruses demonstrated the specificity of the virus-vector interaction, which says that the instances of the biological properties associated with the transmission process may differ significantly (Ng and Falk, 2006). However, the genomics approach to study the virus-vector interaction using NGS technique might contribute to investigate the gene regulation mechanism in the persistent vector insects and the upstream regulation factors, which will be a new approach to study the interaction among the virus, plant, and vector.

Table 3

The genes up-regulated more than 10-fold in RSV viruliferous Laodelphax striatellus

Isotig No.	Ratio (NVLS/RVLS)	Accession No.	Gene description [species] (%identity)
isotig01452	NVLS=0	AEC04824.1	Putative mucin-like protein [Laodelphax striatellus] (100%)
isotig01908	NVLS=0	BAD00045.1	armadillo, mRNA [Culex quinquefasciatus] (98%)
isotig02263	NVLS=0	XP_002425481.1	Chitinase [Apis mellifera] (97%)
isotig00457	NVLS=0	XP_001868736.1	Peptide methionine sulfoxide reductase [Culex quinquefasciatus] (56%)
isotig00475	NVLS=0	BAF75351.1	Vitellogenin [Nilaparvata lugens] (69%)
isotig00928	NVLS=0	EFN70829.1	E3 ubiquitin-protein ligase MARCH6 [Camponotus floridanus] (78%)
isotig00942	NVLS=0	XP_001945803.2	Seryl-tRNA synthetase, cytoplasmic-like [Acyrthosiphon pisum] (79%)
isotig00979	NVLS=0	ABN11925.1	Putative vacuolar ATP synthase subunit B [Maconellicoccus hirsutus] (98%)
isotig01149	NVLS=0	EFN62899.1	Vesicle-associated membrane protein 2 [Camponotus floridanus] (91%)
isotig01191	NVLS=0	EFY90146.1	Cell surface protein, putative [Metarhizium acridum] (85%)
isotig01302	NVLS=0	XP_001948516.2	Serine/threonine-protein kinase 24-like [Acyrthosiphon pisum] (97%)
isotig01406	NVLS=0	XP_002427193.1	Pre-mRNA-processing-splicing factor, putative [Pediculus humanus] (93%)
isotig01764	NVLS=0	AAP59456.1	Cathepsin B precursor [Araneus ventricosus] (88%)
isotig01775	NVLS=0	XP_001859032.1	Importin beta-3 [Culex quinquefasciatus] (88%)
isotig01955	NVLS=0	ADD20429.1	Putative BBOX Zn-finger protein [Glossina morsitans] (81%)
isotig01390	0.02	XP_001600318.1	Zinc finger protein [Nasonia vitripennis] (83%)
isotig01822	0.03	EFN84196.1	Pancreatic lipase-related protein 2 [Harpegnathos saltator] (63%)
isotig01477	0.04	XP_002426267.1	60S ribosomal protein L17, putative [Pediculus humanus] (89%)
isotig00827	0.04	EFZ13522.1	Hypothetical protein SINV_04974 [Solenopsis invicta] (85%)
isotig00861	0.05	XP_001605359.1	Protein disulfide isomerase [Nasonia vitripennis] (67%)
isotig00784	0.06	XP_001968815.1	Hel25E [Drosophila yakuba] (97%)
isotig01188	0.08	EFN87947.1	26S protease regulatory subunit 8 [Harpegnathos saltator] (98%)
isotig00802	0.08	P32122.1	Arrestin homolog [Locusta migratoria] (79%)