The cytochrome P450 genes of channel catfish: their involvement in disease defense responses as revealed by meta-analysis of RNA-Seq data sets

Genome sequencing, comparative analysis, and gene expression responses of cytochrome P450 genes in Oryzias curvinotus provide insights into environmental adaptation

The mangrove fish (Oryzias curvinotus) serves as a model for researching environmental adaptation and sexual development. To further such research, we sequenced and assembled a high-quality 842 Mb reference genome for O. curvinotus. Comparative genomic analysis revealed 891 expanded gene families, including significantly expanded cytochrome P450 (CYP) detoxification genes known to be involved in xenobiotic defense. We identified 69 O. curvinotus CYPs (OcuCYPs) across 18 families and 10 clans using multiple methods. Extensive RNA-seq and qPCR analysis demonstrated diverse spatiotemporal expression patterns of OcuCYPs by developmental stage, tissue type, sex, and pollutant exposure (17β-estradiol (E2) and testosterone (MT)). Many OcuCYPs exhibited sexual dimorphism in gonads, suggesting reproductive roles in steroidogenesis, while their responsiveness to model toxicants indicates their importance in environmental adaptation through enhanced detoxification. Pathway analysis highlighted expanded CYP genes in arachidonic acid metabolism, drug metabolism, and steroid hormone biosynthesis. This chromosome-level genomic resource provides crucial biological insights to elucidate the functional roles of expanded CYPs in environmental adaptation, sexual development, early life history, and conservation in the anthropogenically impacted mangrove habitats of O. curvinotus. It also enables future ecotoxicology research leveraging O. curvinotus as a pollution sentinel species.

Oxidative stress, gene expression and histopathology of cultured gilthead sea bream (Sparus aurata) naturally co-infected with Ergasilus sieboldi and Vibrio alginolyticus

BACKGROUND

Parasitic and bacterial co-infections have been associated with increasing fish mortalities and severe economic losses in aquaculture through the past three decades. The aim of this study was to evaluate the oxidative stress, histopathology, and immune gene expression profile of gilthead sea bream (Sparus aurata) co-infected with Ergasilus sieboldi and Vibrio alginolyticus.

RESULTS

Vibrio alginolyticus and Ergasilus sieboldi were identified using 16 S rRNA and 28 S rRNA sequencing, respectively. The collagenase virulence gene was found in all Vibrio alginolyticus isolates, and the multiple antimicrobial resistance index ranged from 0.286 to 0.857. Oxidant-antioxidant parameters in the gills, skin, and muscles of naturally infected fish revealed increased lipid peroxidation levels and a decrease in catalase and glutathione antioxidant activities. Moreover, naturally co-infected gilthead sea bream exhibited substantial up-regulation of il-1β, tnf-α, and cyp1a1. Ergasilus sieboldi encircled gill lamellae with its second antennae, exhibited severe gill architectural deformation with extensive eosinophilic granular cell infiltration. Vibrio alginolyticus infection caused skin and muscle necrosis in gilthead sea bream.

CONCLUSION

This study described some details about the gill, skin and muscle tissue defense mechanisms of gilthead sea bream against Ergasilus sieboldi and Vibrio alginolyticus co-infections. The prevalence of co-infections was 100%, and no resistant fish were detected. These co-infections imbalance the health status of the fish by hampering the oxidant-antioxidant mechanisms and proinflammatory/inflammatory immune genes to a more detrimental side. Our results suggest that simultaneous screening for bacterial and parasitic pathogens should be considered.

Functional and molecular characterization of the Atlantic salmon gill epithelium cell line ASG-10; a tool for in vitro gill research

