The oyster genome reveals stress adaptation and complexity of shell formation

The Pacific oyster Crassostrea gigas belongs to one of the most species-rich but genomically poorly explored phyla, the Mollusca. Here we report the sequencing and assembly of the oyster genome using short reads and a fosmid-pooling strategy, along with transcriptomes of development and stress response and the proteome of the shell. The oyster genome is highly polymorphic and rich in repetitive sequences, with some transposable elements still actively shaping variation. Transcriptome studies reveal an extensive set of genes responding to environmental stress. The expansion of genes coding for heat shock protein 70 and inhibitors of apoptosis is probably central to the oyster's adaptation to sessile life in the highly stressful intertidal zone. Our analyses also show that shell formation in molluscs is more complex than currently understood and involves extensive participation of cells and their exosomes. The oyster genome sequence fills a void in our understanding of the Lophotrochozoa. © 2012 Macmillan Publishers Limited. All rights reserved

Full-Length RNA Sequencing Provides Insights into Goldfish Evolution under Artificial Selection

Goldfish Carassius auratus is an ideal model for exploring fish morphology evolution. Although genes underlying several ornamental traits have been identified, little is known about the effects of artificial selection on embryo gene expression. In the present study, hybrid transcriptome sequencing was conducted to reveal gene expression profiles of Celestial-Eye (CE) and Ryukin (RK) goldfish embryos. Full-length transcriptome sequencing on the PacBio platform identified 54,218 and 54,106 transcript isoforms in CE and RK goldfish, respectively. Of particular note was that thousands of alternative splicing (AS) and alternative polyadenylation (APA) events were identified in both goldfish breeds, and most of them were inter-breed specific. RT-PCR and Sanger sequencing showed that most of the predicted AS and APA were correct. Moreover, abundant long non-coding RNA and fusion genes were detected, and again most of them were inter-breed specific. Through RNA-seq, we detected thousands of differentially expressed genes (DEGs) in each embryonic stage between the two goldfish breeds. KEGG enrichment analysis on DEGs showed extensive differences between CE and RK goldfish in gene expression. Taken together, our results demonstrated that artificial selection has led to far-reaching influences on goldfish gene expression, which probably laid the genetic basis for hundreds of goldfish variations

RestrictionDigest: A powerful Perl module for simulating genomic restriction digests

Background: Reduced-representation sequencing technology is widely used in genotyping for its economical and efficient features. A popular way to construct the reduced-representation sequencing libraries is to digest the genomic DNA with restriction enzymes. A key factor of this method is to determine the restriction enzyme(s). But there are few computer programs which can evaluate the usability of restriction enzymes in reduced-representation sequencing. SimRAD is an R package which can simulate the digestion of DNA sequence by restriction enzymes and return enzyme loci number as well as fragment number. But for linkage mapping analysis, enzyme loci distribution is also an important factor to evaluate the enzyme. For phylogenetic studies, comparison of the enzyme performance across multiple genomes is important. It is strongly needed to develop a simulation tool to implement these functions. Results: Here, we introduce a Perl module named RestrictionDigest with more functions and improved performance. It can analyze multiple genomes at one run and generate concise comparison of enzyme performance across the genomes. It can simulate single-enzyme digestion, double-enzyme digestion and size selection process and generate comprehensive information of the simulation including enzyme loci number, fragment number, sequences of the fragments, positions of restriction sites on the genome, the coverage of digested fragments on different genome regions and detailed fragment length distribution. Conclusions: RestrictionDigest is an easy-to-use Perl module with flexible parameter settings. With the help of the information produced by the module, researchers can easily determine the most appropriate enzymes to construct the reduced-representation libraries to meet their experimental requirements

Classification and statistical analysis of SNPs in the two subspecies.

<p>A cell with a different background color represents different types of SNPs between the subspecies, and the number in each cell represents the amount of each type of SNP. Lilac: inter-subspecific SNPs, pink: northern subspecies specific SNPs, reseda: southern subspecies specific SNPs, azure: shared SNPs with the same heterozygosis type between the two subspecies and light brown: shared SNPs with a different heterozygosis type between the two subspecies. Abbreviations: N, the northern subspecies of bay scallop; S, the southern subspecies of bay scallop; A, Adenine; T, Thymine; G, Guanine; C, Cytosine; K, Thymine or Guanine; M, Adenine or Cytosine; R, Adenine or Guanine; S, Guanine or Cytosine; W, Adenine or Thymine; Y, Thymine or Cytosine.</p

Seawater Culture Increases Omega-3 Long-Chain Polyunsaturated Fatty Acids (N-3 LC-PUFA) Levels in Japanese Sea Bass (Lateolabrax japonicus), Probably by Upregulating Elovl5

