The Y chromosome sequence of the channel catfish suggests novel sex determination mechanisms in teleost fish.

BACKGROUND: Sex determination mechanisms in teleost fish broadly differ from mammals and birds, with sex chromosomes that are far less differentiated and recombination often occurring along the length of the X and Y chromosomes, posing major challenges for the identification of specific sex determination genes. Here, we take an innovative approach of comparative genome analysis of the genomic sequences of the X chromosome and newly sequenced Y chromosome in the channel catfish. RESULTS: Using a YY channel catfish as the sequencing template, we generated, assembled, and annotated the Y genome sequence of channel catfish. The genome sequence assembly had a contig N50 size of 2.7 Mb and a scaffold N50 size of 26.7 Mb. Genetic linkage and GWAS analyses placed the sex determination locus within a genetic distance less than 0.5 cM and physical distance of 8.9 Mb. However, comparison of the channel catfish X and Y chromosome sequences showed no sex-specific genes. Instead, comparative RNA-Seq analysis between females and males revealed exclusive sex-specific expression of an isoform of the breast cancer anti-resistance 1 (BCAR1) gene in the male during early sex differentiation. Experimental knockout of BCAR1 gene converted genetic males (XY) to phenotypic females, suggesting BCAR1 as a putative sex determination gene. CONCLUSIONS: We present the first Y chromosome sequence among teleost fish, and one of the few whole Y chromosome sequences among vertebrate species. Comparative analyses suggest that sex-specific isoform expression through alternative splicing may underlie sex determination processes in the channel catfish, and we identify BCAR1 as a potential sex determination gene.

Characterization, expression signatures and microbial binding analysis of cathepsin A in turbot, Scophthalmus maximus L.(SmCTSA).

Mucosal immune system is one of the most vital components in the innate immunity and constitutes the first line of host defense against bacterial infections, especially for the teleost, which live in the pathogen-rich aquatic environment. Cathepsins, a superfamily of hydrolytic enzymes produced and enclosed within lysosomes, play multiple roles at physiological and pathological states. In this regard, we sought here to identify Cathepsin A in turbot (SmCTSA), characterize its mucosal expression patterns following Vibrio anguillarum and Streptococcus iniae infections in mucosal tissues, and explore its binding ability with three microbial ligands for the first time. The SmCTSA was 2631 bp long containing a 1422 bp open reading frame (ORF) that encoded 473 amino acids. Phylogenetic analysis revealed that SmCTSA showed the closest relationship to half-smooth tongue sole (Cynoglossus semilaevis). In addition, SmCTSA was ubiquitously expressed in all examined healthy tissues, with high expression levels in head kidney (HK) and intestine, while the lowest expression level in blood. Moreover, SmCTSA was significantly differentially expressed at least two timepoints in each mucosal tissue, suggesting its potential important roles in innate immune responses of turbot. Finally, in vitro assays showed that recombinant SmCTSA bound Lipopolysaccharide (LPS) with high affinity, and lipoteichoic acid (LTA) and peptidoglycan (PGN) with relatively low affinity. This study provides valuable data for understanding the roles of ctsa in the host defense against bacterial infections.

A genome-wide association study in catfish reveals the presence of functional hubs of related genes within QTLs for columnaris disease resistance.

