Molecular characterization and expression patterns of two hormone-sensitive lipase genes in common carp Cyprinus carpio.

The hormone-sensitive lipase (HSL) gene plays an important role in mammals' lipid metabolism. Therefore, its function in fish is capturing increasing attention. In this study, two distinct cDNAs, designated HSL1 and HSL2, are firstly identified from common carp Cyprinus carpio. The full-length cDNA of HSL1 and HSL2 consists of 3379 bp and 2732 bp, encoding polypeptide of 693 and 847 amino acids, respectively, and shares 60.6% amino acid identity. Phylogenetic analysis suggests that HSL1 and HSL2 are derived from paralogous genes, which might have arisen during a teleost-specific genome duplication event. The two HSL mRNAs are differentially expressed, both in terms of distribution among tissues and in terms of abundance during embryogenesis. Moreover, both HSL mRNAs are expressed in various tissues, the highest in abdominal fat. Meanwhile, the two HSLs are detected at all stages of embryonic development, suggesting that they could be functional and involved in embryogenesis. In addition, the results show that the mRNA expression level of HSL2 in the high group of intramuscular fat content is significantly higher than that in the low group (P < 0.01). The research provides basic data for developing a further understanding of the function of HSL as well as molecular regulation mechanism in fat metabolism of common carp.

Correction to: L_RNA_scaffolder: scaffolding genomes with transcripts.

Following the publication of this article [1], the authors reported that the link to the software described in the article is no longer valid.

QTL variations for growth-related traits in eight distinct families of common carp (Cyprinus carpio).

BACKGROUND: Comparing QTL analyses of multiple pair-mating families can provide a better understanding of important allelic variations and distributions. However, most QTL mapping studies in common carp have been based on analyses of individual families. In order to improve our understanding of heredity and variation of QTLs in different families and identify important QTLs, we performed QTL analysis of growth-related traits in multiple segregating families. RESULTS: We completed a genome scan for QTLs that affect body weight (BW), total length (TL), and body thickness (BT) of 522 individuals from eight full-sib families using 250 microsatellites evenly distributed across 50 chromosomes. Sib-pair and half-sib model mapping identified 165 QTLs on 30 linkage groups. Among them, 10 (genome-wide P <0.01 or P < 0.05) and 28 (chromosome-wide P < 0.01) QTLs exhibited significant evidence of linkage, while the remaining 127 exhibited a suggestive effect on the above three traits at a chromosome-wide (P < 0.05) level. Multiple QTLs obtained from different families affect BW, TL, and BT and locate at close or identical positions. It suggests that same genetic factors may control variability in these traits. Furthermore, the results of the comparative QTL analysis of multiple families showed that one QTL was common in four of the eight families, nine QTLs were detected in three of the eight families, and 26 QTLs were found common to two of the eight families. These common QTLs are valuable candidates in marker-assisted selection. CONCLUSION: A large number of QTLs were detected in the common carp genome and associated with growth-related traits. Some of the QTLs of different growth-related traits were identified at similar chromosomal regions, suggesting a role for pleiotropy and/or tight linkage and demonstrating a common genetic basis of growth trait variations. The results have set up an example for comparing QTLs in common carp and provided insights into variations in the identified QTLs affecting body growth. Discovery of these common QTLs between families and growth-related traits represents an important step towards understanding of quantitative genetic variation in common carp.

Development and evaluation of the first high-throughput SNP array for common carp (Cyprinus carpio).

