Jingjie Hu

Transcriptome Sequencing and De Novo Analysis for Yesso Scallop (Patinopecten yessoensis) Using 454 GS FLX

Background Bivalves comprise 30,000 extant species, constituting the second largest group of mollusks. However, limited genetic research has focused on this group of animals so far, which is, in part, due to the lack of genomic resources. The advent of high-throughput sequencing technologies enables generation of genomic resources in a short time and at a minimal cost, and therefore provides a turning point for bivalve research. In the present study, we performed de novo transcriptome sequencing to first produce a comprehensive expressed sequence tag (EST) dataset for the Yesso scallop (Patinopecten yessoensis). Results In a single 454 sequencing run, 805,330 reads were produced and then assembled into 32,590 contigs, with about six-fold sequencing coverage. A total of 25,237 unique protein-coding genes were identified from a variety of developmental stages and adult tissues based on sequence similarities with known proteins. As determined by GO annotation and KEGG pathway mapping, functional annotation of the unigenes recovered diverse biological functions and processes. Transcripts putatively involved in growth, reproduction and stress/immune-response were identified. More than 49,000 single nucleotide polymorphisms (SNPs) and 2,700 simple sequence repeats (SSRs) were also detected. Conclusion Our data provide the most comprehensive transcriptomic resource currently available for P. yessoensis. Candidate genes potentially involved in growth, reproduction, and stress/immunity-response were identified, and are worthy of further investigation. A large number of SNPs and SSRs were also identified and ready for marker development. This resource should lay an important foundation for future genetic or genomic studies on this species.

Transcriptome Sequencing and Characterization for the Sea Cucumber Apostichopus japonicus (Selenka, 1867)

Background Sea cucumbers are a special group of marine invertebrates. They occupy a taxonomic position that is believed to be important for understanding the origin and evolution of deuterostomes. Some of them such as Apostichopus japonicus represent commercially important aquaculture species in Asian countries. Many efforts have been devoted to increasing the number of expressed sequence tags (ESTs) for A. japonicus, but a comprehensive characterization of its transcriptome remains lacking. Here, we performed the large-scale transcriptome profiling and characterization by pyrosequencing diverse cDNA libraries from A. japonicus. Results In total, 1,061,078 reads were obtained by 454 sequencing of eight cDNA libraries representing different developmental stages and adult tissues in A. japonicus. These reads were assembled into 29,666 isotigs, which were further clustered into 21,071 isogroups. Nearly 40% of the isogroups showed significant matches to known proteins based on sequence similarity. Gene ontology (GO) and KEGG pathway analyses recovered diverse biological functions and processes. Candidate genes that were potentially involved in aestivation were identified. Transcriptome comparison with the sea urchin Strongylocentrotus purpuratus revealed similar patterns of GO term representation. In addition, 4,882 putative orthologous genes were identified, of which 202 were not present in the non-echinoderm organisms. More than 700 simple sequence repeats (SSRs) and 54,000 single nucleotide polymorphisms (SNPs) were detected in the A. japonicus transcriptome. Conclusion Pyrosequencing was proven to be efficient in rapidly identifying a large set of genes for the sea cucumber A. japonicus. Through the large-scale transcriptome sequencing as well as public EST data integration, we performed a comprehensive characterization of the A. japonicus transcriptome and identified candidate aestivation-related genes. A large number of potential genetic markers were also identified from the A. japonicus transcriptome. This transcriptome resource would lay an important foundation for future genetic or genomic studies on this species.

Fine-Scale Population Genetic Structure of Zhikong Scallop ( Chlamys farreri ): Do Local Marine Currents Drive Geographical Differentiation?

Marine scallops, with extended planktonic larval stages which can potentially disperse over large distances when advected by marine currents, are expected to possess low geographical differentiation. However, the sessile lifestyle as adult tends to form discrete “sea beds” with unique population dynamics and structure. The narrow distribution of Zhikong scallop (Chlamys farreri), its long planktonic larval stage, and the extremely hydrographic complexity in its distribution range provide an interesting case to elucidate the impact of marine currents on geographical differentiation for marine bivalves at a fine geographical scale. In this study, we analyzed genetic variation at nine microsatellite DNA loci in six locations throughout the distribution of Zhikong scallop in the Northern China. Very high genetic diversity was present in all six populations. Two populations sampled from the same marine gyre had no detectable genetic differentiation (FST = 0.0013); however, the remaining four populations collected from different marine gyres or separated by strong marine currents showed low but significant genetic differentiation (FST range 0.0184–0.0602). Genetic differentiation was further analyzed using the Monmonier algorithm to identify genetic barriers and using the assignment test conducted by software GeneClass2 to ascertain population membership of individuals. The genetic barriers fitting the orientation of marine gyres/currents were clearly identified, and the individual assignment analysis indicated that 95.6% of specimens were correctly allocated to one of the six populations sampled. The results support the hypothesis that significant population structure is present in Zhikong scallop at a fine geographical scale, and marine currents can be responsible for the genetic differentiation.

