Reference genomes of channel catfish and blue catfish reveal multiple pericentric chromosome inversions.

BACKGROUND: Channel catfish and blue catfish are the most important aquacultured species in the USA. The species do not readily intermate naturally but F1 hybrids can be produced through artificial spawning. F1 hybrids produced by mating channel catfish female with blue catfish male exhibit heterosis and provide an ideal system to study reproductive isolation and hybrid vigor. The purpose of the study was to generate high-quality chromosome level reference genome sequences and to determine their genomic similarities and differences. RESULTS: We present high-quality reference genome sequences for both channel catfish and blue catfish, containing only 67 and 139 total gaps, respectively. We also report three pericentric chromosome inversions between the two genomes, as evidenced by long reads across the inversion junctions from distinct individuals, genetic linkage mapping, and PCR amplicons across the inversion junctions. Recombination rates within the inversional segments, detected as double crossovers, are extremely low among backcross progenies (progenies of channel catfish female × F1 hybrid male), suggesting that the pericentric inversions interrupt postzygotic recombination or survival of recombinants. Identification of channel catfish- and blue catfish-specific genes, along with expansions of immunoglobulin genes and centromeric Xba elements, provides insights into genomic hallmarks of these species. CONCLUSIONS: We generated high-quality reference genome sequences for both blue catfish and channel catfish and identified major chromosomal inversions on chromosomes 6, 11, and 24. These perimetric inversions were validated by additional sequencing analysis, genetic linkage mapping, and PCR analysis across the inversion junctions. The reference genome sequences, as well as the contrasted chromosomal architecture should provide guidance for the interspecific breeding programs.

The Effects of Larval Cryopreservation on the Epigenetics of the Pacific Oyster Crassostrea gigas.

High mortalities and highly variable results during the subsequent development of post-thaw larvae have been widely considered as key issues restricting the application of cryopreservation techniques to support genetic improvement programs and hatchery production in farmed marine bivalve species. To date, few studies have been undertaken to investigate the effects of cryodamage at the molecular level in bivalves. This study is the first to evaluate the effect of larval cryopreservation on the epigenetics of the resultant progenies of the Pacific oyster Crassostrea gigas. The results show that the level of DNA methylation was significantly (p < 0.05) higher and lower than that of the control when the trochophore larvae were revived and when they developed to D-stage larvae (day 1 post-fertilization), respectively, but the level returned to the control level from day 8 post-fertilization onwards. The expression of the epigenetic regulator genes DNMT3b, MeCP2, JmjCA, KDM2 and OSA changed significantly (p < 0.05) when the trochophore larvae were thawed, and then they reverted to the control levels at the D- and later larval developmental stages. However, the expression of other epigenetic regulator genes, namely, MBD2, DNMT1, CXXC1 and JmjD6, did not change at any post-thaw larval developmental stage. For the newly thawed trochophore larvae, the amount of methylated H3K4Me1 and H3K27Me1 significantly changed, and the expression of all Jumonji orthologs, except that of Jumonji5, significantly (p < 0.05) decreased. These epigenetic results agree with the data collected on larval performances (e.g., survival rate), suggesting that the effect period of the published cryopreservation technique on post-thaw larvae is short in C. gigas.

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.

Chemokine C-C motif ligand 33 is a key regulator of teleost fish barbel development.

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.

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

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

Comparative transcriptome analysis reveals conserved branching morphogenesis related genes involved in chamber formation of catfish swimbladder.

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.

Comparative transcriptome analysis of the swimbladder reveals expression signatures in response to low oxygen stress in channel catfish, Ictalurus punctatus.

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.

Comparative genome analysis of 52 fish species suggests differential associations of repetitive elements with their living aquatic environments.