BACKGROUND: Columnaris causes severe mortalities among many different wild and cultured freshwater fish species, but understanding of host resistance is lacking. Catfish, the primary aquaculture species in the United States, serves as a great model for the analysis of host resistance against columnaris disease. Channel catfish in general is highly resistant to the disease while blue catfish is highly susceptible. F2 generation of hybrids can be produced where phenotypes and genotypes are segregating, providing a useful system for QTL analysis. To identify genes associated with columnaris resistance, we performed a genome-wide association study (GWAS) using the catfish 250 K SNP array with 340 backcross progenies derived from crossing female channel catfish (Ictalurus punctatus) with male F1 hybrid catfish (female channel catfish I. punctatus × male blue catfish I. furcatus). RESULTS: A genomic region on linkage group 7 was found to be significantly associated with columnaris resistance. Within this region, five have known functions in immunity, including pik3r3b, cyld-like, adcyap1r1, adcyap1r1-like, and mast2. In addition, 3 additional suggestively associated QTL regions were identified on linkage groups 7, 12, and 14. The resistant genotypes on the QTLs of linkage groups 7 and 12 were found to be homozygous with both alleles being derived from channel catfish. The paralogs of the candidate genes in the suggestively associated QTL of linkage group 12 were found on the QTLs of linkage group 7. Many candidate genes on the four associated regions are involved in PI3K pathway that is known to be required by many bacteria for efficient entry into the host. CONCLUSION: The GWAS revealed four QTLs associated with columnaris resistance in catfish. Strikingly, the candidate genes may be arranged as functional hubs; the candidate genes within the associated QTLs on linkage groups 7 and 12 are not only co-localized, but also functionally related, with many of them being involved in the PI3K signal transduction pathway, suggesting its importance for columnaris resistance.

RNA-Seq reveals expression signatures of genes involved in oxygen transport, protein synthesis, folding, and degradation in response to heat stress in catfish.

Temperature is one of the most prominent abiotic factors affecting ectotherms. Most fish species, as ectotherms, have extraordinary ability to deal with a wide range of temperature changes. While the molecular mechanism underlying temperature adaptation has long been of interest, it is still largely unexplored with fish. Understanding of the fundamental mechanisms conferring tolerance to temperature fluctuations is a topic of increasing interest as temperature may continue to rise as a result of global climate change. Catfish have a wide natural habitat and possess great plasticity in dealing with environmental variations in temperature. However, no studies have been conducted at the transcriptomic level to determine heat stress-induced gene expression. In the present study, we conducted an RNA-Seq analysis to identify heat stress-induced genes in catfish at the transcriptome level. Expression analysis identified a total of 2,260 differentially expressed genes with a cutoff of twofold change. qRT-PCR validation suggested the high reliability of the RNA-Seq results. Gene ontology, enrichment, and pathway analyses were conducted to gain insight into physiological and gene pathways. Specifically, genes involved in oxygen transport, protein folding and degradation, and metabolic process were highly induced, while general protein synthesis was dramatically repressed in response to the lethal temperature stress. This is the first RNA-Seq-based expression study in catfish in response to heat stress. The candidate genes identified should be valuable for further targeted studies on heat tolerance, thereby assisting the development of heat-tolerant catfish lines for aquaculture.

Whole genome comparative analysis of channel catfish (Ictalurus punctatus) with four model fish species.

BACKGROUND: Comparative mapping is a powerful tool to study evolution of genomes. It allows transfer of genome information from the well-studied model species to non-model species. Catfish is an economically important aquaculture species in United States. A large amount of genome resources have been developed from catfish including genetic linkage maps, physical maps, BAC end sequences (BES), integrated linkage and physical maps using BES-derived markers, physical map contig-specific sequences, and draft genome sequences. Application of such genome resources should allow comparative analysis at the genome scale with several other model fish species. RESULTS: In this study, we conducted whole genome comparative analysis between channel catfish and four model fish species with fully sequenced genomes, zebrafish, medaka, stickleback and Tetraodon. A total of 517 Mb draft genome sequences of catfish were anchored to its genetic linkage map, which accounted for 62% of the total draft genome sequences. Based on the location of homologous genes, homologous chromosomes were determined among catfish and the four model fish species. A large number of conserved syntenic blocks were identified. Analysis of the syntenic relationships between catfish and the four model fishes supported that the catfish genome is most similar to the genome of zebrafish. CONCLUSION: The organization of the catfish genome is similar to that of the four teleost species, zebrafish, medaka, stickleback, and Tetraodon such that homologous chromosomes can be identified. Within each chromosome, extended syntenic blocks were evident, but the conserved syntenies at the chromosome level involve extensive inter-chromosomal and intra-chromosomal rearrangements. This whole genome comparative map should facilitate the whole genome assembly and annotation in catfish, and will be useful for genomic studies of various other fish species.

