Establishment and identification of a skin cell line from Chinese tongue sole (Cynoglossus semilaevis) and analysis of the changes in its transcriptome upon LPS stimulation.

The Chinese tongue sole (Cynoglossus semilaevis) holds significant economic importance within the fishing industry along the eastern coasts of China. In recent years, the frequent outbreaks of bacterial diseases have become a common concern as the aquaculture scale expands. The majority of the diseased fish exhibit symptoms such as skin congestion, damage and skin ulceration. As the skin serves as the first line of defense against bacterial infections, establishing a skin cell line for immunological research on Chinese tongue sole's response to bacterial infection is of utmost importance. In this study, a cell line named CSS (derived from the skin of the Chinese tongue sole) was successfully established. The cells have demonstrated stability during passages and exhibit a multipolar fibroblast-like morphology. They were cultured in L-15 medium with 20% serum and have been successfully passed through 60 passages over a period of 20 months. The identification of the mitochondrial CO1 gene confirmed that the cell originated from Chinese tongue sole. The karyotype detection revealed that the cell had a chromosome number of 2n = 42. After being stored in liquid nitrogen for 15 months, the cells can maintain more than 75% viability upon recovery. After transfecting with cy3-labeled scramble siRNA and pEGFP-N3 plasmid, clear fluorescence was observed in the transfected cells. We observed that lipopolysaccharide (LPS) from Escherichia coli significantly upregulate the gene expression of various immune-related pathways at 2 h in the CSS cell line. Additionally, the differentially expressed genes showed a higher enrichment in immune-related pathways at 2 and 6 h after stimulation compared to the 24 h point. Moreover, we identified 347 genes that exhibited a gradual increase in expression during the 0-24 h stimulation period. These genes were primarily enriched in pathways related to Autophagy, GABAergic synapse, Apelin signaling and Ferroptosis. In general, the CSS cell line established in this study exhibits stable growth and can serve as a valuable tool for in vitro studies of immunology and other basic biologies of Chinese tongue sole.

Genome-wide DNA methylation mediates the resistance to vibriosis in Cynoglossus semilaevis.

Chinese tongue sole (Cynoglossus semilaevis) is an economically important marine fish in China. However, vibriosis has caused huge mortality and economic losses in its culturing industry. To reveal the effect of DNA methylation on the resistance to vibriosis in tongue sole, we conducted RNA sequencing and whole genome bisulfite sequencing (WGBS), and compared the gene expressions and DNA methylation patterns between the resistant and susceptible families. We identified a total of 741 significantly differentially expressed genes (DEGs) in kidney and 17460 differentially methylated genes (DMGs), which were both enriched in immune-related pathways, such as "cAMP signaling pathway" and "NOD-like receptor signaling pathway". Through the correlation analysis of DEGs and DMGs, we identified two important immune pathways, including "complement and coagulation cascades", and "cytokine-cytokine receptor interaction", which played important roles in regulating the inflammation level and immune homeostasis. For example, the expression of proinflammatory cytokine il17c was down-regulated under the regulation of DNA methylation; in addition, the expression of protease-activated receptor 3 (par3) was up-regulated, which could induce the up-expressionof il8. These results demonstrated that the regulation of DNA methylation on the genes involved in immune responses might contribute to the resistance to vibriosis in tongue sole, and provided a basis for the control of diseases in fish aquaculture.

A novel miRNA, Cse-miR-33, functions as an immune regulator by targeting CsTRAF6 in Chinese tongue sole (Cynoglossus semilaevis).

The tumor necrosis factor receptor-associated factor 6 (TRAF6) can act as a fundamental adaptor protein in a chain reaction of signal transduction and cascade events to finish off immune defenses. However, immunomodulatory research on TRAF6 gene is still limited in fish. In this study, a novel miRNA, Cse-miR-33 was identified from the whole genome of Chinese tongue sole (Cynoglossus semilaevis). After separate infections with three different Vibrio strains (V. harveyi, V. anguillarum, V. parahemolyticus) and one virus (nervous necrosis virus, NNV), the expressions of CsTRAF6 and Cse-miR-33 displayed significant time-dependent changes in immune related tissues and the trends were opposite in general. Through target gene prediction and dual luciferase reporter assay, Cse-miR-33 was proven to regulate CsTRAF6 by combining with 3'-UTR sequence of the gene. The results of qRT-PCR and western blotting (WB) analyses showed that Cse-miR-33 blocked the translation of CsTRAF6 protein at post-transcriptional level, rather than degrading the target mRNA. Further experiment indicated that Cse-miR-33 inhibitor largely reduced the death rate of Chinese tongue sole caused by V. harveyi and NNV. The expressions of CsTRAF6-associated immune genes (such as CsIL-1R, CsMYD88, CsIRAK1, CsTNFα, CsIL6 and CsIL8) were also significantly changed in response to Cse-miR-33 agomir and inhibitor. The study suggested that Cse-miR-33 affected the immune response via targeting CsTRAF6 in C. semilaevis, which would provide us deep insights into miRNA-mediated regulatory network and help improve the immunity in fish.