The fatty acid compositions of the fish muscle and liver are substantially affected by rearing environment. However, the mechanisms underlying this effect have not been thoroughly described. In this study, we investigated the effects of different culture patterns, i.e., marine cage culture and freshwater pond culture, on long-chain polyunsaturated fatty acids (LC-PUFA) biosynthesis in an aquaculturally important fish, the Japanese sea bass (Lateolabrax japonicus). Fish were obtained from two commercial farms in the Guangdong province, one of which raises Japanese sea bass in freshwater, while the other cultures sea bass in marine cages. Fish were fed the same commercial diet. We found that omega-3 long-chain polyunsaturated fatty acids (n-3 LC-PUFA) levels in the livers and muscles of the marine cage cultured fish were significantly higher than those in the livers and muscles of the freshwater pond cultured fish. Quantitative real-time PCRs indicated that fatty acid desaturase 2 (FADS2) transcript abundance was significantly lower in the livers of the marine cage reared fish as compared to the freshwater pond reared fish, but that fatty acid elongase 5 (Elovl5) transcript abundance was significantly higher. Consistent with this, two of the 28 CpG loci in the FADS2 promoter region were heavily methylated in the marine cage cultured fish, but were only slightly methylated in freshwater pond cultured fish (n = 5 per group). Although the Elovl5 promoter was less methylated in the marine cage reared fish as compared to the freshwater pond reared fish, this difference was not significant. Thus, our results might indicate that Elovl5, not FADS2, plays an important role in the enhancing LC-PUFA synthesis in marine cage cultures

Betaine Alleviates High-Fat Diet Induced Excessive Lipid Deposition in Gibel Carp Hepatopancreas and L8824 Cells by Enhancing VLDL Secretion through HNF4α/MTTP Pathway

Betaine, a methyl donor, plays a crucial role in lipid metabolism. Previous studies have shown that appropriate betaine supplementation in a high-fat diet reduces triglycerides (TG) of serum and hepatopancreas in fish. However, the underlying mechanism remains unclear. This study examined whether betaine can enhance the secretion of very low-density lipoprotein (VLDL) and sought to identify the specific mechanisms through which this enhancement occurs. A lipid accumulation model was established in gibel carp and L8824 cells using a high-fat diet and oleic acid, respectively. Different doses of betaine (1, 4, and 16 g/kg in the diet; 400 μmol in cell culture) were administered, and measurements were taken for lipid deposition, gene expression of HNF4α, MTTP, and ApoB, as well as the regulation of Mttp and Apob promoters by HNF4α. The results showed that betaine supplementation mitigated lipid droplet accumulation, TG levels, and VLDL production induced by the high-fat diet in gibel carp hepatopancreas and L8824 cells. Moreover, betaine not only increased VLDL content in the cell culture supernatant but also reversed the inhibitory effects of the high-fat diet on protein expression of MTTP, ApoB, and HNF4α in both gibel carp hepatopancreas and L8824 cells. Additionally, HNF4α exhibits transactivating activity on the promoter of Mttp in gibel carp. These findings suggest that betaine supplementation exerts its effects through the HNF4α/MTTP/ApoB pathway, promoting the assembly and secretion of VLDL and effectively reducing lipid accumulation in the hepatopancreas of farmed gibel carp fed a high-fat diet

SNP Identification by Transcriptome Sequencing and Candidate Gene-Based Association Analysis for Heat Tolerance in the Bay Scallop <i>Argopecten irradians</i>

<div><p>The northern bay scallop <i>Argopecten irradians irradians</i> (Lamarck) and the southern bay scallop <i>Argopecten irradians concentricus</i> (Say) were introduced into China in the 1980s and 1990s, and are now major aquaculture molluscs in China. Here, we report the transcriptome sequencing of the two subspecies and the subsequent association analysis on candidate gene on the trait of heat tolerance. In total, RNA from six tissues of 67 and 42 individuals of northern and southern bay scallops, respectively, were used and 55.5 and 34.9 million raw reads were generated, respectively. There were 82,267 unigenes produced in total, of which 32,595 were annotated. Altogether, 32,206 and 23,312 high-quality SNPs were identified for northern and southern bay scallops, respectively. For case-control analysis, two intercrossed populations were heat stress treated, and both heat-susceptible and heat-resistant individuals were collected. According to annotation and SNP allele frequency analysis, 476 unigenes were selected, and 399 pairs of primers were designed. Genotyping was conducted using the high-resolution melting method, and Fisher’s exact test was performed for allele frequency comparison between the heat-susceptible and heat-resistant groups. SNP all-53308-760 T/C showed a significant difference in allele frequency between the heat-susceptible and heat-resistant groups. Notably, considerable difference in allele frequency at this locus was also observed between the sequenced natural populations. These results suggest that SNP all-53308-760 T/C may be related to the heat tolerance of the bay scallop. Moreover, quantitative expression analysis revealed that the expression level of all-53308 was negatively correlated with heat tolerance of the bay scallop.</p></div

Association results of two populations, assessed using Fisher’s exact test.

<p>*Please see the “Hybridization of the two subspecies and the heat stress experiment” section for ZZ96 and ZN96.</p

Summary of basic information after transcriptome sequencing.

<p>*Assembly using clean reads of the two subspecies combined.</p><p>**Unigene quantity and length information after clustering.</p

General flowchart for SNP identification and association analysis.

<p>Abbreviations: N, the northern subspecies of bay scallop; S, the southern subspecies of bay scallop; QHD, Qinhuangdao, Hebei Province, China; QD, Qingdao, Shandong Province, China; YT, Yantai, Shandong Province, China; ZJ, Zhanjiang, Guangdong Province, China; Z, hybrids of the two subspecies; ZZ, a descendant of the cross between Z (♂) and Z (♀); ZN, a descendant of the cross between Z (♂) and N (♀).</p