Efficient assembly and annotation of the transcriptome of catfish by RNA-Seq analysis of a doubled haploid homozygote.

BACKGROUND: Upon the completion of whole genome sequencing, thorough genome annotation that associates genome sequences with biological meanings is essential. Genome annotation depends on the availability of transcript information as well as orthology information. In teleost fish, genome annotation is seriously hindered by genome duplication. Because of gene duplications, one cannot establish orthologies simply by homology comparisons. Rather intense phylogenetic analysis or structural analysis of orthologies is required for the identification of genes. To conduct phylogenetic analysis and orthology analysis, full-length transcripts are essential. Generation of large numbers of full-length transcripts using traditional transcript sequencing is very difficult and extremely costly. RESULTS: In this work, we took advantage of a doubled haploid catfish, which has two sets of identical chromosomes and in theory there should be no allelic variations. As such, transcript sequences generated from next-generation sequencing can be favorably assembled into full-length transcripts. Deep sequencing of the doubled haploid channel catfish transcriptome was performed using Illumina HiSeq 2000 platform, yielding over 300 million high-quality trimmed reads totaling 27 Gbp. Assembly of these reads generated 370,798 non-redundant transcript-derived contigs. Functional annotation of the assembly allowed identification of 25,144 unique protein-encoding genes. A total of 2,659 unique genes were identified as putative duplicated genes in the catfish genome because the assembly of the corresponding transcripts harbored PSVs or MSVs (in the form of pseudo-SNPs in the assembly). Of the 25,144 contigs with unique protein hits, around 20,000 contigs matched 50% length of reference proteins, and over 14,000 transcripts were identified as full-length with complete open reading frames. The characterization of consensus sequences surrounding start codon and the stop codon confirmed the correct assembly of the full-length transcripts. CONCLUSIONS: The large set of transcripts assembled in this study is the most comprehensive set of genome resources ever developed from catfish, which will provide the much needed resources for functional genome research in catfish, serving as a reference transcriptome for genome annotation, analysis of gene duplication, gene family structures, and digital gene expression analysis. The putative set of duplicated genes provide a starting point for genome scale analysis of gene duplication in the catfish genome, and should be a valuable resource for comparative genome analysis, genome evolution, and genome function studies.

Alternative complement pathway of channel catfish (Ictalurus punctatus): molecular characterization, mapping and expression analysis of factors Bf/C2 and Df.

The complement system is important in both innate and adaptive host defense against microbial infection in vertebrates. It contains three pathways: the classical, alternative, and lectin pathways. Complement component factors B and D are two crucial proteases in the alternative pathway. In this study, the genes of complement factors Bf/C2 and Df from channel catfish, Ictalurus punctatus were identified and characterized. Two complement factor B-related genes, Bf/C2A and Bf/C2B, and factor D gene Df were identified. Phylogenetic analysis suggested that Bf/C2A and Bf/C2B is likely orthologous to factor B and factor C2, respectively. Southern blot results suggested that these three genes are all single-copy genes in the catfish genome. The catfish Bf/C2A, Bf/C2B and Df genes were genetically mapped on linkage group 3, 20 and 29, respectively. Bf/C2A and Bf/C2B are highly expressed in liver and kidney, while Df is highly expressed in gill and spleen. After infection with Edwardsiella ictaluri, the expression of Bf/C2A, Bf/C2B and Df genes were found to be remarkably induced in the gill, liver, spleen and kidney at some sampling times, indicating that these three complement factors play a pivotal role in immune responses after the bacterial infection in catfish.

A pilot study for channel catfish whole genome sequencing and de novo assembly.

