Molecular characterization of Myostatin gene from Zhikong scallop Chlamys farreri (Jones et Preston 1904)

Myocyte enhancer factor 2 upregulates expression of myostatin promoter in Yesso scallop, Patinopecten yessoensis

Myostatin (MSTN) plays an important role in muscle development in animals, especially for mammals and fishes. However, little information has been reported on the regulation of MSTN in marine invertebrates, such as bivalves. In the present study, we cloned the MSTN promoter sequence of Yesso scallop Patinopecten yessoensis, identifying 4 transcription start sites, eleven TATA boxes and one E-box. Additionally, transcription factor binding sites, including myocyte enhancer factor 2 (MEF2) and POU homeodomain protein, were identified. The interaction between the MSTN promoter and MEF2 was analyzed to reveal the transcriptional activity of different fragment sizes of promoters through the dual-luciferase reporter assays. The highest transcriptional activity was found in recombinant plasmids with the most MEF2 binding sites, indicating that this transcription factor upregulates MSTN in Yesso scallop. This study provides new insight into the regulation of muscle growth and development in this species.

A Potential Negative Regulatory Function of Myostatin in the Growth of the Pacific Abalone, Haliotis discus hannai

Myostatin, also known as GDF8, is a member of the transforming growth factor-β (TGF-β) superfamily. In vertebrates, myostatin negatively regulates the growth of skeletal muscle. In invertebrates, it has been reported to be closely related to animal growth. However, knowledge concerning the molecular mechanisms involved in the myostatin regulation of molluscan growth is limited. In this study, we found that the hdh-myostatin open reading frame (ORF) comprised 1470 base pairs that encoded 489 amino acids and contained structural characteristics typical of the TGF-β superfamily, including a C-terminal signal peptide, a propeptide domain, and TGF-β region. Gene expression analysis revealed that hdh-myostatin mRNA was widely expressed at different levels in all of the examined tissues of Haliotis discus hannai. Nine single nucleotide polymorphisms (SNPs) were associated with the growth traits. RNA interference (RNAi) against hdh-myostatin mRNA significantly downregulated hdh-myostatin at days 1, 15, and 30 post injection, and the pattern was correlated with downregulation of the genes TGF-β receptor type-I (hdh-TβR I), activin receptor type-IIB (hdh-ActR IIB), and mothers against decapentaplegic 3 (hdh-Smad3). After one month of the RNAi experiment, the shell lengths and total weights increased in the abalone, Haliotis discus hannai. The results of qRT-PCR showed that the hdh-myostatin mRNA level was higher in the slow-growing group than in the fast-growing group. These results suggest that hdh-myostatin is involved in the regulation of growth, and that these SNPs would be informative for further studies on selective breeding in abalone.

Hyperaccumulation of cadmium by scallop Chlamys farreri revealed by comparative transcriptome analysis

Cadmium (Cd) is a hazardous environmental contaminant, which has a serious effect on the ecosystem, food safety and human health. Scallop could accumulate high concentration of Cd from the environment and has been regarded as a Cd hyper-accumulator. In this work, we investigated the antioxidative defense, detoxification and transport of Cd in the kidneys of scallops by transcriptome analysis. A total of 598 differentially expressed genes including 387 up-regulated and 211 down-regulated ones were obtained during Cd exposure, and 46 up-regulated and 260 down-regulated ones were obtained during depuration. Cadmium exposure could cause oxidative stress in the kidneys, which was particularly shown in the pathways involved in proteasome and oxidative phosphorylation. The mRNA expression of 5 metallothionein (MT) genes were overexpressed under Cd exposure and significantly decreased during Cd depuration, which played a vital role in Cd chelation and detoxification. The expression of divalent metal transporter (DMT) genes were down-regulated insignificantly during accumulation and depuration of Cd, which suggested that the DMT played little roles in Cd transport in scallops. A positive relationship in the expression of the zinc transporter (ZIP6 and ZIP1) genes with Cd exposure and depuration was observed, which confirmed its important role for Cd uptake in the kidneys of scallops. 26S proteasome activities and MT expression were Cd-dependent. This study supplied the important reference on the hyperaccumulation of Cd by scallops and identified some effective bioindicators for the environmental risk assessment.

Molecular characterization, expression analysis of myostatin gene and its negative regulation by miR-29b-3p in Chinese concave-eared frogs (Odorrana tormota)

