Effects of miR-722 on gene expression and alternative splicing in the liver of half-smooth tongue sole after infection with Vibrio anguillarum.

MicroRNAs play crucial roles in various biological processes, including but not limited to differentiation, development, disease, and immunity. However, their immunoregulatory roles in half-smooth tongue sole are lacking. Our previous studies indicated that miR-722 could target C5aR1 to modulate the complement pathway to alleviate inflammatory response and even affect the mortality after the bacterial infection with Vibrio anguillarum. Driven by the purpose of revealing the underlying mechanisms, in this study, we investigated the effects of miR-722 on the gene expression and alternative splicing (AS) in the liver of half-smooth tongue sole after Vibrio anguillarum infection, with the approach of miR-722 overexpression/silencing and subsequent RNA-seq. Among the different comparisons, the I group (miR-722 inhibitor and V. anguillarum) versus blank control (PBS) exhibited the highest number of differentially expressed genes (DEGs), suggesting that the immune response was overactivated after inhibiting the miR-722. In addition, enrichment analyses were performed to reveal the functions of DEGs and differential AS (DAS) genes, reflecting the enrichment of RNA splicing and immune-related pathways including NF-κB and T cell receptor signaling pathway. Comparing the M group (miR-722 mimic and V. anguillarum) with the negative control (random sequence and V. anguillarum), two immune-related genes, cd48 and mapk8, were differentially expressed, of which mapk8 was also differentially spliced, indicating their importance in the immune response. Furthermore, representative gene analysis was performed, suggesting their corresponding functional changes due to AS. To verify the RNA-seq data, quantitative real-time PCR was employed with twenty pairs of primers for DEGs and DAS events. Overall, our results demonstrated that miR-722 could mediate the transcriptome-wide changes of gene expression and AS in half-smooth tongue sole, and provided insights into the regulatory role of miR-722 in immune responses, laying the foundation for further functional analyses and practical applications in aquaculture.

Manipulated C5aR1 over/down-expression associates with IL-6 expression during bacterial inflammation in half-smooth tongue sole (Cynoglossus semilaevis).

The complement component 5a/complement component 5 receptor 1 (C5a/C5aR1) pathway plays a crucial role in the onset and development of inflammation, but relevant studies in fish are lacking. In this study, we successfully characterized the relationship between half-smooth tongue sole (Cynoglossus semilaevis) C5aR1 (CsC5aR1) and bacterial inflammation. First, we showed that the overexpression of CsC5aR1 significantly increased bacterial pathological damage in the liver and intestine, whereas inhibition attenuated the damage. The in vitro experiments suggested that CsC5aR1 was able to positively regulate the phagocytic activity and respiratory burst of tongue sole macrophages. In terms of both transcriptional and translational levels, overexpression/inhibition of CsC5aR1 was followed by a highly consistent up-regulation/decrease of its downstream canonical inflammatory factor interleukin-6 (CsIL-6). Furthermore, we stimulated macrophages by lipopolysaccharide (LPS) and lipoteichoic acid (LTA) and found a broad-spectrum response to bacterial infections by the C5a/C5aR1 complement pathway together with the downstream inflammatory factor CsIL-6. Subsequently, we directly elucidated that CsIL-6 is an indicator of C5a/C5aR1-mediated inflammation at different infection concentrations, different infectious bacteria (Vibrio anguillarum and Mycobacterium marinum), and different detection levels. These results might provide a new inflammation bio-marker for early warning of bacteria-induced hyperinflammation leading to fish mortality and a promising target for the treatment of bacterial inflammation in teleost.

A comparative analysis of alternative splicing patterns in Atlantic salmon (Salmo salar) in response to Moritella viscosa and sea lice (Lepeophtheirus salmonis) infection.