Effects of hypoxia and reoxygenation on oxidative stress, histological structure, and apoptosis in a new hypoxia-tolerant variety of blunt snout bream (Megalobrama amblycephala).

A new hypoxia-tolerant variety of blunt snout bream was obtained by successive breeding of the wild population, which markedly improved hypoxia tolerance. In this study, the hypoxia-tolerant variety was exposed to hypoxia (2.0 mg O2·L-1) for 4, 7 days. The contents of blood biochemical indicators including the number of red blood cells (RBC), total cholesterol (T-CHO), total protein (TP), triglyceride (TG), glucose (GLU), and lactic acid (LD) increased significantly (P < 0.05) under hypoxia. The glycogen content in the liver and muscle decreased significantly (P < 0.05) and the LD content in the brain, muscle and liver increased significantly (P < 0.05) under hypoxia. The levels of oxidative stress-related indicators i.e., superoxide dismutase (SOD), malondialdehyde (MDA), glutathione (GSH), catalase (CAT), and total antioxidant capacity (T-AOC) also changed significantly (P < 0.05) in the heart, liver, and intestine of the new variety under hypoxia. Additionally, hypoxia has caused injuries to the heart, liver, and intestine, but it shows amazing repair ability during reoxygenation. The apoptotic cells and apoptosis rate in the heart, liver, and intestine increased under hypoxia. Under hypoxia, the expression of the B-cell lymphomas 2 (Bcl-2) gene in the heart, liver, and intestine was significantly (P < 0.05) down-regulated, while the expression of the BCL2-associated agonist of cell death (Bad) gene was significantly (P < 0.05) up-regulated. These results are of great significance for enriching the basic data of blunt snout bream new variety in response to hypoxia and promoting the healthy development of its culture industry.

Identification of crucial factors involved in Cynoglossus semilaevis sexual size dimorphism by GWAS and demonstration of zbed1 regulatory network by DAP-seq.

Sexual size dimorphism (SSD), whereby females and males exhibit different body sizes, are widely documented in animals. To explore crucial regulators implicated in female-biased SSD of Chinese tongue sole (Cynoglossus semilaevis), GWAS was conducted on 350 females and 59 males. Twenty SNPs and 25 genes including zbed1, nsd3, cdc45, klhl29, and smad4 with -log(p) > 7 were screened, mainly mapping to sex chromosome. The chromosome W-linked gene zbed1 attracted particular attention because it is a master key for cell proliferation. Thus, the regulatory network of zbed1 in C. semilaevis was explored by DAP-seq and 1352 peaks were discovered in the female brain. Moreover, zbed1 potentially regulated hippo signaling pathway, cell cycle, translation, and PI3k-Akt signaling pathway in C. semilaevis. These findings identify crucial SNPs and genes associated with female-biased SSD in C. semilaevis, also provide the first genome-wide survey for the zbed1 regulatory network in fish species.

Genome-Wide Identification and Involvement in Response to Biotic and Abiotic Stresses of lncRNAs in Turbot (Scophthalmus maximus).

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

Genome-Wide Identification, Molecular Characterization, and Involvement in Response to Abiotic and Biotic Stresses of the HSP70 Gene Family in Turbot (Scophthalmus maximus).