Fish gills are not only the respiratory organ, but also essential for ion-regulation, acid-base control, detoxification, waste excretion and host defense. Multifactorial gill diseases are common in farmed Atlantic salmon, and still poorly understood. Understanding gill pathophysiology is of paramount importance, but the sacrifice of large numbers of experimental animals for this purpose should be avoided. Therefore, in vitro models, such as cell lines, are urgently required to replace fish trials. An Atlantic salmon gill epithelial cell line, ASG-10, was established at the Norwegian Veterinary institute in 2018. This cell line forms a monolayer expressing cytokeratin, e-cadherin and desmosomes, hallmarks of a functional epithelial barrier. To determine the value of ASG-10 for comparative studies of gill functions, the characterization of ASG-10 was taken one step further by performing functional assays and comparing the cell proteome and transcriptome with those of gills from juvenile freshwater Atlantic salmon. The ASG-10 cell line appear to be a homogenous cell line consisting of epithelial cells, which express tight junction proteins. We demonstrated that ASG-10 forms a barrier, both alone and in co-culture with the Atlantic salmon gill fibroblast cell line ASG-13. ASG-10 cells can phagocytose and express several ATP-binding cassette transport proteins. Additionally, ASG-10 expresses genes involved in biotransformation of xenobiotics and immune responses. Taken together, this study provides an overview of functions that can be studied using ASG-10, which will be an important contribution to in vitro gill epithelial research of Atlantic salmon.

Cytochrome P450 superfamily in spotted sea bass: Genome-wide identification and expression profiles under trichlorfon and environmental stresses

Cytochrome P450s (CYPs), as one of the most diverse enzyme superfamilies in nature, play critical functions in antioxidant reactions against endogenous and exogenous compounds. In this study, we performed genome-wide characterization of CYP superfamily members and analyzed their expression patterns under several abiotic stresses in spotted sea bass, which is known as an economically important fish species in the Chinese aquaculture industry. A total of 55 CYP genes were identified and divided into 17 families within 10 clans. The analysis of phylogeny, gene structure, and syntenic relationships provided evidence for the evolution of CYP genes and confirmed their annotation and orthology. The expression of CYP genes was examined in the liver during trichlorfon stress using quantitative real-time PCR. The results showed that 20 tested CYP genes displayed significant mRNA expression changes, indicating that they may play crucial roles in the metabolism of trichlorfon and can be potential biomarkers for trichlorfon pollution. Moreover, by screening transcriptomic databases, 10, 3 and 19 CYP genes exhibited differential expression patterns in response to hypoxia, alkalinity and heat stress, respectively. Taken together, this study provided insights into the regulation of CYP genes by toxicological and environmental stresses, laid basis for extensive functional studies of the CYP superfamily in spotted sea bass and other teleost species.

Proteomic analysis of liver tissue reveals Aeromonas hydrophila infection mediated modulation of host metabolic pathways in Labeo rohita

Aeromonas hydrophila (Ah) is a Gram-negative bacterium and a serious global pathogen causing Motile Aeromonas Septicaemia (MAS) in fish leading to global loss in aquaculture. Investigation of the molecular alterations of host tissues such as liver could be a powerful approach to identify mechanistic and diagnostic immune signatures of disease pathogenesis. We performed a proteomic analysis of Labeo rohita liver tissue to examine the protein dynamics in the host cells during Ah infection. The proteomic data was acquired using two strategies; discovery and targeted proteomics. Label-free quantification was performed between Control and challenged group (AH) to identify the differentially expressed proteins (DEPs). A total of 2525 proteins were identified and 157 were DEPs. DEPs include metabolic enzymes (CS, SUCLG2), antioxidative proteins, cytoskeletal proteins and immune related proteins (TLR3, CLEC4E). Pathways like lysosome pathway, apoptosis, metabolism of xenobiotics by cytochrome P450 were enriched by downregulated proteins. However, upregulated proteins majorly mapped to innate immune system, signaling of B cell receptor, proteosome pathway, ribosome, carbon metabolism and protein processing in ER. Our study would help in exploring the role of Toll-like receptors, C-type lectins and, metabolic intermediates like citrate and succinate in Ah pathogenesis to understand the Ah infection in fish. SIGNIFICANCE: Bacterial diseases such as motile aeromonas septicaemia (MAS) are among the most serious problems in aquaculture industry. Small molecules that target the metabolism of the host have recently emerged as potential treatment possibilities in infectious diseases. However, the ability to develop new therapies is hampered due to lack of knowledge about pathogenesis mechanisms and host-pathogen interactions. We examined alterations in the host proteome during MAS caused by Aeromonas hydrophila (Ah) infection, in Labeo rohita liver tissue to find cellular proteins and processes affected by Ah infection. Upregulated proteins belong to innate immune system, signaling of B cell receptor, proteosome pathway, ribosome, carbon metabolism and protein processing. Our work is an important step towards leveraging host metabolism in targeting the disease by providing a bigger picture on proteome pathology correlation during Ah infection.

