Zexia Gao

Transcriptome analysis and SSR/SNP markers information of the blunt snout bream (Megalobrama amblycephala).

Background Blunt snout bream (Megalobrama amblycephala) is an herbivorous freshwater fish species native to China and has been recognized as a main aquaculture species in the Chinese freshwater polyculture system with high economic value. Right now, only limited EST resources were available for M. amblycephala. Recent advances in large-scale RNA sequencing provide a fast, cost-effective, and reliable approach to generate large expression datasets for functional genomic analysis, which is especially suitable for non-model species with un-sequenced genomes. Methodology and Principal Findings Using 454 pyrosequencing, a total of 1,409,706 high quality reads (total length 577 Mbp) were generated from the normalized cDNA of pooled M. amblycephala individuals. These sequences were assembled into 26,802 contigs and 73,675 singletons. After BLAST searches against the NCBI non-redundant (NR) and UniProt databases with an arbitrary expectation value of E−10, over 40,000 unigenes were functionally annotated and classified using the FunCat functional annotation scheme. A comparative genomics approach revealed a substantial proportion of genes expressed in M. amblycephala tanscriptome to be shared across the genomes of zebrafish, medaka, tetraodon, fugu, stickleback, human, mouse, and chicken, and identified a substantial number of potentially novel M. amblycephala genes. A total number of 4,952 SSRs were found and 116 polymorphic loci have been characterized. A significant number of SNPs (25,697) and indels (23,287) were identified based on specific filter criteria in the M. amblycephala. Conclusions This study is the first comprehensive transcriptome analysis for a fish species belonging to the genus Megalobrama. These large EST resources are expected to be valuable for the development of molecular markers, construction of gene-based linkage map, and large-scale expression analysis of M. amblycephala, as well as comparative genome analysis for the genus Megalobrama fish species. The identified SSR and SNP markers will greatly benefit its breeding program and whole genome association studies.

Identification and characterization of microRNAs involved in growth of blunt snout bream (Megalobrama amblycephala) by Solexa sequencing

BackgroundBlunt snout bream (Megalobrama amblycephala) is an economically important fish species in the Chinese freshwater polyculture system for its delicacy and high economic value. MicroRNAs (miRNAs) play important roles in regulation of almost all biological processes in eukaryotes. Although previous studies have identified thousands of miRNAs from many species, little information is known for miRNAs of M. amblycephala. To investigate functions of miRNAs associated with growth of M. amblycephala, we adopted the Solexa sequencing technology to sequence two small RNA libraries prepared from four growth related tissues (brain, pituitary, liver and muscle) of M. amblycephala using individuals with relatively high and low growth rates.ResultsIn this study, we have identified 347 conserved miRNAs (belonging to 123 families) and 22 novel miRNAs in M. amblycephala. Moreover, we observed sequence variants and seed edits of the miRNAs. Of the 5,166 single nucleotide substitutions observed in two libraries, the most abundant were G-to-U (15.9%), followed by U-to-C (12.1%), G-to-A (11.2%), and A to G (11.2%). Subsequently, we compared the expression patterns of miRNAs in the two libraries (big-size group with high growth rate versus small-size group with low growth rate). Results indicated that 27 miRNAs displayed significant differential expressions between the two libraries (p < 0.05). Of these, 16 were significantly up-regulated and 11 were significantly down-regulated in the big-size group compared to the small-size group. Furthermore, stem-loop RT-PCR was applied to validate and profile the expression of the differentially expressed miRNAs in ten tissues, and the result revealed that the conserved miRNAs expressed at higher levels than the novel miRNAs, especially in brain, liver and muscle. Also, targets prediction of differentially expressed miRNAs and KEGG pathway analysis suggested that differentially expressed miRNAs are involved in growth and metabolism, signal transduction, cell cycle, neural development and functions.ConclusionsThe present study provides the first large-scale characterization of miRNAs in M. amblycephala and miRNA profile related to different growth performances. The discovery of miRNA resource from this study is expected to contribute to a better understanding of the miRNAs roles playing in regulating the growth biological processes and the study of miRNA function and phenotype-associated miRNA identification in fish.

The seahorse genome and the evolution of its specialized morphology

Seahorses have a specialized morphology that includes a toothless tubular mouth, a body covered with bony plates, a male brood pouch, and the absence of caudal and pelvic fins. Here we report the sequencing and de novo assembly of the genome of the tiger tail seahorse, Hippocampus comes. Comparative genomic analysis identifies higher protein and nucleotide evolutionary rates in H. comes compared with other teleost fish genomes. We identified an astacin metalloprotease gene family that has undergone expansion and is highly expressed in the male brood pouch. We also find that the H. comes genome lacks enamel matrix protein-coding proline/glutamine-rich secretory calcium-binding phosphoprotein genes, which might have led to the loss of mineralized teeth. tbx4, a regulator of hindlimb development, is also not found in H. comes genome. Knockout of tbx4 in zebrafish showed a ‘pelvic fin-loss’ phenotype similar to that of seahorses.

