Characterization of two paralogous myostatin genes and evidence for positive selection in Tibet fish: Gymnocypris przewalskii

Identification of two MEF2s and their role in inhibiting the transcription of the mstn2a gene in the yellowfin seabream, Acanthopagrus latus (Hottuyn, 1782)

Myocyte-specific enhancer binding factor 2 (MEF2), which belongs to the MADS superfamily, is a pivotal and conserved transcription factor that combines with the E-box motif to control the expression of muscle genes. Myostatin (mstn), a muscle growth inhibitor, is a vital member of the TGF-β superfamily. Currently, an understanding of the mechanisms of A. latus mstn (Almstn) transcriptional regulation mediated by MEF2 in fish muscle development is lacking. In the present study, two AlMEF2s (AlMEF2A and AlMEF2B) and Almstn2a were characterized from Acanthopagrus latus. AlMEF2A and AlMEF2B had 456 and 315 amino acid (aa) residues, respectively. Two typical regions, a MADS-box, MEF2, and transcriptionally activated (TAD) domains, are present in both AlMEF2s. The expression profiles of the two AlMEF2 genes were similar. The AlMEF2 genes were mainly expressed in the brain, white muscle, and liver, while Almstn2a expression was higher in the brain than in other tissues. Moreover, the expression trends of AlMEF2s and Almstn2a were significantly changed after starvation and refeeding in the five groups. Additionally, truncation experiments showed that -987 to +168 and -105 to +168 were core promoters of Almstn2a that responded to AlMEF2A and AlMEF2B, respectively. The point mutation experiment confirmed that Almstn2a transcription relies on the mutation binding sites 1 or 5 (M1/5) and mutation binding sites 4 or 5 (M4/5) for AlMEF2A and AlMEF2B regulation, respectively. The electrophoretic mobile shift assay (EMSA) further verified that M1 (-527 to -512) was a pivotal site where AlMEF2A acted on the Almstn2a gene. Furthermore, a siRNA interference gene expression experiment showed that reduced levels of AlMEF2A or AlMEF2B could prominently increase Almstn2a transcription. These results provide new information about the regulation of Almstn2a transcriptional activity by AlMEF2s and a theoretical basis for the regulatory mechanisms involved in muscle development in fish.

Functional Characterization of the Almstn2 Gene and Its Association with Growth Traits in the Yellowfin Seabream Acanthopagrus latus (Hottuyn, 1782)

Myostatin (mstn), also known as GDF8, is a growth and differentiation factor of the transforming growth factor-β (TGF-β) superfamily and plays a key inhibitory effect in the regulation of skeletal muscle development and growth in vertebrates. In the present study, to comprehend the role of the mstn2 gene of the yellowfin seabream Acanthopagrus latus (Almstn2b), the genomic sequence of Almstn2b is 2359 bp, which encodes 360 amino acids and is composed of three exons and two introns, was obtained. Two typical regions, a TGF-β propeptide and TGF-β domain, constitute Almstn2b. The topology indicated that Almstn2 was grouped together with other Perciformes, such as the gilthead seabream Sparus aurata. Moreover, Almstn2b was mainly expressed in the brain, fins, and spleen. Furthermore, five SNPs, one in the exons and four in the introns, were identified in the Almstn2b gene. The allele and genotype frequencies of SNP-Almstn2b +1885 A/G were significantly related to the total weight, interorbital distance, stem length, tail length, caudal length, caudal height, body length, and total length (p < 0.05). The allele and genotype frequencies of SNP-Almstn2b +1888 A/G were significantly related to the weight, interorbital distance, long head behind the eyes, body height, tail length, caudal length, and body length. Additionally, the relationship between the SNP-Almstn2b +1915 A/G locus and weight and long head behind the eyes was significant (p < 0.05). Furthermore, the other two SNPs were not significantly associated with any traits. Thus, the SNPs identified in this study could be utilized as candidate SNPs for breeding and marker-assisted selection in A. latus.