Regions of genetic divergence in depth-separated Sebastes rockfish species pairs: Depth as a potential driver of speciation

Depth separation is a proposed driver of speciation in marine fishes, with marine rockfish (genus Sebastes) providing a potentially informative study system. Sebastes rockfishes are commercially and ecologically important. This genus encompasses more than one hundred species and the ecological and morphological variance between these species provides opportunity for identifying speciation-driving adaptations, particularly along a depth gradient. A reduced-representation sequencing method (ddRADseq) was used to compare 95 individuals encompassing six Sebastes species. In this study, we sought to identify regions of divergence between species that were indicative of divergent adaptation and reproductive barriers leading to speciation. A pairwise comparison of S. chrysomelas (black-and-yellow rockfish) and S. carnatus (gopher rockfish) F_ST values revealed three major regions of elevated genomic divergence, two of which were also present in the S. miniatus (vermilion rockfish) and S. crocotulus (sunset rockfish) comparison. These corresponded with regions of both elevated D_XY values and reduced nucleotide diversity in two cases, suggesting a speciation-with-gene-flow evolutionary model followed by post-speciation selective sweeps within each species. Limited whole-genome resequencing was also performed to identify mutations with predicted effects between S. chrysomelas and S. carnatus. Within these islands, we identified important SNPs in genes involved in immune function and vision. This supports their potential role in speciation, as these are adaptive vectors noted in other organisms. Additionally, changes to genes involved in pigment expression and mate recognition shed light on how S. chrysomelas and S. carnatus may have become reproductively isolated.

The genome of walking catfish Clarias magur (Hamilton, 1822) unveils the genetic basis that may have facilitated the development of environmental and terrestrial adaptation systems in air-breathing catfishes

The walking catfish Clarias magur (Hamilton, 1822) (magur) is an important catfish species inhabiting the Indian subcontinent. It is considered as a highly nutritious food fish and has the capability to walk to some distance, and survive a considerable period without water. Assembly, scaffolding and several rounds of iterations resulted in 3,484 scaffolds covering ∼94% of estimated genome with 9.88 Mb largest scaffold, and N50 1.31 Mb. The genome possessed 23,748 predicted protein encoding genes with annotation of 19,279 orthologous genes. A total of 166 orthologous groups represented by 222 genes were found to be unique for this species. The Computational Analysis of gene Family Evolution (CAFE) analysis revealed expansion of 207 gene families and 100 gene families have rapidly evolved. Genes specific to important environmental and terrestrial adaptation, viz. urea cycle, vision, locomotion, olfactory and vomeronasal receptors, immune system, anti-microbial properties, mucus, thermoregulation, osmoregulation, air-breathing, detoxification, etc. were identified and critically analysed. The analysis clearly indicated that C. magur genome possessed several unique and duplicate genes similar to that of terrestrial or amphibians’ counterparts in comparison to other teleostean species. The genome information will be useful in conservation genetics, not only for this species but will also be very helpful in such studies in other catfishes.

Cloning and characterization of estrogen hydroxylase (cyp1a1 and cyp1b1) genes in the stinging catfish Heteropneustes fossilis and induction of mRNA expression during final oocyte maturation

Estrogen hydroxylases (EHs) are cytochrome P450 Family 1 (Cyp1, Clan 2) proteins involved in estrogen hydroxylations at 2-, 4- or 16- carbon positions to form catecholestrogens. EHs are encoded by CYP1A1, CYP1A2 and CYP1B1 in mammals. In the catfish Heteropneustes fossilis, cyp1a1 and cyp1b1 cDNAs were cloned and characterized from liver and ovary. The cyp1a1 cDNA is 2071 bp long and codes for a 518 amino acids (aa) long protein. The cloned cyp1b1 cDNA is 1927 bp long and codes for a 509 residue protein. The deduced proteins clustered distinctly into teleost Cyp1a1 and Cyp1b1 clades, distinct from the tetrapod clusters and featured common function domains and homology with other teleost proteins. In the qPCR assay, the transcripts were the most abundant in the liver, followed by brain and ovary, and moderate in gill, kidney and muscle. Evidence was presented to show the involvement of the genes in reproduction. Expression of brain and ovarian transcripts showed significant seasonal variations with the highest abundance in the spawning phase. In situ hybridization showed the transcripts in the follicular layer (theca and granulosa) of the ovarian follicles. Periovulatory changes in the expression cyp1a1 and cyp1b1 were obtained during final oocyte maturation (FOM) and ovulation induced by human chorionic gonadotropin (hCG), both in vivo and in vitro, and by 2-hydroxyestradiol-17β (catecholestrogen) in vitro. In the brain, the transcript levels increased with time but in the ovary, the increase was maximal at 16 h and decreased at 24 h. The periovulatory activation of the cyp1 genes was reported in this study and discussed on the basis of complex regulation of FOM and ovulation.