Heat shock proteins 70 (HSP70s) are known to play essential roles in organisms' response mechanisms to various environmental stresses. However, no systematic identification and functional analysis has been conducted for HSP70s in the turbot (Scophthalmus maximus), a commercially important worldwide flatfish. Herein, 16 HSP70 genes unevenly distributed on nine chromosomes were identified in the turbot at the genome-wide level. Analyses of gene structure, motif composition, and phylogenetic relationships provided valuable data on the HSP70s regarding their evolution, classification, and functional diversity. Expression profiles of the HSP70 genes under five different stresses were investigated by examining multiple RNA-seq datasets. Results showed that 10, 6, 8, 10, and 9 HSP70 genes showed significantly up- or downregulated expression after heat-induced, salinity-induced, and Enteromyxum scophthalmi, Vibrio anguillarum, and Megalocytivirus infection-induced stress, respectively. Among them, hsp70 (hspa1a), hspa1b, and hspa5 showed significant responses to each kind of induced stress, and qPCR analyses further validated their involvement in comprehensive anti-stress, indicating their involvement in organisms' anti-stress mechanisms. These findings not only provide new insights into the biological function of HSP70s in turbot adapting to various environmental stresses, but also contribute to the development of molecular-based selective breeding programs for the production of stress-resistant turbot strains in the aquaculture industry.

Identification and Functional Analysis of foxo Genes in Chinese Tongue Sole (Cynoglossus semilaevis).

The Chinese tongue sole (Cynoglossus semilaevis) is a traditional, precious fish in China. Due to the large growth difference between males and females, the investigation of their sex determination and differentiation mechanisms receives a great deal of attention. Forkhead Box O (FoxO) plays versatile roles in the regulation of sex differentiation and reproduction. Our recent transcriptomic analysis has shown that foxo genes may participate in the male differentiation and spermatogenesis of Chinese tongue sole. In this study, six Csfoxo members (Csfoxo1a, Csfoxo3a, Csfoxo3b, Csfoxo4, Csfoxo6-like, and Csfoxo1a-like) were identified. Phylogenetic analysis indicated that these six members were clustered into four groups corresponding to their denomination. The expression patterns of the gonads at different developmental stages were further analyzed. All members showed high levels of expression in the early stages (before 6 months post-hatching), and this expression was male-biased. In addition, promoter analysis found that the addition of C/EBPα and c-Jun transcription factors enhanced the transcriptional activities of Csfoxo1a, Csfoxo3a, Csfoxo3b, and Csfoxo4. The siRNA-mediated knockdown of the Csfoxo1a, Csfoxo3a, and Csfoxo3b genes in the testicular cell line of Chinese tongue sole affected the expression of genes related to sex differentiation and spermatogenesis. These results have broadened the understanding of foxo's function and provide valuable data for studying the male differentiation of tongue sole.

Reconstruction of the Origin of a Neo-Y Sex Chromosome and Its Evolution in the Spotted Knifejaw, Oplegnathus punctatus.

Sex chromosomes are a peculiar constituent of the genome because the evolutionary forces that fix the primary sex-determining gene cause genic degeneration and accumulation of junk DNA in the heterogametic partner. One of the most spectacular phenomena in sex chromosome evolution is the occurrence of neo-Y chromosomes, which lead to X1X2Y sex-determining systems. Such neo-sex chromosomes are critical for understanding the processes of sex chromosome evolution because they rejuvenate their total gene content. We assembled the male and female genomes at the chromosome level of the spotted knifejaw (Oplegnathus punctatus), which has a cytogenetically recognized neo-Y chromosome. The full assembly and annotation of all three sex chromosomes allowed us to reconstruct their evolutionary history. Contrary to other neo-Y chromosomes, the fusion to X2 is quite ancient, estimated at 48 Ma. Despite its old age and being even older in the X1 homologous region which carries a huge inversion that occurred as early as 55-48 Ma, genetic degeneration of the neo-Y appears to be only moderate. Transcriptomic analysis showed that sex chromosomes harbor 87 genes, which may serve important functions in the testis. The accumulation of such male-beneficial genes, a large inversion on the X1 homologous region and fusion to X2 appear to be the main drivers of neo-Y evolution in the spotted knifejaw. The availability of high-quality assemblies of the neo-Y and both X chromosomes make this fish an ideal model for a better understanding of the variability of sex determination mechanisms and of sex chromosome evolution.

Dynamic changes of rhizosphere soil bacterial community and nutrients in cadmium polluted soils with soybean-corn intercropping.