BACKGROUND: Repetitive elements make up significant proportions of genomes. However, their roles in evolution remain largely unknown. To provide insights into the roles of repetitive elements in fish genomes, we conducted a comparative analysis of repetitive elements of 52 fish species in 22 orders in relation to their living aquatic environments. RESULTS: The proportions of repetitive elements in various genomes were found to be positively correlated with genome sizes, with a few exceptions. More importantly, there appeared to be specific enrichment between some repetitive element categories with species habitat. Specifically, class II transposons appear to be more abundant in freshwater bony fish than in marine bony fish when phylogenetic relationship is not considered. In contrast, marine bony fish harbor more tandem repeats than freshwater species. In addition, class I transposons appear to be more abundant in primitive species such as cartilaginous fish and lamprey than in bony fish. CONCLUSIONS: The enriched association of specific categories of repetitive elements with fish habitats suggests the importance of repetitive elements in genome evolution and their potential roles in fish adaptation to their living environments. However, due to the restriction of the limited sequenced species, further analysis needs to be done to alleviate the phylogenetic biases.

Genome sequence of walking catfish (Clarias batrachus) provides insights into terrestrial adaptation.

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.

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

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

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

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

A deletion in the Hermansky-Pudlak syndrome 4 (Hps4) gene appears to be responsible for albinism in channel catfish.

Albinism is caused by a series of genetic abnormalities leading to reduction of melanin production. Albinism is quite frequent in catfish, but the causative gene and the molecular basis were unknown. In this study, we conducted a genome-wide association study (GWAS) using the 250 K SNP array. The GWAS analysis allowed mapping of the albino phenotype in the Hermansky-Pudlak syndrome 4 (Hps4) gene, which is known to be involved in melanosome biosynthesis. Sequencing analysis revealed that a 99-bp deletion was present in all analyzed albino catfish at the intron 2 and exon 3 junction. This deletion led to the skipping of the entire exon 3 which was confirmed by RT-PCR. Therefore, Hps4 was determined to be the candidate gene of the catfish albinism.

Multiple across-strain and within-strain QTLs suggest highly complex genetic architecture for hypoxia tolerance in channel catfish.

The ability to survive hypoxic conditions is important for various organisms, especially for aquatic animals. Teleost fish, representing more than 50 % of vertebrate species, are extremely efficient in utilizing low levels of dissolved oxygen in water. However, huge variations exist among various taxa of fish in their ability to tolerate hypoxia. In aquaculture, hypoxia tolerance is among the most important traits because hypoxia can cause major economic losses. Genetic enhancement for hypoxia tolerance in catfish is of great interest, but little was done with analysis of the genetic architecture of hypoxia tolerance. The objective of this study was to conduct a genome-wide association study to identify QTLs for hypoxia tolerance using the catfish 250K SNP array with channel catfish families from six strains. Multiple significant and suggestive QTLs were identified across and within strains. One significant QTL and four suggestive QTLs were identified across strains. Six significant QTLs and many suggestive QTLs were identified within strains. There were rare overlaps among the QTLs identified within the six strains, suggesting a complex genetic architecture of hypoxia tolerance. Overall, within-strain QTLs explained larger proportion of phenotypic variation than across-strain QTLs. Many of genes within these identified QTLs have known functions for regulation of oxygen metabolism and involvement in hypoxia responses. Pathway analysis indicated that most of these genes were involved in MAPK or PI3K/AKT/mTOR signaling pathways that were known to be important for hypoxia-mediated angiogenesis, cell proliferation, apoptosis and survival.

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

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

Suppression of adipocyte differentiation and lipid accumulation by stearidonic acid (SDA) in 3T3-L1 cells.

