Haigang Qi

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.

Molecular Basis for Adaptation of Oysters to Stressful Marine Intertidal Environments

Oysters that occupy estuarine and intertidal habitats have well-developed stress tolerance mechanisms to tolerate harsh and dynamically changing environments. In this review, we summarize common pathways and genomic features in oyster that are responsive to environmental stressors such as temperature, salinity, hypoxia, air exposure, pathogens, and anthropogenic pollutions. We first introduce the key genes involved in several pathways, which constitute the molecular basis for adaptation to stress. We use genome analysis to highlight the strong cellular homeostasis system, a unique adaptive characteristic of oysters. Next, we provide a global view of features of the oyster genome that contribute to stress adaptation, including oyster-specific gene expansion, highly inducible expression, and functional divergence. Finally, we review the consequences of interactions between oysters and the environment from ecological and evolutionary perspectives by discussing mass mortality and adaptive divergence among populations and related species of the genus Crassostrea. We conclude with prospects for future study.

Discovery and validation of genic single nucleotide polymorphisms in the Pacific oyster Crassostrea gigas

The economic and ecological importance of the oyster necessitates further research on the molecular mechanisms, which both regulate the commercially important traits of the oyster and help it to survive in the variable marine environment. Single nucleotide polymorphisms (SNPs) have been widely used to assess genetic variation and identify genes underlying target traits. In addition, high‐resolution melting (HRM) analysis is a potentially powerful method for validating candidate SNPs. In this study, we adopted a rapid and efficient pipeline for the screening and validation of SNPs in the genic region of Crassostrea gigas based on transcriptome sequencing and HRM analysis. Transcriptomes of three wild oyster populations were sequenced using Illumina sequencing technology. In total, 50–60 million short reads, corresponding to 4.5–5.4 Gbp, from each population were aligned to the oyster genome, and 5.8 × 105 SNPs were putatively identified, resulting in a predicted SNP every 47 nucleotides on average. The putative SNPs were unevenly distributed in the genome and high‐density (≥2%), nonsynonymous coding SNPs were enriched in genes related to apoptosis and responses to biotic stimuli. Subsequently, 1,671 loci were detected by HRM analysis, accounting for 64.7% of the total selected candidate primers, and finally, 1,301 polymorphic SNP markers were developed based on HRM analysis. All of the validated SNPs were distributed into 897 genes and located in 672 scaffolds, and 275 of these genes were stress inducible under unfavourable salinity, temperature, and exposure to air and heavy metals. The validated SNPs in this study provide valuable molecular markers for genetic mapping and characterization of important traits in oysters.

Characterization of 12 single nucleotide polymorphisms (SNPs) in Pacific abalone, Haliotis discus hannai

We report here for the first time 12 polymorphic single nucleotide polymorphisms (SNPs) in a commercially important gastropod, Pacific abalone (Haliotis discus hannai) that were identified by searching expressed sequence tag database. These SNP loci (seven nuclear and five mitochondrial SNPs) were polymorphic among 37 wild abalone individuals, based on a four‐primer allele‐specific polymerase chain reaction analysis. All loci had two alleles and the minor allele frequency ranged from 0.027 to 0.473. For the seven nuclear SNPs, the expected and observed heterozygosities ranged from 0.053 to 0.499 and from 0.054 to 0.811, respectively.

Development and characterization of microsatellite markers for the Pacific Oyster Crassostrea gigas

Thirty-three polymorphic microsatellite DNA markers were developed and characterized for the Pacific Oyster Crassostrea gigas, a widespread oyster species with great economic and ecological importance, by the end-sequencing of an in-house fosmid library. The polymorphism of each locus was analyzed by screening 30 Pacific oysters from a natural population. The size of repeat motif ranged from 2 to 6 and the number of alleles per locus ranged from 3 to 9. The observed and expected heterozygosity ranged from 0.080 to 0.880 and from 0.174 to 0.876, respectively. These microsatellite markers will contribute to the increasing genetic studies in C. gigas.

Candidate Gene Polymorphisms and their Association with Glycogen Content in the Pacific Oyster Crassostrea gigas.

Background The Pacific oyster Crassostrea gigas is an important cultivated shellfish that is rich in nutrients. It contains high levels of glycogen, which is of high nutritional value. To investigate the genetic basis of this high glycogen content and its variation, we conducted a candidate gene association analysis using a wild population, and confirmed our results using an independent population, via targeted gene resequencing and mRNA expression analysis. Results We validated 295 SNPs in the 90 candidate genes surveyed for association with glycogen content, 86 of were ultimately genotyped in all 144 experimental individuals from Jiaonan (JN). In addition, 732 SNPs were genotyped via targeted gene resequencing. Two SNPs (Cg_SNP_TY202 and Cg_SNP_3021) in Cg_GD1 (glycogen debranching enzyme) and one SNP (Cg_SNP_4) in Cg_GP1 (glycogen phosphorylase) were identified as being associated with glycogen content. The glycogen content of individuals with genotypes TT and TC in Cg_SNP_TY202 was higher than that of individuals with genotype CC. The transcript abundance of both glycogen-associated genes was differentially expressed in high glycogen content and low glycogen content individuals. Conclusions This study identified three polymorphisms in two genes associated with oyster glycogen content, via candidate gene association analysis. The transcript abundance differences in Cg_GD1 and Cg_GP1 between low- and the high-glycogen content individuals suggests that it is possible that transcript regulation is mediated by variations of Cg_SNP_TY202, Cg_SNP_3021, and Cg_SNP_4. These findings will not only provide insights into the genetic basis of oyster quality, but also promote research into the molecular breeding of oysters.