Molecular cloning and expression analysis of major histocompatibility complex class I, IIA and IIB genes of blunt snout bream (Megalobrama amblycephala).

Major histocompatibility complex (MHC) plays an important role in the immune response of vertebrates. In this study, we isolated MHC class IIA and IIB genes from blunt snout bream (Megalobrama amblycephala) by rapid amplification of cDNA ends polymerase chain reaction (RACE-PCR). In order to study the function of the MHC genes in M. amblycephala, tissue distribution and immune response of the MHC genes to bacterial challenge were analyzed. All the characteristic features of MHC class II chain structure could be identified in the deduced amino sequences of MHC IIA and IIB, including the leader peptide, α1/β1 and α2/β2 domains, connecting peptide and transmembrane and cytoplasmic regions, as well as conserved cysteines and N-glycosylation site. The deduced amino acid sequence of the MHC IIA and IIB molecules shared from 48% to 88% and from 65% to 77% similarity with those of other teleosts, respectively. Quantitative real-time PCR (qRT-PCR) demonstrated that MHC I and II genes were ubiquitously expressed in ten tissues, with high level in immune related tissues, including kidney, intestine, gill and spleen. Challenge of M. amblycephala with the extracellular pathogen, Aeromonas hydrophila, resulted in a significant increase in the expression of MHC I, MHC IIA and IIB mRNA within 72 h after infection in gill, kidney, intestine and liver, followed by a recovery to normal level after 120 h. The changes of expression levels for MHC IIA and IIB in most tissues were significantly higher than that of MHC I in the corresponding tissues at most time points (P<0.05). These results demonstrated the MHC genes played an important role in response to bacterial infection in M. amblycephala; however, MHC class I and II genes showed different functional activity, which need be further investigated in teleost.

Identification of MicroRNA for Intermuscular Bone Development in Blunt Snout Bream (Megalobrama amblycephala)

Intermuscular bone (IB), which occurs only in the myosepta of the lower teleosts, is attracting more attention of researchers due to its particular development and lack of genetic information. MicroRNAs (miRNAs) are emerging as important regulators for biological processes. In the present study, miRNAs from IBs and connective tissue (CT; encircled IBs) from six-month-old Megalobrama amblycephala were characterized and compared. The results revealed the sequences and expression levels of 218 known miRNA genes (belonging to 97 families). Of these miRNAs, 44 known microRNA sequences exhibited significant expression differences between the two libraries, with 24 and 20 differentially-expressed miRNAs exhibiting higher expression in the CT and IBs libraries, respectively. The expressions of 11 miRNAs were selected to validate in nine tissues. Among the high-ranked predicted gene targets, differentiation, cell cycle, metabolism, signal transduction and transcriptional regulation were implicated. The pathway analysis of differentially-expressed miRNAs indicated that they were abundantly involved in regulating the development and differentiation of IBs and CT. This study characterized the miRNA for IBs of teleosts for the first time, which provides an opportunity for further understanding of miRNA function in the regulation of IB development.

Molecular characterization of trypsinogens and development of trypsinogen gene expression and tryptic activities in grass carp (Ctenopharyngodon idellus) and topmouth culter (Culter alburnus)

This study examined the gene structures and expression of trypsinogens, as well as the trypsin activities of the grass carp Ctenopharyngodon idellus (herbivorous) and the topmouth culter Culter alburnus (carnivorous), which are commercially important freshwater species of the family Cyprinidae in China. Isolated full-length trypsinogen cDNA clones were 869 bp and 857 bp. The deduced amino acid sequences were 242aa and 247aa long, both containing the highly conserved residues essential for serine protease catalytic and conformational maintenance. The results from isoelectric and phylogenetic analyses suggest that grass carp trypsinogen is grouped with teleost trypsinogen group I, while topmouth culter trypsinogen is grouped with group II. The expression pattern of trypsinogen mRNA was similar between these two species, appearing 2 days post-hatching (dph) and reaching peaks at 11 and 23 dph. The trypsin-specific activities in both species were detected 2 dph and reached the major peaks at 8 dph, however the minor peaks were observed at 20 dph in the grass carp and 17 dph in the topmouth culter. The trypsin-specific activity was significantly higher in the grass carp than in the topmouth culter, which may be attributed to the nature of their different nutritional habits.

The draft genome of blunt snout bream (Megalobrama amblycephala) reveals the development of intermuscular bone and adaptation to herbivorous diet