BACKGROUND: Soybean-corn intercropping is widely practised by farmers in Southwest China. Although rhizosphere microorganisms are important in nutrient cycling processes, the differences in rhizosphere microbial communities between intercropped soybean and corn and their monoculture are poorly known. Additionally, the effects of cadmium (Cd) pollution on these differences have not been examined. Therefore, a field experiment was conducted in Cd-polluted soil to determine the effects of monocultures and soybean-corn intercropping systems on Cd concentrations in plants, on rhizosphere bacterial communities, soil nutrients and Cd availability. Plants and soils were examined five times in the growing season, and Illumina sequencing of 16S rRNA genes was used to analyze the rhizosphere bacterial communities. RESULTS: Intercropping did not alter Cd concentrations in corn and soybean, but changed soil available Cd (ACd) concentrations and caused different effects in the rhizosphere soils of the two crop species. However, there was little difference in bacterial community diversity for the same crop species under the two planting modes. Proteobacteria, Chloroflexi, Acidobacteria, Actinobacteria and Firmicutes were the dominant phyla in the soybean and corn rhizospheres. In ecological networks of bacterial communities, intercropping soybean (IS) had more module hubs and connectors, whereas intercropped corn (IC) had fewer module hubs and connectors than those of corresponding monoculture crops. Soil organic matter (SOM) was the key factor affecting soybean rhizosphere bacterial communities, whereas available nutrients (N, P, K) were the key factors affecting those in corn rhizosphere. During the cropping season, the concentration of soil available phosphorus (AP) in the intercropped soybean-corn was significantly higher than that in corresponding monocultures. In addition, the soil available potassium (AK) concentration was higher in intercropped soybean than that in monocropped soybean. CONCLUSIONS: Intercropped soybean-corn lead to an increase in the AP concentration during the growing season, and although crop absorption of Cd was not affected in the Cd-contaminated soil, soil ACd concentration was affected. Intercropped soybean-corn also affected the soil physicochemical properties and rhizosphere bacterial community structure. Thus, intercropped soybean-corn was a key factor in determining changes in microbial community composition and networks. These results provide a basic ecological framework for soil microbial function in Cd-contaminated soil.

Clinical characteristics and postoperative outcomes of systemic artery-to-pulmonary vessel fistula in hemoptysis patients.

OBJECTIVES: To investigate the clinical characteristics and outcomes on the success of bronchial arterial embolization (BAE) in patients with and without systemic artery-to-pulmonary vessel fistula (SA-PF) and to evaluate the feasibility of CTA in the assessment of SA-PF. METHODS: We retrospectively enrolled 420 consecutive patients that underwent BAE for hemoptysis control in our hospital from September 2011 to May 2019. The clinical characteristics, preprocedural CTA findings, BAE procedural findings, and follow-up outcomes were collected. Patients were divided into two groups according to DSA findings: patients with SA-PF and those without. RESULTS: A total of 184 (43.7%) patients presented with SA-PF. Pneumonia was less likely to be the concomitant condition in patients with SA-PF (p < 0.001). The mean number of culprit arteries per patient was significantly higher in patients with SA-PF compared to that in patients without SA-PF (p = 0.017). The SA-PF patients saw a greater probability of recurrence (HR: 2.782, 95% CI: 1.617-4.784, p < 0.001). SA-pulmonary venous fistula (SA-PVF) favored lower hemoptysis recurrence rate (HR: 0.199, 95%CI: 0.052-0.765, p = 0.019). SA-pulmonary artery fistula (SA-PAF) can be detected by optimized CTA protocol with a detection rate of 65.3% (49/75). CONCLUSIONS: The presence of SA-PF is an independent risk factor predicting early recurrence of hemoptysis after BAE. SA-PVF seems to be a protective factor for longer hemoptysis control compared to SA-PAF. Optimized preprocedural CTA is a reliable examination to identify SA-PAF. KEY POINTS: • The appearance of SA-PF is associated with a greater probability of early recurrent hemoptysis after bronchial artery embolization. • The presence of SA-PVF seems to be a protective factor for longer hemoptysis control after BAE compared to SA-PAF. • Optimized CTA protocol seems to be a promising auxiliary examination to detect SA-PAF.

Genome-wide identification, immune response profile and functional characterization of IL-10 from spotted knifejaw (Oplegnathus punctatus) during host defense against bacterial and viral infection.

