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.

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.

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.

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.

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.

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.

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.

The SWine IMputation (SWIM) haplotype reference panel enables nucleotide resolution genetic mapping in pigs.

Genetic mapping to identify genes and alleles associated with or causing economically important quantitative trait variation in livestock animals such as pigs is a major goal in animal genetic improvement. Despite recent advances in high-throughput genotyping technologies, the resolution of genetic mapping in pigs remains poor due in part to the low density of genotyped variant sites. In this study, we overcame this limitation by developing a reference haplotype panel for pigs based on 2259 whole genome-sequenced animals representing 44 pig breeds. We evaluated software combinations and breed composition to optimize the imputation procedure and achieved an average concordance rate in excess of 96%, a non-reference concordance rate of 88%, and an r2 of 0.85. We demonstrated in two case studies that genotype imputation using this resource can dramatically improve the resolution of genetic mapping. A public web server has been developed to allow the pig genetics community to fully utilize this resource. We expect this resource to facilitate genetic mapping and accelerate genetic improvement in pigs.

Genomic imprinting-like monoallelic paternal expression determines sex of channel catfish.

The X and Y chromosomes of channel catfish have the same gene contents. Here, we report allelic hypermethylation of the X chromosome within the sex determination region (SDR). Accordingly, the X-borne hydin-1 gene was silenced, whereas the Y-borne hydin-1 gene was expressed, making monoallelic expression of hydin-1 responsible for sex determination, much like genomic imprinting. Treatment with a methylation inhibitor, 5-aza-dC, erased the epigenetic marks within the SDR and caused sex reversal of genetic females into phenotypic males. After the treatment, hydin-1 and six other genes related to cell cycle control and proliferative growth were up-regulated, while three genes related to female sex differentiation were down-regulated in genetic females, providing additional support for epigenetic sex determination in catfish. This mechanism of sex determination provides insights into the plasticity of genetic sex determination in lower vertebrates and its connection with temperature sex determination where DNA methylation is broadly involved.

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.

FishExp: A comprehensive database and analysis platform for gene expression and alternative splicing of fish species.

The publicly archived RNA-seq data has grown exponentially, while its valuable information has not yet been fully discovered and utilized, such as alternative splicing and its integration with gene expression. This is especially true for fish species which play important roles in ecology, research and the food industry. Furthermore, there is a lack of online platform to analyze users' new data individually and jointly with existing data for the comprehensive analysis of alternative splicing and gene expression. Here, we present FishExp, a web-based data platform covering gene expression and alternative splicing in 26,081 RNA-seq experiments from 44 fishes. It allows users to query the data in a variety of ways, including gene identifier/symbol, functional term, and BLAST alignment. Moreover, users can customize experiments and tools to perform differential/specific expression and alternative splicing analysis, co-expression and cross-species analysis. In addition, functional enrichment is provided to confer biological significance. Notably, users are allowed to submit their own data and perform various analyses using the new data alone or alongside existing data in FishExp. Results of retrieval and analysis can be visualized on the gene-, transcript- and splicing event-level webpage in a highly interactive and intuitive manner. All data in FishExp can be downloaded for more in-depth analysis. The manually curated sample information, uniform data processing and various tools make it efficient for users to gain new insights from these large data sets, facilitating scientific hypothesis generation. FishExp is freely accessible at https://bioinfo.njau.edu.cn/fishExp.

Feminization of channel catfish with 17ß-oestradiol involves methylation and expression of a specific set of genes independent of the sex determination region.

Exogenous oestrogen 17ß-oestradiol (E2) has been shown to effectively induce feminization in teleosts. However, the molecular mechanisms underlying the process remain unclear. Here, we determined global DNA methylation and gene expression profiles of channel catfish (Ictalurus punctatus) during early sex differentiation after E2 treatment. Overall, the levels of global DNA methylation after E2 treatment were not significantly different from those of controls. However, a specific set of genes were differentially methylated, which included many sex differentiation-related pathways, such as MARK signalling, adrenergic signalling, Wnt signalling, GnRH signalling, ErbB signalling, and ECM-receptor interactions. Many genes involved in these pathways were also differentially expressed after E2 treatment. Specifically, E2 treatments resulted in upregulation of female-related genes and downregulation of male-related genes in genetic males during sex reversal. However, E2-induced sex reversal did not cause sex-specific changes in methylation profiles or gene expression within the sex determination region (SDR) on chromosome 4, suggesting that E2-induced sex reversal was a downstream process independent of the sex determination process that was regulated by sex-specific methylation within the SDR.

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.

Allelically and Differentially Expressed Genes After Infection of Edwardsiella ictaluri in Channel Catfish as Determined by Bulk Segregant RNA-Seq.

Identification of genetic markers associated with resistance against enteric septicemia of catfish (ESC) is of great interest for genetic enhancement programs of catfish. In the present study, bulk segregant RNA-Seq analysis was applied to determine differentially expressed genes and alleles after ESC infection. Here we report three genomic regions on LG1, LG12, and LG26, containing significant single-nucleotide polymorphisms (SNPs). These genomic regions aligned well with quantitative trait loci (QTL) previously identified. Within the QTL regions, eleven genes were found to be differentially regulated between phenotypic bulks. Importantly, the QTL on linkage group 1 (LG1) were found to be expressed in the liver, whereas the QTL on LG12 and LG26 were expressed in the intestine, suggesting multiple mechanisms of ESC resistance. It is apparent that apolipoproteins may be important for ESC resistance as the QTL on LG1 included the 14-kDa apolipoprotein genes that are both allelically expressed and differentially expressed between the resistant and susceptible bulks. Traf2 and NCK-interacting protein kinase (TNIK) were found in the QTL on LG12, and it was downregulated in resistant fish, suggesting the importance of NCK downregulation in ESC resistance, as previously reported. In addition, we observed divergent gene expression patterns between the liver and intestine after infection. Immune/inflammatory-related processes were overrepresented from liver DEGs, while those DEGs identified from intestine were enriched for proteolysis and wounding processes. Taken together, the BSR-Seq analysis presented here advanced the knowledge of ESC resistance, providing information of not only positions of QTL but also genes and their differential expression between resistant and susceptible fish, making it one step closer to the identification of the causal genes for ESC resistance.

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.

A bioenergetic shift is required for spermatogonial differentiation.

A bioenergetic balance between glycolysis and mitochondrial respiration is particularly important for stem cell fate specification. It however remains to be determined whether undifferentiated spermatogonia switch their preference for bioenergy production during differentiation. In this study, we found that ATP generation in spermatogonia was gradually increased upon retinoic acid (RA)-induced differentiation. To accommodate this elevated energy demand, RA signaling concomitantly switched ATP production in spermatogonia from glycolysis to mitochondrial respiration, accompanied by increased levels of reactive oxygen species. Disrupting mitochondrial respiration significantly blocked spermatogonial differentiation. Inhibition of glucose conversion to glucose-6-phosphate or pentose phosphate pathway also repressed the formation of c-Kit+ differentiating germ cells, suggesting that metabolites produced from glycolysis are required for spermatogonial differentiation. We further demonstrated that the expression levels of several metabolic regulators and enzymes were significantly altered upon RA-induced differentiation, with both RNA-seq and quantitative proteomic analyses. Taken together, our data unveil a critically regulated bioenergetic balance between glycolysis and mitochondrial respiration that is required for spermatogonial proliferation and differentiation.

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.

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.