BACKGROUND: A large number of single nucleotide polymorphisms (SNPs) have been identified in common carp (Cyprinus carpio) but, as yet, no high-throughput genotyping platform is available for this species. C. carpio is an important aquaculture species that accounts for nearly 14% of freshwater aquaculture production worldwide. We have developed an array for C. carpio with 250,000 SNPs and evaluated its performance using samples from various strains of C. carpio. RESULTS: The SNPs used on the array were selected from two resources: the transcribed sequences from RNA-seq data of four strains of C. carpio, and the genome re-sequencing data of five strains of C. carpio. The 250,000 SNPs on the resulting array are distributed evenly across the reference C.carpio genome with an average spacing of 6.6 kb. To evaluate the SNP array, 1,072 C. carpio samples were collected and tested. Of the 250,000 SNPs on the array, 185,150 (74.06%) were found to be polymorphic sites. Genotyping accuracy was checked using genotyping data from a group of full-siblings and their parents, and over 99.8% of the qualified SNPs were found to be reliable. Analysis of the linkage disequilibrium on all samples and on three domestic C.carpio strains revealed that the latter had the longer haplotype blocks. We also evaluated our SNP array on 80 samples from eight species related to C. carpio, with from 53,526 to 71,984 polymorphic SNPs. An identity by state analysis divided all the samples into three clusters; most of the C. carpio strains formed the largest cluster. CONCLUSIONS: The Carp SNP array described here is the first high-throughput genotyping platform for C. carpio. Our evaluation of this array indicates that it will be valuable for farmed carp and for genetic and population biology studies in C. carpio and related species.

The transcriptomes of the crucian carp complex (Carassius auratus) provide insights into the distinction between unisexual triploids and sexual diploids.

Both sexual reproduction and unisexual reproduction are adaptive strategies for species survival and evolution. Unisexual animals have originated largely by hybridization, which tends to elevate their heterozygosity. However, the extent of genetic diversity resulting from hybridization and the genomic differences that determine the type of reproduction are poorly understood. In Carassius auratus, sexual diploids and unisexual triploids coexist. These two forms are similar morphologically but differ markedly in their modes of reproduction. Investigation of their genomic differences will be useful to study genome diversity and the development of reproductive mode. We generated transcriptomes for the unisexual and sexual populations. Genes were identified using homology searches and an ab initio method. Estimation of the synonymous substitution rate in the orthologous pairs indicated that the hybridization of gibel carp occurred 2.2 million years ago. Microsatellite genotyping in each individual from the gibel carp population indicated that most gibel carp genes were not tri-allelic. Molecular function and pathway comparisons suggested few gene expansions between them, except for the progesterone-mediated oocyte maturation pathway, which is enriched in gibel carp. Differential expression analysis identified highly expressed genes in gibel carp. The transcriptomes provide information on genetic diversity and genomic differences, which should assist future studies in functional genomics.

L_RNA_scaffolder: scaffolding genomes with transcripts.

BACKGROUND: Generation of large mate-pair libraries is necessary for de novo genome assembly but the procedure is complex and time-consuming. Furthermore, in some complex genomes, it is hard to increase the N50 length even with large mate-pair libraries, which leads to low transcript coverage. Thus, it is necessary to develop other simple scaffolding approaches, to at least solve the elongation of transcribed fragments. RESULTS: We describe L_RNA_scaffolder, a novel genome scaffolding method that uses long transcriptome reads to order, orient and combine genomic fragments into larger sequences. To demonstrate the accuracy of the method, the zebrafish genome was scaffolded. With expanded human transcriptome data, the N50 of human genome was doubled and L_RNA_scaffolder out-performed most scaffolding results by existing scaffolders which employ mate-pair libraries. In these two examples, the transcript coverage was almost complete, especially for long transcripts. We applied L_RNA_scaffolder to the highly polymorphic pearl oyster draft genome and the gene model length significantly increased. CONCLUSIONS: The simplicity and high-throughput of RNA-seq data makes this approach suitable for genome scaffolding. L_RNA_scaffolder is available at http://www.fishbrowser.org/software/L_RNA_scaffolder.

Dynamic RBM47 ISGylation confers broad immunoprotection against lung injury and tumorigenesis via TSC22D3 downregulation.