Interleukin 10 (IL-10), a pleiotropic cytokine, plays an essential role in multiple immunity responses. In the current study, the sequences of IL-10 family were identified from spotted knifejaw (Oplegnathus punctatus) whole genome, and O. punctatus IL-10 (OpIL-10) was cloned and characterized. OpIL-10 encodes 187 amino acids with a typical IL-10 family signature motif and predicted α-helices. It shared high identities with Notolabrus celidotus IL-10 and Epinephelus Lanceolatus IL-10. OpIL-10 was widely detected in healthy tissues, with the abundant expression in liver and skin. It was significantly up-regulated in the six immune-related tissues (liver, spleen, kidney, intestine, gill and skin) after infection against Vibrio harveyi and spotted knifejaw iridovirus (SKIV). Dual-luciferase analysis showed that OpIL-10 overexpression could suppress the activity of NF-κB. Meanwhile, OpIL-10 knockdown caused the down-regulation of five immune-related genes in JAK2/STAT3 signaling pathway and NF-κB signaling pathway, including IL-10R2, TYK2, STAT3, NOD2, and IκB. In addition, LPS and poly I:C stimulated expression of pro-inflammatory cytokines, including IL-6, IL-1ß, IL-8, and IL-12, were lower with recombinant OpIL-10 (rOp IL-10) than the control group, indicating the anti-inflammatory roles of rOpIL-10. Taken together, these results indicated OpIL-10 as a negative regulator in the inflammatory responses of spotted knifejaw against bacterial and viral infection, which would help us better understand the role of IL-10 in teleost immunity.

Transcriptome analysis and protein-protein interaction in resistant and susceptible families of Japanese flounder (Paralichthys olivaceus) to understand the mechanism against Edwardsiella tarda.

Edwardsiella tarda is one of the most harmful bacterial pathogens for aquaculture flatfish. After artificial infection of 47 Japanese flounder (Paralichthys olivaceus) families, resistant and susceptible families were identified in this study. High-throughput sequencing was performed on the liver transcriptome of uninfected groups (PoRU and PoSU) and infected groups (PoRC and PoSC). Through assembly and annotation, a total of 3012 and 1386 differentially expressed genes (DEGs) were identified in PoRU vs. PoSU and PoRC vs. PoSC. The significant enrichment pathways between PoRU and PoSU were mainly in metabolic and biosynthesis pathways. A total of thirty dominant enrichment pathways between PoRC and PoSC mainly focused on some immune-related pathways, including the hematopoietic cell lineage, cytokine-cytokine receptor interaction, complement and coagulation cascades, antigen processing and presentation, the intestinal immune network for immunoglobulin A (IgA) production and T/B cell receptor signaling pathway. Under the protein-protein interaction (PPI) analysis, hub genes, including CD molecules, complement component factors and chemokines, were identified in the network, and 16 core genes were differentially expressed in resistant and sustainable families in quantitative polymerase chain reaction (qPCR) validation. This study represents the first transcriptome analysis based on resistant and susceptible families and provides resistant genes to understand the potential molecular mechanisms of antibacterial function in marine fish. The results obtained in this study provide crucial information on gene markers for resistant breeding of Japanese flounder.

Molecular cloning, expression analysis of the IgT gene and detection of IgT+ B cells in the half-smooth tongue sole (Cynoglossus semilaevis).

IgT is a specific Ig isotype in teleosts, which plays extremely important roles in the mucosal immunity of fish. In the present study, the membrane-bound and secretory IgT of the half-smooth tongue sole (Cynoglossus semilaevis) were identified for the first time. The V-D-J-C structure of two forms of csIgT are translated by the same Cτ gene, and the secretory tail and transmembrane domain were encoded through alternative splicing at the 3' end of the Cτ4. The CH regions of csIgT had high similarity with that of other flatfish (P. olivaceus and S. maximus). In healthy C. semilaevis, sIgT and mIgT were mainly expressed in mucus related tissues such as skin, intestine and gill. The transcript levels of sIgT and mIgT mRNA showed a significant induction in the immune-related tissues upon Vibrio Harveyi infection. A polyclonal rabbit anti-csIgT was successfully prepared using the csIgT heavy chain recombinant protein. Using this antibody, we detected the native IgT with the molecular mass at 220 kDa in skin total protein under non-reducing SDS-PAGE condition. Immunofluorescence analysis indicated that IgT+ B lymphocytes were intensively located in the skin, gill, intestine, and head kidney of C. semilaevis. These results suggest that IgT may participate in the immune response of C. semilaevis, which will facilitate the investigations of the immunoglobulins of marine fish.

Gill remodeling increases the respiratory surface area of grass carp (Ctenopharyngodon idella) under hypoxic stress.