The molecular characteristics, expression patterns and functions of the amphibian myostatin (MSTN) gene are unknown. Here, we isolated a full-length Odorrana tormota MSTN cDNA sequence of 1701 bp (Ot-MSTN), containing a putative N-terminal signal peptide, a TGF-β propeptide domain and an active peptide. Ot-MSTN was expressed in 9 selected tissues examined, and the highest level of expression was in thigh muscle, followed by brain and female gonadal tissue. The expression of Ot-MSTN in multiple O. tormota tissues supported that the activities of MSTN may be not limited to skeletal muscle. Ot-MSTN expression was decreased from stage 31 to stage 40, while the growth rate was increased. The expression of Ot-MSTN in adult male frogs increased with age, indicating that adult male frogs may inhibit the continued hypertrophy of thigh muscle fibers and decrease the growth rate of thigh muscle to ensure muscles do not grow too large. Luciferase reporter assays showed that miR-29b-3p directly targeted the 3'-UTR of Ot-MSTN. miR-29b-3p expression in the thigh muscle of 2 yrs. females who grew faster was significantly lower than that of the slow-growing 2 yrs. male individuals, which showed an opposite trend with Ot-MSTN expression. In addition，miR-29b-3p expression reversed trends of Ot-MSTN expression at different developmental stages in thigh muscle. Therefore, these data indicate that miR-29-3p may negatively regulate the expression of MSTN and regulate thigh muscle growth and development in O. tormota.

Myostatin gene promoter: structure, conservation and importance as a target for muscle modulation

Myostatin (MSTN) is one of the key factors regulating myogenesis. Because of its role as a negative regulator of muscle mass deposition, much interest has been given to its protein and, in recent years, several studies have analysed MSTN gene regulation. This review discusses the MSTN gene promoter, focusing on its structure in several animal species, both vertebrate and invertebrate. We report the important binding sites considering their degree of phylogenetic conservation and roles they play in the promoter activity. Finally, we discuss recent studies focusing on MSTN gene regulation via promoter manipulation and the potential applications they have both in medicine and agriculture.

Dynamics of DNA Methylation and DNMT Expression During Gametogenesis and Early Development of Scallop Patinopecten yessoensis

DNA methylation reprograms during gametogenesis and embryo development, which is essential for germ cell specification and genomic imprinting in mammals. Corresponding process remains poorly investigated in molluscs. Here, we examined global DNA methylation level in the gonads of scallop Patinopecten yessoensis during gametogenesis and in embryos/larvae at different stages. DNA methylation level fluctuates during gametogenesis and early development, peaking at proliferative stage of ovary, growing stage of testis, and in blastulae. To understand the mechanisms underlying these changes, we conducted genome-wide characterization of DNMT family and investigated their expression profiles based on transcriptomes and in situ hybridization. Three genes were identified, namely PyDNMT1, PyDNMT2, and PyDNMT3. Expression of PyDnmt3 agrees with DNA methylation level during oogenesis and early development, suggesting PyDNMT3 may participate in de novo DNA methylation that occurs mainly at proliferative stage of ovary and testis, and in blastulae and gastrulae. PyDnmt1 expression is positively correlated with DNA methylation level during spermatogenesis, and is higher at maturation stage of ovary and in 2-8 cell embryos than other stages, implying possible involvement of PyDNMT1 in DNA methylation maintenance during meiosis and embryonic development. This study will facilitate better understanding of the developmental epigenetic reprogramming in bivalve molluscs.

Molecular cloning and characteristics analysis of Pmtgfbr1 from Pinctada fucata martensii

•

This study obtains the full length of Pmtgfbr1 of the pearl oyster P. fucata martensii.

•

Pmtgfbr1 possesses the conserved domain of Tgfbr1.

•

Pmtgfbr1 holds negatively effect on the growth of P. fucata martensii.

Pinctada fucata martensii is cultured for pearl production. Growth improvement has received considerable research interest. Transforming growth factor β type Ⅰ receptor (TβR-I), which is involved in signals transmission of transforming growth factor beta (TGF-β), participates in cell proliferation and growth. In this study, we characterized a Tgfbr1 gene which encoded TβR-I from P. fucata martensii (Pmtgfbr1). Pmtgfbr1 cDNA contains an open reading frame of 1569 bp and encodes a polypeptide of 522 amino acids (aa). Pmtgfbr1 possesses a typical TβR-I structure (extracellular receptor ligand domain, transmembrane domain, and cytoplasmic tyrosine kinase catalytic domain). Pmtgfbr1 is expressed in all the studied tissues and exhibited the highest expression level in the adductor muscle. Moreover, Pmtgfbr1 exhibited the lower expression level in the larger group (L) than that in the smaller group (S) and is negatively correlated with growth traits (P < 0.01). Our results indicated that Pmtgfbr1 is a candidate functional gene associated with growth traits.