The reference genome of the selfing fish Kryptolebias hermaphroditus: Identification of phases I and II detoxification genes

The selfing fish Kryptolebias hermaphroditus has unique reproductive system for self-fertilization, making genetically homozygous offsprings. Here, we report on high density genetic map-based genome assembly for the K. hermaphroditus Panama line (PanRS). The numbers of scaffolds were 5212 and the genome was 683,992,224 bp (N50 = 27.45 Mb). The length of anchored scaffold onto 24 linkage groups was 652,231,070 bp (95.3% of genome) with 0.01% of the gap and 39.33% of GC content and complete Benchmarking Universal Single-Copy Orthologs value was 96.6%. The numbers of annotated genes were 36,756 (average gene length 1368 bp) with the GC content of 54.1%. To examine the difference between the two sister species in the genus Kryptolebias, we compared the genomes of K. hermaphroditus PanRS and Kryptolebias marmoratus PAN line on the composition of transposable elements. To demonstrate applications of genome library, phase I and II detoxification related gene families have been analyzed, and compared the syntenies containing loci of CYP and GST genes on linkage groups. This K. hermaphroditus genome information will be helpful for a better understanding on genome-wide mechanistic view of detoxification and antioxidant-related genes over evolution in the view of fish environmental ecotoxicology.

Effects of the sulfated polysaccharide ulvan from Ulva ohnoi on the modulation of the immune response in Senegalese sole (Solea senegalensis)

Sulfated polysaccharides derived from green seaweeds exhibit many beneficial biological activities and have great potential to be used as nutraceutical in aquaculture. In this work, we evaluated the effects of the sulfated polysaccharide ulvan from Ulva ohnoi on Senegalese sole (Solea senegalensis) juveniles at the transcriptomic level. Cytotoxicity assay performed in liver primary cell cultures from sole determined that the different ulvan concentrations assayed did not impair cell viability. Juveniles were intraperitoneally (IP) injected with ulvan (0.5 mg/fish) followed by a challenge with Photobacterium damselae subsp. piscicida (Phdp) at 7 days. RNASeq analyses at 2 days post injection (dpi) revealed that 402 transcripts were differentially expressed in liver between ulvan IP injected and control groups before the challenge. Genes related to bacterial and antiviral defence, complement system, chemokines, proteasomes and antigen presentation were upregulated in ulvan treated groups. A detailed expression analysis of sixteen genes related to innate and adaptive immune system was performed in two systemic tissues: liver and spleen. Ulvan injection provoked the upregulation of tlr22 and a transient inflammatory response was initiated in both liver and spleen at 2 dpi. As consequence, expression of acute phase proteins, antimicrobial peptides and complement genes was induced. Moreover, expression of mhcI, mhcII, psmb10 and bcl6 was also induced 2 dpi. At 2 dpi with Phdp, inflammatory cytokines and genes related to bacterial and antiviral defense, iron metabolism, complement system and antigen presentation were differentially modulated in survival juveniles previously IP injected with ulvan. Moreover, mortality was retarded in ulvan treated juveniles. These results provide new evidence about the role of ulvan as a bioactive compound with immunomodulatory activity in Senegalese sole as well as its possible use as vaccine adjuvant against Phdp. This is the first published study that evaluates the transcriptomic response of Senegalese sole IP injected with ulvan.

Integrated analysis of mRNA-miRNA expression in Tilapia infected with Tilapia lake virus (TiLV) and identifies primarily immuneresponse genes