BACKGROUND: Increased consumption of omega-3 (ω-3) fatty acids found in cold-water fish and fish oil has been reported to protect against obesity. A potential mechanism may be through reduction in adipocyte differentiation. Stearidonic acid (SDA), a plant-based ω-3 fatty acid, has been targeted as a potential surrogate for fish-based fatty acids; however, its role in adipocyte differentiation is unknown. This study was designed to evaluate the effects of SDA on adipocyte differentiation in 3T3-L1 cells. METHODS: 3T3-L1 preadipocytes were differentiated in the presence of SDA or vehicle-control. Cell viability assay was conducted to determine potential toxicity of SDA. Lipid accumulation was measured by Oil Red O staining and triglyceride (TG) quantification in differentiated 3T3-L1 adipocytes. Adipocyte differentiation was evaluated by adipogenic transcription factors and lipid accumulation gene expression by quantitative real-time polymerase chain reaction (qRT-PCR). Fatty acid analysis was conducted by liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS). RESULTS: 3T3-L1 cells treated with SDA were viable at concentrations used for all studies. SDA treatment reduced lipid accumulation in 3T3-L1 adipocytes. This anti-adipogenic effect by SDA was a result of down-regulation of mRNA levels of the adipogenic transcription factors CCAAT/enhancer-binding proteins alpha and beta (C/EBPα, C/EBPß), peroxisome proliferator-activated receptor gamma (PPARγ), and sterol-regulatory element binding protein-1c (SREBP-1c). SDA treatment resulted in decreased expression of the lipid accumulation genes adipocyte fatty-acid binding protein (AP2), fatty acid synthase (FAS), stearoyl-CoA desaturase (SCD-1), lipoprotein lipase (LPL), glucose transporter 4 (GLUT4) and phosphoenolpyruvate carboxykinase (PEPCK). The transcriptional activity of PPARγ was found to be decreased with SDA treatment. SDA treatment led to significant EPA enrichment in 3T3-L1 adipocytes compared to vehicle-control. CONCLUSION: These results demonstrated that SDA can suppress adipocyte differentiation and lipid accumulation in 3T3-L1 cells through down-regulation of adipogenic transcription factors and genes associated with lipid accumulation. This study suggests the use of SDA as a dietary treatment for obesity.

Ribosomal protein genes are highly enriched among genes with allele-specific expression in the interspecific F1 hybrid catfish.

Interspecific hybrids provide a rich source for the analysis of allele-specific expression (ASE). In this work, we analyzed ASE in F1 hybrid catfish using RNA-Seq datasets. While the vast majority of genes were expressed with both alleles, 7-8 % SNPs exhibited significant differences in allele ratios of expression. Of the 66,251 and 177,841 SNPs identified from the datasets of the liver and gill, 5420 (8.2 %) and 13,390 (7.5 %) SNPs were identified as significant ASE-SNPs, respectively. With these SNPs, a total of 1519 and 3075 ASE-genes were identified. Gene Ontology analysis revealed that genes encoding cytoplasmic ribosomal proteins (RP) were highly enriched among ASE genes. Parent-of-origin was determined for 27 and 30 ASE RP genes in the liver and gill, respectively. The results indicated that genes from both channel catfish and blue catfish were involved in ASE. However, each RP gene appeared to be almost exclusively expressed from only one parent, indicating that ribosomes in the hybrid catfish were in the "hybrid" form. Overall representation of RP transcripts among the transcriptome appeared lower in the F1 hybrid catfish than in channel catfish or blue catfish, suggesting that the "hybrid" ribosomes may work more efficiently for translation in the F1 hybrid catfish.

A genome-wide association study in catfish reveals the presence of functional hubs of related genes within QTLs for columnaris disease resistance.