Moritella viscosa (M. viscosa) and sea lice (Lepeophtheirus salmonis) are severe pathogens that primarily infect the skin of Atlantic salmon (Salmo salar), which cause significant economic losses in the farming industry. However, the pathogenesis and molecular mechanisms underlying the host's immune defence at the post-transcriptional level remain unclear. Alternative splicing (AS) is an evolutionarily conserved post-transcriptional mechanism that can greatly increase the richness of the transcriptome and proteome. In this study, transcriptomic data derived from skin tissues of Atlantic salmon after M. viscosa and sea lice infections were used to examine the AS profiles and their differential expression patterns. In total, we identified 33,044 AS events (involving 13,718 genes) in the control (CON) group, 35,147 AS events (involving 14,340 genes) in the M. viscosa infection (MV) group, and 30,364 AS events (involving 13,142 genes) in the sea lice infection (LC) group, respectively. Among the five types of AS identified in our study (i.e., SE, A5SS, A3SS, MXE, and RI), SE was the most prevalent type in all three groups (i.e., CON, MV, and LC groups). Decreased percent-spliced-in (PSI) levels were observed in SE events under both MV- and LC-infected conditions, suggesting that MV or LC infection elevated exon-skipping isoforms and promoted the selection of shorter transcripts in numerous DAS genes. In addition, most of the differential AS genes were found to be associated with pathways related to mRNA regulation, epithelial or muscle development, and immune response. These findings provide novel insights into the role of AS in host-pathogen interactions and represent the first comparative analysis of AS in response to bacterial and parasitic infections in fish.

MetazExp: a database for gene expression and alternative splicing profiles and their analyses based on 53 615 public RNA-seq samples in 72 metazoan species.

