# Subtracted Transcriptome Profile of Tiger Shrimp (Penaeus monodon) That Survived WSSV Challenge

Authors: Benedict A. Maralit; Minerva Fatimae H. Ventolero; Mary Beth B. Maningas; Edgar C. Amar; Mudjekeewis D. Santos
Journal: Dataset Papers in Science
 (2014)

---

## Abstract

There is increased interest in the development of virus-resistant or improved shrimp stock because production is currently hindered by outbreaks and limited understanding of shrimp defense. Recent advancement now allows for high-throughput molecular studies on shrimp immunity. We used next-generation sequencing (NGS) coupled with suppression subtractive hybridization (SSH) to generate a transcriptome database of genes from tiger shrimp that survived White spot syndrome virus (WSSV) challenge. A total of 9,597 unique sequences were uploaded to NCBI Sequence Read Archive with accession number SRR577080. Sixty-five unique sequences, 6% of the total, were homologous to genes ofPenaeus monodon. Genes that were initially related to bacterial infection and environmental stress such as 14-3-3 gene, heat shock protein 90, and calreticulin were also found including a few full-length gene sequences. Initial analysis of the expression of some genes was done. Hemocyanin, ferritin, and fortilin-binding protein exhibited differential expression between survivor and control tiger shrimps. Furthermore, candidate microsatellite markers for brood stock selection were mined and tested. Four trinucleotide and one dinucleotide microsatellites were successfully amplified. The study highlights the advantage of the NGS platform coupled with SSH in terms of gene discovery and marker generation.


---

## Body

## 1. Introduction

Although a lot of best management practices, detection methods, and disease intervention have been developed, shrimp production is still continuously affected by viral infections. Because of this, greater attention has now been placed on molecular studies for application on genetic breeding programs. As seen in recent years, there is an increase in the number of studies that aim to elucidate the function of several previously unknown genes and their viral response pathways. Expression and characterization studies of one or a group of similar genes are currently gaining pace. Likewise, large databases from which immune-related genes may be determined are now publicly available and are continuously generated.Some notable genes that have been recently described from Expressed Sequence Tags (EST) databases are hemocyanin, ferritin, and fortilin-binding protein [1–5]. Many subtraction libraries have also been used to associate mined genes from EST to specific roles in immunity. For instance, cathepsins [1, 6] and penaeidins [7–9] were proven to be immune related based on several separate cDNA library studies [10]. Using the same large-scale sequencing approach, genes currently related to bacterial infection and environmental stress, such as heat shock proteins [11], hemocyte homeostasis-associated protein [12], 14-3-3 gene, and calreticulin, are also suggested to play important roles in White spot syndrome virus (WSSV) infection or resistance.As seen, large-scale transcriptome analysis and EST databases have established many immune related genes, which were further characterized by expression or functional analyses [1, 8, 9, 12, 13]. Indeed, its usefulness has been proven in many aspects of shrimp health. Recent transcriptome approaches such as microarrays and next-generation sequencing (NGS) have also provided efficient applications for rapid determination of quantitative trait loci and genetic markers from immune-related genes [1, 14–16] because these approaches provide longer average sequence reads of 600 base pairs (bp) than traditional ones. The markers generated from these platforms may be used for a lot of genetic studies integral to taxonomy or aquaculture. In fact, the latest breeding technologies in plants make use of next-generation sequencing data for rapid marker development [17]. In shrimp, traditional selective breeding programs have already been useful for producing WSSV resistantPenaeus vannamei [1, 18–21]. However, this classical approach entails high-input requirements of capital, time, and research effort. As seen in plants, marker-assisted selection [22], a modified approach in traditional selective breeding, proves successful. It is, thus, the goal of this study to create a database of genes from which candidate markers may be mined and used in marker-assisted selection (MAS) of WSSV-resistant or tolerantP. monodon. Suppression subtractive hybridization (SSH) coupled with NGS was used in this study to generate the desired database. We show that this platform can be used to further study genes that may have relevant function in WSSV immunity. Likewise, we identified resistance-related markers such as microsatellites that may be amplified and analyzed for downstream genetic applications such as semiquantitative expression analysis.