BACKGROUND: Columnaris causes severe mortalities among many different wild and cultured freshwater fish species, but understanding of host resistance is lacking. Catfish, the primary aquaculture species in the United States, serves as a great model for the analysis of host resistance against columnaris disease. Channel catfish in general is highly resistant to the disease while blue catfish is highly susceptible. F2 generation of hybrids can be produced where phenotypes and genotypes are segregating, providing a useful system for QTL analysis. To identify genes associated with columnaris resistance, we performed a genome-wide association study (GWAS) using the catfish 250 K SNP array with 340 backcross progenies derived from crossing female channel catfish (Ictalurus punctatus) with male F1 hybrid catfish (female channel catfish I. punctatus × male blue catfish I. furcatus). RESULTS: A genomic region on linkage group 7 was found to be significantly associated with columnaris resistance. Within this region, five have known functions in immunity, including pik3r3b, cyld-like, adcyap1r1, adcyap1r1-like, and mast2. In addition, 3 additional suggestively associated QTL regions were identified on linkage groups 7, 12, and 14. The resistant genotypes on the QTLs of linkage groups 7 and 12 were found to be homozygous with both alleles being derived from channel catfish. The paralogs of the candidate genes in the suggestively associated QTL of linkage group 12 were found on the QTLs of linkage group 7. Many candidate genes on the four associated regions are involved in PI3K pathway that is known to be required by many bacteria for efficient entry into the host. CONCLUSION: The GWAS revealed four QTLs associated with columnaris resistance in catfish. Strikingly, the candidate genes may be arranged as functional hubs; the candidate genes within the associated QTLs on linkage groups 7 and 12 are not only co-localized, but also functionally related, with many of them being involved in the PI3K signal transduction pathway, suggesting its importance for columnaris resistance.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach.

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Whole genome comparative analysis of channel catfish (Ictalurus punctatus) with four model fish species.

BACKGROUND: Comparative mapping is a powerful tool to study evolution of genomes. It allows transfer of genome information from the well-studied model species to non-model species. Catfish is an economically important aquaculture species in United States. A large amount of genome resources have been developed from catfish including genetic linkage maps, physical maps, BAC end sequences (BES), integrated linkage and physical maps using BES-derived markers, physical map contig-specific sequences, and draft genome sequences. Application of such genome resources should allow comparative analysis at the genome scale with several other model fish species. RESULTS: In this study, we conducted whole genome comparative analysis between channel catfish and four model fish species with fully sequenced genomes, zebrafish, medaka, stickleback and Tetraodon. A total of 517 Mb draft genome sequences of catfish were anchored to its genetic linkage map, which accounted for 62% of the total draft genome sequences. Based on the location of homologous genes, homologous chromosomes were determined among catfish and the four model fish species. A large number of conserved syntenic blocks were identified. Analysis of the syntenic relationships between catfish and the four model fishes supported that the catfish genome is most similar to the genome of zebrafish. CONCLUSION: The organization of the catfish genome is similar to that of the four teleost species, zebrafish, medaka, stickleback, and Tetraodon such that homologous chromosomes can be identified. Within each chromosome, extended syntenic blocks were evident, but the conserved syntenies at the chromosome level involve extensive inter-chromosomal and intra-chromosomal rearrangements. This whole genome comparative map should facilitate the whole genome assembly and annotation in catfish, and will be useful for genomic studies of various other fish species.

Bulk segregant RNA-seq reveals expression and positional candidate genes and allele-specific expression for disease resistance against enteric septicemia of catfish.

BACKGROUND: The application of RNA-seq has accelerated gene expression profiling and identification of gene-associated SNPs in many species. However, the integrated studies of gene expression along with SNP mapping have been lacking. Coupling of RNA-seq with bulked segregant analysis (BSA) should allow correlation of expression patterns and associated SNPs with the phenotypes. RESULTS: In this study, we demonstrated the use of bulked segregant RNA-seq (BSR-Seq) for the analysis of differentially expressed genes and associated SNPs with disease resistance against enteric septicemia of catfish (ESC). A total of 1,255 differentially expressed genes were found between resistant and susceptible fish. In addition, 56,419 SNPs residing on 4,304 unique genes were identified as significant SNPs between susceptible and resistant fish. Detailed analysis of these significant SNPs allowed differentiation of significant SNPs caused by genetic segregation and those caused by allele-specific expression. Mapping of the significant SNPs, along with analysis of differentially expressed genes, allowed identification of candidate genes underlining disease resistance against ESC disease. CONCLUSIONS: This study demonstrated the use of BSR-Seq for the identification of genes involved in disease resistance against ESC through expression profiling and mapping of significantly associated SNPs. BSR-Seq is applicable to analysis of genes underlining various performance and production traits without significant investment in the development of large genotyping platforms such as SNP arrays.