Seasonal changes, diurnal variations, and eutrophication result in periodic hypoxia in fish habitats, thus affecting the success of commercial aquaculture. In this study, the grass carp (Ctenopharyngodon idella) presented moderate hypoxia tolerance; they showed a medium critical oxygen tension during the loss of equilibrium. In response to 7 d of hypoxic exposure, the erythrocyte count and hemoglobin (Hb) concentration significantly increased (p < 0.01). To cope with the hypoxic environment, the grass carp underwent gill remodeling marked by reduction in the interlamellar cell mass (ILCM) and an increase in respiratory surface area. The gill remodeling under hypoxia was enabled by apoptosis induction. Although apoptotic signals were not found on ILCM cells, transferase dUTP nick end labeling (TUNEL) assay results indicated that after 1 d of hypoxic exposure, the number of TUNEL-positive cells per lamella increased until 4 d and then began to decrease. Consistent with the results of the TUNEL assay, the mRNA expression of apoptosis-related genes, caspase-3, Bax, and Bcl-2, increased at 1, 4, and 7 d of the hypoxia treatment. In addition, gill remodeling significantly (p < 0.01) decreased the concentration of sodium and chloride ions in the fish serum. These findings provide evidence that grass carps increase their respiratory surface area through gill remodeling by apoptosis in the gill filaments to acclimate to a hypoxic environment. This study expands our understanding of the morphological and physiological changes in grass carp in response to a hypoxic environment; therefore, it could be useful for maintaining grass carp production.

Development of a 38 K single nucleotide polymorphism array and application in genomic selection for resistance against Vibrio harveyi in Chinese tongue sole, Cynoglossus semilaevis.

Based on 1572 re-sequenced Chinese tongue sole (Cynoglossus semilaevis), we investigated the accuracy of four genomic methods at predicting genomic estimated breeding values (GEBVs) of Vibrio harveyi resistance in C. semilaevis when SNPs varying from 500 to 500 k. All methods outperformed the pedigree-based best linear unbiased prediction when SNPs reached 50 k or more. Then, we developed an SNP array "Solechip No.1" for C. semilaevis breeding using the Affymetrix Axiom technology. This array contains 38,295 SNPs with an average of 10.5 kb inter-spacing between two adjacent SNPs. We selected 44 candidates as the parents of 23 families and genotyped them by the array. The challenge survival rates of offspring families had a correlation of 0.706 with the mid-parental GEBVs. This SNP array is a convenient and reliable tool in genotyping, which could be used for improving V. harveyi resistance in C. semilaevis coupled with the genomic selection methods.

Effects of long-term sex steroid hormones (estradiol and testosterone)-supplemented feeds on the growth performance of Chinese tongue sole (Cynoglossus semilaevis).

The phenomenon of sexual size dimorphism (SSD), existing in mammals, birds, reptiles, spiders, amphibians, insects, and fishes, is generally related to feeding efficiency, energy allocation, sex steroids, and somatotropic and reproductive endocrine axes. Recently, positive and negative regulations of sex steroids have been reported on SSD in various species. Chinese tongue soles (Cynoglossus semilaevis) at 4 months were fed with 17ß-estradiol (E2) and testosterone (T) supplemented feeds for 8 months to assess the effect of sex steroids on growth traits in different sexes. The potential genetic regulation was examined using several growth-related genes. The results showed that two sex steroid hormones had inhibitory effects on the growth performance of different sexes of C. semilaevis. At the age of 8 months, the expression of insulin-like growth factor 2 gene (igf2), 24-dehydrocholesterol reductase (dhcr24), leptin, and estrogen receptor 2 (esr2) in the liver showed an overall downward trend. The expression of insulin-like growth factor 1 (igf1) was reduced, while thyroid hormone receptor-associated protein 3 (thrap3) expression tended to increase in the gonad after T and E2 treatments. In the brain, somatostatin 1, tandem duplicate 2 (sst1.2) expression increased with the treatment of T and E2 (P < 0.05), while growth hormone-releasing hormone (ghrh) expression decreased. E2 and T had different effects on growth differentiation factor 8 (gdf8) and insulin-like growth factor-binding protein 7 (igfbp7) expression in the muscle. Expression of gdf8 increased in the treated fishes in contrast to the reduction expression of igfbp7. This study provided important clues for understanding the role of sex steroids in flatfish SSD.

Footprints of global change in marine life: Inferring past environment based on DNA methylation and gene expression marks.