Expression of the myostatin gene in the adductor muscle of the Pacific lion-paw scallop Nodipecten subnodosus in association with growth and environmental conditions

The cDNA sequence of the myostatin gene in the Pacific lion-paw Nodipecten subnodosus (Ns-mstn) was characterized, and the temporal expression during grow-out was analyzed for the first time in a scallop. Ns-mstn encodes a 459-amino-acid protein in which two propeptide proteolytic sites were identified, the previously recognized (RSKR) and a second one at position 266-269 aa (RRKR). The alternative furin cleavage site could be related with post-translational processing, or it could be a tissue-specific mechanism for signaling activity. The Ns-mstn transcript was located by in situ hybridization in sarcomeres and around the nucleus of muscle fibers. The temporal expression analysis by qPCR in the adductor muscle showed that Ns-mstn expression was significantly different (P < 0.05) between months during the grow-out period, increasing largely during the summer months when both biomass and muscle weight did not increase or even decreased; muscle fiber size and number were found to decrease significantly. Exogenous and endogenous factors such as high temperature and low food availability, as well as gametogenesis and reproduction, can be associated with the growth pattern and Ns-mstn expression changes. Our results indicate that MSTN is involved in adductor muscle growth regulation in N. subnodosus as it occurs in vertebrate skeletal muscle although Ns-mstn expression in non-muscle organs/tissues suggests additional functions.

The myostatin gene of Mytilus chilensis evidences a high level of polymorphism and ubiquitous transcript expression