ISGylation is a well-established antiviral mechanism, but its specific function in immune and tissue homeostasis regulation remains elusive. Here, we reveal that the RNA-binding protein RBM47 undergoes phosphorylation-dependent ISGylation at lysine 329 to regulate immune activation and maintain lung homeostasis. K329R knockin (KI) mice with defective RBM47-ISGylation display heightened susceptibility to LPS-induced acute lung injury and lung tumorigenesis, accompanied with multifaceted immunosuppression characterized by elevated pro-inflammatory factors, reduced IFNs/related chemokines, increased myeloid-derived suppressor cells, and impaired tertiary lymphoid structures. Mechanistically, RBM47-ISGylation regulation of the expression of TSC22D3 mRNA, a glucocorticoid-inducible transcription factor, partially accounts for the effects of RBM47-ISGylation deficiency due to its broad immunosuppressive activity. We further demonstrate the direct inhibitory effect of RBM47-ISGylation on TSC22D3 expression in human cells using a nanobody-targeted E3 ligase to induce site-specific ISGylation. Furthermore, epinephrine-induced S309 phosphorylation primes RBM47-ISGylation, with epinephrine treatment exacerbating dysregulated cytokine expression and ALI induction in K329R KI mice. Our findings provide mechanistic insights into the dynamic regulation of RBM47-ISGylation in supporting immune activation and maintaining lung homeostasis.

Evaluation of qPCR reference genes for taimen (Hucho taimen) under heat stress.

As a powerful and attractive method for detecting gene expression, qRT-PCR has been broadly used in aquaculture research. Understanding the biology of taimen (Hucho taimen) has drawn increasing interest because of its ecological and economic value. Stable reference genes are required for the reliable quantification of gene expression, but such genes have not yet been optimized for taimen. In this study, the stability levels of 10 commonly used candidate reference genes were evaluated using geNorm, NormFinder, BestKeeper, and RefFinder. The expression levels of the 10 genes were detected using 240 samples from 48 experimental groups consisting of 40 individuals treated under four heat-stress conditions (18, 20, 22, and 24 °C) for 24 h and 26 °C for 4, 24, 48, and 72 h. Six tissues (blood, heart, brain, gill, skin, and liver) were collected from each individual. Ribosomal protein S29 (RPS29) and ribosomal protein L19 (RPL19) were the most stable genes among all of the samples, whereas 28S ribosomal RNA (28S rRNA), attachment region binding protein (ARBP), and 18S ribosomal RNA (18S rRNA) were the least stable. These results were verified by an expression analysis of taimen heat-stress genes (heat shock protein 60, hsp60, and heat shock protein 70, hsp70). In conclusion, RPS29 and RPL19 are the optimal reference genes for qRT-PCR analyses of taimen, irrespective of the tissue and experimental conditions. These results allow the reliable study of gene expression in taimen.

In Silico Screening and Development of Microsatellite Markers for Genetic Analysis in Perca fluviatilis.

Perca fluviatilis is an economically important species of freshwater fish. To understand the genetic structure of P. fluviatilis in China, 268 samples were collected from Wulungu Lake (WL), Jili Lake (JL), the Wulungu River (WR), and the Kalaeerqisi River (KR). These samples were then analyzed using microsatellite markers. A total of 98,425 microsatellite markers were developed based on the genomic data, and 29 polymorphic microsatellite markers were selected to analyze genetic diversity in this study. The number of alleles (Na) and observed heterozygosity (Ho) per population ranged from 4.621 (KR) to 11.172 (WL) and from 0.510 (KR) to 0.716 (JL), respectively. The results of the polymorphic information content (PIC) showed that the WL, JL, and WR populations were highly polymorphic (PIC≥ 0.5) and that the KR population was moderately polymorphic (0.25 ≤ PIC < 0.5). The genetic differentiation coefficient (Fst) among the four P. fluviatilis populations was 0.074, indicating moderate genetic differentiation among the populations in Xinjiang. The reason for the significant difference between the rivers and lakes could be the presence of a dam blocking the flow of P. fluviatilis. The development of microsatellite markers provides support for population genetics in the future. The evaluation of the genetic structure of P. fluviatilis in Xinjiang provides a reference for the reproduction and conservation of P. fluviatilis.

Transcriptomic Analysis Provides Insights to Reveal the bmp6 Function Related to the Development of Intermuscular Bones in Zebrafish.