Genetic Diversity and Population Differentiation of Naked Carp (Gymnocypris przewalskii) Revealed by Cytochrome Oxidase Subunit I and D-Loop

The genetic diversity in the naked carp (Gymnocypris przewalskii) of China is threatened by climate change, human activities, as well as natural factors, eliciting conservation concerns. To explore the genetic aspects of G. przewalskii, the genetic diversity, genetic structure, population differentiation, and historical demography of 566 representative individuals from seven geographically distinct ranges of Qinghai Lake were evaluated by mitochondrial DNA cytochrome oxidase subunit I (COI) and D-loop sequences. Estimates of genetic parameters showed that the seven populations of G. przewalskii had high levels of haplotype diversity (0.50243–0.94620) and low levels of nucleotide diversity (0.00079–0.00624). Haplotype genealogy indicated there was no obvious phylogenetic pattern between haplotypes. Both markers denoted the absence of population genetic structure [the genetic differentiation coefficient F-statistics (Fst) &lt; 0] and the presence of high genetic flow (COI: 0.9731–1.0441; D-loop: 0.9480–1.0398). The mismatch between the distribution and neutrality tests supported the evidence of population expansion, which occurred during the late middle Pleistocene [COI: 0.36–0.108 MYA (Million Years Ago); D-loop: 0.497–0.165 MYA]. Furthermore, this work illustrated two simple, reliable, and inexpensive molecular markers for analysis of genetic diversity, while the sensitivity of the mitochondrial D-loop region as a reflection of genetic diversity in G. przewalskii is higher than that of the COI gene.

Riverbed Substrate Requirements for Natural Reproduction of Gymnocypris przewalskii

Gymnocypris przewalskii (i.e., Qinghai Lake naked carp) is a migratory fish species that lives in highland brackish water. It is important to understand the abiotic environment required by this fish to reproduce naturally so that its habitat can be protected and the wild population can be conserved. Here, artificial simulation and spawning ground substrate transformation experiments were conducted to examine the riverbed substrate requirements for G. przewalskii to naturally reproduce. Using various techniques (in vitro markers, videography, and Ethovision XT behavior tracking), this study systematically investigated the riverbed substrate preferences of G. przewalskii as well as the characteristics and effectiveness of natural reproduction induced by pebble riverbed substrate. The findings can be summarized as follows: (1) the habitat preferences of G. przewalskii differed significantly between various riverbed substrate, with pebble substrate being preferred during natural reproduction, and sand substrate being preferred pre- and post-spawning, and (2) the natural reproduction of G. przewalskii was heavily reliant on pebble riverbed substrate. Specifically, pebble substrate significantly improved spawn quantity and fertilization rate. These findings provide scientific evidence for the improvement and restoration of G. przewalskii spawning grounds, and insights regarding the artificial bionic reproduction of G. przewalskii.

Characterization of two myostatin genes in pufferfish Takifugu bimaculatus: sequence, genomic structure, and expression

Myostatin (MSTN) is a negative regulator of muscle growth, which restrains the proliferation and differentiation of myoblasts. To understand the role of two mstn genes of Takifugu bimaculatus, the full-length cDNAs of 1131 bp Tbmstn1 and 1,080 bp Tbmstn2 were obtained from the T. bimaculatus' genomic database, which encodes 376 and 359 amino acids, respectively. The results of qRT-PCR showed that Tbmstn1 was expressed in the eye, kidney, spleen, skeletal muscle, gill, and brain, and the expression level in the skeletal muscle was extremely significantly higher than in other examined tissues. Tbmstn2 was expressed in the skin, skeletal muscle, gill, and brain, and had the highest expression in the skeletal muscle, followed by expression in the brain. Meanwhile, in different stages of embryonic development, the expression of Tbmstn1 started from the gastrula stage. Its expression in the eye-pigment formation stage and hatching stage was significantly higher than that in other stages. The Tbmstn2 was expressed in all examined embryonic stages with different levels, and the highest expression was detected in the eye-pigment formation stage. These results suggested that Tbmstn1 and Tbmstn2 may involve in the development of skeletal muscle, and Tbmstn2 may be related to the formation of nervous system.