BACKGROUND: Recent advances in next-generation sequencing technologies have drastically increased throughput and significantly reduced sequencing costs. However, the average read lengths in next-generation sequencing technologies are short as compared with that of traditional Sanger sequencing. The short sequence reads pose great challenges for de novo sequence assembly. As a pilot project for whole genome sequencing of the catfish genome, here we attempt to determine the proper sequence coverage, the proper software for assembly, and various parameters used for the assembly of a BAC physical map contig spanning approximately a million of base pairs. RESULTS: A combination of low sequence coverage of 454 and Illumina sequencing appeared to provide effective assembly as reflected by a high N50 value. Using 454 sequencing alone, a sequencing depth of 18 X was sufficient to obtain the good quality assembly, whereas a 70 X Illumina appeared to be sufficient for a good quality assembly. Additional sequencing coverage after 18 X of 454 or after 70 X of Illumina sequencing does not provide significant improvement of the assembly. Considering the cost of sequencing, a 2 X 454 sequencing, when coupled to 70 X Illumina sequencing, provided an assembly of reasonably good quality. With several software tested, Newbler with a seed length of 16 and ABySS with a K-value of 60 appear to be appropriate for the assembly of 454 reads alone and Illumina paired-end reads alone, respectively. Using both 454 and Illumina paired-end reads, a hybrid assembly strategy using Newbler for initial 454 sequence assembly, Velvet for initial Illumina sequence assembly, followed by a second step assembly using MIRA provided the best assembly of the physical map contig, resulting in 193 contigs with a N50 value of 13,123 bp. CONCLUSIONS: A hybrid sequencing strategy using low sequencing depth of 454 and high sequencing depth of Illumina provided the good quality assembly with high N50 value and relatively low cost. A combination of Newbler, Velvet, and MIRA can be used to assemble the 454 sequence reads and the Illumina reads effectively. The assembled sequence can serve as a resource for comparative genome analysis. Additional long reads using the third generation sequencing platforms are needed to sequence through repetitive genome regions that should further enhance the sequence assembly.

Generation of genome-scale gene-associated SNPs in catfish for the construction of a high-density SNP array.

BACKGROUND: Single nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies. In order to provide the best genome coverage for the analysis of performance and production traits, a large number of relatively evenly distributed SNPs are needed. Gene-associated SNPs may fulfill these requirements of large numbers and genome wide distribution. In addition, gene-associated SNPs could themselves be causative SNPs for traits. The objective of this project was to identify large numbers of gene-associated SNPs using high-throughput next generation sequencing. RESULTS: Transcriptome sequencing was conducted for channel catfish and blue catfish using Illumina next generation sequencing technology. Approximately 220 million reads (15.6 Gb) for channel catfish and 280 million reads (19.6 Gb) for blue catfish were obtained by sequencing gene transcripts derived from various tissues of multiple individuals from a diverse genetic background. A total of over 35 billion base pairs of expressed short read sequences were generated. Over two million putative SNPs were identified from channel catfish and almost 2.5 million putative SNPs were identified from blue catfish. Of these putative SNPs, a set of filtered SNPs were identified including 342,104 intra-specific SNPs for channel catfish, 366,269 intra-specific SNPs for blue catfish, and 420,727 inter-specific SNPs between channel catfish and blue catfish. These filtered SNPs are distributed within 16,562 unique genes in channel catfish and 17,423 unique genes in blue catfish. CONCLUSIONS: For aquaculture species, transcriptome analysis of pooled RNA samples from multiple individuals using Illumina sequencing technology is both technically efficient and cost-effective for generating expressed sequences. Such an approach is most effective when coupled to existing EST resources generated using traditional sequencing approaches because the reference ESTs facilitate effective assembly of the expressed short reads. When multiple individuals with different genetic backgrounds are used, RNA-Seq is very effective for the identification of SNPs. The SNPs identified in this report will provide a much needed resource for genetic studies in catfish and will contribute to the development of a high-density SNP array. Validation and testing of these SNPs using SNP arrays will form the material basis for genome association studies and whole genome-based selection in catfish.

Oncogenic extrachromosomal DNA functions as mobile enhancers to globally amplify chromosomal transcription.