Intermuscular bones (IBs) are small, hard-boned spicules located in the muscle tissue that mainly exist in the myosepta of lower teleosts, which hurt the edibleness and economic value of fish. The study of the development of IBs is very important for freshwater aquaculture fish, but the molecular mechanism of its formation and the key regulatory genes remain unclear. In this study, we first constructed two types of zebrafish mutants (the mutants losing IBs and the mutants with partial deletion of IBs) by knocking out bmp6. We then carried out a transcriptomic analysis to reveal the role of bmp6 in the developmental mechanism of IBs; we used the caudal musculoskeletal tissues of these mutants and wild-type zebrafish at three development stages (20, 45, and 60 dph) to perform transcriptomic analysis. The results showed that the deficiency of bmp6 upregulated sik1 and activated the TNF-A signaling via the NF-KB pathway, which inhibited the development of osteoblasts and promoted osteoclast formation, thereby inhibiting the formation of IBs. These results provided insights to understand the role of bmp6 in the development of IBs in zebrafish and are useful for selective breeding of IBs in cyprinids.

A genetic linkage map and comparative genome analysis of common carp (Cyprinus carpio L.) using microsatellites and SNPs.

A genetic linkage map is a powerful research tool for mapping traits of interest and is essential to understanding genome evolution. The aim of this study is to provide an expanded genetic linkage map of common carp to effectively carry out quantitative trait loci analysis and conduct comparative mapping analysis between lineages. Here, we constructed a genetic linkage map of common carp (Cyprinus carpio L.) using microsatellite and single-nucleotide polymorphism (SNP) markers in a 159 sibling family. A total of 246 microsatellites and 306 SNP polymorphic markers were genotyped in this family. Linkage analysis using JoinMap 4.0 organized 427 markers (186 microsatellites and 241 SNPs) to 50 linkage groups, ranging in size from 1.4 to 130.1 cM. Each group contained 2-30 markers. The linkage map covered a genetic distance of 2,039.2 cM and the average interval for markers within the linkage groups was approximately 6.4 cM. In addition, comparative genome analysis within five model teleost fish revealed a high percentage (74.7%) of conserved loci corresponding to zebrafish chromosomes. In most cases, each zebrafish chromosome comprised two common carp linkage groups. The comparative analysis also revealed independent chromosome rearrangements in common carp and zebrafish. The linkage map will be of great assistance in mapping genes of interest and serve as a reference to approach comparative mapping and enable further insights into the comprehensive investigations of genome evolution of common carp.

[Quantitative trait locus analysis of standard length, body depth and body thickness in mirror carp (Cyprinus carpio L.)].

Based on a full-sib family, the genetic linkage map was constructed with 246 microsatellite and 306 SNP markers, which was used to detect the QTLs for standard length (SL), body depth (H), body thickness (BT), and the ratio of standard length and body depth (SLH) in mirror carp by GridQTL software. The results indicated that a total of 14 related QTLs distributed on the 7 linkage groups were obtained. Seven QTLs were related to standard length, of which the linkage groups of LG6, LG17, LG21, LG23, and LG35 were at 5% significant level, and linkage group LG1 and LG28 were at 1% significant level, which explained 6.6%-12.6% of the phenotypic variance. Three QTLs were identified for body depth on the linkage groups of LG17, LG23 and LG28 (P amp; 0.01), accounting for 11.6%, 12.7%, and 15.6% of the phenotypic variance, respectively. Two QTLs were associated with body thickness on the linkage of LG23 and LG28 (P amp; 0.05), which explained 8.6% and 7.2% of the phenotypic variation, respectively. Two QTLs were responsible for the ratio of standard length and body depth on the linkage of LG21 and LG35 (P amp; 0.05), both of which explained 8.2% of the phenotypic variance. The results provide a useful reference for further candidate gene research and molecular marker assisted selection in mirror carp.

Resequencing and SNP discovery of Amur ide (Leuciscus waleckii) provides insights into local adaptations to extreme environments.

