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.

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.

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.

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.

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.

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.

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.

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.

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.

A zebrafish model of myelodysplastic syndrome produced through tet2 genomic editing.

The ten-eleven translocation 2 gene (TET2) encodes a member of the TET family of DNA methylcytosine oxidases that converts 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) to initiate the demethylation of DNA within genomic CpG islands. Somatic loss-of-function mutations of TET2 are frequently observed in human myelodysplastic syndrome (MDS), which is a clonal malignancy characterized by dysplastic changes of developing blood cell progenitors, leading to ineffective hematopoiesis. We used genome-editing technology to disrupt the zebrafish Tet2 catalytic domain. tet2(m/m) (homozygous for the mutation) zebrafish exhibited normal embryonic and larval hematopoiesis but developed progressive clonal myelodysplasia as they aged, culminating in myelodysplastic syndromes (MDS) at 24 months of age, with dysplasia of myeloid progenitor cells and anemia with abnormal circulating erythrocytes. The resultant tet2(m/m) mutant zebrafish lines show decreased levels of 5hmC in hematopoietic cells of the kidney marrow but not in other cell types, most likely reflecting the ability of other Tet family members to provide this enzymatic activity in nonhematopoietic tissues but not in hematopoietic cells. tet2(m/m) zebrafish are viable and fertile, providing an ideal model to dissect altered pathways in hematopoietic cells and, for small-molecule screens in embryos, to identify compounds with specific activity against tet2 mutant cells.

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.

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 pike perch (Sander canadensis).

Pike perch (Sander canadensis) is a member of the largest order of Osteichthyes, Perciformes, and is an important ecological and economic freshwater species, which distributes in Ili River and Ergis River of Xinjiang Province, China. In this study, we sequenced the whole mitochondrial genome of pike perch, and analyzed the similarity with its related species. The mitochondrial genome of S. canadensis is 16,542 bp in length with 55.05% AT content, contained 13 protein coding genes, 22 tRNA genes, 2 ribosomal genes and an 892 bp non-coding region. In control region, 6 CSBs (CSB-1, CSB-2, CSB-3, CSB-D, CSB-E and CSB-F), one potential TAS and one poly-T region were identified. Comparing all protein-coding genes and whole genome sequence with 4 species of Perciformes (three species of Percidae, Perca flavescens. Percina macrolepida. Etheostoma radiosum and one outgroup Oreochromis sp. red tilapia), ND3 gene has the highest mutation rate, and S. canadensis has higher similarity with Perca flavescens than others. The mitochondrial genomic sequence will help us to study the conservation genetic and evolution of Percidae.

Genome sequence and genetic diversity of the common carp, Cyprinus carpio.

The common carp, Cyprinus carpio, is one of the most important cyprinid species and globally accounts for 10% of freshwater aquaculture production. Here we present a draft genome of domesticated C. carpio (strain Songpu), whose current assembly contains 52,610 protein-coding genes and approximately 92.3% coverage of its paleotetraploidized genome (2n = 100). The latest round of whole-genome duplication has been estimated to have occurred approximately 8.2 million years ago. Genome resequencing of 33 representative individuals from worldwide populations demonstrates a single origin for C. carpio in 2 subspecies (C. carpio Haematopterus and C. carpio carpio). Integrative genomic and transcriptomic analyses were used to identify loci potentially associated with traits including scaling patterns and skin color. In combination with the high-resolution genetic map, the draft genome paves the way for better molecular studies and improved genome-assisted breeding of C. carpio and other closely related species.

Duplication and differentiation of common carp (Cyprinus carpio) myoglobin genes revealed by BAC analysis.

Two distinct myoglobin (mb) transcripts have been reported in common carp, Cyprinus carpio, which is a hypoxia-tolerant fish living in habitats with greatly fluctuant dissolved oxygen levels. Recombinant protein analysis has shown functional specialization of the two mb transcripts. In this work, analysis for mb-containing bacterial artificial chromosome (BAC) clones indicated different genome loci for common carp myoglobin-1 (mb-1) and myoglobin-2 (mb-2) genes. Fluorescence in situ hybridization (FISH) revealed that mb-1 and mb-2 are located on separate chromosomes. In both of the mb-1 and mb-2 containing BAC clones, gene synteny was well conserved with the homologous region on zebrafish chromosome 1, supporting that the common carp specific mb-2 gene originated from the recent whole genome duplication event in cyprinid lineage. Transcription factor binding sites search indicated that both common carp mb genes lacked specificity Protein 1 (Sp1) and myocyte enhancer factor-2 (MEF2) binding sites, which mediated muscle-specific and calcium-dependent expression in the well-studied mb promoters. Potential hypoxia response elements (HREs) were predicted in the regulatory region of common carp mb genes. These characteristics of common carp mb gene regulatory region well interpreted the hypoxia-inducible, non-muscle expression pattern of mb-1. In the case of mb-2, a 10 bp insertion to the binding site of upstream stimulatory factor (USF), which was a co-factor of hypoxia inducible factor (HIF), might cause the non-response to hypoxia treatment of mb-2. The case of common carp mb gene duplication and subsequent differentiation in expression pattern and protein function provided an example for adaptive evolution toward aquatic hypoxia tolerance.

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.