We investigated differential gene expression in Tilapia infected with the Tilapia Lake virus (TiLV).We used high-throughput sequencing to identify mRNAs and miRNAs involved in TiLV infection progression We identified 25,359 differentially expressed genes that included 863 new genes. We identified 1770, 4142 and 4947 differently expressed genes comparing non-infected controls with 24 and 120 h infections and between the infected groups, respectively. These genes were enriched to 291 GO terms and 62 KEGG pathways and included immune system progress and virion genes. High-throughput miRNA sequencing identified 316 conserved miRNAs, 525 known miRNAs and 592 novel miRNAs. Furthermore, 138, 198 and 153 differently expressed miRNAs were found between the 3 groups listed above, respectively. Target prediction revealed numerous genes including erythropoietin isoform X2, double-stranded RNA-specific adenosine deaminase isoform X1, bone morphogenetic protein 4 and tapasin-related protein that are involved in immune responsiveness. Moreover, these target genes overlapped with differentially expressed mRNAs obtained from RNA-seq. These target genes were significantly enriched to GO terms and KEGG pathways including immune system progress, virion and Wnt signaling pathways. Expression patterns of differentially expressed mRNA and miRNAs were validated in 20 mRNA and 19 miRNAs by qRT-PCR. We also were able to construct a miRNA-mRNA target network that can further understand the molecular mechanisms on the pathogenesis of TiLV and guide future research in developing effective agents and strategies to combat TiLV infections in Tilapia.

The round goby genome provides insights into mechanisms that may facilitate biological invasions

Background

The invasive benthic round goby (Neogobius melanostomus) is the most successful temperate invasive fish and has spread in aquatic ecosystems on both sides of the Atlantic. Invasive species constitute powerful in situ experimental systems to study fast adaptation and directional selection on short ecological timescales and present promising case studies to understand factors involved the impressive ability of some species to colonize novel environments. We seize the unique opportunity presented by the round goby invasion to study genomic substrates potentially involved in colonization success.

Results

We report a highly contiguous long-read-based genome and analyze gene families that we hypothesize to relate to the ability of these fish to deal with novel environments. The analyses provide novel insights from the large evolutionary scale to the small species-specific scale. We describe expansions in specific cytochrome P450 enzymes, a remarkably diverse innate immune system, an ancient duplication in red light vision accompanied by red skin fluorescence, evolutionary patterns of epigenetic regulators, and the presence of osmoregulatory genes that may have contributed to the round goby’s capacity to invade cold and salty waters. A recurring theme across all analyzed gene families is gene expansions.

Conclusions

The expanded innate immune system of round goby may potentially contribute to its ability to colonize novel areas. Since other gene families also feature copy number expansions in the round goby, and since other Gobiidae also feature fascinating environmental adaptations and are excellent colonizers, further long-read genome approaches across the goby family may reveal whether gene copy number expansions are more generally related to the ability to conquer new habitats in Gobiidae or in fish.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0731-8) contains supplementary material, which is available to authorized users.

Benzothiazole inhibits the growth of Phytophthora capsici through inducing apoptosis and suppressing stress responses and metabolic detoxification

Benzothiazole (BZO) is an antimicrobial secondary metabolite volatilized by many plants and microbes. However, the mechanism of BZO against phytopathogens is still unclear. Here, we found that BZO has antimicrobial activity against the oomycete pathogen Phytophthora capsici. Transcriptome and proteome analyses demonstrated that BZO significantly suppressed the expression of genes and proteins involved in morphology, abiotic stress defense and detoxification, but induced the activity of apoptosis. Annexin V-FITC/PI staining confirmed that the process of apoptosis was significantly induced by BZO at concentration of 150 mg L^-1. FITC-phalloidin actin-cytoskeleton staining combined with hyphal cell wall staining and hyphal ultrastructure studies further confirmed that BZO disrupted the cell membrane and hyphal morphology through disrupting the cytoskeleton, eventually inhibiting the growth of hyphae. These data demonstrated that BZO has multiple modes of action and may act as potential leading compound for the development of new oomycete fungicides. These results also showed that the combination of transcriptomic and proteomic approaches was a useful method for exploring the novel antifungal mechanisms of natural compounds.

Genome-wide identification of HSP70/110 genes in sea cucumber Apostichopus japonicus and comparative analysis of their involvement in aestivation