Construction of microsatellite-based linkage maps and identification of size-related quantitative trait loci for Zhikong scallop (Chlamys farreri).

We constructed the microsatellite-based linkage maps using 318 markers typed in two F(1) outbred families of Zhikong scallop (Chlamys farreri). The results showed an extremely high proportion (56.2%) of non-amplifying null alleles and a high ratio (30%) of segregation distortion. By aligning different individual-based linkage maps, 19 linkage groups were identified, which are consistent with the haploid chromosome number of Zhikong scallop. The integrated linkage map contains 154 markers covering 1561.8 cM with an average intermarker spacing of 12.3 cM and 77.0% of genome coverage. We found that the heterogeneity in recombination rate was not determined by sexes but by different individuals on 18 linkage regions. The phenotypic marker of general shell colour was placed on LG4, which was flanked by microsatellite markers CFLD064 and CFBD055. Four size-related traits including shell length (SL), shell width (SW), shell height (SH) and gross weight (GW) were analysed to identify the putative quantitative trait loci (QTL). Under the half-sib model, using dam as common parent, three, two, two and one QTL affecting SL, SW, SH and GW exceeded the genome-wide thresholds respectively. While using sir as common parent, a larger number of QTL were detected for these four traits: four, five, three and two for SL, SW, SH and GW respectively. The single QTL explained 3.7-19.2% of the phenotypic variation. The linkage map and the QTL associated with economic traits will provide useful information for marker-assisted selection of Zhikong scallop.

Identification of Reference Genes for qRT-PCR Analysis in Yesso Scallop Patinopecten yessoensis

Background Bivalves comprise around 30,000 extant species and have received much attention for their importance in ecosystems, aquaculture and evolutionary studies. Despite the increasing application of real-time quantitative reverse transcription PCR (qRT-PCR) in gene expression studies on bivalve species, little research has been conducted on reference gene selection which is critical for reliable and accurate qRT-PCR analysis. For scallops, systematic evaluation of reference genes that can be used among tissues or embryo/larva stages is lacking, and β-actin (ACT) is most frequently used as qRT-PCR reference gene without validation. Results In this study, 12 commonly used candidate reference genes were selected from the transcriptome data of Yesso scallop ( Patinopecten yessoensis ) for suitable qRT-PCR reference genes identification. The expression of these genes in 36 tissue samples and 15 embryo/larva samples under normal physiological conditions was examined by qRT-PCR, and their expression stabilities were evaluated using three statistic algorithms, geNorm, NormFinder, and comparative ∆Ct method. Similar results were obtained by the three approaches for the most and the least stably expressed genes. Final comprehensive ranking for the 12 genes combing the results from the three programs showed that, for different tissues, DEAD-box RNA helicase (HELI), ubiquitin (UBQ), and 60S ribosomal protein L16 (RPL16) were the optimal reference genes combination, while for different embryo/larva stages, gene set containing Cytochrome B (CB), Cytochrome C (CC), Histone H3.3 (His3.3), and Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were recommended for qRT-PCR normalization. ACT was among the least stable genes for both adult tissues and embryos/larvae. Conclusions This work constitutes the first systematic analysis on reference genes selection for qRT-PCR normalization in scallop under normal conditions. The suitable reference genes we recommended will be useful for the identification of genes related to biological processes in Yesso scallop, and also in the reference gene selection for other scallop or bivalve species.

Solving the SAT problem using a DNA computing algorithm based on ligase chain reaction

A new DNA computing algorithm based on a ligase chain reaction is demonstrated to solve an SAT problem. The proposed DNA algorithm can solve an n-variable m-clause SAT problem in m steps and the computation time required is O (3m+n). Instead of generating the full-solution DNA library, we start with an empty test tube and then generate solutions that partially satisfy the SAT formula. These partial solutions are then extended step by step by the ligation of new variables using Taq DNA ligase. Correct strands are amplified and false strands are pruned by a ligase chain reaction (LCR) as soon as they fail to satisfy the conditions. If we score and sort the clauses, we can use this algorithm to markedly reduce the number of DNA strands required throughout the computing process. In a computer simulation, the maximum number of DNA strands required was 2(0.48n) when n=50, and the exponent ratio varied inversely with the number of variables n and the clause/variable ratio m/n. This algorithm is highly space-efficient and error-tolerant compared to conventional brute-force searching, and thus can be scaled-up to solve large and hard SAT problems.

