The scallop IGF2 mRNA-binding protein gene PyIMP and association of a synonymous mutation with growth traits

Mutations in Growth-Related Genes Induced by EMS Treatment in Scallops

Background: The goal of genetic breeding is to select variants with mutations that are related to expected traits, such as fast growth. Artificial induction has been widely used to obtain strains with more mutations for further selection. Ethylmethylsulfone (EMS) is one of the most commonly used chemical mutagens in plant and microorganism breeding. However, the application of EMS mutagenesis in shellfish has not been reported. The aim of this study is to evaluate the potential use of EMS as a mutagen in scallop breeding, especially in characterization of mutations in growth-related genes. Results: Our results indicated that hatching of about 50% of fertilized eggs was blocked by treatment with 20 mM EMS for 3 h and the resulted larvae developed normally into adult stages. We then evaluated the mutagenic effects of EMS by sequencing the genomes of 4 adult scallops from the control group and 12 from the treatment group at 8 months after fertilization. On average, after removing shared types of mutations, there were 1,151,380 ± 258,188 SNPs (Single Nucleotide Polymorphisms) and 229,256 ± 51,714 InDels (insertion-deletion) in each animal in the EMS treatment group, while there were only134841 ± 10,115 SNPs and 42,605 ± 5,136 InDels in the control group. The average mutation rate in the genome of the EMS treatment group (0.0137 ± 0.0013%) was about 9 times that of the control group (0.0015 ± 0.0002%). GO (Gene Ontology) annotation and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analyses revealed that mutations induced by EMS occurred evenly in most biological processes, cellular components and functions, as well in most pathways. However, significant lower percentage of mutations were found in the exonic region, in non-synonymous or Stopgain/Stoploss SNPs and in coding domains, suggesting apparent DNA repair or selection during grow-out stage. Analyses of the growth-related genes with mutations indicated that mutations in MFS (Major Facilitator Superfamily) and Tubulin were only found in the large-sized group (Five largest scallops: Treated-1, Treated-2, Treated-3, Treated-4, and Treated-5) and Homeobox and Socs (Suppressor of cytokine signaling) only in the small group (Two smallest scallops: Treated-11 and Treated-12). These results suggested that these genes may be involved in the regulation of growth in these animals, although further verification is certainly warranted. Conclusion: Treatment of fertilized eggs with 20 mM EMS for 3 h induced 9 times more mutations in scallop genomes. We found that mutations in MFS and Tubulin may be related to fast growth in the large-sized group and those mutations in Homeobox and SOCs may be involved in the slow growth in the small-sized scallops. EMS can be used to accelerate selection of economically important traits in molluscs.

Mantle Transcriptome Provides Insights into Biomineralization and Growth Regulation in the Eastern Oyster (Crassostrea virginica)

Growth of the eastern oyster Crassostrea virginica, a major aquaculture species in the USA, is highly variable and not well understood at molecular levels. As growth of mollusks is confined in shells constructed by the mantle, mantle transcriptomes of large (fast-growing) and small (slow-growing) eastern oysters were sequenced and compared in this study. Transcription was observed for 31,186 genes, among which 104 genes were differentially expressed between the large and small oysters, including 48 upregulated and 56 downregulated in large oysters. Differentially expressed genes (DEGs) included genes from diverse pathways highlighting the complexity of shell formation and growth regulations. Seventeen of the 48 upregulated DEGs were related to shell matrix formation, most of which were upregulated in large oysters, indicating that large oysters are more active in biomineralization and shell formation. Genomic and transcriptomic analyses identified 22 genes encoding novel polyalanine containing proteins (Pacps) with characteristic motifs for matrix function that are tandemly duplicated on one chromosome, all specifically expressed in mantle and at higher levels in large oysters, suggesting that these expanded Pacps play important roles in shell formation and growth. Analysis of sequence variation identified 244,964 SNPs with 328 associated with growth. This study provides novel candidate genes and markers for shell formation and growth, and suggests that genes related to shell formation are important for the complex regulation of growth in the eastern oyster and possibly other bivalve mollusks. Results of this study show that both transcriptional modulation and functional polymorphism are important in determining growth.