Amur ide (Leuciscus waleckii), a Cyprinid species, is broadly distributed in Northeast Asia. Different from its freshwater counterparts, the population in Lake Dali Nor has a strong alkalinity tolerance and can adapt to extremely alkali-saline water with bicarbonate over 50 mmol/L. To uncover the genetic basis of its alkaline adaptation, three populations, including one alkali form from Lake Dali Nor (DL), one freshwater form from its adjacent sister Lake Ganggeng Nor (GG), and one freshwater form from its historical origin, namely, the Songhua River (SH), were analyzed using genome resequencing technology. A total of 679.82 Gb clean data and 38,091,163 high-quality single-nucleotide polymorphism (SNP) loci were detected in the three populations. Nucleotide diversity and population structure analysis revealed that the DL and GG populations have lower nucleotide diversities and different genetic structures than those of the SH population. Selective sweeping showed 21 genes involved in osmoregulatory regulation (DLG1, VIPR1, AKT1, and GNAI1), inflammation and immune responses (DLG1, BRINP1, CTSL, TRAF6, AKT1, STAT3, GNAI1, SEC22b, and PSME4b), and cardiorespiratory development (TRAF6, PSME4b, STAT3, AKT1, and COL9A1) to be associated with alkaline adaption of the DL population. Interestingly, selective pressure (CodeML, MEME, and FEL) methods identified two functional codon sites of VIPR1 to be under positive selection in the DL population. The subsequent 3D protein modeling confirmed that these selected sites will incur changes in protein structure and function in the DL population. In brief, this study provides molecular evidence of population divergence and alkaline adaptation, which will be very useful for revealing the genetic basis of alkaline adaptation in Amur ide.

De novo assembly and characterization of the Hucho taimen transcriptome.

Taimen (Hucho taimen) is an important ecological and economic species that is classified as vulnerable by the IUCN Red List of Threatened Species; however, limited genomic information is available on this species. RNA-Seq is a useful tool for obtaining genetic information and developing genetic markers for nonmodel species in addition to its application in gene expression profiling. In this study, we performed a comprehensive RNA-Seq analysis of taimen. We obtained 157 M clean reads (14.7 Gb) and used them to de novo assemble a high-quality transcriptome with a N50 size of 1,060 bp. In the assembly, 82% of the transcripts were annotated using several databases, and 14,666 of the transcripts contained a full open reading frame. The assembly covered 75% of the transcripts of Atlantic salmon and 57.3% of the protein-coding genes of rainbow trout. To learn about the genome evolution, we performed a systematic comparative analysis across 11 teleosts including eight salmonids and found 313 unique gene families in taimen. Using Atlantic salmon and rainbow trout transcriptomes as the background, we identified 250 positive selection transcripts. The pathway enrichment analysis revealed a unique characteristic of taimen: It possesses more immune-related genes than Atlantic salmon and rainbow trout; moreover, some genes have undergone strong positive selection. We also developed a pipeline for identifying microsatellite marker genotypes in samples and successfully identified 24 polymorphic microsatellite markers for taimen. These data and tools are useful for studying conservation genetics, phylogenetics, evolution among salmonids, and selective breeding for threatened taimen.

JDP2: An oncogenic bZIP transcription factor in T cell acute lymphoblastic leukemia.

A substantial subset of patients with T cell acute lymphoblastic leukemia (T-ALL) develops resistance to steroids and succumbs to their disease. JDP2 encodes a bZIP protein that has been implicated as a T-ALL oncogene from insertional mutagenesis studies in mice, but its role in human T-ALL pathogenesis has remained obscure. Here we show that JDP2 is aberrantly expressed in a subset of T-ALL patients and is associated with poor survival. JDP2 is required for T-ALL cell survival, as its depletion by short hairpin RNA knockdown leads to apoptosis. Mechanistically, JDP2 regulates prosurvival signaling through direct transcriptional regulation of MCL1. Furthermore, JDP2 is one of few oncogenes capable of initiating T-ALL in transgenic zebrafish. Notably, thymocytes from rag2:jdp2 transgenic zebrafish express high levels of mcl1 and demonstrate resistance to steroids in vivo. These studies establish JDP2 as a novel oncogene in high-risk T-ALL and implicate overexpression of MCL1 as a mechanism of steroid resistance in JDP2-overexpressing cells.

Quantitative trait loci for morphometric traits in multiple families of common carp (Cyprinus carpio).