Transcriptomic signature of rapidly evolving immune genes in a highland fish

Recent genome-wide studies have begun to elucidate the genomic basis of hypoxia, long-term cold and high saline and alkaline adaptation in highland fish, and a number of key genes contributed to its highland adaptation were identified. An increasing number of studies indicated that immune genes of Tibetan endemic fish species underwent positive selection towards functional shift, while the insight into immune gene repertoire of Tibetan highland fishes from genome-wide studies has largely lagged behind. In this study, we performed one of the first comparative genomics study in particular focusing on the signatures of immune genes in a highland fish, Gymnocypris przewalskii based on immune-relevant tissue transcriptome assemblies. We identified seven putative rapidly evolving immune genes with elevated molecular evolutionary rate (dN/dS) relative to lowland fish species. Using tissue-transcriptome data, we found most of rapidly evolving immune genes were broadly expressed in head-kidney, spleen, gills and skin tissues, which significantly enriched for complement activation and inflammatory response processes. In addition, we found a set of complement activation related genes underwent accelerated evolution and showed consistently repressed expression patterns in response to parasite Ichthyophthirius multifiliis infection. Moreover, we detected a number of immune genes involved in adaptive immune system exhibited distinct signature of upregulated expression patterns after parasite infection. Taken together, this study provided putative transcriptomic signatures of rapidly evolving immune genes, and will gain the insight into Schizothoracine fish adaptation to high-altitude extreme aquatic environments including diversified pathogen challenge.

Total substitution of dietary fish oil by vegetable oils stimulates muscle hypertrophic growth in Senegalese sole and the upregulation of fgf6

The long term effects of fish oil (FO) substitution by increasing the levels of vegetable oils (VO), 0% (CTR), 50% (VO50) and 100% (VO100), in diets for Senegalese sole were evaluated in terms of skeletal muscle cellularity and expression of related genes. After 140 days of feeding, all fish had similar body weight and length. The inclusion of 50% VO did not result in differences in muscle cellularity, but dorsal muscle cross-sectional area and fast-twitch fibre diameter increased in fish fed total FO substitution, whilst fibre density was reduced (P < 0.05) in relation to CTR. The total number of fibres was similar in all treatments. FO substitution did not affect the transcript levels of myogenic genes (myf5, mrf4, myog, myod1, myod2), but resulted in a two-fold increase of fgf6 transcript levels compared to CTR (P < 0.05). The relative expression of igf-I was higher in VO100 than in VO50, but was similar to CTR. FO substitution resulted in cellularity changes related to the stimulation of muscle hypertrophic growth, but not hyperplastic growth, and associated with a nutritional modulation of fgf6 by dietary VO. This study indicates that 50% VO does not affect the muscle phenotype, but total FO substitution stimulates muscle hypertrophy.

The microRNA repertoire of Tibetan naked carp Gymnocypris przewalskii: A case study in Schizothoracinae fish on the Tibetan Plateau