HSP70/110s are a subgroup of heat shock proteins and play crucial roles in protein homeostasis. HSP70/110s can enhance cell survival in response to a multitude of stressful stimuli, of which the most studied one is heat stress. To perform a systematic study of HSP70/110s in sea cucumber Apostichopus japonicus, 15 HSP70/110 genes, including 13 HSP70s and two HSP110s, were identified and characterized from the transcriptome and genome of sea cucumber. Moderate expansion and conserved structure were found by the phylogenetic and syntenic analysis. Differential expression patterns of HSP70/110s were observed in adult individuals during aestivation, with the comparison of juvenile individuals without aestivation in chronic heat stress. Tissue-specific expression profiles were found both in adult and juvenile individuals, which might indicate that the functional tissues (intestine and respiratory tree) could be restored to normal physiological activity prior to protecting and sporting tissues (body wall and muscle). Differential expression profiles were also observed between the adult and juvenile individuals, which was mainly due to the hypometabolism in aestivation. Taken together, tissue-specific pattern and individual-specific pattern were observed in the HSP70/110 expression profiles in sea cucumber during aestivation. These findings could provide early insight into the involvement of HSP70/110s in the aestivation of marine invertebrate.

The effects of sewage treatment plant effluents on hepatic and intestinal biomarkers in common carp (Cyprinus carpio)

Sewage treatment plants (STPs) are one of the major source of pharmaceuticals and personal care products in the aquatic environment. Generally, the effects of individual chemicals on fish are studied under laboratory conditions, which leads to results that are potentially not realistic regarding the effects of these chemicals under environmental conditions. Therefore, in this study, common carps were held in exposed pond that receive water from STP effluents for 360 days under natural conditions. Elimination of xenobiotics starts in the fish intestine, in which the microbial community strongly influences its function. Moreover, the fish intestine functions as crucial organ for absorbing lipids and fatty acids (FA), with consequent transport to the liver where their metabolism occurs. The liver is the primary organ performing xenobiotic metabolism in fish, and therefore, the presence of pollutants may interact with the metabolism of FA. The catalytic activity of CYP1A and CYP3A-like enzymes, their gene expression, FA composition and intestinal microbiome consortia were measured. The catalytic activity of enzymes and their gene and protein expression, were induced in hepatic and intestinal tissues of fish from the exposed pond. Also, fish from the exposed pond had different compositions of FA than those from the control pond: concentration of 18:1 n-9 and 18:2 n-6 were significantly elevated and the longer chain n-3 FA 20:5 n-3, 22:5 n-3 and 22:6 n-3 were significantly lowered. There were clear differences among microbiome consortia in fish intestines across control and exposed groups. Microbiome taxa measured in exposed fish were also associated with those found in STP activated sludge. This study reveals that treated STP water, which is assumed to be clean, affected measured biomarkers in common carp.

A Review of Molecular Responses of Catfish to Bacterial Diseases and Abiotic Stresses

Catfish is one of the major aquaculture species in the United States. However, the catfish industry is threatened by several bacterial diseases such as enteric septicemia of catfish (ESC), columnaris disease and Aeromonas disease, as well as by abiotic stresses such as high temperature and low oxygen. Research has been conducted for several decades to understand the host responses to these diseases and abiotic stresses. With the development of sequencing technologies, and the application of genome-wide association studies in aquaculture species, significant progress has been made. This review article summarizes recent progress in understanding the molecular responses of catfish after bacterial infection and stress challenges, and in understanding of genomic and genetic basis for disease resistance and stress tolerance.

Applications of transcriptomics and proteomics in understanding fish immunity

With the development of intensive aquaculture, economic losses increasingly result from fish mortality due to pathogen infection. In recent years, a growing number of researchers have used transcriptomic and proteomic analyses to study fish immune responses to exogenous pathogen infection. Integrating transcriptomic and proteomic analyses provides a better understanding of the fish immune system including gene expression, regulation, and the intricate biological processes underlying immune responses against infection. This review focuses on the recent advances in the fields of transcriptomics and proteomics, which have contributed to our understanding of fish immunity to exogenous pathogens.

Identification of 74 cytochrome P450 genes and co-localized cytochrome P450 genes of the CYP2K, CYP5A, and CYP46A subfamilies in the mangrove killifish Kryptolebias marmoratus

BACKGROUND

The mangrove killifish Kryptolebias marmoratus is the only vertebrate that reproduces by self-fertilizing and is an important model species in genetics and marine ecotoxicology. Using whole-genome and transcriptome sequences, we identified all members of the cytochrome P450 (CYP) family in this model teleost and compared them with those of other teleosts.