Quantitative trait locus (QTL) mapping is frequently used to understand the genetic architecture of quantitative traits. Herein, we performed a genome scan for QTL affecting the morphometric characters in eight full-sib families containing 522 individuals using different statistical methods (Sib-pair and half-sib model). A total of 194 QTLs were detected in 25 different regions on 10 linkage groups (LGs). Among them, 37 QTLs on five LGs (eight, 13, 24, 40 and 45) were significant (5% genome-wide level), while the remaining 40 (1% chromosome-wide level) and 117 (5% chromosome-wide level) indicated suggestive effect on those traits. Heritabilities for most morphometric traits were moderate to high, ranging from 0.21 to 0.66, with generally strong phenotypic and genetic correlations between the traits. A large number of QTLs for morphometric traits were co-located, consistent with their high correlations, and may reflect pleiotropic effect on the same genes. Biological pathways were mapped for possible candidate genes on QTL regions. One significantly enriched pathway was identified on LG45, which had a P-value of 0.04 and corresponded to the "regulation of actin cytoskeleton pathway". The results are expected to be useful in marker-assisted selection (MAS) and provide valuable information for the study of gene pathway for morphometric and growth traits of the common carp.

Genome-Wide Association Study for Muscle Fat Content and Abdominal Fat Traits in Common Carp (Cyprinus carpio).

Muscle fat content is an important phenotypic trait in fish, as it affects the nutritional, technical and sensory qualities of flesh. To identify loci and candidate genes associated with muscle fat content and abdominal fat traits, we performed a genome-wide association study (GWAS) using the common carp 250 K SNP assay in a common carp F2 resource population. A total of 18 loci surpassing the genome-wide suggestive significance level were detected for 4 traits: fat content in dorsal muscle (MFdo), fat content in abdominal muscle (MFab), abdominal fat weight (AbFW), and AbFW as a percentage of eviscerated weight (AbFP). Among them, one SNP (carp089419) affecting both AbFW and AbFP reached the genome-wide significance level. Ten of those loci were harbored in or near known genes. Furthermore, relative expressions of 5 genes related to MFdo were compared using dorsal muscle samples with high and low phenotypic values. The results showed that 4 genes were differentially expressed between the high and low phenotypic groups. These genes are, therefore, prospective candidate genes for muscle fat content: ankyrin repeat domain 10a (ankrd10a), tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (tanc2), and four jointed box 1 (fjx1) and choline kinase alpha (chka). These results offer valuable insights into the complex genetic basis of fat metabolism and deposition.

The genetic map of goldfish (Carassius auratus) provided insights to the divergent genome evolutions in the Cyprinidae family.

A high-density linkage map of goldfish (Carassius auratus) was constructed using RNA-sequencing. This map consists of 50 linkage groups with 8,521 SNP markers and an average resolution of 0.62 cM. Approximately 84% of markers are in protein-coding genes orthologous to zebrafish proteins. We performed comparative genome analysis between zebrafish and medaka, common carp, grass carp, and goldfish to study the genome evolution events in the Cyprinidae family. The comparison revealed large synteny blocks among Cyprinidae fish and we hypothesized that the Cyprinidae ancestor undergone many inter-chromosome rearrangements after speciation from teleost ancestor. The study also showed that goldfish genome had one more round of whole genome duplication (WGD) than zebrafish. Our results illustrated that most goldfish markers were orthologous to genes in common carp, which had four rounds of WGD. Growth-related regions and genes were identified by QTL analysis and association study. Function annotations of the associated genes suggested that they might regulate development and growth in goldfish. This first genetic map enables us to study the goldfish genome evolution and provides an important resource for selective breeding of goldfish.

Complete mitochondrial genome sequence of Romanogobio tenuicorpus (Amur whitefin gudgeon).

Amur whitefin gudgeon (Romanogobio tenuicorpus) belongs to the family Cyprinidae, it is freshwater aquaculture species in China. In the report, we determined the complete mitochondrial genome sequence of Romanogobio tenuicorpus, which is 16,600 bp long circular molecule with 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a control region, the conserved sequence blocks, CSB1, CSB2 and CSB3 were also detected.

The complete mitochondrial genome sequence of Hemiculter leucisculus.

The complete mitochondrial genome of Hemiculter leucisculus was determined to be 16,617 bp. It contains the structure of 22 transfer RNA genes, 13 protein-coding genes, 2 ribosomal RNA genes, and non-coding control region (D-loop). The critical central conserved sequences (CSB-D, CSB-E, and CSB-F) were also detected. The determination of H. leucisculus mitogenome would play an important role in genetic diversity and population vitality in Cyprinidae.