Tibetan naked carp Gymnocypris przewalskii is an ideal model system to study highland adaptation of fish, because it evolved specific genetic and phenotypic characteristics to adapt to chronic cold and alkaline environments in Lake Qinghai. MicroRNAs (miRNAs) are small noncoding RNAs that regulating gene expression post-transcriptionally in a wide range of biological processes. In this study, we focus on the role of miRNAs in adaptation of G. przewalskii to extreme conditions in Lake Qinghai. We generate the first miRNAome of G. przewalskii in Schizothoracinae fish. Using several genomic resources, we inferred 341 conserved miRNAs belonged to 152 miRNA families and 43 novel miRNAs in G. przewalskii, and also identified 15 teleost-specific miRNAs. Using a large scale of conserved miRNAs, we constructed a high-confidence phylogenetic tree between teleost and mammals than mitochondria and nuclear genes. In addition, we found that several miRNA family (e.g. miR-10 and let-7) members highly expressed in G. przewalskii, which may function in multiple biological processes. Finally, we predicted a total of 34,258 miRNA targets genes. Conserved miRNAs target genes participating in signal transduction, cell differentiation and biosynthetic process, and showed signature of functional constraint. While novel miRNAs in a species displayed species-specific targets and involved in ion binding, transport and oxidoreductase activity, may affect the expression patterns of targets with signature of gene family expansion or positive selection under extreme environment. Taken together, this study demonstrated that miRNAs may involve into roles of adaptation of G. przewalskii to highland aquatic environment, and also provide insights into miRNA regulatory network in Schizothoracinae fish as a case study.

Positive Darwinian selection within interferon regulatory factor genes of Gymnocypris przewalskii (Cyprinidae) on the Tibetan Plateau

Tibetan Plateau (TP) had experienced phased uplift, resulting in inhospitable environment of low temperature, hypoxia and high ultraviolet radiation for Tibetan wildlife. Many organisms can well adapt to TP, it is of ecological and evolutionary interest to untangle how organisms adapt to extreme environment on TP through evolution. Previous studies mainly focused on hypoxia and metabolism related genes, but we know little about the evolutionary history of immune genes in Tibetan wildlife. In this study, we first identified 10 interferon regulatory factor (IRF) genes from Tibetan naked carp Gymnocypris przewalskii. Within this gene family, IRF3, IRF5, IRF7 and IRF8 contained positive selection sites. Evidences indicated that positive selection may lead to IRF genes functional alternations, presumably driving genes towards adaptation to the environmental changes. Taken together, our results suggested 4 candidate genes as interesting targets for further experimental confirmation of their functional variations and contributions to high altitude adaptation in Tibet fish.

Transcriptome profiling analysis of naked carp (Gymnocypris przewalskii) provides insights into the immune-related genes in highland fish

The naked carp, Gymnocypris przewalskii, is one of the dominant aquaculture fish species in Qinghai Province, China. Its wild stocks have severely suffered from overfishing, and the farming species are vulnerable to various pathogens infections. Here we report the first immune-related tissues transcriptome of a wild naked carp using a deep sequencing approach. A total of 158,087 unigenes are generated, 2687 gill-specific gene and 3215 kidney-specific genes are identified, respectively. Gene ontology analysis shows that 51,671 unigenes are involved in three major functional categories: biological process, cellular component, and molecular function. Further analysis shows that numerous consensus sequences are homologous to known immune-related genes. Pathways mapping annotate 56,270 unigenes and identify a large number of immune-related pathways. In addition, we focus on the immune-related genes and gene family in Toll-like receptor signaling pathway involved in innate immunity, including toll-like receptors (TLRs), interferon regulatory factors (IRFs), interleukins (ILs) and tumor necrosis factors (TNFs). Eventually, we identify 5 TLRs, 4 IRFs, 3 ILs and 2 TNFs with a completed coding sequence though mining the transcriptome data. Phylogeny analysis shows these genes of naked carp are mostly close to zebrafish. Protein domain and selection pressure analyses together show that all these genes are highly conserved in gene sequence and protein domain structure with other species, and purifying selection underwent in these genes, implied functionally important features are conserved in the genes above. Intriguingly, we detect positive selection signals in naked carp TLR4, and significant divergence occurred among tested species TLR4, suggested that naked carp TLR4 function may be affected. Finally, we identify 23,867 simple sequence repeat (SSR) marks in this transcriptome. Taken together, this study not only contributes a large number of candidate genes in naked carp immunity, and also helps improve current understanding of immunogenetics basis and evolutionary history of immune related genes and gene family in highland fish species.