Extrachromosomal, circular DNA (ecDNA) is emerging as a prevalent yet less characterized oncogenic alteration in cancer genomes. We leverage ChIA-PET and ChIA-Drop chromatin interaction assays to characterize genome-wide ecDNA-mediated chromatin contacts that impact transcriptional programs in cancers. ecDNAs in glioblastoma patient-derived neurosphere and prostate cancer cell cultures are marked by widespread intra-ecDNA and genome-wide chromosomal interactions. ecDNA-chromatin contact foci are characterized by broad and high-level H3K27ac signals converging predominantly on chromosomal genes of increased expression levels. Prostate cancer cells harboring synthetic ecDNA circles composed of characterized enhancers result in the genome-wide activation of chromosomal gene transcription. Deciphering the chromosomal targets of ecDNAs at single-molecule resolution reveals an association with actively expressed oncogenes spatially clustered within ecDNA-directed interaction networks. Our results suggest that ecDNA can function as mobile transcriptional enhancers to promote tumor progression and manifest a potential synthetic aneuploidy mechanism of transcription control in cancer.

Identification and characterisation of a novel small galectin in razor clam (Sinonovacula constricta) with multiple innate immune functions.

Galectins are lectins possessing an evolutionarily conserved carbohydrate recognition domain (CRD) with affinity for ß-galactoside. The key role played by innate immunity in invertebrates has recently become apparent. Herein, a full-length galectin (ScGal) was identified in razor clam (Sinonovacula constricta). The 528 bp open reading frame encodes a polypeptide of 176 amino acids with a single CRD and no signal peptide. ScGal mRNA transcripts were mainly expressed in hemolymph and gill, and were significantly up-regulated following bacterial challenge. Recombinant rScGal protein binds to and aggregates various bacteria, and has affinity for peptidoglycan, lipoteichoic acid and d-galactose. The protein also stimulates hemocytes to phagocytose invading bacterial pathogens. ScGal is an important immune factor in innate immunity, and a small protein with multiple important functions.

Three novel ANO5 missense mutations in Caucasian and Chinese families and sporadic cases with gnathodiaphyseal dysplasia.