Ocean global warming affects the distribution, life history and physiology of marine life. Extreme events, like marine heatwaves, are increasing in frequency and intensity. During sensitive stages of early fish development, the consequences may be long-lasting and mediated by epigenetic mechanisms. Here, we used European sea bass as a model to study the possible long-lasting effects of a marine heatwave during early development. We measured DNA methylation and gene expression in four tissues (brain, muscle, liver and testis) and detected differentially methylated regions (DMRs). Six genes were differentially expressed and contained DMRs three years after exposure to increased temperature, indicating direct phenotypic consequences and representing persistent changes. Interestingly, nine genes contained DMRs around the same genomic regions across tissues, therefore consisting of common footprints of developmental temperature in environmentally responsive loci. These loci are, to our knowledge, the first metastable epialleles (MEs) described in fish. MEs may serve as biomarkers to infer past life history events linked with persistent consequences. These results highlight the importance of subtle phenotypic changes mediated by epigenetics to extreme weather events during sensitive life stages. Also, to our knowledge, it is the first time the molecular effects of a marine heatwave during the lifetime of individuals are assessed. MEs could be used in surveillance programs aimed at determining the footprints of climate change on marine life. Our study paves the way for the identification of conserved MEs that respond equally to environmental perturbations across species. Conserved MEs would constitute a tool of assessment of global change effects in marine life at a large scale.

Characterization of Chinese tongue sole (Cynoglossus semilaevis) 24-dehydrocholesterol reductase: Expression profile, epigenetic modification, and its knock-down effect.

The sexual size dimorphism of the Chinese tongue sole (Cynoglossus semilaevis) has greatly obstructed its sustainable development; however, the underlying mechanism remains unclear. Based on C. semilaevis transcriptomic information, 24-dehydrocholesterol reductase (dhcr24) was identified in steroid biosynthesis, showing female-liver-biased expression. Dhcr24 has been reported to participate in various processes, such as cholesterol synthesis, oxidative stress response, neuroprotection, and cell survival. The present study assessed its role in the sexual size dimorphism in fish. First, detailed expression pattern analysis showed that dhcr24 mRNAs were extensively expressed in tissues and the highest levels were found in the liver and gonads of females. Analysis of the dhcr24 promoter region demonstrated different DNA methylation statuses in female, male, and pseudomale gonads with higher epigenetic modification in males. The confirmation of transcription activity of the dhcr24 promoter and putative transcription factors (e.g., ER, AR, SREBP, and POU1F1a) provides the foundation for studying its regulatory mechanism. Finally, dhcr24-siRNA mediated knock-down assay using C. semilaevis liver cells showed that steroid biosynthesis related genes (e.g., ebp, dhcr7, and sc5d), core component of PI3K/Akt pathway (e.g., pi3k), and igf1r exhibited different expression patterns. Further investigation on the interplay between steroid hormones, dhcr24, PI3K/Akt, and IGF-1 systems will be valuable to better understand the mechanism underlying the sexual size dimorphism in C. semilaevis.

A novel and efficient method to induce allospecific CD8+ memory T lymphocytes.

The aim of the current study was to establish a simple method for effectively inducing memory T lymphocytes by the intraperitoneal injection of spleen lymphocytes into mice. In total, 75 mice were divided into the following groups: an injection group administered three doses of spleen lymphocytes (1 × 106 , 5 × 106 , and 1 × 107 cells), a transplantation group in which a 0.25-cm2 skin section from C57BL/6 mice was transplanted onto the back of the recipient, and a control group in which an equal volume of phosphate-buffered saline was injected. At 1, 2, or 3 months following transplantation, the following parameters were evaluated: quantity of T lymphocytes, percentage of cluster of differentiation 8+ (CD8+ ) memory T cells, and proliferation index of purified CD8+ memory T cells. No significant differences among groups were detected at 1 month (p > .05). However, the injection group administered 1 × 106 cells exhibited the highest proportion of CD8+ memory T cells among all groups at 2 months, and the proportions of CD8+ T cells were higher in the three injection groups than in the skin transplantation and control groups at 3 months. The proportions of memory T cells were higher in the injection groups administered 5 × 106 or 1 × 107 cells than in the skin transplantation and control groups at 3 months. The newly established method effectively induces memory T lymphocytes via the intraperitoneal injection of spleen lymphocytes in vivo and has potential applications in the field of immunotherapy.