ABSTRACT
Target enrichment sequencing techniques are gaining widespread use in the field of genomics, prized for their economic efficiency and swift processing times. However, their success depends on the performance of probes and the evenness of sequencing depth among each probe. To accurately predict probe coverage depth, a model called Deqformer is proposed in this study. Deqformer utilizes the oligonucleotides sequence of each probe, drawing inspiration from Watson-Crick base pairing and incorporating two BERT encoders to capture the underlying information from the forward and reverse probe strands, respectively. The encoded data are combined with a feed-forward network to make precise predictions of sequencing depth. The performance of Deqformer is evaluated on four different datasets: SNP panel with 38 200 probes, lncRNA panel with 2000 probes, synthetic panel with 5899 probes and HD-Marker panel for Yesso scallop with 11 000 probes. The SNP and synthetic panels achieve impressive factor 3 of accuracy (F3acc) of 96.24% and 99.66% in 5-fold cross-validation. F3acc rates of over 87.33% and 72.56% are obtained when training on the SNP panel and evaluating performance on the lncRNA and HD-Marker datasets, respectively. Our analysis reveals that Deqformer effectively captures hybridization patterns, making it robust for accurate predictions in various scenarios. Deqformer leads to a novel perspective for probe design pipeline, aiming to enhance efficiency and effectiveness in probe design tasks.
Subject(s)
Deep Learning , RNA, Long Noncoding , DNA Probes/genetics , Nucleic Acid Hybridization , GenomicsABSTRACT
Recently, elevated seawater temperatures have resulted numerous adverse effects, including significant mortality among bivalves. The dwarf surf clam, Mulinia lateralis, is considered a valuable model species for bivalve research due to its rapid growth and short generation time. The successful cultivation in laboratory setting throughout its entire life cycle makes it an ideal candidate for exploring the potential mechanisms underlying bivalve responses to thermal stress. In this study, a total of 600 clams were subjected to a 17-day thermal stress experiment at a temperature of 30 °C which is the semi-lethal temperature for this species. Ninety individuals who perished initially were classified as heat-sensitive populations (HSP), while 89 individuals who survived the experiment were classified as heat-tolerant populations (HTP). Subsequently, 179 individuals were then sequenced, and 21,292 single nucleotide polymorphisms (SNPs) were genotyped for downstream analysis. The heritability estimate for survival status was found to be 0.375 ± 0.127 suggesting a genetic basis for thermal tolerance trait. Furthermore, a genome-wide association study (GWAS) identified three SNPs and 10 candidate genes associated with thermal tolerance trait in M. lateralis. These candidate genes were involved in the ETHR/EHF signaling pathway and played pivotal role in signal sensory, cell adhesion, oxidative stress, DNA damage repair, etc. Additionally, qPCR results indicated that, excluding MGAT4A, ZAN, and RFC1 genes, all others exhibited significantly higher expression in the HTP (p < 0.05), underscoring the critical involvement of the ETHR/EHF signaling pathway in M. lateralis' thermal tolerance. These results unveil the presence of standing genetic variations associated with thermal tolerance in M. lateralis, highlighting the regulatory role of the ETHR/EHF signaling pathway in the bivalve's response to thermal stress, which contribute to comprehension of the genetic basis of thermal tolerance in bivalves.
ABSTRACT
It is of vital importance to understand the population structure, dissect the genetic bases of performance traits, and make proper strategies for selection in breeding programs. However, there is no single webserver covering the specific needs in aquaculture. We present Aquaculture Molecular Breeding Platform (AMBP), the first web server for genetic data analysis in aquatic species of farming interest. AMBP integrates the haplotype reference panels of 18 aquaculture species, which greatly improves the accuracy of genotype imputation. It also supports multiple tools to infer genetic structures, dissect the genetic architecture of performance traits, estimate breeding values, and predict optimum contribution. All the tools are coherently linked in a web-interface for users to generate interpretable results and evaluate statistical appropriateness. The webserver supports standard VCF and PLINK (PED, MAP) files, and implements automated pipelines for format transformation and visualization to simplify the process of analysis. As a demonstration, we applied the webserver to Pacific white shrimp and Atlantic salmon datasets. In summary, AMBP constitutes comprehensive resources and analytical tools for exploring genetic data and guiding practical breeding programs. AMBP is available at http://mgb.qnlm.ac.
Subject(s)
Aquaculture , DNA Shuffling , Polymorphism, Single Nucleotide , Software , Aquaculture/methods , Genotype , Phenotype , Animals , Breeding , InternetABSTRACT
A specific strain of Vibrio parahaemolyticus (VpAHPND) causes acute hepatopancreatic necrosis disease (AHPND), leading to significant losses in shrimp aquaculture. Outer membrane vesicles (OMVs) are naturally secreted by Gram-negative bacteria, and their significant roles in host-pathogen interactions and pathogenicity have been recognized. In the present study, OMVs were isolated from VpAHPND by differential-ultracentrifugation and used for proteomics analysis. In the Nano-HPLC-MS/MS analysis, totally 645 proteins were determined, including virulence factors, immunogenic proteins, outer membrane protein, bacterial secretory proteins, ribosomal proteins, protease, and iron regulation proteins. Furthermore, GO and KEGG annotations indicated that proteins identified in VpAHPND-OMVs are involved in metabolism, regulation of multiple biological processes, genetic information processes, immunity and more. Meanwhile, toxin proteins PirAvp and PirBvp, associated with VpAHPND pathogenicity, were also identified in the proteome of VpAHPND-OMVs. Our objective is to identify the protein composition of OMVs released by VpAHPND, analyzing the potential for cytotoxicity and immunomodulatory activity of these granule hosts. This study is crucial for understanding the roles played by bacterial-derived vesicles in the disease process, given that these vesicles carry relevant activities inherent to the bacteria that produce them.
Subject(s)
Penaeidae , Proteome , Vibrio parahaemolyticus , Vibrio parahaemolyticus/pathogenicity , Proteome/analysis , Animals , Penaeidae/microbiology , Hepatopancreas/microbiology , Hepatopancreas/pathology , Bacterial Outer Membrane Proteins/metabolism , Proteomics , Vibrio Infections/veterinary , Vibrio Infections/microbiology , Extracellular Vesicles/metabolismABSTRACT
Nitrite exposure has become a significant concern in the aquaculture industry, posing a severe threat to aquatic animals such as shrimp. While studies have reported the adverse effects of nitrite on shrimp growth, the part played by the gut microbiota in shrimp mortality resulting from nitrite exposure is poorly understood. Here, the effects of nitrite on shrimp gut bacterial community were investigated using 16S rRNA amplicon sequencing, bacterial isolation, genomic analysis, and infection experiments. Compared to the control_healthy group, changes in the bacterial composition of the nitrite_dead group were associated with reduced abundance of specific beneficial bacteria and increased abundance of certain pathogenic bacteria. Notably, members of the Photobacterium genus were found to be significantly enriched in the nitrite_dead group. Genomic analysis of a representative Photobacterium strain (LvS-8n3) revealed a variety of genes encoding bacterial toxins, including hemolysin, adhesin, and phospholipase. Furthermore, it was also found that LvS-8n3 exhibits strong pathogenicity, probably due to its high production of pathogenic factors and the ability to utilize nitrite for proliferation. Therefore, the proliferation of pathogenic Photobacterium species appears pivotal for driving shrimp mortality caused by nitrite exposure. These findings provide novel insights into the disease mechanism in shrimp under conditions of environmental change.
ABSTRACT
BACKGROUND: Mytilidae, also known as marine mussels, are widely distributed in the oceans worldwide. Members of Mytilidae show a tremendous range of ecological adaptions, from the species distributed in freshwater to those that inhabit in deep-sea. Mitochondria play an important role in energy metabolism, which might contribute to the adaptation of Mytilidae to different environments. In addition, some bivalve species are thought to lack the mitochondrial protein-coding gene ATP synthase F0 subunit 8. Increasing studies indicated that the absence of atp8 may be caused by annotation difficulties for atp8 gene is characterized by highly divergent, variable length. RESULTS: In this study, the complete mitochondrial genomes of three marine mussels (Xenostrobus securis, Bathymodiolus puteoserpentis, Gigantidas vrijenhoeki) were newly assembled, with the lengths of 14,972 bp, 20,482, and 17,786 bp, respectively. We annotated atp8 in the sequences that we assembled and the sequences lacking atp8. The newly annotated atp8 sequences all have one predicted transmembrane domain, a similar hydropathy profile, as well as the C-terminal region with positively charged amino acids. Furthermore, we reconstructed the phylogenetic trees and performed positive selection analysis. The results showed that the deep-sea bathymodiolines experienced more relaxed evolutionary constraints. And signatures of positive selection were detected in nad4 of Limnoperna fortunei, which may contribute to the survival and/or thriving of this species in freshwater. CONCLUSIONS: Our analysis supported that atp8 may not be missing in the Mytilidae. And our results provided evidence that the mitochondrial genes may contribute to the adaptation of Mytilidae to different environments.
Subject(s)
Genome, Mitochondrial , Mytilidae , Animals , Mytilidae/genetics , Phylogeny , Genes, Mitochondrial , Mitochondrial Proton-Translocating ATPases/genetics , Genomics/methodsABSTRACT
Blood clams differ from their molluscan kins by exhibiting a unique red-blood (RB) phenotype; however, the genetic basis and biochemical machinery subserving this evolutionary innovation remain unclear. As a fundamental step toward resolving this mystery, we presented the first chromosome-level genome and comprehensive transcriptomes of the blood clam Tegillarca granosa for an integrated genomic, evolutionary, and functional analyses of clam RB phenotype. We identified blood clam-specific and expanded gene families, as well as gene pathways that are of RB relevant. Clam-specific RB-related hemoglobins (Hbs) showed close phylogenetic relationships with myoglobins (Mbs) of blood clam and other molluscs without the RB phenotype, indicating that clam-specific Hbs were likely evolutionarily derived from the Mb lineage. Strikingly, similar to vertebrate Hbs, blood clam Hbs were present in a form of gene cluster. Despite the convergent evolution of Hb clusters in blood clam and vertebrates, their Hb clusters may have originated from a single ancestral Mb-like gene as evidenced by gene phylogeny and synteny analysis. A full suite of enzyme-encoding genes for heme synthesis was identified in blood clam, with prominent expression in hemolymph and resembling those in vertebrates, suggesting a convergence of both RB-related Hb and heme functions in vertebrates and blood clam. RNA interference experiments confirmed the functional roles of Hbs and key enzyme of heme synthesis in the maintenance of clam RB phenotype. The high-quality genome assembly and comprehensive transcriptomes presented herein serve new genomic resources for the super-diverse phylum Mollusca, and provide deep insights into the origin and evolution of invertebrate RB.
Subject(s)
Arcidae/genetics , Biological Evolution , Hemoglobins/genetics , Animals , Arcidae/metabolism , Chromosomes , Genome , Heme/biosynthesis , Hemolymph/metabolism , Humans , Multigene Family , TranscriptomeABSTRACT
Targeted genotyping of transcriptome-scale genetic markers is highly attractive for genetic, ecological, and evolutionary studies, but achieving this goal in a cost-effective manner remains a major challenge, especially for laboratories working on nonmodel organisms. Here, we develop a high-throughput, sequencing-based GoldenGate approach (called HD-Marker), which addresses the array-related issues of original GoldenGate methodology and allows for highly multiplexed and flexible targeted genotyping of more than 12,000 loci in a single-tube assay (in contrast to fewer than 3100 in the original GoldenGate assay). We perform extensive analyses to demonstrate the power and performance of HD-Marker on various multiplex levels (296, 795, 1293, and 12,472 genic SNPs) across two sequencing platforms in two nonmodel species (the scallops Chlamys farreri and Patinopecten yessoensis), with extremely high capture rate (98%-99%) and genotyping accuracy (97%-99%). We also demonstrate the potential of HD-Marker for high-throughput targeted genotyping of alternative marker types (e.g., microsatellites and indels). With its remarkable cost-effectiveness (as low as $0.002 per genotype) and high flexibility in choice of multiplex levels and marker types, HD-Marker provides a highly attractive tool over array-based platforms for fulfilling genome/transcriptome-wide targeted genotyping applications, especially in nonmodel organisms.
Subject(s)
Genetic Markers , Genotyping Techniques , High-Throughput Nucleotide Sequencing , Alleles , Genotype , INDEL Mutation , Microsatellite Repeats , Polymerase Chain Reaction , Polymorphism, Single Nucleotide , Reproducibility of Results , Sensitivity and SpecificityABSTRACT
Barbels are important sensory organs in teleosts, reptiles, and amphibians. The majority of â¼4,000 catfish species, such as the channel catfish (Ictalurus punctatus), possess abundant whisker-like barbels. However, barbel-less catfish, such as the bottlenose catfish (Ageneiosus marmoratus), do exist. Barbeled catfish and barbel-less catfish are ideal natural models for determination of the genomic basis for barbel development. In this work, we generated and annotated the genome sequences of the bottlenose catfish, conducted comparative and subtractive analyses using genome and transcriptome datasets, and identified differentially expressed genes during barbel regeneration. Here, we report that chemokine C-C motif ligand 33 (ccl33), as a key regulator of barbel development and regeneration. It is present in barbeled fish but absent in barbel-less fish. The ccl33 genes are differentially expressed during barbel regeneration in a timing concordant with the timing of barbel regeneration. Knockout of ccl33 genes in the zebrafish (Danio rerio) resulted in various phenotypes, including complete loss of barbels, reduced barbel sizes, and curly barbels, suggesting that ccl33 is a key regulator of barbel development. Expression analysis indicated that paralogs of the ccl33 gene have both shared and specific expression patterns, most notably expressed highly in various parts of the head, such as the eye, brain, and mouth areas, supporting its role for barbel development.
Subject(s)
Chemokines/metabolism , Fish Proteins/metabolism , Sense Organs/growth & development , Animals , Catfishes/genetics , Catfishes/growth & development , Catfishes/metabolism , Chemokines/genetics , Chemokines/physiology , Fish Proteins/genetics , Fish Proteins/physiology , Gene Expression Profiling , Genome/genetics , Male , Sense Organs/metabolism , Zebrafish/genetics , Zebrafish/growth & development , Zebrafish/metabolismABSTRACT
The oriental river prawn Macrobrachium nipponense is a commercially important freshwater shrimp that is widely farmed in China. Aeromonas veronii is a conditional pathogen of farmed shrimp, which has caused huge economic losses to the industry. Therefore, there is urgency to study the host-pathogen interactions between M. nipponense and A. veronii to screen individuals with antimicrobial resistance. In this study, we examined the hepatopancreas of moribund M. nipponense infected with A. veronii and healthy individuals at both the histopathological and transcriptomic levels. We showed that A. veronii infection resulted in tubular necrosis of the M. nipponense hepatopancreas. Such changes likely affect assimilation, storage, and excretion by the hepatopancreas, which could ultimately affect the survival and growth of infected individuals. Among the 61,345 unigenes obtained through RNA sequencing and de novo transcriptome assembly, 232 were differentially expressed between the two groups. KEGG and GO analyses revealed that these differentially expressed genes were implicated in pathways, including PPAR, PI3K/AKT, and AMPK signaling. The results of this study will contribute to an analysis of the immune response of M. nipponense to A. veronii infection at the transcriptomic level. Furthermore, the RNA-seq data generated here provide an important genomic resource for research on M. nipponense in the absence of a reference genome.
Subject(s)
Aeromonas veronii/physiology , Hepatopancreas/immunology , Palaemonidae/microbiology , Seafood/microbiology , Transcriptome/immunology , Animals , China , Hepatopancreas/pathology , Host-Pathogen Interactions , Necrosis , Palaemonidae/immunology , Signal TransductionABSTRACT
Spiroplasma eriocheiris (S. eriocheiris) infection causes a significant economic loss in Penaeus vannamei (P. vannamei) culture industry. However, the response of P. vannamei hemocytes to S. eriocheiris infection has not been extensively studied. In this study, we conducted full-length transcriptome and long non-coding RNA (lncRNA) analyses of P. vannamei hemocytes by a challenge test with S. eriocheiris. Following assembly and annotation, there were 8077 high-quality unigenes. A total of 1168 differentially expressed genes (DEGs) were obtained, including 792 up-regulated and 376 down-regulated genes by differential expression analysis. Gene ontology (GO) enrichment analysis showed that the up-regulated DEGs were mainly clustered into immune system process, defense response, cell cycle and organelle organization. On the other hand, the down-regulated DEGs included that genes that were mainly clustered into metabolic processes related to organic compounds, metabolic process and cellular metabolic process. Protein-protein interaction (PPI) network analysis of DEGs indicated that the pivotal gene interactions were connected to stress response, immune system process and cell cycle. The lncRNA analysis identified multiple lncRNAs, which were highly co-expressed with the immune-related genes, such as lncRNA transcript-12631 and transcript-12631, suggesting that lncRNAs may be involved in the regulation of immune defense in shrimp hemocytes. Additionally, 20 hub unigenes and putative lncRNAs related to immune system were validated by quantitative real-time PCR (qRT-PCR), validating the reliability of RNA-Seq. This study revealed a close connection between the immune and metabolic systems of S. eriocheiris infected P. vannamei.
Subject(s)
Arthropod Proteins/genetics , Arthropod Proteins/immunology , Gram-Negative Bacterial Infections/immunology , Penaeidae/genetics , Penaeidae/immunology , RNA, Long Noncoding/immunology , Spiroplasma , Animals , Gram-Negative Bacterial Infections/veterinary , Hemocytes/immunology , Penaeidae/microbiology , TranscriptomeABSTRACT
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.
Subject(s)
Ictaluridae/genetics , Sex Determination Processes/genetics , Y Chromosome , Animals , Chromosome Mapping , Female , Genetic Linkage , Genome , Male , Sequence Analysis, DNAABSTRACT
The oriental river prawn, Macrobrachium nipponense, is an economically and nutritionally important species of the Palaemonidae family of decapod crustaceans. Ammonia is a major aquatic environmental pollutant that negatively affects the health of prawns and their associated commercial productivity. Here, we used high-throughput sequencing techniques for detecting the effects of ammonia stress (22.1â¯mg/L ammonia-N for 48â¯h) on gene expression in the hepatopancreas of M. nipponense. We generated 176,228,782 high-quality reads after eliminating adapter sequences and filtering out low-quality reads, which were assembled into 63453 unigenes. Comparative analysis of the expression profiles of the ammonia-treated and control groups identified 887 differentially expressed genes (Pâ¯<â¯0.05), including 481 upregulated genes and 406 downregulated genes. Analyses of the GO and KEGG databases revealed significant differences between the two groups in 32 pathways. Immune-related pathways under ammonia stress included Complement and coagulation cascades, Platelet activation, B cell receptor signaling pathway, Antigen processing and presentation, Chemokine signaling pathway, NOD-like receptor signaling pathway, RIG-I-like receptor signaling pathway, T cell receptor signaling pathway and Toll-like receptor signaling pathway. Remarkably, ammonia stress altered the expression patterns of key immune genes (lectin3, syntenin, alpha-2-macroglobulin, cathepsin L, PIM3, serine protease inhibitor, suppressor of cytokine signaling-2 like protein), indicating that ammonia-stress induce immune response. These data provide new insights into the immune response of M. nipponense and pave a new way for fighting ammonia stress. The genes and pathways identified here represent valuable genetic resources for development of molecular markers and genetic breeding studies.
Subject(s)
Ammonia/adverse effects , Arthropod Proteins/genetics , Gene Expression Regulation/immunology , Immunity, Innate/genetics , Palaemonidae/drug effects , Transcriptome/drug effects , Water Pollutants, Chemical/adverse effects , Animals , Arthropod Proteins/metabolism , Gene Expression Profiling , Hepatopancreas/drug effects , Hepatopancreas/metabolism , Palaemonidae/genetics , Palaemonidae/metabolism , Transcriptome/geneticsABSTRACT
The oriental river prawn Macrobrachium nipponense is a highly adaptable, tolerant, and fecund freshwater prawn that inhabits a wide range of aquatic environments. The hepatopancreas of crustaceans is not only a site for secretion of digestive enzymes, and also plays important roles in several metabolic processes, such as lipid and carbohydrate metabolism. It is the main organ for the detoxification and immunity. In this study, high-throughput sequencing techniques were used to detect the effect of nitrite stress (10â¯mg/L nitrite-N for 48â¯h) on gene expression in the hepatopancreas of M. nipponense. A total of 13,769 million reads were harvested, and 94,534 transcripts were de novo assembled using Trinity software and produced 56,054 non-redundant transcripts. A total of 825 differentially expressed genes were obtained comparing 48â¯h nitrite stress with control group. In the analysis of GO and KEGG database, significant differences were found in 49 pathways. Immune-related pathways under nitrite stress included arginine and proline metabolism, glutamate metabolism, Jak-Stat signaling pathway, endocytosis, wnt signaling pathway, RIG-I-like receptor signaling pathway, TGF-beta signaling pathway, GnRH signaling pathway and phagosome. Apoptosis-related pathway was also significantly altered, such as lysosome and apoptosis. Remarkably, nitrite stress altered the expression patterns of key apoptosis genes (tetraspanins-like protein, LAMP, CD63, caspase 3C and Caspase 1) and immune genes (Serine proteinase-like protein, C-type lectin, daf-36, SOCS-2, alpha-2-macroglobulin), confirmed that nitrite-stress induce immune response and eventually even apoptosis. This study provided a new insight into the role of hepatopancreas in crustaceans, and further investigation will continue.
Subject(s)
Hepatopancreas/drug effects , Nitrites/toxicity , Palaemonidae/drug effects , Water Pollutants, Chemical/toxicity , Animals , Apoptosis/drug effects , Arthropod Proteins/metabolism , Gene Expression Profiling , Hepatopancreas/metabolism , Immunity, Innate/drug effects , Palaemonidae/genetics , Palaemonidae/metabolism , Signal Transduction , Stress, PhysiologicalABSTRACT
The swimbladder is an internal gas-filled organ in teleosts. Its major function is to regulate buoyancy. The swimbladder exhibits great variation in size, shape, and number of compartments or chambers among teleosts. However, genomic control of swimbladder variation is unknown. Channel catfish ( Ictalurus punctatus), blue catfish ( Ictalurus furcatus), and their F1 hybrids of female channel catfish × male blue catfish (C × B hybrid catfish) provide a good model in which to investigate the swimbladder morphology, because channel catfish possess a single-chambered swimbladder, whereas blue catfish possess a bichambered swimbladder; C × B hybrid catfish possess a bichambered swimbladder but with a significantly reduced posterior chamber. Here we determined the transcriptional profiles of swimbladder from channel catfish, blue catfish, and C × B hybrid catfish. We examined their transcriptomes at both the fingerling and adult stages. Through comparative transcriptome analysis, ~4,000 differentially expressed genes (DEGs) were identified. Among these DEGs, members of the Wnt signaling pathway ( wnt1, wnt2, nfatc1, rac2), Hedgehog signaling pathway ( shh), and growth factors ( fgf10, igf-1) were identified. As these genes were known to be important for branching morphogenesis of mammalian lung and of mammary glands, their association with budding of the posterior chamber primordium and progressive development of bichambered swimbladder in fish suggest that these branching morphogenesis-related genes and their functions in branching are evolutionarily conserved across a broad spectrum of species.
Subject(s)
Air Sacs/metabolism , Catfishes/genetics , Transcriptome/genetics , Animals , Female , Gene Expression Profiling , Male , Morphogenesis/genetics , Morphogenesis/physiologyABSTRACT
Channel catfish is the leading aquaculture species in the US, and one of the reasons for its application in aquaculture is its relatively high tolerance against hypoxia. However, hypoxia can still cause huge economic losses to the catfish industry. Studies on hypoxia tolerance, therefore, are important for aquaculture. Fish swimbladder has been considered as an accessory respiration organ surrounded by a dense capillary countercurrent exchange system. In this regard, we conducted RNA-Seq analysis with swimbladder samples of catfish under hypoxic and normal conditions to determine if swimbladder was responsive to low oxygen treatment and to reveal genes, their expression patterns, and pathways involved in hypoxia responses in catfish. A total of 155 differentially expressed genes (DEGs) were identified from swimbladder of adult catfish, whereas a total of 2,127 DEGs were identified from swimbladder of fingerling catfish under hypoxic condition as compared with untreated controls. Subsequent pathway analysis revealed that many DEGs under hypoxia were involved in HIF signaling pathway ( nos2, eno2, camk2d2, prkcb, cdkn1a, eno1, and tfrc), MAPK signaling pathway (voltage-dependent calcium channel subunit genes), PI3K/Akt/mTOR signaling pathway ( itga6, g6pc, and cdkn1a), Ras signaling pathway ( efna3 and ksr2), and signaling by VEGF ( fn1, wasf3, and hspb1) in catfish swimbladder. This study provided insights into regulation of gene expression and their involved gene pathways in catfish swimbladder in response to low oxygen stresses.
Subject(s)
Air Sacs/metabolism , Gene Expression Profiling/methods , Ictaluridae/genetics , Oxygen/metabolism , Transcriptome , Animals , Fish Proteins/genetics , Hypoxia , Signal Transduction/genetics , Stress, PhysiologicalABSTRACT
BACKGROUND: Walking catfish (Clarias batrachus) is a freshwater fish capable of air-breathing and locomotion on land. It usually inhabits various low-oxygen habitats, burrows inside the mudflat, and sometimes "walks" to search for suitable environments during summer. It has evolved accessory air-breathing organs for respiring air and corresponding mechanisms to survive in such challenging environments. Thereby, it serves as a great model for understanding adaptations to terrestrial life. RESULTS: Comparative genomics with channel catfish (Ictalurus punctatus) revealed specific adaptations of C. batrachus in DNA repair, enzyme activator activity, and small GTPase regulator activity. Comparative analysis with 11 non-air-breathing fish species suggested adaptive evolution in gene expression and nitrogenous waste metabolic processes. Further, myoglobin, olfactory receptor related to class A G protein-coupled receptor 1, and sulfotransferase 6b1 genes were found to be expanded in the air-breathing walking catfish genome, with 15, 15, and 12 copies, respectively, compared to non-air-breathing fishes that possess only 1-2 copies of these genes. Additionally, we sequenced and compared the transcriptomes of the gill and the air-breathing organ to characterize the mechanism of aerial respiration involved in elastic fiber formation, oxygen binding and transport, angiogenesis, ion homeostasis and acid-base balance. The hemoglobin genes were expressed dramatically higher in the air-breathing organ than in the gill of walking catfish. CONCLUSIONS: This study provides an important genomic resource for understanding the adaptive mechanisms of walking catfish to terrestrial environments. It is possible that the coupling of enhanced abilities for oxygen storage and oxygen transport through genomic expansion of myoglobin genes and transcriptomic up-regulation of hemoglobin and angiogenesis-related genes are important components of the molecular basis for adaptation of this aquatic species to terrestrial life.
Subject(s)
Catfishes/genetics , Gene Expression Profiling/veterinary , Genome , Gills/metabolism , Sequence Analysis, DNA/veterinary , Adaptation, Physiological , Animals , Catfishes/physiology , Gills/physiology , Respiration , TranscriptomeABSTRACT
Rhamnose-binding lectins (RBLs) are crucial elements associated with innate immune responses to infections and have been characterized from a variety of teleost fishes. Given the importance of RBL in teleost fishes, we sought to study the diversity and expression profiles of RBLs in an important cultured fish, Nile tilapia (Oreochromis niloticus) following experimental infection with Streptococcus agalactiae, a major cause of streptococcosis in farmed tilapia. In this study, four predicted RBL genes were identified from Nile tilapia and were designated as OnRBL3a, OnRBL3b, OnRBL3c, and OnRBL3d. These OnRBLs were composed of two tandem-repeated type five carbohydrate recognition domains (CRDs), classified as type IIIc, and all clustered together phylogenetically. OnRBL-CRDs shared conserved topology of eight cysteine residues, characteristic peptide motifs of -YGR- and -DPC- (or -FGR- and -DTC-), and similar exon/intron organization. OnRBLs had the highest expression in immune-related tissues, gills, intestine or liver. However, the changes of OnRBL expression in the gills and intestine at 2 h, 4 h and 24 h post S. agalactiae challenge were modest, suggesting that tilapia may not mediate the entry or confront the infection of S. agalactiae through induction of RBL genes. The observed expression pattern may be related to the RBL type and CRD composition, S. agalactiae pathogenesis, the accessibility of ligands on the bacterial surface, and/or the species of fish. OnRBLs characterized in this study were the first RBL members identified in Nile tilapia and their characterization will expand our knowledge of RBLs in immunity.
Subject(s)
Cichlids/genetics , Cichlids/immunology , Fish Diseases/immunology , Gene Expression Regulation/immunology , Immunity, Innate/genetics , Lectins/genetics , Lectins/immunology , Amino Acid Sequence , Animals , Fish Proteins/chemistry , Fish Proteins/genetics , Fish Proteins/immunology , Gene Expression Profiling , Mucous Membrane , Phylogeny , Rhamnose , Sequence Alignment/veterinaryABSTRACT
Chemokines are a superfamily of structurally related chemotactic cytokines exerting significant roles in regulating cell migration and activation. They are defined by the presence of four conserved cysteine residues and are divided into four subfamilies depending on the arrangement of the first two conserved cysteines residues: CXC, CC, C and CX3C. In this study, a complete set of 17 CXC chemokine ligand (CXCL) genes was systematically identified and characterized from channel catfish genome through data mining of existing genomic resources. Phylogenetic analysis allowed annotation of the 17 CXC chemokines. Extensive comparative genomic analyses supported their annotations and orthologies, revealing the existence of fish-specific CXC chemokines and the expansion of CXC chemokines in the teleost genomes. The analysis of gene expression after bacterial infection indicated the CXC chemokines were expressed in a gene-specific manner. CXCL11.3 and CXCL20.3 were expressed significantly higher in resistant fish than in susceptible fish after ESC infection, while CXCL20.2 were expressed significantly higher in resistant fish than in susceptible fish after columnaris infection. The expression of those CXC chemokines, therefore can be a useful indicator of disease resistance. A similar pattern of expression was observed between resistant and susceptible fish with biotic and abiotic stresses, ESC, columnaris and hypoxia, suggesting that high levels of expression of the majority of CXC chemokines, with exception of CXC11 and CXC20, are detrimental to the host.