RESULTS

A total of 74 cytochrome P450 genes and one pseudogene were identified in K. marmoratus. Phylogenetic analysis indicated that the CYP genes in clan 2 were most expanded, while synteny analysis with other species showed orthologous relationships of CYP subfamilies among teleosts. In addition to the CYP2K expansions, five tandem duplicated gene copies of CYP5A were observed. These features were unique to K. marmoratus.

CONCLUSIONS

These results shed a light on CYP gene evolution, particularly the co-localized CYP2K, CYP5A, and CYP46A subfamilies in fish. Future studies of CYP expression could identify specific endogenous and exogenous environmental factors that triggered the evolution of tandem CYP duplication in K. marmoratus.

Training the salmon's genes: influence of aerobic exercise, swimming performance and selection on gene expression in Atlantic salmon

BACKGROUND

Farmed and wild Atlantic salmon are exposed to many infectious and non-infectious challenges that can cause mortality when they enter the sea. Exercise before transfer promotes growth, health and survival in the sea. Swimming performance in juveniles at the freshwater parr stage is positively associated with resistance to some diseases. Genetic variation is likely to affect response to exercise. In this study we map genetic differences associated with aerobic exercise, swimming performance and genetic origin. Eggs from the selectively bred Bolaks salmon and wild Lærdal River salmon strains were reared until parr in a common environment. Swimming performance was assessed by subjecting the fish to either continuous hard exercise or control conditions for 18 days. Heart was sampled for examination of gene expression using RNA-seq (~60 fish/treatment).

RESULTS

Lower expression of genes affecting immune function was found in domesticated than wild parr. Among wild parr under control exercise the expression of a large number of genes involved in general metabolism, stress and immune response was lower in superior swimmers suggesting that minimisation of energy expenditure during periods of low activity makes parr better able to sustain bursts of swimming for predator avoidance. A similar set of genes were down-regulated with training among wild parr with inferior swimming performance. These parr react to training in a way that their cardiac expression patterns become like the superior performing wild parr under control exercise conditions. Diversifying selection caused by breeding of domesticated stock, and adaptive pressures in wild stock, has affected the expression and frequency of single nucleotide polymorphisms (SNPs) for multiple functional groups of genes affecting diverse processes. SNPs associated with swimming performance in wild parr map to genes involved in energetic processes, coding for contractile filaments in the muscle and controlling cell proliferation.

CONCLUSIONS

Domesticated parr have less phenotypic plasticity in response to training and lower expression of genes with functions affecting immune response. The genetic response to training is complex and depends on the background of parr and their swimming ability. Exercise should be tailored to the genetics and swimming performance of fish.

Omics in fish mucosal immunity

The mucosal immune system of fish is a complex network of immune cells and molecules that are constantly surveilling the environment and protecting the host from infection. A number of "omics" tools are now available and utilized to understand the complexity of mucosal immune systems in non-traditional animal models. This review summarizes recent advances in the implementation of "omics" tools pertaining to the four mucosa-associated lymphoid tissues in teleosts. Genomics, transcriptomics, proteomics, and "omics" in microbiome research require interdisciplinary collaboration and careful experimental design. The data-rich datasets generated are proving really useful at discovering new innate immune players in fish mucosal secretions, identifying novel markers of specific mucosal immune responses, unraveling the diversity of the B and T cell repertoires and characterizing the diversity of the microbial communities present in teleost mucosal surfaces. Bioinformatics, data analysis and storage platforms should be developed to facilitate rapid processing of large datasets, especially when mammalian tools such as bioinformatics analysis software are not available in fishes.

Estimation of random effects and identifying heterogeneous genes in meta-analysis of gene expression studies

Combining effect sizes from individual studies using random-effects meta-analysis models are commonly applied in high-dimensional gene expression data. However, unknown study heterogeneity can arise from inconsistencies in sample quality and experimental conditions. High heterogeneity of effect sizes can reduce statistical power of the models. In this study, we describe three hypothesis-testing frameworks for meta-analysis of microarray data, and review several existing meta-analytic techniques that have been used in the genomic setting. These include P-value-based methods, rank-based methods and effect-size-based methods. We then discuss limitations of some of these methods and describe random-effects-based methods in detail. We introduce two methods for estimating the inter-study variance in random-effects meta-analytic models and another method for identifying heterogeneous genes for gene expression data. We compared various methods with the standard and existing meta-analytic techniques in the genomic framework. We demonstrate our results through a series of simulations and application in Alzheimer's gene expression data.