FISH Mapping and Identification of Zhikong Scallop ( Chlamys farreri ) Chromosomes

Chromosome identification is the first step in genomic research of a species, but it remains a challenge in scallops. In the present study, fluorescence in situ hybridization (FISH) mapping of 19 fosmid clones was attempted and used for chromosome identification in Zhikong scallop (Chlamys farreri Jones et Preston, 1904). Data showed that 10 clones were successfully mapped, including 7 without and 3 with C0t-1 DNA. Among them, 2 represented multiple signals and made no contribution to chromosome identification. Karyotypic analysis and cohybridization indicated that the remaining 8 clones realized the identification of 8 chromosomes. All 10 clones were sequenced at both ends, which could be developed as sequence-tagged sites and used for the unification of the cytological and genetic linkage maps. This study shows that fosmid clones can benefit chromosome identification and will undoubtedly be useful for cytogenetic research in Zhikong scallop.

Characterization of 95 novel microsatellite markers for Zhikong scallop Chlamys farreri using FIASCO-colony hybridization and EST database mining

In order to construct a simple sequence repeat (SSR)-based genetic linkage map and to promote molecular marker-assisted selection (MAS) in scallop breeding, the methods of Fast Isolation by AFLP of Sequences COntaining repeats (FIASCO)-colony hybridization and expressed sequence tag (EST) database mining were modified and used to develop 95 novel microsatellite markers for Zhikong scallop. The SSR-enriched library constructed by the FIASCO method consisted of 830 clones, and 295 (35.5%) positive clones were identified after colony hybridization. One hundred and fifty clones were randomly sequenced and the results showed all clones contained at least one microsatellite. Of 91 primer pairs designed, 72 were amplified scorable polymerase chain reaction (PCR) products and 70 were polymorphic with the allele number range of 3–16 alleles/locus (average 7.0 alleles/locus). When EST database mining was performed, 66 microsatellites containing ESTs were identified from 3467 sequences deposited in GenBank. Based on cluster analysis of length and GC content of the flanking regions, 47 primer pairs were designed and 23 scorable EST SSRs were obtained. Compared with genomic SSRs developed in this study, EST SSRs showed lower genetic variability with an average of 4.2 alleles/locus. The results in the present study demonstrate that modified FIASCO-colony hybridization is an efficient and low-cost method for the isolation of large numbers of microsatellite markers for scallop species.

Analysis of the Secondary Structure of ITS1 in Pectinidae: Implications for Phylogenetic Reconstruction and Structural Evolution

It is at present difficult to accurately position gaps in sequence alignment and to determine substructural homology in structure alignment when reconstructing phylogenies based on highly divergent sequences. Therefore, we have developed a new strategy for inferring phylogenies based on highly divergent sequences. In this new strategy, the whole secondary structure presented as a string in bracket notation is used as phylogenetic characters to infer phylogenetic relationships. It is no longer necessary to decompose the secondary structure into homologous substructural components. In this study, reliable phylogenetic relationships of eight species in Pectinidae were inferred from the structure alignment, but not from sequence alignment, even with the aid of structural information. The results suggest that this new strategy should be useful for inferring phylogenetic relationships based on highly divergent sequences. Moreover, the structural evolution of ITS1 in Pectinidae was also investigated. The whole ITS1 structure could be divided into four structural domains. Compensatory changes were found in all four structural domains. Structural motifs in these domains were identified further. These motifs, especially those in D2 and D3, may have important functions in the maturation of rRNAs.

Initial analysis of tandemly repetitive sequences in the genome of Zhikong scallop (Chlamys farreri Jones et Preston).

Tandemly repetitive sequences are widespread in all eukaryotic genomes, but data on tandem repeats are limited in Zhikong scallop (Chlamys farreri). In the present study, paired-end sequencing of 2016 individual fosmid clones resulted in 3646 sequences. A total of 2,286,986 bp of genomic sequences were generated, representing approximately 1.84‰ of the Zhikong scallop genome. Using tandem repeats finder (TRF) software, a total of 2500 tandem repeats were found, including 313 satellites, 1816 minisatellites and 371 microsatellites. The cumulative length of tandem repeats was 552,558 bp, accounting for 24.16% of total length. Specifically, the length of microsatellites, minisatellites and satellites was 9425, 336,001 and 207,132 bp, accounting for 1.71, 60.81 and 37.49% of the length of tandem repeats, and 0.41, 14.69 and 9.06% of total length, respectively. The detailed information on the characteristic of all repeat units was also represented, which will provide a useful resource for physical mapping and better utilization of the existing genomic information in Zhikong scallop.