Myostatin (MSTN) is a protein of the Transforming Growth Factor-β (TGF-β) superfamily and plays a crucial role in muscular development for higher vertebrates. However, its biological function in marine invertebrates remains undiscovered. This study characterizes the full-length sequence of the Mytilus chilensis myostatin gene (Mc-MSTN). Furthermore, tissue transcription patterns and putative single nucleotide polymorphisms (SNPs) were also identified. The Mc-MSTN cDNA sequence showed 3528 base pairs (bp), consisting of 161 bp of 5' UTR, 2,110 bp of 3' UTR, and an open reading frame of 1,257 bp encoding for 418 amino acids and with an RXXR proteolytic site and nine cysteine-conserved residues. Gene transcription analysis revealed that the Mc-MSTN has ubiquitous expression among several tissues, with higher expression in the gonads and mantle than in the digestive gland, gills, and hemolymph. Furthermore, high levels of polymorphisms were detected (28 SNPs in 3'-UTR and 9 SNPs in the coding region). Two SNPs were non-synonymous and involved amino acid changes between Glu/Asp and Thr/Ile. Until now, the MSTN gene has been mainly related to muscle growth in marine bivalves. However, the present study suggests a putative biological function not entirely associated to muscle tissue and contributes molecular evidence to the current debate about the function of the MSTN gene in marine invertebrates.

A SNP in the 5' flanking region of the myostatin-1b gene is associated with harvest traits in Atlantic salmon (Salmo salar)

Background

Myostatin (MSTN) belongs to the transforming growth factor-β superfamily and is a potent negative regulator of skeletal muscle development and growth in mammals. Most teleost fish possess two MSTN paralogues. However, as a consequence of a recent whole genome-duplication event, salmonids have four: MSTN-1 (−1a and -1b) and MSTN-2 (−2a and -2b). Evidence suggests that teleost MSTN plays a role in the regulation of muscle growth. In the current study, the MSTN-1b gene was re-sequenced and screened for SNP markers in a commercial population of Atlantic salmon. After genotyping 4,800 progeny for the discovered SNPs, we investigated their association with eight harvest traits - four body-weight traits, two ratios of weight traits, flesh colour and fat percentage - using a mixed model association analysis.

Results

Three novel SNPs were discovered in the MSTN-1b gene of Atlantic salmon. One of the SNPs, located within the 5′ flanking region (g.1086C > T), had a significant association with harvest traits (p < 0.05), specifically for: Harvest Weight (kg), Gutted Weight (kg), Deheaded Weight (kg) and Fillet Weight (kg). The haplotype-based association analysis was consistent with this result because the two haplotypes that showed a significant association with body-weight traits, hap4 and hap5 (p < 0.05 and p < 0.01, respectively), differ by a single substitution at the g.1086C > T locus. The alleles at g.1086C > T act in an additive manner and explain a small percentage of the genetic variation of these phenotypes.

Conclusions

The association analysis revealed that g.1086C > T had a significant association with all body-weight traits under study. Although the SNP explains a small percentage of the variance, our results indicate that a variation in the 5′ flanking region of the myostatin gene is associated with the genetic regulation of growth in Atlantic salmon.

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 (Patinopectenyessoensis) 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.

Transcriptome sequencing of Zhikong scallop (Chlamys farreri) and comparative transcriptomic analysis with Yesso scallop (Patinopecten yessoensis)

Background

Bivalves play an important role in the ecosystems they inhabit and represent an important food source all over the world. So far limited genetic research has focused on this group of animals largely due to the lack of sufficient genetic or genomic resources. Here, we performed de novo transcriptome sequencing to produce the most comprehensive expressed sequence tag resource for Zhikong scallop (Chlamys farreri), and conducted the first transcriptome comparison for scallops.

Results

In a single 454 sequencing run, 1,033,636 reads were produced and then assembled into 26,165 contigs. These contigs were then clustered into 24,437 isotigs and further grouped into 20,056 isogroups. About 47% of the isogroups showed significant matches to known proteins based on sequence similarity. Transcripts putatively involved in growth, reproduction and stress/immune-response were identified through Gene ontology (GO) and KEGG pathway analyses. Transcriptome comparison with Yesso scallop (Patinopecten yessoensis) revealed similar patterns of GO representation. Moreover, 38 putative fast-evolving genes were identified through analyzing the orthologous gene pairs between the two scallop species. More than 46,000 single nucleotide polymorphisms (SNPs) and 350 simple sequence repeats (SSRs) were also detected.

Conclusion

Our study provides the most comprehensive transcriptomic resource currently available for C. farreri. Based on this resource, we performed the first large-scale transcriptome comparison between the two scallop species, C. farreri and P. yessoensis, and identified a number of putative fast-evolving genes, which may play an important role in scallop speciation and/or local adaptation. A large set of single nucleotide polymorphisms and simple sequence repeats were identified, which are ready for downstream marker development. This transcriptomic resource should lay an important foundation for future genetic or genomic studies on C. farreri.

Molecular characterization of TGF-β type I receptor gene (Tgfbr1) in Chlamys farreri, and the association of allelic variants with growth traits

BACKGROUND

Scallops are an economically important aquaculture species in Asian countries, and growth-rate improvement is one of the main focuses of scallop breeding. Investigating the genetic regulation of scallop growth could benefit scallop breeding, as such research is currently limited. The transforming growth factor beta (TGF-β) signaling through type I and type II receptors, plays critical roles in regulating cell proliferation and growth, and is thus a plausible candidate growth regulator in scallops.

RESULTS

We cloned and characterized the TGF-β type I receptor (Tgfbr1) gene from Zhikong scallops (Chlamys farreri). The deduced amino acid sequence contains characteristic residues and exhibits the conserved structure of Tgfbr1 proteins. A high expression level of scallop Tgfbr1 was detected during early embryonic stages, whereas Tgfbr1 expression was enriched in the gonad and striated muscle in adults. A single nucleotide polymorphism (SNP, c. 1815C>T) in the 3' UTR was identified. Scallops with genotype TT had higher growth traits values than those with genotype CC or CT in a full-sib family, and significant differences were found between genotypes CC and TT for shell length, shell height, and striated muscle weight. An expression analysis detected significantly more Tgfbr1 transcripts in the striated muscle of scallops with genotype CC compared to those with genotype TT or CT. Further evaluation in a population also revealed higher striated muscle weight in scallops with genotype TT than those with the other two genotypes. The inverse correlation between striated muscle mass and Tgfbr1 expression is consistent with TGF-β signaling having a negative effect on cell growth.

CONCLUSION

The scallop Tgfbr1 gene was cloned and characterized, and an SNP potentially associated with both scallop growth and Tgfbr1 expression was identified. Our results suggest the negative regulation of Tgfbr1 in scallop growth and provide a candidate marker for Zhikong scallop breeding.

Rheb, an activator of target of rapamycin, in the blackback land crab, Gecarcinus lateralis: cloning and effects of molting and unweighting on expression in skeletal muscle

Molt-induced claw muscle atrophy in decapod crustaceans facilitates exuviation and is coordinated by ecdysteroid hormones. There is a 4-fold reduction in mass accompanied by remodeling of the contractile apparatus, which is associated with an 11-fold increase in myofibrillar protein synthesis by the end of the premolt period. Loss of a walking limb or claw causes a loss of mass in the associated thoracic musculature; this unweighting atrophy occurs in intermolt and is ecdysteroid independent. Myostatin (Mstn) is a negative regulator of muscle growth in mammals; it suppresses protein synthesis, in part, by inhibiting the insulin/metazoan target of rapamycin (mTOR) signaling pathway. Signaling via mTOR activates translation by phosphorylating ribosomal S6 kinase (s6k) and 4E-binding protein 1. Rheb (Ras homolog enriched in brain), a GTP-binding protein, is a key activator of mTOR and is inhibited by Rheb-GTPase-activating protein (GAP). Akt protein kinase inactivates Rheb-GAP, thus slowing Rheb-GTPase activity and maintaining mTOR in the active state. We hypothesized that the large increase in global protein synthesis in claw muscle was due to regulation of mTOR activity by ecdysteroids, caused either directly or indirectly via Mstn. In the blackback land crab, Gecarcinus lateralis, a Mstn-like gene (Gl-Mstn) is downregulated as much as 17-fold in claw muscle during premolt and upregulated 3-fold in unweighted thoracic muscle during intermolt. Gl-Mstn expression in claw muscle is negatively correlated with hemolymph ecdysteroid level. Full-length cDNAs encoding Rheb orthologs from three crustacean species (G. lateralis, Carcinus maenas and Homarus americanus), as well as partial cDNAs encoding Akt (Gl-Akt), mTOR (Gl-mTOR) and s6k (Gl-s6k) from G. lateralis, were cloned. The effects of molting on insulin/mTOR signaling components were quantified in claw closer, weighted thoracic and unweighted thoracic muscles using quantitative polymerase chain reaction. Gl-Rheb mRNA levels increased 3.4-fold and 3.9-fold during premolt in claw muscles from animals induced to molt by eyestalk ablation (ESA) and multiple leg autotomy (MLA), respectively, and mRNA levels were positively correlated with hemolymph ecdysteroids. There was little or no effect of molting on Gl-Rheb expression in weighted thoracic muscle and no correlation of Gl-Rheb mRNA with ecdysteroid titer. There were significant changes in Gl-Akt, Gl-mTOR and Gl-s6k expression with molt stage. These changes were transient and were not correlated with hemolymph ecdysteroids. The two muscles differed in terms of the relationship between Gl-Rheb and Gl-Mstn expression. In thoracic muscle, Gl-Rheb mRNA was positively correlated with Gl-Mstn mRNA in both ESA and MLA animals. By contrast, Gl-Rheb mRNA in claw muscle was negatively correlated with Gl-Mstn mRNA in ESA animals, and no correlation was observed in MLA animals. Unweighting increased Gl-Rheb expression in thoracic muscle at all molt stages; the greatest difference (2.2-fold) was observed in intermolt animals. There was also a 1.3-fold increase in Gl-s6k mRNA level in unweighted thoracic muscle. These data indicate that the mTOR pathway is upregulated in atrophic muscles. Gl-Rheb, in particular, appears to play a role in the molt-induced increase in protein synthesis in the claw muscle.

Sequencing and analysis of four BAC clones containing innate immune genes from the Zhikong scallop (Chlamys farreri)

The sequencing of BAC clones (~100 kb) can reveal some characteristics of a genome that are challenging to obtain based on short sequences. Additionally, although the immune genes of the Zhikong scallop (Chlamys farreri) have been studied widely, few analyses have been conducted at the DNA level. In this study, four C. farreri BAC clones containing innate immune genes, including hsp70, l gbp (lipopolysaccharide and beta-1,3-glucan binding protein), serine protease and a gene with an immunoglobulin-like domain, were sequenced and analyzed both to explore the genomic characteristics of C. farreri based on long DNA sequences and to promote the study of C. farreri immune genes at the DNA level. The total length of the four BACs was 389.98 kb. A total of 34 genes were predicted in these sequences, and several features of protein-coding regions in the C. farreri genome were inferred based on this information. Two LGBP genes were located close together in a 22-kb region in one BAC clone, indicating the physical linkage of some immune genes in C. farreri. A cluster of membrane transport genes was also observed; these genes might play important roles in eliminating toxins in C. farreri, which lives as a filter feeder. Further analysis showed 15.43% of the BAC sequence was repetitive. Tandem repeats were the most abundant repeat type, followed by transposable elements. A total of 31 SSRs were predicted in the four BACs. An IS10 family transposon was identified, and a suspected regulatory non-coding RNA gene for this transposon (RNA-OUT) was observed to overlap with it complementarily. This work will promote future studies on the genomics, immune system and non-coding regions of C. farreri.