Identification and expression of cytochrome P450 genes in the Chinese giant salamander Andrias davidianus

Cytochrome P450 (cyp) genes play vital roles in biological processes, including in metabolism of endogenous and exogenous compounds. Exogenous hormone influences on gene expression that leads to disruption of gonad development have been studied in several species, but whether exogenous hormones affect cyp genes that show sexually dimorphic expression remains to be determined. Here, we identified and characterized cyp genes from the widely-cultured Chinese giant salamander Andrias davidianus. We obtained 20 cyp genes including 11 genes with complete sequences. Phylogenetic analyses supported the classification of cyp genes similar to other vertebrates. Expression profile of female and male salamanders showed multiple cyp genes to exhibit higher expression in ovary than in testis, including cyp26a, cyp19a, cyp1a1, cyp4v2, cyp3a24, cyp2c20, cyp2d14, cyp2d15, and cyp4b; while cyp11a, cyp2b11, cyp11b1, cyp4f22, cyp2j6, cyp2k1, cyp2k4, cyp1a5 was higher in testis than in ovary. Seven sex-biased genes were detected after aromatase inhibitor injection and with exposure to high water temperatures. Cyp2k1, cyp11b1, and cyp2j6 expression were increased, while cyp26a, cyp2c20, cyp4b genes were decreased with aromatase inhibition. In ovary, cyp26a and cyp2c20 expression were significantly reduced; whereas cyp2k1, cyp2j6, and cyp2k4 expression were significantly elevated with no changes in cyp11a and cyp11b1 expression after temperature treatment. These findings provide valuable information for further study of sex differentiation mechanisms and cyp gene evolution.

Transcriptomic responses in the fish intestine

The intestine, being a multifunctional organ central to both nutrient uptake, pathogen recognition and regulating the intestinal microbiome, has been subjected to intense research. This review will focus on the recent studies carried out using high-throughput gene expression approaches, such as microarray and RNA sequencing (RNA-seq). These techniques have advanced greatly in recent years, mainly as a result of the massive changes in sequencing methodologies. At the time of writing, there is a transition between relatively well characterised microarray platforms and the developing RNA-seq, with the prediction that within a few years as costs decrease and computation power increase, RNA-seq related approaches will supersede the microarrays. Comparisons between the approaches are made and specific examples of how the techniques have been used to examine intestinal responses to pathogens, dietary manipulations and osmoregulatory challenges are given.

Characterization and sexual dimorphic expression of Cytochrome P450 genes in the hypothalamic-pituitary-gonad axis of yellow catfish

Yellow catfish (Pelteobagrus fulvidraco) is an important freshwater fish species in China. In particular, an all-male population has been commercially produced for the males grow faster than females. However, the molecular mechanisms underlying sexual dimorphism of body size and sex differentiation are still unclear in yellow catfish. This study attempts to characterize and analyze the expression of Cytochrome P450 (CYP) family members that have been shown to play an important role in sex differentiation and metabolism in teleosts. A total of 25 CYP genes were identified from our transcriptomes by 454 pyrosequencing and Solexa sequencing, including 17 genes with complete open reading frame (ORF). Phylogenetic analyses were conducted to compare these genes with their counterparts from other teleosts. In the tissues of hypothalamic-pituitary-gonad (HPG) axis, most of the genes were expressed at uniform level in both sexes. However, multiple CYP genes displayed sexual dimorphic expression, such as cyp2AD, cyp4b, cyp8a, cyp11b2, cyp17a and cyp27a expressed at higher level in testis than in ovary, whereas cyp2g, cyp7a, cyp8b, cyp19a1a and cyp26a expressed at higher level in ovary than in testis. The expression response of six CYP genes in ovary was also assessed after 17α-methyltestosterone (MT) treatment. Testis-biased expressed cyp11b2 and cyp17a were significantly up-regulated, while cyp11a and cyp19a1a were reduced in ovary after MT treatment. Our work is helpful for understanding molecular evolution of CYP genes in vertebrates and the mechanism of sexual dimorphism in teleosts.