Gnathodiaphyseal dysplasia (GDD; MIM#166260) is an autosomal dominant syndrome with characteristic cemento-osseous lesions of jawbones, bone fragility, and diaphyseal sclerosis of tubular bones. To date, only five mutations in the proposed calcium-activated chloride channel ANO5/TMEM16E gene have been identified. In this study, we describe two families and two singular patients with three new mutations. One Caucasian family with seven affected members exhibited frequent bone fractures and florid osseous dysplasia (p.Cys356Tyr), while one Chinese family with two affected members suffered from cementoma and purulent osteomyelitis (p.Cys360Tyr). In addition, two different novel mutations (p.Gly518Glu and p.Arg215Gly) were identified in sporadic patients without family history. In vitro studies overexpressing GDD mutations (p.Cys356Tyr and p.Cys360Tyr) showed significantly reduced ANO5 protein. It appears that all GDD mutations known so far locate in an extracellular domain following the first transmembrane domain or in the 4th putative transmembrane domain. Both wild-type and mutant ANO5 protein localize to the endoplasmic reticulum. After Ano5 gene knock-down with shRNA in MC3T3-E1 osteoblast precursors we saw elevated expression of osteoblast-related genes such as Col1a1, osteocalcin, osterix and Runx2 as well as increased mineral nodule formation in differentiating cells. Our data suggest that ANO5 plays a role in osteoblast differentiation.

Identification of ASAH1 as a susceptibility gene for familial keloids.

Keloids result from abnormal proliferative scar formation with scar tissue expanding beyond the margin of the original wound and are mostly found in individuals of sub-Saharan African descent. The etiology of keloids has not been resolved but previous studies suggest that keloids are a genetically heterogeneous disorder. Although possible candidate genes have been suggested by genome-wide association studies using common variants, by upregulation in keloids or their involvement in syndromes that include keloid formation, rare coding variants that contribute to susceptibility in non-syndromic keloid formation have not been previously identified. Through analysis of whole-genome data we mapped a locus to chromosome 8p23.3-p21.3 with a statistically significant maximum multipoint LOD score of 4.48. This finding was followed up using exome sequencing and led to the identification of a c.1202T>C (p.(Leu401Pro)) variant in the N-acylsphingosine amidohydrolase (ASAH1) gene that co-segregates with the keloid phenotype in a large Yoruba family. ASAH1 is an acid ceramidase known to be involved in tumor formation by controlling the ratio of ceramide and sphingosine. ASAH1 is also involved in cell proliferation and inflammation, and may affect the development of keloids via multiple mechanisms. Functional studies need to clarify the role of the ASAH1 variant in wound healing.

Developmental Transcriptome Analysis and Identification of Genes Involved in Larval Metamorphosis of the Razor Clam, Sinonovacula constricta.

The razor clam Sinonovacula constricta is an important commercial species. The deficiency of developmental transcriptomic data is becoming the bottleneck of further researches on the mechanisms underlying settlement and metamorphosis in early development. In this study, de novo transcriptome sequencing was performed for S. constricta at different early developmental stages by using Illumina HiSeq 2000 paired-end (PE) sequencing technology. A total of 112,209,077 PE clean reads were generated. De novo assembly generated 249,795 contigs with an average length of 585 bp. Gene annotation resulted in the identification of 22,870 unigene hits against the NCBI database. Eight unique sequences related to metamorphosis were identified and analyzed using real-time PCR. The razor clam reference transcriptome would provide useful information on early developmental and metamorphosis mechanisms and could be used in the genetic breeding of shellfish.

The channel catfish genome sequence provides insights into the evolution of scale formation in teleosts.

Catfish represent 12% of teleost or 6.3% of all vertebrate species, and are of enormous economic value. Here we report a high-quality reference genome sequence of channel catfish (Ictalurus punctatus), the major aquaculture species in the US. The reference genome sequence was validated by genetic mapping of 54,000 SNPs, and annotated with 26,661 predicted protein-coding genes. Through comparative analysis of genomes and transcriptomes of scaled and scaleless fish and scale regeneration experiments, we address the genomic basis for the most striking physical characteristic of catfish, the evolutionary loss of scales and provide evidence that lack of secretory calcium-binding phosphoproteins accounts for the evolutionary loss of scales in catfish. The channel catfish reference genome sequence, along with two additional genome sequences and transcriptomes of scaled catfishes, provide crucial resources for evolutionary and biological studies. This work also demonstrates the power of comparative subtraction of candidate genes for traits of structural significance.

Saliva as a source of genomic DNA for genetic studies: review of current methods and applications.

DNA isolated from saliva has become an attractive alternative to the use of blood-derived DNA in genetic studies and is now extensively used in many applications. This review addresses advantages of saliva DNA for large population-based studies and points out caveats for use in certain techniques. We compiled data comparing saliva-derived genomic DNA to blood-derived DNA with regards to quality, quantity and convenience. Special attention is given to the usefulness of saliva-derived DNA for PCR-based methods and genome-wide analyses.

Comprehensive transcriptome study to develop molecular resources of the copepod Calanus sinicus for their potential ecological applications.

Calanus sinicus Brodsky (Copepoda, Crustacea) is a dominant zooplanktonic species widely distributed in the margin seas of the Northwest Pacific Ocean. In this study, we utilized an RNA-Seq-based approach to develop molecular resources for C. sinicus. Adult samples were sequenced using the Illumina HiSeq 2000 platform. The sequencing data generated 69,751 contigs from 58.9 million filtered reads. The assembled contigs had an average length of 928.8 bp. Gene annotation allowed the identification of 43,417 unigene hits against the NCBI database. Gene ontology (GO) and KEGG pathway mapping analysis revealed various functional genes related to diverse biological functions and processes. Transcripts potentially involved in stress response and lipid metabolism were identified among these genes. Furthermore, 4,871 microsatellites and 110,137 single nucleotide polymorphisms (SNPs) were identified in the C. sinicus transcriptome sequences. SNP validation by the melting temperature (T m )-shift method suggested that 16 primer pairs amplified target products and showed biallelic polymorphism among 30 individuals. The present work demonstrates the power of Illumina-based RNA-Seq for the rapid development of molecular resources in nonmodel species. The validated SNP set from our study is currently being utilized in an ongoing ecological analysis to support a future study of C. sinicus population genetics.

Development of the catfish 250K SNP array for genome-wide association studies.

BACKGROUND: Quantitative traits, such as disease resistance, are most often controlled by a set of genes involving a complex array of regulation. The dissection of genetic basis of quantitative traits requires large numbers of genetic markers with good genome coverage. The application of next-generation sequencing technologies has allowed discovery of over eight million SNPs in catfish, but the challenge remains as to how to efficiently and economically use such SNP resources for genetic analysis. RESULTS: In this work, we developed a catfish 250K SNP array using Affymetrix Axiom genotyping technology. The SNPs were obtained from multiple sources including gene-associated SNPs, anonymous genomic SNPs, and inter-specific SNPs. A set of 640K high-quality SNPs obtained following specific requirements of array design were submitted. A panel of 250,113 SNPs was finalized for inclusion on the array. The performance evaluated by genotyping individuals from wild populations and backcross families suggested the good utility of the catfish 250K SNP array. CONCLUSIONS: This is the first high-density SNP array for catfish. The array should be a valuable resource for genome-wide association studies (GWAS), fine QTL mapping, high-density linkage map construction, haplotype analysis, and whole genome-based selection.

Genome sequence and genetic diversity of the common carp, Cyprinus carpio.

The common carp, Cyprinus carpio, is one of the most important cyprinid species and globally accounts for 10% of freshwater aquaculture production. Here we present a draft genome of domesticated C. carpio (strain Songpu), whose current assembly contains 52,610 protein-coding genes and approximately 92.3% coverage of its paleotetraploidized genome (2n = 100). The latest round of whole-genome duplication has been estimated to have occurred approximately 8.2 million years ago. Genome resequencing of 33 representative individuals from worldwide populations demonstrates a single origin for C. carpio in 2 subspecies (C. carpio Haematopterus and C. carpio carpio). Integrative genomic and transcriptomic analyses were used to identify loci potentially associated with traits including scaling patterns and skin color. In combination with the high-resolution genetic map, the draft genome paves the way for better molecular studies and improved genome-assisted breeding of C. carpio and other closely related species.

Development of molecular resources for an intertidal clam, Sinonovacula constricta, using 454 transcriptome sequencing.

BACKGROUND: The razor clam Sinonovacula constricta is a benthic intertidal bivalve species with important commercial value. Despite its economic importance, knowledge of its transcriptome is scarce. Next generation sequencing technologies offer rapid and efficient tools for generating large numbers of sequences, which can be used to characterize the transcriptome, to develop effective molecular markers and to identify genes associated with growth, a key breeding trait. RESULTS: Total RNA was isolated from the mantle, gill, liver, siphon, gonad and muscular foot tissues. High-throughput deep sequencing of S. constricta using 454 pyrosequencing technology yielded 859,313 high-quality reads with an average read length of 489 bp. Clustering and assembly of these reads produced 16,323 contigs and 131,346 singletons with average lengths of 1,376 bp and 458 bp, respectively. Based on transcriptome sequencing, 14,615 sequences had significant matches with known genes encoding 147,669 predicted proteins. Subsequently, previously unknown growth-related genes were identified. A total of 13,563 microsatellites (SSRs) and 13,634 high-confidence single nucleotide polymorphism loci (SNPs) were discovered, of which almost half were validated. CONCLUSION: De novo sequencing of the razor clam S. constricta transcriptome on the 454 GS FLX platform generated a large number of ESTs. Candidate growth factors and a large number of SSRs and SNPs were identified. These results will impact genetic studies of S. constricta.