RNA-seq has been widely used in experimental studies and produced a massive amount of data deposited in public databases. New biological insights can be obtained by retrospective analyses of previously published data. However, the barrier to efficiently utilize these data remains high, especially for those who lack bioinformatics skills and computational resources. We present MetazExp (https://bioinfo.njau.edu.cn/metazExp), a database for gene expression and alternative splicing profiles based on 53 615 uniformly processed publicly available RNA-seq samples from 72 metazoan species. The gene expression and alternative splicing profiles can be conveniently queried by gene IDs, symbols, functional terms and sequence similarity. Users can flexibly customize experimental groups to perform differential and specific expression and alternative splicing analyses. A suite of data visualization tools and comprehensive links with external databases allow users to efficiently explore the results and gain insights. In conclusion, MetazExp is a valuable resource for the research community to efficiently utilize the vast public RNA-seq datasets.

Inflammatory reaction and immune response of half-smooth tongue sole (Cynoglossus semilaevis) after infection with Vibrio anguillarum.

Frequently occurred bacterial diseases have seriously affected the aquaculture industry of half-smooth tongue sole (Cynoglossus semilaevis). Notably, vibriosis, with Vibrio anguillarum as one of the causative pathogens, is the most severe bacterial disease with severe inflammatory response of the host, leading to high mortality rates. In the present study, we explored the relationship between bacterial concentrations and host mortality, inflammatory reaction, and immune response in half-smooth tongue sole after infection with V. anguillarum at different concentrations (Treatment 1, 6.4 × 105 CFU/mL; Treatment 2, 6.4 × 106 CFU/mL). The mortality of Treatment 2 (77.5%) was significantly higher than that of Treatment 1 (10%), corresponding with bacterial concentrations. Although the number of deaths varies, intensive deaths were observed within 24 h post infection (hpi) in both bacterial concentration groups. Histopathological analyses revealed that fish tissues were most severely damaged at 24 or 48 hpi, and Treatment 2 was more severe than Treatment 1. A qRT-PCR-based detection method with virulence factor gene empA was established to quantify the bacterial loads in various tissues, and the bacterial loads were the highest at 24 hpi in Treatment 2, and at 48 hpi in Treatment 1. Additionally, the expression levels of complement genes (C5a, C3, C5, and C6), inflammatory factors (IL-1ß, TNF-α, and IL-10), and other immune-related genes (jak2, NF-κB1, stat3, and tlr3) were increased in various tissues after infection in both treatment groups, with most genes being most expressed at 24 or 48 hpi, and expression levels of inflammatory factors in Treatment 2 were higher than those in Treatment 1. Moreover, the expression of C5a was positively correlated with that of proinflammatory cytokines in both bacterial concentration groups. According to the results of this study, 24-48 hpi was a key node for early vibriosis detection and intervention. Compared with the low mortality of Treatment 1, the mass death of fish in Treatment 2 was suggested to be caused by uncontrolled excessive inflammatory reaction induced by the overactivation of complement system, especially C5a. We believe these results could provide theoretical basis for prevention, evaluation, and treatment of vibrio disease in tongue sole aquaculture, and lay a solid foundation for future functional analyses.

C5a drives the inflammatory response with bacterial dose effect by binding to C5aR1 in zebrafish infected with Aeromonas hydrophila.

The complement system is essential to host defense, but its excessive activation caused by severe pathogen invasion is a driving force in adverse inflammatory. The binding of complement component 5a (C5a) and complement component 5a receptor 1 (C5aR1) is the key to trigger complement-mediated inflammatory response in mammals. However, the role of C5a-C5aR1 axis in fish immune response remains obscure. In this study, the role of C5a-C5aR1 axis of zebrafish (Danio rerio) after serious infection with Aeromonas hydrophila was investigated. C5a and C5aR1 of zebrafish were cloned, with CDS sequences of 228 and 1041 bp, respectively, and they were widely expressed in various tissues with the highest expression in the liver and spleen, respectively. The survival of zebrafish was closely correlated to the dose of A. hydrophila. The cytokine storm occurred at high concentrations of A. hydrophila infection. At 24 h post infection (hpi), the expression of C5a and C5aR1 in the spleen increased 26.8-fold and 9.9-fold in treatment group 1 (TG1, 3.0 × 107 CFU/mL) (P < 0.01), and 4.7-fold and 3.4-fold in treatment group 2 (TG2, 1.0 × 107 CFU/mL) (P < 0.05), respectively. Correspondingly, proinflammatory cytokines interleukin-1ß (IL-1ß), interleukin-8 (IL-8), and interleukin-17 (IL-17) were positively correlated to C5a and C5aR1 at mRNA and protein expression levels. The expression of IL-1ß was significantly increased in the spleen at 6 hpi, with a 599.2-fold and 203.2-fold upregulation in TG1 and TG2 (P < 0.001), respectively. Moreover, after inhibition of C5a-C5aR1 binding treated with C5aR1 antagonist (W-54011), zebrafish showed lower expression of C5a, C5aR1, and cytokines, less intestinal damage, and significantly enhancement of survival (P < 0.05) after A. hydrophila challenge. This study revealed that the inflammatory effect of C5a was achieved by binding to C5aR1 in zebrafish, providing novel insights into using C5a-C5aR1 axis as an effective target to reduce bacterial inflammation and disease in fish.

A web-based database server using 43,710 public RNA-seq samples for the analysis of gene expression and alternative splicing in livestock animals.

BACKGROUND: Livestock animals is of great significance in agricultural production. However, the role of specific gene expression, especially alternative splicing in determining phenotype, is not well understood. The livestock research community needs a gene expression and alternative splicing database contributing to livestock genetic improvement. DESCRIPTION: We report the construction of LivestockExp ( https://bioinfo.njau.edu.cn/livestockExp ), a web-based database server for the exploration of gene expression and alternative splicing using 43,710 uniformly processed RNA-seq samples from livestock animals and several relative species across six orders. The database is equipped with basic querying functions and multiple online analysis modules including differential/specific expression analysis, co-expression network analysis, and cross-species gene expression conservation analysis. In addition to the re-analysis of public datasets, users can upload personal datasets to perform co-analysis with public datasets. The database also offers a wide range of visualization tools and diverse links to external databases enabling users to efficiently explore the results and to gain additional insights. CONCLUSION: LivestockExp covers by far the largest number of livestock animal species and RNA-seq samples and provides a valuable data resource and analysis platform for the convenient utilization of public RNA-seq datasets.

Identification of an Epigenetically Marked Locus within the Sex Determination Region of Channel Catfish.

Channel catfish has an XY sex determination system. However, the X and Y chromosomes harbor an identical gene content of 950 genes each. In this study, we conducted comparative analyses of methylome and transcriptome of genetic males and genetic females before gonadal differentiation to provide insights into the mechanisms of sex determination. Differentially methylated CpG sites (DMCs) were predominantly identified on the sex chromosome, most notably within the sex determination region (SDR), although the overall methylation profiles across the entire genome were similar between genetic males and females. The drastic differences in methylation were located within the SDR at nucleotide position 14.0-20.3 Mb of the sex chromosome, making this region an epigenetically marked locus within the sex determination region. Most of the differentially methylated CpG sites were hypermethylated in females and hypomethylated in males, suggesting potential involvement of methylation modification in sex determination in channel catfish. Along with the differential methylation in the SDR, a number of differentially expressed genes within the SDR were also identified between genetic males and females, making them potential candidate genes for sex determination and differentiation in channel catfish.

ASlive: a database for alternative splicing atlas in livestock animals.

BACKGROUND: Alternative splicing is an important biological process whose precision must be tightly regulated during growth and development. Although there are species, disease (e.g. cancers), or study specific databases available in many organisms, no database exists in livestock animals specifically tailored for alternative splicing. DESCRIPTION: We present in this study the development and implementation of a database for alternative splicing atlas in livestock animals (ASlive.org). Using publicly available RNASeq data sets across many tissues, cell types, and biological conditions totaling 28.6 T bases, we built a database of alternative splicing events in five major livestock and poultry animal species (cattle, sheep, pigs, horses, and chickens). The database contains many types of information on alternative splicing events, including basic information such as genomic locations, genes, and event types, quantitative measurements of alternative splicing in the form of percent spliced in (PSI), overlap with known DNA variants, as well as orthologous events across different lineage groups. CONCLUSIONS: This database, the first of its kind in livestock animals, will provide a useful exploratory tool to assist functional annotation of animal genomes.

The Y chromosome sequence of the channel catfish suggests novel sex determination mechanisms in teleost fish.

BACKGROUND: Sex determination mechanisms in teleost fish broadly differ from mammals and birds, with sex chromosomes that are far less differentiated and recombination often occurring along the length of the X and Y chromosomes, posing major challenges for the identification of specific sex determination genes. Here, we take an innovative approach of comparative genome analysis of the genomic sequences of the X chromosome and newly sequenced Y chromosome in the channel catfish. RESULTS: Using a YY channel catfish as the sequencing template, we generated, assembled, and annotated the Y genome sequence of channel catfish. The genome sequence assembly had a contig N50 size of 2.7 Mb and a scaffold N50 size of 26.7 Mb. Genetic linkage and GWAS analyses placed the sex determination locus within a genetic distance less than 0.5 cM and physical distance of 8.9 Mb. However, comparison of the channel catfish X and Y chromosome sequences showed no sex-specific genes. Instead, comparative RNA-Seq analysis between females and males revealed exclusive sex-specific expression of an isoform of the breast cancer anti-resistance 1 (BCAR1) gene in the male during early sex differentiation. Experimental knockout of BCAR1 gene converted genetic males (XY) to phenotypic females, suggesting BCAR1 as a putative sex determination gene. CONCLUSIONS: We present the first Y chromosome sequence among teleost fish, and one of the few whole Y chromosome sequences among vertebrate species. Comparative analyses suggest that sex-specific isoform expression through alternative splicing may underlie sex determination processes in the channel catfish, and we identify BCAR1 as a potential sex determination gene.

Post-transcriptional regulation through alternative splicing after infection with Flavobacterium columnare in channel catfish (Ictalurus punctatus).

Columnaris disease has long been recognized as a serious problem worldwide which affects both wild and cultured freshwater fish including the commercially important channel catfish (Ictalurus punctatus). The fundamental molecular mechanisms of the host immune response to the causative agent Flavobacterium columnare remain unclear, though gene expression analysis after the bacterial infection has been conducted. Alternative splicing, a post-transcriptional regulation process to modulate gene expression and increase the proteomic diversity, has not yet been studied in channel catfish following infection with F. columnare. In this study, genomic information and RNA-Seq datasets of channel catfish were used to characterize the changes of alternative splicing after the infection. Alternative splicing was shown to be induced by F. columnare infection, with 8.0% increase in alternative splicing event at early infection stage. Intriguingly, genes involved in RNA binding and RNA splicing themselves were significantly enriched in differentially alternatively spliced (DAS) gene sets after infection. This finding was consistent with our previous study in channel catfish following infection with Edwardsiella ictaluri. It was suggested to be a universal mechanism that genes involved in RNA binding and splicing were regulated to undergo differential alternative splicing after stresses in channel catfish. Moreover, many immune genes were observed to be differentially alternatively spliced after infection. Further studies need to be performed to get a deeper view of molecular regulation on alternative splicing after stresses, setting a foundation for developing catfish broodstocks with enhanced disease resistance.

GWAS analysis using interspecific backcross progenies reveals superior blue catfish alleles responsible for strong resistance against enteric septicemia of catfish.

Infectious diseases pose significant threats to the catfish industry. Enteric septicemia of catfish (ESC) caused by Edwardsiella ictaluri is the most devastating disease for catfish aquaculture, causing huge economic losses annually. Channel catfish and blue catfish exhibit great contrast in resistance against ESC, with channel catfish being highly susceptible and blue catfish being highly resistant. As such, the interspecific backcross progenies provide an ideal system for the identification of quantitative trait locus (QTL). We previously reported one significant QTL on linkage group (LG) 1 using the third-generation backcrosses, but the number of founders used to make the second- and third-generation backcross progenies was very small. Although the third-generation backcross progenies provided a greater power for fine mapping than the first-generation backcrosses, some major QTL for disease resistance may have been missing due to the small numbers of founders used to produce the higher generation backcrosses. In this study, we performed a genome-wide association study using first-generation backcrosses with the catfish 690 K SNP arrays to identify additional ESC disease resistance QTL, especially those at the species level. Two genomic regions on LG1 and LG23 were determined to be significantly associated with ESC resistance as revealed by a mixed linear model and family-based association test. Examination of the resistance alleles indicated their origin from blue catfish, indicating that at least two major disease resistance loci exist among blue catfish populations. Upon further validation, markers linked with major ESC disease resistance QTL should be useful for marker-assisted introgression, allowing development of highly ESC resistant breeds of catfish.

Genome-wide association analysis of intra-specific QTL associated with the resistance for enteric septicemia of catfish.

Disease resistance is one of the most important traits for aquaculture industry. For catfish industry, enteric septicemia of catfish (ESC), caused by the bacterial pathogen Edwardsiella ictaluri, is the most severe disease, causing enormous economic losses every year. In this study, we used three channel catfish families with 900 individuals (300 fish per family) and the 690K catfish SNP array, and conducted a genome-wide association study to detect the quantitative trait loci (QTL) associated with ESC resistance. Three significant QTL, with two of located on LG1 and one on LG26, and three suggestive QTL located on LG1, LG3, and LG21, respectively, were identified to be associated with ESC resistance. With a well-assembled- and -annotated reference genome sequence, genes around the involved QTL regions were identified. Among these genes, 37 genes had known functions in immunity, which may be involved in ESC resistance. Notably, nlrc3 and nlrp12 identified here were also found in QTL regions of ESC resistance in the channel catfish × blue catfish interspecific hybrid system, suggesting this QTL was operating within both intra-specific channel catfish populations and interspecific hybrid backcross populations. Many of the genes of the Class I MHC pathway, for mediated antigen processing and presentation, were found in the QTL regions. The positional correlation found in this study and the expressional correlation found in previous studies indicated that Class I MHC pathway was significantly associated with ESC resistance. This study validated one QTL previously identified using the second and fourth generation of the interspecific hybrid backcross progenies, and identified five additional QTL among channel catfish families. Taken together, it appears that there are only a few major QTL for ESC disease resistance, making marker-assisted selection an effective approach for genetic improvements of ESC resistance.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research.

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

GWAS analysis of QTL for enteric septicemia of catfish and their involved genes suggest evolutionary conservation of a molecular mechanism of disease resistance.

Disease problems cause major economic losses for the aquaculture industries. In catfish, enteric septicemia of catfish (ESC), caused by the bacterial pathogen Edwardsiella ictaluri, is the leading disease problem, causing tens of millions of dollars of annual economic losses. In this study, we conducted a genome-wide association study to determine quantitative trait loci (QTL) for resistance against ESC using an interspecific hybrid system. Five hundred fish were used in the analysis and 192 phenotypic extremes were used for genotyping with the catfish 250K SNP arrays. A genomic region on linkage group (LG) 1 was found significantly associated with ESC disease resistance. In addition, two suggestively associated QTL for ESC resistance were identified on LG 12 and LG 16. The nlrc3 duplicates were identified within all the three QTL, suggesting their importance in association with the QTL. Within the significant QTL on LG 1, 16 genes with known functions in immunity were identified. Of particular interest is the nck1 gene nearby the most significantly associated SNP. Nck1 was known to function as an adaptor to facilitating the pathogenesis of enteropathogenic Escherichia coli (EPEC) in humans. E. ictaluri and EPEC pathogens belong to the same bacterial family and share many common characteristics. The fact that nck1 is mapped in the QTL and that it was significantly upregulated in channel catfish intestine after ESC challenge suggested its candidacy of being involved in resistance/susceptibility of ESC.

The survival, gene expression, and DNA methylation of Paralichthys olivaceus impacted by the decay of green tide and bacterial infection in both laboratory and field simulation experiments.

The recurring appearance of Ulva prolifera green tides has become a pressing environmental issue, especially for marine transportation, tourism, and aquaculture in the stage of decomposition. An abundance of decaying U. prolifera leads to water acidification, hypoxia and pathogenic microorganism proliferation, threatening marine germplasm resources, particularly benthic organisms with weak escape ability. Epigenetic modification is considered to be one of the molecular mechanisms involved in the plastic adaptive response to environmental changes. However, few studies concerning the specific impact of decaying green tide on benthic animals at the epigenetic level. In this study, decomposing algal effluents of U. prolifera, sediments containing uncorrupted U. prolifera, pathogenic microorganism were considered as impact factors, to reveal the effect of decaying U. prolifera on marine economic benthic species, Paralichthys olivaceus, using both field and laboratory simulation experiments. Field simulation experiment showed higher mortality rates and serious histopathological damage than the laboratory simulation experiment. And both the decaying U. prolifera and the sediment containing U. prolifera were harmful to P. olivaceus. Genome-wide DNA methylation and transcription correlation analyses showed that the response of P. olivaceus to green tide stress and bacterial infection was mainly mediated by immune signaling pathways such as PI3K-Akt signaling pathway. DNA methylation regulates the expression of immune-related genes involved in the PI3K-Akt signaling pathway, which enables P. olivaceus to adapt to the adverse environmental stresses by resisting apoptosis. In summary, this research analyzed the potential role of P. olivaceus in decaying U. prolifera, which is of great significance for understanding the impact of decaying green tide on marine commercial fish and also provides some theoretical guidance for the proliferation and release of fish seedlings.

Multi-omic characterization of allele-specific regulatory variation in hybrid pigs.

Hybrid mapping is a powerful approach to efficiently identify and characterize genes regulated through mechanisms in cis. In this study, using reciprocal crosses of the phenotypically divergent Duroc and Lulai pig breeds, we perform a comprehensive multi-omic characterization of regulatory variation across the brain, liver, muscle, and placenta through four developmental stages. We produce one of the largest multi-omic datasets in pigs to date, including 16 whole genome sequenced individuals, as well as 48 whole genome bisulfite sequencing, 168 ATAC-Seq and 168 RNA-Seq samples. We develop a read count-based method to reliably assess allele-specific methylation, chromatin accessibility, and RNA expression. We show that tissue specificity was much stronger than developmental stage specificity in all of DNA methylation, chromatin accessibility, and gene expression. We identify 573 genes showing allele specific expression, including those influenced by parent-of-origin as well as allele genotype effects. We integrate methylation, chromatin accessibility, and gene expression data to show that allele specific expression can be explained in great part by allele specific methylation and/or chromatin accessibility. This study provides a comprehensive characterization of regulatory variation across multiple tissues and developmental stages in pigs.

Exploring the effect of Ulva prolifera decay on the immune tissue of Paralichthys olivaceus based on transcriptomics and histopathological analysis.

For 17 consecutive years, the outbreak of Ulva prolifera in the South Yellow Sea area of China has caused significant negative impacts on coastal ecological environment. However, its specific influence on fish immunity is rare. In this study, the juvenile Paralichthys olivaceus was exposed to fresh U. prolifera algae (FU) and decomposing algal effluent (DU). After short-term stress for 14 days, the histopathological and transcriptome analysis were performed to study the effect of U. prolifera decay on P. olivaceus. Histopathological analysis found that the liver, spleen and head kidneys of P. olivaceus were damaged after the short-term stress. The transcriptome results showed that the steroid biosynthesis signaling pathway and the PI3K-Akt signaling pathway were significantly enriched. Some immune related genes, including c1qc-like, dusp1, dusp16, HSP90 and metabolic related genes serotransferrin, were differentially expressed. These results highlighted the harmfulness of U. prolifera on marine fish, setting a solid foundation for further analyses.

Genome-wide characterization of integrin (ITG) gene family and their expression profiling in half-smooth tongue sole (Cynoglossus semilaevis) upon Vibrio anguillarum infection.

Integrins (ITGs) are transmembrane heterodimer receptors with ITGα subunit and ITGß subunit, participating in various physiological processes, including immunity. At present, systematic research on ITGs in teleost is scarce, especially in half-smooth tongue sole (Cynoglossus semilaevis). In this study, a set of 28 ITG genes in half-smooth tongue sole have been identified and characterized. The phylogenetic analysis showed that ITGα and ITGß subunits were respectively classified into five and two clusters, consistent with previous studies. The selection pressure analysis indicated that most of ITG genes were under purifying selection, except for ITGα11b and ITGαL with positive selection. The expression profiles of eight selected ITG genes, including ITGα1, ITGα5, ITGα8, ITGα11, ITGß1, ITGß2, ITGß3, and ITGß8, were analyzed in healthy tissues and after infection with Vibrio anguillarum, revealed their implications in immune response. The study provided a comprehensive characterization and expression analysis of ITG genes in half-smooth tongue sole, setting a solid foundation for further functional studies and promising potential in disease control.

Regulatory mechanism of miR-722 on C5aR1 and its functions against bacterial inflammation in half-smooth tongue sole (Cynoglossus semilaevis).

MicroRNAs (miRNAs) are small non-coding RNAs involved in various biological processes, including immunity. Previously, we investigated the miRNAs of half-smooth tongue sole (Cynoglossus semilaevis) and found that miR-722 (designated Cse-miR-722) was significantly differentially expressed after infection with Vibrio anguillarum, reflecting its importance in immune response. Our preliminary bioinformatic analysis suggested that Cse-miR-722 could target C5aR1 (designated CsC5aR1), which was known to play crucial roles in complement activation and inflammatory response, as a receptor of C5a. However, the underlying mechanisms of their interactions and specific functions in inflammatory and immune response are still enigmas. In this study, we successfully cloned the precursor sequence of Cse-miR-722 (94 bp) and the full length of CsC5aR1 (1541 bp, protein molecular weight 39 kDa). The target gene of Cse-miR-722 was verified as CsC5aR1 by a dual luciferase reporter assay, and Cse-miR-722 was confirmed to regulate CsC5aR1 at the protein level using quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting. The expression of CsC5aR1 and Cse-miR-722 in liver cells and four immune tissues of half-smooth tongue sole changed significantly after LPS stimulation and infection with V. anguillarum. To explore the functional role of Cse-miR-722 in half-smooth tongue sole, we performed both in vitro and in vivo experiments. Cse-miR-722 was observed to affect phagocytosis and respiratory burst activity of macrophages by regulating CsC5aR1 in half-smooth tongue sole. Furthermore, we found that Cse-miR-722 regulated the expression of CsC5aR1, CsC5a, and the inflammatory factors CsIL1-ß, CsIL6, CsIL8, and CsTNF-α both in vitro and in vivo. In addition, Cse-miR-722 reduced mortality and pathological damage. This study clarified the regulatory mechanism of Cse-miR-722 on CsC5aR1 and provided insight into the regulatory roles of Cse-miR-722 in immune responses, laying a theoretical foundation for the feasibility of using miR-722 to prevent and control bacterial diseases in teleost.