Identification and functional characterization of two Bcl-2 family protein genes in Zhikong scallop Chlamys farreri

Apoptosis plays significant roles in maintenance of homeostasis, immune defense and development. The Bcl-2 family proteins are important regulators of the intrinsic apoptosis. In the study, we have characterized a Bcl-2-like gene (named CfBcl-2) and a Bax-like gene (named CfBax) from the Zhikong scallop Chlamys farreri. The full-length of the CfBcl-2 cDNA is 944 nucleotides (nt) encoding a putative protein of 225 amino acid residues (aa) that contains four Bcl-2 homology (BH) domains, and the CfBax cDNA is 505 nt encoding a putative protein of 115 aa that contains three Bcl-2 BH domains. Sequence and phylogenetic analysis demonstrate that CfBcl-2 and CfBax present typical domain organization of the corresponding Bcl-2 related proteins and are more similar and clustered with their homologues of other molluscs. The two genes are ubiquitously expressed in six tissues of C. farreri, with the highest expression level of CfBcl-2 in adductor muscle and highest expression level of CfBax in gill. The expressions of CfBcl-2 and CfBax in hemocytes were both significantly up-regulated after an in vivo exposure of scallops to air, injection with lipopolysaccharide and infection with acute viral necrobiotic disease virus, and the expression patterns of the two genes after the three treatments vary in different change magnitude and up-regulation timespan. Yeast two-hybrid assay reveals a direct interaction between the CfBcl-2 and CfBax proteins. These results indicate that the CfBcl-2 and CfBax may participate in the apoptosis-based stress and immune responses against noxious stimulation.

Development of gene-associated intronic TR markers for the Pacific abalone Haliotis discus hannai

Source/description: The Pacific abalone, Haliotis discus hannai, is one of the most important shellfish cultured in China. The large-scale hatchery production of this marine mollusc calls for the development and implementation of sound genetic management practices and selective breeding programmes. Genetic markers such as RAPD, AFLP and microsatellites have been developed and used in H. discus hannai. Useful and novel DNA markers, especially type I markers associated with genes of interest, are particularly useful for association studies to elucidate possible gene function. We developed intronic tandem repeat (TR) markers associated with genes of known function using exon-primed intron-crossing polymerase chain reaction (EPIC-PCR) sequencing.

Genome-wide survey and analysis of microsatellites in the Pacific oyster genome: abundance, distribution, and potential for marker development

Microsatellites are a ubiquitous component of the eukaryote genome and constitute one of the most popular sources of molecular markers for genetic studies. However, no data are currently available regarding microsatellites across the entire genome in oysters, despite their importance to the aquaculture industry. We present the first genome-wide investigation of microsatellites in the Pacific oyster Crassostrea gigas by analysis of the complete genome, resequencing, and expression data. The Pacific oyster genome is rich in microsatellites. A total of 604 653 repeats were identified, in average of one locus per 815 base pairs (bp). A total of 12 836 genes had coding repeats, and 7 332 were expressed normally, including genes with a wide range of molecular functions. Compared with 20 different species of animals, microsatellites in the oyster genome typically exhibited 1) an intermediate overall frequency; 2) relatively uniform contents of (A)n and (C)n repeats and abundant long (C)n repeats (≥24 bp); 3) large average length of (AG)n repeats; and 4) scarcity of trinucleotide repeats. The microsatellite-flanking regions exhibited a high degree of polymorphism with a heterozygosity rate of around 2.0%, but there was no correlation between heterozygosity and microsatellite abundance. A total of 19 462 polymorphic microsatellites were discovered, and dinucleotide repeats were the most active, with over 26% of loci found to harbor allelic variations. In all, 7 451 loci with high potential for marker development were identified. Better knowledge of the microsatellites in the oyster genome will provide information for the future design of a wide range of molecular markers and contribute to further advancements in the field of oyster genetics, particularly for molecular-based selection and breeding.

Fosmid library construction and end sequences analysis of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is globally distributed and is one of the most commercially and ecologically important marine organisms. However, little is known about the genome of this species. In this study, a C. gigas fosmid library was constructed that contains 459,936 clones with an average insert size of approximately 40 kb, representing 22.34-fold haploid genome equivalents. End sequencing generated 90,240 fosmid end sequences (FESs) with an average length of 384.27 base pairs (bp), covering approximately 2.58% of the Pacific oyster genome. The FESs were subsequently assembled and annotated, resulting in 6332 sequences with predicted open reading frames≥300 and 1,189,100 bp repeats. Furthermore, a total of 3200 microsatellite repeats were identified, and dinucleotide repeats were found to occur most abundantly, with AG and AAT being the most abundant repeat class of dinucleotides and trinucleotides. We also found that the repeat number was generally negatively proportional to the repeat element length. Microsatellites composition between the transcribed sequences and genomic sequences was shown to be different. Point mutations of microsatellite were non-random and underwent strong selection stress. Overall, a comprehensive sequence resource for the Pacific oyster was created, including annotated transposable elements, tandem repeats, protein coding sequences and microsatellites. These initial findings will serve as resources for further in-depth studies of physical mapping, gene discovery, microsatellite marker developing and evolution studies.