## 2. Methodology

Sample collection was as follows. Shrimps for challenge experiments were collected from a farm located in Negros Occidental, Philippines. Shrimp fry used to stock the farm came from a hatchery that practices biosecurity monitoring. Before they were transported to a biosecure facility in SEAFDEC AQD, samples were tested for WSSV and other viruses. Subsequently, they were quarantined in this facility for one month as they were being acclimated to tank conditions.For WSSV challenge, anin vivo titration test by intramuscular injection of the viral inoculum was conducted to determine the dose that resulted in 50% mortality (SID50) of test shrimps (155 days old). To do this, gills from moribund shrimp were collected and homogenized in PBS. The homogenate was centrifuged at 5000 rpm. The resulting supernatant was filtered in 0.45 μm sterile membrane filter, serially diluted (10−3 to 10−9) in PBS, and 100 μL was injected into the abdominal muscle region of each shrimp. Mortality was observed for 15 days. The computed SID50 was 10−6.45 mL−0.1 and the titer of the undiluted supernatant was 107.45 SID50 mL−1. The challenge test was conducted by injecting each shrimp with 1 SID50 of the inoculum following the same protocol used in thein vivo titration test. Surviving shrimps from this challenge were kept for genetic analysis. Tissue sampling from surviving shrimps (survivor) was done after 30 days of challenge test termination.DNA extraction and WSSV detection were as follows. A total of 10 representativeP. monodon samples were collected. Five shrimps that survived viral challenge experiments (30 days after infection) were labeled as “survivor” while five apparently healthy shrimps (no gross signs of disease) were labeled as “control.” Extraction of DNA from gills was performed using CTAB extraction buffer as described by Santos et al. [23]. Spectrophotometry of DNA extracts was performed before proceeding to amplification (PCR) experiments. A total of 10 μL of reaction mixture was used containing 2 μL of DNA template, 0.2 mM dNTP, 2 μM of each primer [24, 25], and 1x buffer already containing 1.5 mM MgCl2 and amplified under the following conditions: initial denaturation at 95°C for 5 minutes, followed by 30 cycles of denaturation at 95°C for 30 seconds, annealing at 52°C for 30 seconds, and extension at 72°C for another 30 seconds and a final extension for 5 minutes and storage at 4°C. All DNA extracts were confirmed to contain DNA by PCR amplification using universal 16Sar primers [26]. Based on the PCR detection method described, none of the samples was positive for WSSV (Figure 1).Figure 1
PCR-based detection of White spot syndrome virus in allPenaeus monodon samples. Top: internal controls; bottom lane 1: ladder; lane 2: negative control; lane 3: positive control; lanes 4–8: survivor samples; lanes 9–13: control shrimp samples.Library preparation was as follows. Suppression subtractive hybridization (SSH) is one of the most powerful and popular methods for generating subtracted cDNA libraries because the probability of obtaining low-abundance differentially expressed cDNA fragments is very high. It is a suppression PCR technique that uses normalization and subtraction in a single procedure. It was used in this study because downstream analyses, such as filtering after next-generation sequencing of the subtracted library, are simplified in this platform.RNA extracted from gills using Nucleospin RNA II kit (Machery-Nagel) was subjected to spectrophotometry for the determination of purity and concentration. cDNA was synthesized from the RNA extracts using Clonetech SMARTer cDNA synthesis kit. Immediately after, cDNA was amplified using Clonetech Advantage II PCR kit through long distance-PCR (ld-PCR) following the recommended protocol. Incubation time was increased from 3 minutes to 6 minutes to ensure synthesis and amplification of longer cDNA transcripts. The cDNA templates with highest quality were used in the subtraction experiment.cDNA subtraction was carried out using Clonetech PCR-Select cDNA Subtraction Kit in accordance with the manufacturer’s protocol. For forward subtraction, cDNA synthesized from shrimp sample that survived viral challenge (survivor) was used as driver while cDNA from control shrimp was used as tester. Forward and reverse SSH reactions were performed simultaneously. Tester and driver cDNA populations were briefly digested with RsaI and ligated each with one of the two adaptors (adaptor 1 and adaptor 2R) included in the kit. After denaturation at 98°C for 90 s, tester cDNAs were ligated with the adaptor 1 and adaptor 2R separately and then hybridized with an excess of driver cDNAs at 68°C for 8 h. After the first hybridization, the two products (adaptor 1 and adaptor 2R) were mixed together and immediately hybridized again with an excess of freshly denatured driver cDNA for 20 h at 68°C. Differentially expressed products were then selectively amplified by two rounds of PCR using Clonetech Advantage II PCR kit with primers specific to the two adaptors, resulting in suppression of genes common to the two cDNA populations and enrichment of the differentially expressed genes. The efficiency of the subtraction was finally evaluated by comparing the abundance of the constitutively expressed-actin in the subtracted and unsubtracted products by PCR as described by Zhao et al. [1]. Figure 2 shows that, after subtraction, the constitutive abundance of β-actin was efficiently reduced after SSH.Figure 2
Evaluation of the subtraction efficiency of the subtracted cDNAs. PCRs withβ-actin-specific primers were performed using forward (lanes 2 and 3) and reverse (lanes 4 and 5) SSH products. Lane 1: ladder; lanes 2 and 4: unsubtracted; lanes 3 and 5: subtracted.Next-generation sequencing and bioinformatics were as follows. The forward subtracted product from SSH was sent to Macrogen Inc., Korea (http://www.macrogen.com/eng/), for 454 sequencing of 1/4 plate. De novo assembly of the acquired reads was conducted using GS De Novo Assembler v 2.6 (Newbler) with cDNA/transcriptome option. Newbler parameters such as base calling, quality calling, and assembly quality filters were set to default. Singletons were cleaned by SeqClean (http://sourceforge.net/projects/seqclean/) and further assessed and controlled by Lucy (http://lucy.sourceforge.net/) with minimum length filter of 100 bp. A complete workflow of the study is exhibited in Figure 3.Figure 3
Experiment workflow.As seen in Table2, a total of 76,889,928 bases, which comprised 240,616 reads, were generated from NGS. Approximately 52% of the reads were assembled to form 125,598 sequences with an average read length of 319 ranging from 100 to 600 base pairs. De novo assembly generated 4,090 contigs with N50 of 504. Likewise, 9,526 singletons were produced leaving out 10,108 outliers and 2,650 too short sequences. Two phases of sequence cleaning using SeqClean and Lucy generated 5,507 singletons, completing a total of 9,597. This number of unique sequences is comparable to those stored inPenaeus monodon Genome Database (http://sysbio.iis.sinica.edu.tw/page/index.php) and thePenaeus monodon gene discovery project [27]. The genome size of tiger shrimp is not yet fully determined but the 1C nuclear DNA content ofP. monodon was estimated to be ~72.2% of the human genome, that is, ~2.53 pg DNA per nucleus or 2.17×109 bp per haploid genome [28]. In Table 2, the first part presents Reads Status where the number of reads counts the reads used in the assembly computation; number of bases, the read’s bases used in the assembly computation; assembled, the number of reads that are fully incorporated into the assembly; partial, the number of reads with only a part being included in the assembly; singleton, the number of reads that did not overlap with any other reads in the input; repeat, the number of reads deemed to be from repeat regions; outlier, the number of reads that were identified by the GS De Novo Assembler as problematic; and too short, the number of reads that were too short to be used in the computation. The second part of the table presents the frequency information on isogroups or collection of contigs containing reads that imply connections between them. The number of isogroups corresponds to the number of isogroups identified; average contig count, the average count of contigs in the isogroups; largest contig count, the largest count of contigs in the isogroups; number with one contig, the number of isogroups assembled by one contig; average isotig count, the average count of isotigs in the isogroups; largest isotig count, the largest count of isotigs in the isogroups; number with one isotig, the number of isogroups assembled by one isotig. The third part presents the frequency information of isotigs (analogous to an individual transcript). The number of isotigs correspond to the number of contigs identified; average contig count, the average count of contigs in the isotigs; largest contig count, the largest count of contigs in the isotigs; number with one contig, the number of isotigs assembled by one contig; number of bases, the total number of bases in the isotig; average isotig size, the average isotig size; N50 isotig size, the N50 isotig size (an N50 means that half of all bases reside in isotigs of this size or longer); largest isotig size, the size of the largest isotig. The last part in the table presents the frequency information of singletons. The number of singletons corresponds to the number of singletons identified; valid, the number of valid singletons after SeqClean; short, the number of sequences less than 100 bp; LowQual, undetermined base >3% in clear range; dust, standalone by low complexity (dust) filter <40 nt; ShortQ, trimmed sequence less than 100 bp; valid (after second cleaning), the number of valid singletons after Lucy; trashed, the number of low-quality sequences.Processed reads and generated contigs were aligned to the sequences publicly available in the NCBI GenBank database. Separate BLAST analyses were done for isotigs (Dataset Item 1 (Table)) and singletons (Dataset Item 2 (Table)). Of the sequences generated from NGS, 598 (14.62%) had significant matches (≤1.0E-3) to the nonredundant (nr) nucleotide database of GenBank after BLAST analysis with default options. Twenty-six (26) sequences (4.25%) of those with BLASTn hits had E-values of ≤1.0E-100 and, thus, were considered to have highly significant homology [29]. About 296 sequences (49.50%) had homology values between 1.0E-99 and 1.0E-20, which were considered as moderately similar to the sequences in GenBank. The remaining 276 isotigs (46.15%) had weak homologies or similarities (E-values between 1.0E-19 and 1.0E-3). Meanwhile, out of the 5,507 cleaned singletons, 414 had significant hits to GenBank sequences. One hundred sixty-one (161) (38.89%) of those annotated had moderate similarity to the gene sequences in GenBank. Out of the combined 1,012 BLASTn homologous sequences, 9% were similar toDaphnia pulex nucleotide sequences (Figure 4).Figure 4
Organisms with most nucleotide hits: the proportion of unique sequences with BLASTn matches (E-value ≤10−3) in the nonredundant nucleotide database classified according to the organism with the “top hit” nucleotide sequence.The remaining majority had significant homologies to several shrimp species such asPenaeus monodon (6%),Litopenaeus vannamei (3%),Marsupenaeus japonicus (3%), andFenneropenaeus chinensis (2%). Dataset Item 3 (Table) lists the database sequences with homology toP. monodon genes in the GenBank nucleotide database. Resulting gene ontology (GO) assignment showed that 88% of the total isotigs (Dataset Item 4 (Table)) and 94% of the singletons (Dataset Item 5 (Table)) had no significant hits to the Genbank GO database. Tables 3 and 4 show the isotigs and singletons GO analysis results, respectively. 488 isotigs and 336 singletons were assigned and distributed to different categories based on gene ontology. Complete lists of isotig and singleton BLASTx or GO assignments are exhibited in Dataset Item 6 (Table) and Dataset Item 7 (Table), respectively.In terms of similarities to model organisms, the majority of the 824 GO assigned sequences (40%) had significant hits to genes ofDrosophila melanogaster. Genes fromDanio rerio also had homologies to 23% of the sequences. The remaining sequences were significantly homologous to the genes ofRattus norvegicus (14%),Mus musculus (10%), and other species (Figure 5).Figure 5
Organisms with most protein hits: the proportion of putative transcripts with BLASTx matches (E-value ≤10−3) in the nonredundant protein database classified according to the organism with the “top hit” protein sequence.Microsatellite mining was as follows. Using msatcommander with search options for trinucleotide repeats or higher, 25 microsatellite markers were identified in Dataset Item 8 (Table). Five (20%) of these have BLASTn matches. The homologous genes, however, were putative or predicted. Some of these were amplified using RT-PCR with the following mixture:μL 1 cDNA template in 1x PCR buffer with 2 mM MgCl2, 0.2 mM dNTPs, 0.5 μM forward and reverse primers, 1 unit Taq, and distilled water. The PCR mix was run in a thermal cycler with the following conditions: 95°C for 5 minutes, 30 cycles of 95°C for 30 seconds, 52°C for 30 seconds, 72°C for 30 seconds, and a final extension of 72°C for 5 minutes. Successfully amplified microsatellite markers are shown in Figure 6.Figure 6
Successful amplification of microsatellites inP. monodon samples. Leftmost: marker lengths; lane 1: ladder; lanes 2 to 4: control shrimps; lanes 5 to 7: “survivor” shrimps; lanes 8: negative control; rightmost: microsatellite identifier.The successful amplification of these microsatellites indicates that repetitive segments may also be found in the transcriptome and they may be used as genetic markers relating to desired phenotype. Once confirmed to be polymorphic and after further field testing using genomic DNA, these microsatellite markers may eventually be established as markers for WSSV resistance in shrimp.Semiquantitative expression analysis was as follows. The concentration and purity of the RNA extracts were determined by spectrophotometry prior to cDNA synthesis and RT-PCR. All extracts were found to contain quality RNA (OD 260/280 of 1.8 to 2.0). Gene expression was analyzed using semiquantitative PCR. PCR mix was prepared using the same concentrations as described in microsatellite amplification. Running conditions, however, were modified into 25 cycles with a 50°C cycle annealing temperature. After the run, 2μL of each amplicon was mixed with 2 μL of loading dye and then electrophoresed in 1% agarose gel. Relative expression was determined based on the fluorescence of 2 μL DNA ladder loaded on the same gel and computed using BioDocAnalyze Software bundled with the Biometra Gel Documentation system. Ferritin and fortilin-binding protein, whose full-length coding regions were found in the library, were successfully amplified. Differential expression of these genes was clearly observed in the semiquantitative PCR profiles (Figure 7).Figure 7
Successful amplification of immune-related genes inP. monodon samples. Lane 1: ladder; lanes 2 to 4: control shrimps; lanes 5 to 7: “survivor” shrimps; lane 8: negative control. Asterisk indicates significant difference between β-actin normalized fluorescence scores of control and “survivor” shrimps based on t-test at 95% confidence interval.β-Actin was used as negative control [1] while hemocyanin, a known differentially expressed gene in virus resistant shrimps [13], was used as a positive control. It can be seen that fortilin-binding protein is significantly upregulated in “survivor” shrimps. Primers and annealing temperatures used were listed in Table 1.Table 1
Primers used in RT-PCR.
Product
Identifier
5′ to 3′ sequence
Melting temp.
170
isotig001639
Forward: TCAAAGTTAGCTCCTAGACCA
52
Reverse: AGCTCCTACTCCTGTTTCTGT
55
146
isotig001684
Forward: GCCTCCTATAAACGCTAAAAT
50
Reverse: ATTACCCAGAATTCACTGGTT
51
325
isotig00854
Forward: CTGCAGGTCATGCTCAACG
60
Reverse: AGCTCCCACACATCAAGGG
60
314
isotig01256
Forward: ACGACTCGAAGGAAGGCTC
60
Reverse: TCCCGAATGTCCGTTCTCC
60
154
HGUVU5Q01ANPJ4
Forward: TCGCGACGAGTACATGGG
60
Reverse: GTCTTCACTGTATTCAATTGGCTG
60
268
Hemocyanin
Forward: ACAGTCTCCCTCCATACACCGTGG
63
Reverse: AGGCCAGGCGAAAACGCGAA
64
279
Ferritin
Forward: GCCCGGGCAGGTACTCTCCA
66
Reverse: TGGCCGCATCGTTCTCCAGC
64
223
Fortilin-binding protein
Forward: GCTCGCAGTTAACCTCCTCTGCC
63
Reverse: GGCGTCCACCGTGCAAGAAGT
64Table 2
A list of database reads, sequences, and other components with their corresponding frequencies.
Number of reads
Number of bases
Assembled
Partial
Singleton
Repeat
Outlier
Too short
240,616
76,099,104
125,598
91,137
9,526
1,597
10,108
2,650
Number of isogroups
Average contig count
Largest contig count
Number with one contig
Average isotig count
Largest isotig count
Number with one isotig
3,724
2
2,813
3,573
1
13
3,577
Number of isotigs
Average contig count
Largest contig count
Number with one contig
Number of bases
Average isotig size
N50 isotig size
Largest isotig size
4,090
1
6
3,732
1,635,969
400
504
1,921
Number of singletons
Valid
Short
LowQual
Dust
ShortQ
Valid (after 2nd cleaning)
Trashed
9,526
5,771
3,489
1
73
192
5,507
264Table 3
Isotigs GO analysis result.
GO analysis result
Biological process
168.6118225
Cellular component
127.6988942
Molecular function
191.6892833
No hits
3602Table 4
Singletons GO analysis result.
GO analysis result
Biological process
101.4409221
Cellular component
93.16145113
Molecular function
141.3976268
No hits
5171


## 3. Dataset Description

The dataset associated with this Dataset Paper consists of 8 items which are described as follows.1
A list of isotigs and BLAST results.
The column Query Name presents the isotig number designated by the assembly program;
Query Length, the length of isotig sequence in base pairs; Query Start, the start of alignment in query;
Query End, the end of alignment in query; Query Coverage, the percentage of query covered by alignment to the database sequence;
Subject Description, the description or title of the matched database sequence; Subject Accession, the accession of the matched database sequence;
Subject Length, the length of the matched database sequence in base pairs; Subject Start, the start of alignment in subject;
Subject End, the end of alignment in subject; Subject Coverage, the percentage of database sequence covered by alignment to the query sequence;
Bit Score, the alignment bit score;E-Value Score, the alignment expect value from the database sequence; Match/Total Identities,
the number of identical matches; Identities Percentage, the percentage of identical matches;
Match/Total Positives, the number of positive scoring matches; Positives Percentage,
the percentage of positive scoring matches; Match/Total Gaps, the number of gaps in the alignment;
Gaps Percentage, the percentage of gaps in the alignment; Query Frame, the reading frame of the query sequence.
4091
20
59350
Query Name
Query Length
Query Start
Query End
Query Coverage
Subject Description
Subject Accession
Subject Length
Subject Start
Subject End
Subject Coverage
Score Bit
E-Value Score
Match/Total Identities
Identities Percentage
Match/Total Positives
Positives Percentage
Match/Total Gaps
Gaps Percentage
Query Frame2
A list of singletons and BLAST results. The column Query Name presents the singleton code designated by the assembly program; Query Length, the length of isotig sequence in base pairs; Query Start, the start of alignment in query; Query End, the end of alignment in query; Query Coverage, the percentage of query covered by alignment to the database sequence; Subject Description, the description or title of the matched database sequence; Subject Accession, the accession of the matched database sequence; Subject Length, the length of the matched database sequence in base pairs; Subject Start, the start of alignment in subject; Subject End, the end of alignment in subject; Subject Coverage, the percentage of database sequence covered by alignment to the query sequence; Bit Score, the alignment bit score;E-Value Score, the alignment expect value from the database sequence; Match/Total Identities, the number of identical matches; Identities Percentage, the percentage of identical matches; Match/Total Positives, the number of positive scoring matches; Positives Percentage, the percentage of positive scoring matches; Match/Total Gaps, the number of gaps in the alignment; Gaps Percentage, the percentage of gaps in the alignment; Query Frame, the reading frame of the query sequence.
5508
20
86583
Query Name
Query Length
Query Start
Query End
Query Coverage
Subject Description
Subject Accession
Subject Length
Subject Start
Subject End
Subject Coverage
Score Bit
E-Value Score
Match/Total Identities
Identities Percentage
Match/Total Positives
Positives Percentage
Match/Total Gaps
Gaps Percentage
Query Frame3
A list of dataset sequence names (Name),
their corresponding annotation according to the Genbank nonredundant nucleotide database (NR Annotation),
and related organisms (Related Organism).
9598
3
17181
Name
NR Annotation
Related Organism4
A list of isotig frequencies of major gene ontology (GO) categories and subcategories.
Since a gene product could be assigned to more than one GO term, the percentages in each main category do not add up to exactly 100 percent. Supplementary information is included where the isotig number, isogroup number, length, and number of contigs comprising the contig are listed per GO term.
2662
4
0
GO Category
GO Subcategory
Frequency Percentage
Supplementary Information5
A list of singleton frequencies of major gene ontology (GO) categories and subcategories.
Since a gene product could be assigned to more than one GO term, the percentages in each main category do not add up to exactly 100 percent.
Supplementary information is also included where the singleton code is listed per GO term.
1724
4
0
GO Category
GO Subcategory
Frequency Percentage
Supplementary Information6
A list of isotigs and BLAST2GO results. The column Query Name presents the isotig number designated by the assembly program;
Query Length, the length of isotig sequence in base pairs; Query Start, the start of alignment in query; Query End, the end of alignment in query;
Query Coverage, the percentage of query covered by alignment to the database sequence; Subject Description,
the description or title of the matched database sequence with specific details that contains hit ID,
sequence symbol, species number and name, GO number, specific GO term, and evidence source;
Subject Length, the length of the matched database sequence in base pairs; Subject Start, the start of alignment in subject;
Subject End, the end of alignment in subject; Subject Coverage, the percentage of database sequence covered by alignment to the query sequence;
Bit Score, the alignment bit score;E-Value Score, the alignment expect value from the database sequence;
Match/Total Identities, the number of identical matches; Identities Percentage, the percentage of identical matches;
Match/Total Positives, the number of positive scoring matches; Positives Percentage, the percentage of positive scoring matches;
Match/Total Gaps, the number of gaps in the alignment; Gaps Percentage, the percentage of gaps in the alignment; Query Frame, the reading frame of the query sequence.
4091
19
57632
Query Name
Query Length
Query Start
Query End
Query Coverage
Subject Description
Subject Length
Subject Start
Subject End
Subject Coverage
Score Bit
E-Value Score
Match/Total Identities
Identities Percentage
Match/Total Positives
Positives Percentage
Match/Total Gaps
Gaps Percentage
Query Frame7
A list of singletons and BLAST2GO results.
The column Query Name presents the singleton code designated by the assembly program;
Query Length, the length of isotig sequence in base pairs; Query Start, the start of alignment in query;
Query End, the end of alignment in query; Query Coverage, the percentage of query covered by alignment to the database sequence;
Subject Description, the description or title of the matched database sequence with specific details that contains hit ID, sequence symbol,
species number and name, GO number, specific GO term, and evidence source; Subject Length, the length of the matched database sequence in base pairs;
Subject Start, the start of alignment in subject; Subject End, the end of alignment in subject; Subject Coverage,
the percentage of database sequence covered by alignment to the query sequence; Bit Score, the alignment bit score;E-Value Score, the alignment expect value from the database sequence; Match/Total Identities,
the number of identical matches; Identities Percentage, the percentage of identical matches; Match/Total Positives,
the number of positive scoring matches; Positives Percentage, the percentage of positive scoring matches; Match/Total Gaps,
the number of gaps in the alignment; Gaps Percentage, the percentage of gaps in the alignment; Query Frame, the reading frame of the query sequence.
5508
19
82736
Query Name
Query Length
Query Start
Query End
Query Coverage
Subject Description
Subject Length
Subject Start
Subject End
Subject Coverage
Score Bit
E-Value Score
Match/Total Identities
Identities Percentage
Match/Total Positives
Positives Percentage
Match/Total Gaps
Gaps Percentage
Query Frame8
A list dataset sequences with microsatellites.
The column Sequence identifier presents the isotig number or singleton code; Length, the sequence length in base pairs;
Repeats, the number and motif of microsatellite repeats; Left Start Primer, the position of forward primer along the sequence;
Right Start Primer, the position of reverse primer along the sequence; Blastn Match, the annotation based on Blastn.
25
6
0
Sequence Identifier
Length
Repeats
Left Start Primer
Right Start Primer
Blastn Match


## 4. Concluding Remarks

The transcriptome assembly described by this paper is to be used in conjunction with Dataset Items 1–8 (Tables) to facilitate easy and rapid identification of any gene of interest from the “survivor” tiger shrimp transcriptome. This dataset represents a large pool of sequences that needs to be further characterized for additional understanding of WSSVP. monodon interaction. Once validated, candidate markers mined from the library including the generated microsatellites may be used for genetic linkage mapping.



---

*Source: 807806-2014-05-08.xml*