Vital Role of Glutamate Dehydrogenase Gene in Ammonia Detoxification and the Association Between its SNPs and Ammonia Tolerance in Sinonovacula constricta

Increasing evidence has revealed accumulated ammonia will cause adverse effects on the growth, reproduction, and survival of aquatic animals. As a marine benthic mollusk, the razor clam Sinonovacula constricta shows better growth and survival under high ammonia nitrogen environment. However, little is known about its adaptation mechanisms to high ammonia stress in an integrated mariculture system. In this study, we analyzed the association between the polymorphism of glutamate dehydrogenase gene (GDH), a key gene involved in ammonia nitrogen detoxification, and ammonia tolerance. The results showed that 26 and 22 single-nucleotide polymorphisms (SNPs) of GDH in S. constricta (denoted as Sc-GDH) were identified from two geographical populations, respectively. Among them, two SNPs (c.323T > C and c.620C > T) exhibited a significant and strong association with ammonia tolerance, suggesting that Sc-GDH gene could serve as a potential genetic marker for molecular marker–assisted selection to increase survival rate and production of S. constricta. To observe the histological morphology and explore the histocellular localization of Sc-GDH, by paraffin section and hematoxylin–eosin staining, the gills were divided into gill filament (contains columnar and flattened cells) and gill cilia, whereas hepatopancreas was made up of individual hepatocytes. The results of immunohistochemistry indicated that the columnar cells of gill filaments and the endothelial cells of hepatocytes were the major sites for Sc-GDH secretion. Under ammonia stress (180 mg/L), the expression levels of Sc-GDH were extremely significantly downregulated at 24, 48, 72, and 96 h (P < 0.01) after RNA interference. Thus, we can speculate that Sc-GDH gene may play an important role in the defense process against ammonia stress. Overall, these findings laid a foundation for further research on the adaptive mechanisms to ammonia–nitrogen tolerance for S. constricta.

Identification of a Growth-Associated Single Nucleotide Polymorphism (SNP) in Cyclin C of the Giant Tiger Shrimp Penaeus monodon

The full-length cDNA of cyclin C of the giant tiger shrimp Penaeus monodon (PmCyC) was isolated by RACE-PCR. It was 1443 bp in length containing an open reading frame (ORF) of 804 bp and 267 deduced amino acids. Tissue distribution analysis indicated that PmCyC was more abundantly expressed in ovaries and testes than other tissues of female and male juveniles (P < 0.05). A pair of primers was designed, and an amplification product of 403 bp containing an intron of 123 bp was obtained. Polymorphism of amplified PmCyC gene segments of the 5th (3-month-old G5, N = 30) and 7th (5-month-old G7, N = 18) generations of domesticated juveniles was analyzed. Four conserved SNPs (T>C₁₃₄, T>C₁₈₈, G>A₃₇₉, and T>C₃₈₂) were found within the examined sequences. A TaqMan genotyping assay was developed for detection of a T>C₁₃₄ SNP. Association analysis indicated that this SNP displayed significant association with body weight (P < 4.2e-10) and total length (P < 2e-09) of the examined G7 P. monodon (N = 419) with an allele substitution effect of 5.02 ± 0.78 g and 1.41 ± 0.19 cm, respectively. Juveniles with C/C₁₃₄ (22.80 ± 2.51 g and 12.97 ± 0.53 cm, N = 19) and T/C₁₃₄ (20.41 ± 0.93 g and 12.77 ± 0.21 cm, N = 129) genotypes exhibited a significantly greater average body weight and total length than those with a T/T₁₃₄ genotype (14.72 ± 0.53 g and 11.37 ± 0.13 cm, N = 271) (P < 0.05).