Abstract The blunt snout bream Megalobrama amblycephala is the economically most important cyprinid fish species. As an herbivore, it can be grown by eco-friendly and resource-conserving aquaculture. However, the large number of intermuscular bones in the trunk musculature is adverse to fish meat processing and consumption. As a first towards optimizing this aquatic livestock, we present a 1.116-Gb draft genome of M. amblycephala, with 779.54 Mb anchored on 24 linkage groups. Integrating spatiotemporal transcriptome analyses, we show that intermuscular bone is formed in the more basal teleosts by intramembranous ossification and may be involved in muscle contractibility and coordinating cellular events. Comparative analysis revealed that olfactory receptor genes, especially of the beta type, underwent an extensive expansion in herbivorous cyprinids, whereas the gene for the umami receptor T1R1 was specifically lost in M. amblycephala. The composition of gut microflora, which contributes to the herbivorous adaptation of M. amblycephala, was found to be similar to that of other herbivores. As a valuable resource for the improvement of M. amblycephala livestock, the draft genome sequence offers new insights into the development of intermuscular bone and herbivorous adaptation.

Dynamic mRNA and miRNA expression analysis in response to intermuscular bone development of blunt snout bream (Megalobrama amblycephala)

Intermuscular bone (IB), which occurs only in the myosepta of lower teleosts, is attracting more attention because they are difficult to remove and make the fish unpleasant to eat. By gaining a better understanding of the genetic regulation of IB development, an integrated analysis of miRNAs and mRNAs expression profiling was performed on Megalobrama amblycephala. Four key development stages were selected for transcriptome and small RNA sequencing. A number of significantly differentially expressed miRNAs/genes associated with bone formation and differentiation were identified and the functional characteristics of these miRNAs/genes were revealed by GO function and KEGG pathway analysis. These were involved in TGF-β, ERK and osteoclast differentiation pathways known in the literature to affect bone formation and differentiation. MiRNA-mRNA interaction pairs were detected from comparison of expression between different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating bone development and differentiation. A negative regulation effect of two miRNAs was verified through dual luciferase reporter assay. As a unique public resource for gene expression and regulation during the IB development, this study is expected to provide forwards ideas and resources for further biological researches to understand the IBs’ development.

Leptin Genes in Blunt Snout Bream: Cloning, Phylogeny and Expression Correlated to Gonads Development

To investigate the leptin related genes expression patterns and their possible function during the gonadal development in fish, the cDNA and genomic sequences of leptin, leptin receptor (leptinR), and leptin receptor overlapping transcript like-1 (leprotl1) were cloned and their expression levels were quantified in the different gonadal development stages of Megalobrama amblycephala. The results showed that the full length cDNA sequences of leptin, leptinR and leprotl1 were 953, 3432 and 1676 bp, coding 168, 1082, and 131 amino acid polypeptides, and the genomic sequences were 1836, 28,528 and 5480 bp, which respectively had 3, 15 and 4 exons, respectively. The phylogenetic analysis revealed that three genes were relatively conserved in fish species. Quantitative real-time PCR results showed that the three genes were ubiquitously expressed in all examined tissues during the different gonadal development stages. The leptin and leptinR took part in the onset of puberty, especially in female M. amblycephala, by increasing the expression levels in brain during the stage I to III of ovary. The expression levels of leptin and leptinR had significant differences between male and female in hypothalamic-pituitary-gonadal (HPG) axis tissues (p < 0.05). The leptinR had the same variation tendency with leptin, but the opposite changes of expression levels were found in leprotl1, which may resist the expression of leptinR for inhibiting the function of leptin in target organ. These findings revealed details about the possible role of these genes in regulating gonadal maturation in fish species.

A simple and affordable method of DNA extraction from fish scales for polymerase chain reaction.

Most fish studies today are performed at the molecular level, which needs reliable, quick, and cost-effective DNA extraction protocols (Liu et al. 2009; Ming et al. 2010; Shen et al. 2010). Generally, the most widely used materials for DNA extraction in fish are fins or red blood cell samples (Yang et al. 2001; Liu et al. 2009; Zhang et al. 2009; Li et al. 2009). As fins are used by fish in locomotion, their removal would adversely impact survival of sampled fish. Similarly, collection of samples of enough red blood for DNA extraction from small fish would require them to be cultured further, and this would also be difficult for the researcher. Thus, simpler methods of sample collection are needed for DNA extraction. Since collection of fish scale samples is easy and noninvasive (Kumar et al. 2007), there is a need for a simple, rapid, and affordable method of DNA extraction from them. Based on a literature search, only five papers have described methods of DNA isolation from fish scales (Estoup et al. 1996; Nelson et al. 1998; Adcock et al. 2000; Yue and Orban 2001; Kumar et al. 2007). The method described by Nelson et al. (1998) is relatively complicated for DNA isolation from a large number of individuals. Adcock et al. (2000) successfully isolated DNA from old dried scales by the traditional phenol–chloroform procedure. Yue and Orban (2001) found that