A genome-wide association study in catfish reveals the presence of functional hubs of related genes within QTLs for columnaris disease resistance
© Geng et al.; licensee BioMed Central. 2015
Received: 7 November 2014
Accepted: 27 February 2015
Published: 17 March 2015
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).
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.
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.
Commercial production of catfish (Ictalurus spp.) accounts for approximately 60% of US aquaculture production (www.ers.usda.gov). However, due to devastating diseases, the catfish industry has been hindered by unprecedented challenges. Diseases can cause losses of up to one third of the industry each year. Flavobacterium columnare, a Gram-negative bacterium, is the causative agent of columnaris disease, which is very common in wild and cultured freshwater fish worldwide . This pathogen can infect a variety of fish species through mucosal attachment points on the gill and skin, causing external erosion and necrosis . The bacterium can also enter the blood stream and invade the internal organs . The economically important foodfish channel catfish (Ictalurus punctatus) and the other members of the family Ictaluridae are extremely susceptible to columnaris disease. Columnaris disease is considered as one of the most important diseases in the catfish industry, causing tens of millions of dollars in losses every year . Therefore, improving disease resistance is a major goal for genetic stock enhancement using various approaches including strain selection, crossbreeding, hybridization, and transgenics .
Heterosis is an important genetic force that contributes to world food production. Interspecific hybrids are particularly effective in generating significant heterosis . Hybrids of catfish have been investigated for about 50 years . Among all possible combinations, only one cross (female channel catfish × male blue catfish) exhibited better performance than its parental species in growth rate and feed conversion efficiency .
Various techniques have been used to dissect the genes responsible for production and performance traits in aquaculture species. For instance, traditional QTL mapping was conducted to locate the region associated with resistance to infectious pancreatic necrosis virus in Atlantic salmon (Salmo salar) [8-13]. Utilizing RNA-seq, Li et al.  characterized the role of catfish intestinal epithelial barrier following enteric septicemia of catfish (ESC, Edwarsiella ictaluri) challenge. Wang et al.  conducted bulked segregant analysis to study candidate gene locations and allele-specific expression associated with ESC resistance in catfish. However, genetic analysis of resistance against bacterial diseases such as columnaris has been limited.
In the last decade, various studies were conducted aiming at elucidating the mechanisms of columnaris entry, immune evasion, and the host response to the disease. Some studies have been conducted on the modes and dynamics of columnaris adhesion [16,17]. In recent years, the application of RNA-seq has allowed a significantly greater level of understanding of the complexities of columnaris induced gene expression . Several central signatures following infection were revealed by gene expression enrichment analysis and gene pathway analysis of differentially expressed genes . For instance, Beck et al.  revealed that rhamnose-binding lectin was induced dramatically after infection, which was correlated with columnaris susceptibility. Peatman et al.  carried out RNA-seq analysis to compare basal and post-challenge differences in expression between susceptible and resistant channel catfish lines. Some genes involved in critical innate immunity, such as iNOS2b, lysozyme C, IL-8, and TNF-alpha were constitutively expressed higher in resistant than in susceptible catfish gill tissues. In contrast, secreted mucin forms, rhamnose-binding lectin, and some mucosal immune factors were found to be expressed at higher levels in the susceptible catfish line than in the resistant line. Despite these efforts, the knowledge of molecular mechanism of columnaris resistance is still limited.
Genome-wide association mapping allows detection of markers closely linked to QTLs. Although QTL mapping is well-suited for family-based samples, association studies, especially genome-wide association studies, can potentially offer higher mapping resolution using markers with higher density. Moreover, recent developments in GWAS methodologies have offered mature software packages for association analysis. Although GWAS has been widely utilized in economically important plants  and animals [22,23], it has been seldom applied in fish. To identify the QTLs in the hybrid catfish related to columnaris resistance, a genome-wide association study using backcross hybrid catfish was conducted, and here we report the identified QTLs and their associated genes within the highly associated genomic regions.
The pedigree information of catfish samples used in this study
Susceptible sample number
Resistant sample number
Analysis of linkage disequilibrium (LD) blocks
The LD block was defined as a set of contiguous SNPs with the minimum pairwise r2 value exceeding 0.50 . The number of independent SNP markers and LD blocks was 14,420. Thus the threshold P-value for genome-wide significance was 0.05/14420=3.47e-6 (−log10 (P-value)=5.46). The threshold P-value for the significance of “suggestive association”, which allows one false positive effect in a genome-wide test, was 1/14420=6.93e-5 (−log10 (P-value)=4.16).
Linkage groups with associated QTLs for columnaris resistance
Genomic region with significantly associated QTL for columnaris resistance
The significantly associated SNPs on linkage group 7
upstream of pik3r3b
downstream of akr1a1b
upstream of pik3r3b
EMMAX was used to investigate the contribution of the significantly associated QTL to the phenotype. Because of high correlation between SNPs in one locus (Figure 4 and Section Correlation of the SNPs associated with columnaris resistance), when analyzing the fraction of variance explained by the QTL, we chose only the most significant SNP (AX-86060479) to represent this region, which could explain 6.6% of the phenotypic variance. Nevertheless, binary phenotype and selective genotyping used in our study may cause potential overestimation of the QTL effect.
Based on the SNPs placed on the catfish 250 K SNP panel , the parental origins of the SNPs could be determined. All the 12 significant SNPs are interspecific, which means on these loci, two species are simply fixed for alternate alleles. The resistant genotypes for the 12 SNP loci are homozygous with both alleles being originated from channel catfish.
Genes located within the significantly associated QTL region for columnaris resistance
The genes within the genomic region harboring the significant SNPs associated with columnaris resistance were annotated by BLAST analysis against the non-redundant protein database . To provide additional supporting evidence for the proper annotation of the region, syntenic analyses were also conducted to compare the gene contents in the genomic neighborhood around the significant SNPs. As shown in Figure 4, the conserved synteny was identified between the catfish and zebrafish genomes. The flanking genes of catfish are all conserved with zebrafish except efcab7 and tspan1 are missing in catfish and the gene orders are slightly different.
Within the 620 Kb region containing the most significant SNPs, a total of 10 genes were identified (Figure 4). Of the 10 genes, five genes were found to have known functions in immunity, and these include phosphatidylinositol 3-kinase regulatory subunit gamma b (pik3r3b), cylindromatosis-like (cyld-like), pituitary adenylate cyclase-activating polypeptide type 1 receptor (adcyap1r1), adcyap1r1-like, and microtubule-associated serine and threonine kinase 2 (mast2). In order to be sure all the candidate genes were included in the analysis, the extended genomic region was examined, and no gene was found with known function in immunity.
Suggestively associated QTLs
The most significant SNPs with the closest and candidate genes in 3 suggestively associated regions
SNP position (bp)
-log 10 (P-value)
Gene position (bp)
cysteinyl leukotriene receptor 2
guanine nucleotide-binding protein (G protein) subunit beta 1*
voltage-dependent calcium channel subunit alpha 2/delta 2#
tumor suppressor candidate 2
diphosphoinositol polyphosphate phosphohydrolase 1
protein kinase C and casein kinase substrate in neurons protein 1
N-terminal EF-hand calcium-binding protein 1
probable G-protein coupled receptor 27
phosphatidylinositol 3-kinase regulatory subunit 5†
phosphatidylinositol 3-kinase regulatory subunit 6†
guanine nucleotide-binding protein subunit beta 3*
voltage-dependent calcium channel gamma 6 subunit#
chondroitin sulfate proteoglycan 4
Spectrin beta chain, non-erythrocytic 2
The first suggestively associated region is on linkage group 7. About 12 Mb downstream of the significantly associated region on linkage group 7, there is another locus suggestively associated with columnaris resistance. A series of SNPs on that region exhibit relatively low P values with the lowest ones listed on Table 3, with –log10 (P-value) ranging from 4.83-5.12. In addition to linkage group 7, there are SNPs on linkage group 12 that reach suggestive genome-wide significance (P-value < 6.93e-5). The most significant SNP could explain 5.5% of the variance. Compared with the region strongly associated with the columnaris resistance on linkage group 7, the region detected on linkage group 12 is larger. According to the linkage map , the recombinant frequency is very low on this region of linkage group 12. The distance of two most significant SNPs is 2.68 Mb on linkage group 12 expanding about 2 centimorgans, while in catfish, on average, 1 cM is equivalent to approximately 250 Kb. According to the suggestively associated interspecific SNPs on linkage group 12, the resistant genotypes of the SNPs are homozygous with two channel catfish alleles, like the SNPs on linkage group 7. On linkage group 14, there is only one SNP reaching suggestive Bonferroni genome-wide significance, explaining 4.6% of the variance. The candidate genes surrounding the most significant SNPs of these 3 loci are listed in Table 3.
Correlation of the SNPs associated with columnaris resistance
Conditioned analyses were conducted to examine the correlation of the SNPs associated with columnaris resistance . Genotypes of the most significant SNP (AX-86060479) on linkage group 7 were included as a covariate in the mixed linear model. After conditioning, the –log10 (P-value) of all the other SNPs on linkage group 7 dropped below 2.0, while the SNP P-values remained generally unchanged on the other linkage groups. On the linkage groups 12 and 14, there were also strong correlations among these associated SNPs within the same linkage group (data not shown).
In this study, for the first time, we identified a significantly associated QTL on linkage group 7 and three additional suggestive QTLs on linkage groups 7, 12, and 14 for columnaris resistance. The significant QTL on linkage group 7 was narrowed down to a small region of 620 Kb. Therefore, in spite of being just an initial quantitative analysis for the complex disease resistance trait, this work is very significant because it has set the foundation for fine mapping of the columnaris resistance genes, providing basis for application of genome technologies for catfish aquaculture through marker- or genome-based selection. Additional fine mapping using larger or more families could narrow down the candidate genes underlining columnaris disease resistance.
A number of genes involved in the PI3K pathway were found to be within the significantly associated genomic region of 620 Kb, suggesting the involvement of PI3K pathway in the resistance against columnaris. Among the 10 genes found in the 620 Kb region, five were involved in PI3K gene pathway. These are pik3r3b, cyld-like, adcyap1r1, adcyap1r1-like, and mast2. The phosphatidylinositol 3-kinase regulatory subunit gamma b (pik3r3b) is the closest gene to the most significantly associated SNP (AX-86060479) (Table 2). Although the causative SNP could be within any one of the 10 genes found in the 620 Kb region, the fact that PI3K pathway was reported to be involved in immunity and resistance makes them particularly interesting as candidate genes . In addition, many genes found in the suggestive QTL regions were also involved in the PI3K pathway (see below), further increasing the likelihood of PI3K pathway involvement in the resistance against columnaris.
PI3 kinases have been known to play important roles in innate and adaptive immunity [28-31]. For example, PI3K activity is important for NF-κB pathway activation in different mechanisms [32,33]. It was also shown that the infectious agents can manipulate the PI3K pathway to create a favorable environment by various mechanisms . PI3K is required for modifying the cytoskeleton dynamics, regulating membrane traffic, coordinating exocytic membrane insertion and pseudopod extension, which could be utilized by some pathogenic bacteria for entry into host cell [34-38]. Ireton et al.  showed that the protein InlB from Listeria monocytogenes is an agonist of mammalian PI3K, which causes rapid increases in cellular amounts of PI (3,4) P2 and PI (3,4,5) P3. Lambotin et al.  found Neisseria meningitidis requires cortactin recruitment by triggering a PI3K/Rac1 signaling to elicit an efficient invasion in non-phagocytic cells. Kierbel et al.  reported Pseudomonas aeruginosa requires the PI3K and Akt pathway for internalization. It was reported that Porphyromonas gingivalis could activate PI3K, which blocks phagocytosis and promotes inflammation . In addition, it was shown that blockade or deficiency of PI3Kγ significantly enhanced resistance against Leishmania mexicana, which revealed the unique role for Class IB PI3K in parasite invasion . The fact that pik3r3b is closest to the most significant SNP makes it an interesting candidate for future studies.
Pituitary adenylate cyclase-activating polypeptide type I receptor (ADCYAP1R1) is the pituitary adenylate cyclase-activating polypeptide (PACAP) specific receptor . PACAP could activate adenylate cyclase and phospholipase C (PLC) through an interaction with ADCYAP1R1 and stimulates cAMP and inositol phosphate formation [47,48]. In fish larvae, Lugo et al.  reported that the PACAP influences immune functions. They observed an elevated level of nitric oxide synthase-derived metabolites and total immunoglobulin M concentration in serum of juvenile catfish and tilapia after intraperitoneal injection of PACAP.
Microtubule-associated serine and threonine kinase 2 (MAST2) can interact with phosphatase and tensin homolog deleted on linkage group 10 (PTEN) which antagonizes PI3K-dependent signaling pathways [50,51]. MAST2 inhibits TNF receptor associated factor 6 (TRAF6) activity , which represents a molecular bridge for innate immunity and adaptive immunity . For example, TRAF6 could activate PI3K-dependent cytoskeletal changes and activate IκB kinase (IKK) in response to proinflammatory cytokines . Taken together, PI3 kinases themselves, or genes involved in PI3K pathway could play important roles in disease resistance.
Interestingly, many genes mapped within the suggestively associated QTL regions on linkage groups 7, 12, and 14 are functionally related to the genes mapped within the significant QTL region on linkage group 7, further supporting the possibility that PI3K pathway may be important for disease resistance against columnaris (Table 3, Figure 5). On linkage group 7, cysteinyl leukotriene receptor 2 is a G protein-coupled receptor (GPCR) with various functions such as modulation of chemokine gene transcription and calcium signaling [55,56]. G beta gamma activates the class IB p110 gamma/p101 PI3K gamma on the stimulation of GPCR . Voltage-dependent calcium channel subunit alpha-2/delta-2 gene was found in this region, and PI3K could enhance native voltage-dependent calcium channel currents . Diphosphoinositol polyphosphate phosphohydrolase 1 could cleave a beta-phosphate from the diphosphate groups in PP-InsP5. PP-InsP5 is similar to Ins (1,3,4,5) P4, the headgroup of PI (3,4,5) P3, which implied their competition relationship . Chockalingam  reported that alpha 1-syntrophin could bind PI (4,5) P2. Hyaluronidase-2, a glycosylphosphatidylinositol-anchored receptor, could hydrolyze hyaluronic acid which could be degraded by F. columnare . On linkage group 12, the genes within the suggestive QTL region seemed to be related with those within the QTL region on linkage group 7, because the paralogs of some genes in this region on linkage group 12 are found on linkage group 7 including phosphatidylinositol 3-kinase regulatory subunit 5 (p101-PI3K), phosphatidylinositol 3-kinase regulatory subunit 6 (p87-PI3K), guanine nucleotide-binding protein subunit beta 3, and voltage-dependent calcium channel gamma 6 subunit. On linkage group 14, spectrin beta chain, non-erythrocytic 2 (βIII spectrin) is located closest to the most significant SNP (AX-85234783). With a PH domain, βIII spectrin can bind PIP2 and get involved in membrane skeleton [61,62]. As presented in Figure 5, clearly many of these genes mapped in the significant QTL region or the suggestive QTL regions are involved in the related gene pathways.
It is notable that genes involved in the PI3K pathway were located together in “hubs” that were significantly associated with disease resistance. Theoretically, genes that are located together could be readily expressed in a coordinated fashion. However, here we do not have any evidence to indicate that the genes mapped within the QTLs involved in PI3K pathway are coordinately expressed. Analysis of RNA-seq data  indicated that some candidate genes indeed exhibited differences in baseline expression or after bacterial infection with columnaris. For instance, hyaluronidase-2 was expressed at a relative higher level in resistant fish than in susceptible fish before infection, and it was induced more in susceptible fish than resistant fish after infection. Guanine nucleotide-binding protein subunit beta-1, another gene that mapped within the suggestive QTL region, is expressed significantly higher in the susceptible fish than in the resistant fish. After infection, its expression in susceptible fish, but not in resistant fish, was drastically induced. However, because the experimental system is quite different, such expression data may not be directly transferable to our results here.
The most striking finding of our study was that the genes closest to the most significant SNPs were both positionally and functionally related, i.e., they are structurally organized as “functional hubs”. Although it is unknown at present whether these genes are expressed in a coordinated fashion, it was reported that neighboring genes tend to have similar expression patterns and get involved in related functions [63-65]. For instance, Schmid et al.  elucidated that genes in close proximity are much more likely to be co-expressed than expected by chance. Future analysis for coordinated expression of genes involved in PI3K pathway is warranted.
The QTLs identified in this study explained a limited fraction of the phenotypic variance of columnaris disease resistance. Firstly, the population specificity of QTLs is the most important reason why our family-based association mapping cannot detect all the QTLs associated with columnaris resistance . Even within the same strain, various families showed drastically different susceptibilities to columnaris disease . Secondly, segregating alleles within one species may lead to decreased power of analysis, especially in the case that one parental species systematically carries resistance alleles while the other one carries susceptibility alleles . Thus the region cannot be detected with strong significance using intraspecific SNPs. The associated SNPs on linkage group 14 are intraspecific, so the effect of the locus may be underestimated and we cannot infer the origins of the resistant alleles. Thirdly, because of the lack of recombination between nearby QTLs, the locus with minor effect cannot be detected if the favorite allele on a close QTL with a major effect have a different origin. Lastly, but not leastly, genome level variations such as allele variations can account for only a fraction of phenotypic variations. Gene expression regulations at various levels such as transcriptional and posttranscriptional levels, as well as environment and genotype-environment interactions can have a profound impact on the final phenotype in performance.
In summary, our GWAS using backcross interspecific hybrid population allowed mapping of associated QTLs and estimation of their effects for columnaris resistance. On linkage groups 7 and 12, the resistant genotypes were homozygous with both alleles from channel catfish. Examination of genes in the mapped QTL regions allowed further analysis of candidate disease resistance genes. It appears that signal transduction pathways involving PI3Ks may play a crucial role for disease resistance against columnaris. This notion is not only supported by the presence of PI3K pathway genes in the significantly associated QTL on linkage group 7, but also by the fact that many genes within the suggestive QTL on linkage group 12 are paralogs of those found on linkage group 7. In addition, many genes found within suggestive QTLs on linkage groups 7, 12 and 14 are also involved in PI3K pathway. Future studies are required for the identification of the causative genes for disease resistance. For example, GWAS using larger or more families can be conducted to increase the power and resolution. Ultimately, gene knockout experiments are needed to demonstrate the candidate genes as the disease resistance genes.
The most interesting discovery of this work is that functionally related genes that may be responsible for columnaris disease resistance are located closely in a limited number of positions, forming “functional hubs”. Future analysis of expression of genes in the PI3K pathway in relation to the resistance phenotype should determine whether the co-localized and functionally related genes are indeed expressed in a coordinated fashion.
All experiments involving the handling and treatment of fish were approved by the Institutional Animal Care and Use Committee (IACUC) at Auburn University. Tissue samples were collected after euthanasia. All animal procedures were carried out according to the Guide for the Care and Use of Laboratory Animals and the Animal Welfare Act in the United States.
Experimental fish, bacteria challenge and sample collection
The study population was the Auburn University one year old catfish generated from back cross of male F1 hybrid catfish (female channel catfish × male blue catfish) with female channel catfish. The backcross progenies were produced by using the F1 as the male parent to avoid possible maternal effects in the early growth stage. The population consisted of six families. Since the offspring were mixed for culture, the genotypes of the samples could be used to assign the offspring into their families (see Table 1). Three reasons make us to choose the backcross family-based population as the study sample. First, the channel catfish × blue catfish interspecific system provides a useful research system for the understanding of columnaris resistance because they exhibit clear contrast in their phenotypes . Segregation of genotypes in F2, along with the highly contrasted phenotypes, provides a good system for QTL analysis. Second, family-based association mapping is usually more powerful in detecting QTLs, since the lack of recombination between a QTL and linked marker increases the power of detection . Third, family-based association mapping is preferred to detect rare markers related with QTLs .
1200 fish (average body weight 55.3 grams) were randomly obtained from Auburn University Fish Genetics Facility and acclimated for one week in the aerated flow-through water [240 × 60 × 45 cm (L × W × H)]. The average water temperature was 25°C. A total of 340 backcross progenies were selected from the extremes of the disease resistance distribution of the 1200 fish based on the selective genotyping method . The first 169 fish that died of columnaris were continuously sampled as the susceptible group. When the fish lost balance, blood samples were collected. 171 fish that survived from the disease and showed no symptoms were selected randomly as the resistant group.
The bacteria challenge procedure was conducted as previously described . The bacteria F. columnare were provided by the Aquatic Microbiology Laboratory, Auburn University. To get a single F. columnare colony, several fish were experimentally infected with a virulent isolate (BGFS-27; genomovar II) , and bacteria were re-isolated from one symptomatic fish after confirmed visually and biochemically. We selected BGFS-27 as representative of F. columnare for the present study, to which channel catfish was more resistant than blue catfish . A single colony was cultured in modified Shieh broth for 24 h in a shaker incubator (100 rpm) at 28°C. The final concentration of the bacteria was determined using colony forming unit (CFU) per mL. Challenge experiments were then conducted by immersion exposure for 2 h at a final concentration of 3 × 106 CFU/mL. Control fish were treated with identical procedures except that they were exposed to sterile modified Shieh broth.
DNA isolation, genotyping, and quality control
DNA was isolated using standard protocols. Briefly, the blood samples were incubated at 55°C overnight and DNA was extracted twice with phenol and once with chloroform. DNA was precipitated by isopropanol and collected by brief centrifugation, washed twice with 70% ethanol, air-dried, and resuspended in reduced EDTA TE buffer (10 mM Tris–HCl, 0.1 mM EDTA, pH 8.0). DNA was quantified using spectroscopy by Nanodrop (Thermo Scientific). The integrity of DNA samples was checked by 1% agarose gel electrophoresis stained with ethidium bromide. DNA was diluted to 50 ng/uL and the quality of genomic DNA satisfied the requirement for the genotyping platform of SNP arrays.
We have developed a catfish 250 K SNP array using Affymetrix Axiom genotyping technology [25,72]. Genotyping using the catfish 250 K SNP array was performed at GeneSeek (Lincoln, Nebraska, USA). The informative SNPs in this GWAS were distributed across the catfish genome at an average interval of 3.6 Kb. No sample was excluded due to low quality or low call rate (<95%). 214,797 SNPs were kept after filtering out SNPs with an inheritance or genotyping error, a minor allele frequency (MAF) <5%, or a call rate < 95%.
Statistical analysis was carried out using the SVS software package (SNP & Variation Suite, Version 8.0). Pairwise linkage disequilibrium for the backcross progeny population was calculated according to r2 value. LD pruning was conducted with a window size of 50 SNPs, a step of 5 SNPs, and r2 threshold of 0.5, resulting in 14,420 independent SNP markers. The population structure was assessed by principal component analysis with the independent SNP markers.
The threshold P-value for genome-wide significance was calculated based on Bonferroni-correction with the estimated number of independent markers and LD blocks . A Manhattan plot of the P-value results was produced using the SVS software. Although channel catfish and blue catfish are two species, apparently their genome architecture is extremely similar according to our former studies and unpublished data [75-77]. Thus, the genetic marker map was constructed according to channel catfish genome sequence.
The upstream and downstream genes of the significant SNPs were determined. GENSCAN program  and FGENESH+  were used to analyze the catfish genome sequences (unpublished data) that surround the SNPs to identify the upstream and downstream genes. The identified genes were annotated by searching against the non-redundant protein database . Genomicus  was utilized to construct the synteny of the counterpart genes from zebrafish to provide evidence for orthology.
Availability of supporting data and requirements
No new SNPs were discovered in this manuscript. The SNPs used in this manuscript are from the catfish 250 K SNP array . The data sets supporting the results of this article are included within the article. The other information can be obtained from the corresponding author upon request.
This project was supported by Agriculture and Food Research Initiative Competitive Grant no. 2015-67015-22975 from the USDA National Institute of Food and Agriculture (NIFA), and partially supported by USDA Aquaculture Research Program Competitive Grant no. 2014-70007-22395. The authors wish to thank Ludmilla Kaltenboeck and Huseyin Kucuktas for their technical assistance. X. Geng was supported by scholarships from the China Scholarship Council (CSC). X. Geng received a Travel Award from Short Course on Statistical Genetics and Genomics at UAB (Grant No. R25GM093044). D. Zhi was supported by an NIH grant (R01 HG008115).
- Plumb JA, Hanson LA, Plumb JA. Health maintenance and principal: microbial diseases of cultured fishes. 3rd ed. Wiley-Blackwell: Ames, Iowa; 2011.Google Scholar
- Declercq AM, Haesebrouck F, Van den Broeck W, Bossier P, Decostere A. Columnaris disease in fish: a review with emphasis on bacterium-host interactions. Vet Res. 2013;44(1):27–44.PubMed CentralPubMedView ArticleGoogle Scholar
- Hawke JP, Thune RL. Systemic isolation and antimicrobial susceptibility of Cytophaga columnaris from commercially reared channel catfish. J Aquat Anim Health. 1992;4(2):109–13.View ArticleGoogle Scholar
- Arias CR, Cai W, Peatman E, Bullard SA. Catfish hybrid Ictalurus punctatus × I. furcatus exhibits higher resistance to columnaris disease than the parental species. Dis Aquat Organ. 2012;100(1):77–81.PubMedView ArticleGoogle Scholar
- Birchler JA, Auger DL, Riddle NC. In search of the molecular basis of heterosis. The Plant Cell Online. 2003;15(10):2236–9.View ArticleGoogle Scholar
- Giudice JJ. Growth of a blue × channel catfish hybrid as compared to its parent species. The Progressive Fish-Culturist. 1966;28(3):142–5.View ArticleGoogle Scholar
- Dunham RA, Umali GM, Beam R, Kristanto AH, Trask M. Comparison of production traits of NWAC103 channel catfish, NWAC103 channel catfish × blue catfish hybrids, Kansas select 21 channel catfish, and blue catfish grown at commercial densities and exposed to natural bacterial epizootics. N Am J Aquac. 2008;70(1):98–106.View ArticleGoogle Scholar
- Moen T, Baranski M, Sonesson AK, Kjøglum S. Confirmation and fine-mapping of a major QTL for resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar): population-level associations between markers and trait. BMC Genomics. 2009;10(1):368–82.PubMed CentralPubMedView ArticleGoogle Scholar
- Gheyas A, Haley C, Guy D, Hamilton A, Tinch A, Mota‐Velasco J, et al. Effect of a major QTL affecting IPN resistance on production traits in Atlantic salmon. Anim Genet. 2010;41(6):666–8.PubMedView ArticleGoogle Scholar
- Gheyas A, Houston R, Mota‐Velasco J, Guy D, Tinch A, Haley C, et al. Segregation of infectious pancreatic necrosis resistance QTL in the early life cycle of Atlantic Salmon (Salmo salar). Anim Genet. 2010;41(5):531–6.PubMedView ArticleGoogle Scholar
- Phillips R, Ventura A, Dekoning J, Nichols K. Mapping rainbow trout immune genes involved in inflammation reveals conserved blocks of immune genes in teleosts. Anim Genet. 2013;44(1):107–13.PubMedView ArticleGoogle Scholar
- Houston RD, Davey JW, Bishop SC, Lowe NR, Mota-Velasco JC, Hamilton A, et al. Characterisation of QTL-linked and genome-wide restriction site-associated DNA (RAD) markers in farmed Atlantic salmon. BMC Genomics. 2012;13(1):244.PubMed CentralPubMedView ArticleGoogle Scholar
- Houston R, Gheyas A, Hamilton A, Guy D, Tinch A, Taggart J, et al. Detection and confirmation of a major QTL affecting resistance to infectious pancreatic necrosis (IPN) in Atlantic salmon (Salmo salar). Dev Biol. 2008;132:199–204.Google Scholar
- Li C, Zhang Y, Wang R, Lu J, Nandi S, Mohanty S, et al. RNA-seq analysis of mucosal immune responses reveals signatures of intestinal barrier disruption and pathogen entry following Edwardsiella ictaluri infection in channel catfish, Ictalurus punctatus. Fish Shellfish Immunol. 2012;32(5):816–27.PubMedView ArticleGoogle Scholar
- Wang R, Sun L, Bao L, Zhang J, Jiang Y, Yao J, et al. Bulk segregant RNA-seq reveals expression and positional candidate genes and allele-specific expression for disease resistance against enteric septicemia of catfish. BMC Genomics. 2013;14(1):929–47.PubMed CentralPubMedView ArticleGoogle Scholar
- Decostere A, Haesebrouck F, Van Driessche E, Charlier G, Ducatelle R. Characterization of the adhesion of Flavobacterium columnare (Flexibacter columnaris) to gill tissue. J Fish Dis. 1999;22(6):465–74.View ArticleGoogle Scholar
- Olivares-Fuster O, Bullard SA, McElwain A, Llosa MJ, Arias CR. Adhesion dynamics of Flavobacterium columnare to channel catfish Ictalurus punctatus and zebrafish Danio rerio after immersion challenge. Dis Aquat Organ. 2011;96(3):221–7.PubMedView ArticleGoogle Scholar
- Sun F, Peatman E, Li C, Liu S, Jiang Y, Zhou Z, et al. Transcriptomic signatures of attachment, NF-kappaB suppression and IFN stimulation in the catfish gill following columnaris bacterial infection. Dev Comp Immunol. 2012;38(1):169–80.PubMedView ArticleGoogle Scholar
- Beck BH, Farmer BD, Straus DL, Li C, Peatman E. Putative roles for a rhamnose binding lectin in Flavobacterium columnare pathogenesis in channel catfish Ictalurus punctatus. Fish Shellfish Immunol. 2012;33(4):1008–15.PubMedView ArticleGoogle Scholar
- Peatman E, Li C, Peterson BC, Straus DL, Farmer BD, Beck BH. Basal polarization of the mucosal compartment in Flavobacterium columnare susceptible and resistant channel catfish (Ictalurus punctatus). Mol Immunol. 2013;56(4):317–27.PubMedView ArticleGoogle Scholar
- Li F, Chen B, Xu K, Wu J, Song W, Bancroft I, et al. Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (Brassica napus L.). DNA Res. 2014;21(4):355–67.PubMed CentralPubMedView ArticleGoogle Scholar
- Gu X, Feng C, Ma L, Song C, Wang Y, Da Y, et al. Genome-wide association study of body weight in chicken F2 resource population. PLoS One. 2011;6(7):e21872.PubMed CentralPubMedView ArticleGoogle Scholar
- Nishimura S, Watanabe T, Mizoshita K, Tatsuda K, Fujita T, Watanabe N, et al. Genome-wide association study identified three major QTL for carcass weight including the PLAG1-CHCHD7 QTN for stature in Japanese Black cattle. BMC Genet. 2012;13(1):40–51.PubMed CentralPubMedView ArticleGoogle Scholar
- Darvasi A, Soller M. Selective genotyping for determination of linkage between a marker locus and a quantitative trait locus. Theor Appl Genet. 1992;85(2–3):353–9.PubMedGoogle Scholar
- Liu S, Sun L, Li Y, Sun F, Jiang Y, Zhang Y, et al. Development of the catfish 250 K SNP array for genome-wide association studies. BMC Research Notes. 2014;7(1):135–47.PubMed CentralPubMedView ArticleGoogle Scholar
- Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35 suppl 1:D61–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Li Y, Liu S, Qin Z, Waldbieser G, Wang R, Sun L, et al. Construction of a high-density, high-resolution genetic map and its integration with BAC-based physical map in channel catfish. DNA Res. 2015;22(1):39–52.PubMedView ArticleGoogle Scholar
- Koyasu S. Role of phosphatidylinositol 3-kinase in the immune system. Tanpakushitsu Kakusan Koso. 2006;51(11):1569–79.PubMedGoogle Scholar
- Jiang K, Zhong B, Gilvary DL, Corliss BC, Hong-Geller E, Wei S, et al. Pivotal role of phosphoinositide-3 kinase in regulation of cytotoxicity in natural killer cells. Nat Immunol. 2000;1(5):419–25.PubMedView ArticleGoogle Scholar
- Vieira OV, Botelho RJ, Rameh L, Brachmann SM, Matsuo T, Davidson HW, et al. Distinct roles of class I and class III phosphatidylinositol 3-kinases in phagosome formation and maturation. J Cell Biol. 2001;155(1):19–25.PubMed CentralPubMedView ArticleGoogle Scholar
- Fukao T, Tanabe M, Terauchi Y, Ota T, Matsuda S, Asano T, et al. PI3K-mediated negative feedback regulation of IL-12 production in DCs. Nat Immunol. 2002;3(9):875–81.PubMedView ArticleGoogle Scholar
- Bone H, Williams NA. Antigen-receptor cross-linking and lipopolysaccharide trigger distinct phosphoinositide 3-kinase-dependent pathways to NF-kappa B activation in primary B cells. Int Immunol. 2001;13(6):807–16.PubMedView ArticleGoogle Scholar
- Reddy SAG, Huang JH, Liao WSL. Phosphatidylinositol 3-kinase in interleukin 1 signaling-Physical interaction with the interleukin 1 receptor and requirement in NF kappa B and AP-1 activation. J Biol Chem. 1997;272(46):29167–73.PubMedView ArticleGoogle Scholar
- Ireton K, Payrastre B, Cossart P. The Listeria monocytogenes protein InlB is an agonist of mammalian phosphoinositide 3-kinase. J Biol Chem. 1999;274(24):17025–32.PubMedView ArticleGoogle Scholar
- Cox D, Tseng CC, Bjekic G, Greenberg S. A requirement for phosphatidylinositol 3-kinase in pseudopod extension. J Biol Chem. 1999;274(3):1240–7.PubMedView ArticleGoogle Scholar
- Ireton K, Payrastre B, Chap H, Ogawa W, Sakaue H, Kasuga M, et al. A role for phosphoinositide 3-kinase in bacterial invasion. Science. 1996;274(5288):780–2.PubMedView ArticleGoogle Scholar
- Pizarro-Cerda J, Cossart P. Bacterial adhesion and entry into host cells. Cell. 2006;124(4):715–27.PubMedView ArticleGoogle Scholar
- Lambotin M, Hoffmann I, Laran-Chich MP, Nassif X, Couraud PO, Bourdoulous S. Invasion of endothelial cells by Neisseria meningitidis requires cortactin recruitment by a phosphoinositide-3-kinase/Rac1 signalling pathway triggered by the lipo-oligosaccharide. J Cell Sci. 2005;118(Pt 16):3805–16.PubMedView ArticleGoogle Scholar
- Kierbel A, Gassama-Diagne A, Mostov K, Engel JN. The phosphoinositol-3-kinase-protein kinase B/Akt pathway is critical for Pseudomonas aeruginosa strain PAK internalization. Mol Biol Cell. 2005;16(5):2577–85.PubMed CentralPubMedView ArticleGoogle Scholar
- Maekawa T, Krauss JL, Abe T, Jotwani R, Triantafilou M, Triantafilou K, et al. Porphyromonas gingivalis Manipulates Complement and TLR Signaling to Uncouple Bacterial Clearance from Inflammation and Promote Dysbiosis. Cell Host Microbe. 2014;15(6):768–78.PubMedView ArticleGoogle Scholar
- Cummings HE, Barbi J, Reville P, Oghumu S, Zorko N, Sarkar A, et al. Critical role for phosphoinositide 3-kinase gamma in parasite invasion and disease progression of cutaneous leishmaniasis. Proc Natl Acad Sci U S A. 2012;109(4):1251–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Massoumi R. Ubiquitin chain cleavage: CYLD at work. Trends Biochem Sci. 2010;35(7):392–9.PubMedView ArticleGoogle Scholar
- Sun SC. CYLD: a tumor suppressor deubiquitinase regulating NF-kappa B activation and diverse biological processes. Cell Death Differ.2010;17(1):25–34.PubMedView ArticleGoogle Scholar
- Yang W-L, Jin G, Li C-F, Jeong YS, Moten A, Xu D, et al. Cycles of ubiquitination and deubiquitination critically regulate growth factor-mediated activation of Akt signaling. Sci Signal. 2013;6(257):ra3.PubMedGoogle Scholar
- Gao J, Huo L, Sun X, Liu M, Li D, Dong J-T, et al. The tumor suppressor CYLD regulates microtubule dynamics and plays a role in cell migration. J Biol Chem. 2008;283(14):8802–9.PubMedView ArticleGoogle Scholar
- Arimura A, Shioda S. Pituitary adenylate cyclase activating polypeptide (PACAP) and its receptors: neuroendocrine and endocrine interaction. Front Neuroendocrinol. 1995;16(1):53–88.PubMedView ArticleGoogle Scholar
- Romanelli F, Fillo S, Isidori A, Conte D. Pituitary adenylate cyclase-activating polypeptide regulates rat Leydig cell function in vitro. Neuropeptides. 1997;31(4):311–7.PubMedView ArticleGoogle Scholar
- Bodart V, Babinski K, Ong H, De Lean A. Comparative effect of pituitary adenylate cyclase-activating polypeptide on aldosterone secretion in normal bovine and human tumorous adrenal cells. Endocrinology. 1997;138(2):566–73.PubMedGoogle Scholar
- Lugo JM, Carpio Y, Oliva A, Morales A, Estrada MP. Pituitary adenylate cyclase-activating polypeptide (PACAP): a regulator of the innate and acquired immune functions in juvenile fish. Fish Shellfish Immunol. 2010;29(3):513–20.PubMedView ArticleGoogle Scholar
- Downes C, Bennett D, McConnachie G, Leslie N, Pass I, MacPhee C, et al. Antagonism of PI 3-kinase-dependent signalling pathways by the tumour suppressor protein, PTEN. Biochem Soc Trans. 2001;29(Pt 6):846–51.PubMedView ArticleGoogle Scholar
- Terrien E, Chaffotte A, Lafage M, Khan Z, Prehaud C, Cordier F, et al. Interference with the PTEN-MAST2 interaction by a viral protein leads to cellular relocalization of PTEN. Sci Signal. 2012;5(237):ra58.PubMedGoogle Scholar
- Xiong H, Li H, Chen Y, Zhao J, Unkeless JC. Interaction of TRAF6 with MAST205 regulates NF-kappaB activation and MAST205 stability. J Biol Chem. 2004;279(42):43675–83.PubMedView ArticleGoogle Scholar
- Wu H, Arron JR. TRAF6, a molecular bridge spanning adaptive immunity, innate immunity and osteoimmunology. Bioessays. 2003;25(11):1096–105.PubMedView ArticleGoogle Scholar
- Wang KZ, Wara-Aswapati N, Boch JA, Yoshida Y, Hu C-D, Galson DL, et al. TRAF6 activation of PI 3-kinase-dependent cytoskeletal changes is cooperative with Ras and is mediated by an interaction with cytoplasmic Src. J Cell Sci. 2006;119(8):1579–91.PubMedView ArticleGoogle Scholar
- Thompson C, Cloutier A, Bossé Y, Poisson C, Larivée P, McDonald PP, et al. Signaling by the Cysteinyl-Leukotriene Receptor 2 involvement in chemokine gene transcription. J Biol Chem. 2008;283(4):1974–84.PubMedView ArticleGoogle Scholar
- Sjöström M, Johansson A-S, Schröder O, Qiu H, Palmblad J, Haeggström JZ. Dominant expression of the CysLT2 receptor accounts for calcium signaling by cysteinyl leukotrienes in human umbilical vein endothelial cells. Arterioscler Thromb Vasc Biol. 2003;23(8):e37–41.PubMedView ArticleGoogle Scholar
- Brock C, Schaefer M, Reusch HP, Czupalla C, Michalke M, Spicher K, et al. Roles of G beta gamma in membrane recruitment and activation of p110 gamma/p101 phosphoinositide 3-kinase gamma. J Cell Biol. 2003;160(1):89–99.PubMed CentralPubMedView ArticleGoogle Scholar
- Viard P, Butcher AJ, Halet G, Davies A, Nürnberg B, Heblich F, et al. PI3K promotes voltage-dependent calcium channel trafficking to the plasma membrane. Nat Neurosci. 2004;7(9):939–46.PubMedView ArticleGoogle Scholar
- Shears SB. Diphosphoinositol polyphosphates: metabolic messengers? Mol Pharmacol. 2009;76(2):236–52.PubMed CentralPubMedView ArticleGoogle Scholar
- Chockalingam PS, Gee SH, Jarrett HW. Pleckstrin homology domain 1 of mouse α1-syntrophin binds phosphatidylinositol 4, 5-bisphosphate. Biochemistry. 1999;38(17):5596–602.PubMedView ArticleGoogle Scholar
- Viel A, Branton D. Spectrin: on the path from structure to function. Curr Opin Cell Biol. 1996;8(1):49–55.PubMedView ArticleGoogle Scholar
- Holleran EA, Ligon LA, Tokito M, Stankewich MC, Morrow JS, Holzbaur EL. beta III spectrin binds to the Arp1 subunit of dynactin. J Biol Chem. 2001;276(39):36598–605.PubMedView ArticleGoogle Scholar
- Williams EJ, Bowles DJ. Coexpression of neighboring genes in the genome of Arabidopsis thaliana. Genome Res. 2004;14(6):1060–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Michalak P. Coexpression, coregulation, and cofunctionality of neighboring genes in eukaryotic genomes. Genomics. 2008;91(3):243–8.PubMedView ArticleGoogle Scholar
- Sémon M, Duret L. Evolutionary origin and maintenance of coexpressed gene clusters in mammals. Mol Biol Evol. 2006;23(9):1715–23.PubMedView ArticleGoogle Scholar
- Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, et al. A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005;37(5):501–6.PubMedView ArticleGoogle Scholar
- Luo C, Qu H, Ma J, Wang J, Li C, Yang C, et al. Genome-wide association study of antibody response to Newcastle disease virus in chicken. BMC Genet. 2013;14(1):42–51.PubMed CentralPubMedView ArticleGoogle Scholar
- LaFrentz BR, Shoemaker CA, Booth NJ, Peterson BC, Ourth DD. Spleen index and mannose-binding lectin levels in four channel catfish families exhibiting different susceptibilities to Flavobacterium columnare and Edwardsiella ictaluri. J Aquat Anim Health. 2012;24(3):141–7.PubMedView ArticleGoogle Scholar
- Ledur M, Navarro N, Pérez-Enciso M. Large-scale SNP genotyping in crosses between outbred lines: how useful is it ? Heredity. 2009;105(2):173–82.PubMedView ArticleGoogle Scholar
- Mackay I, Powell W. Methods for linkage disequilibrium mapping in crops. Trends Plant Sci. 2007;12(2):57–63.PubMedView ArticleGoogle Scholar
- Olivares-Fuster O, Arias CR. Development and characterization of rifampicin-resistant mutants from high virulent strains of Flavobacterium columnare. J Fish Dis. 2011;34(5):385–94.PubMedView ArticleGoogle Scholar
- Liu S, Zhou Z, Lu J, Sun F, Wang S, Liu H, et al. Generation of genome-scale gene-associated SNPs in catfish for the construction of a high-density SNP array. BMC Genomics. 2011;12:53–66.PubMed CentralPubMedView ArticleGoogle Scholar
- Kang HM, Sul JH, Service SK, Zaitlen NA, Kong S-y, Freimer NB, et al. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 2010;42(4):348–54.PubMed CentralPubMedView ArticleGoogle Scholar
- Duggal P, Gillanders EM, Holmes TN, Bailey-Wilson JE. Establishing an adjusted p-value threshold to control the family-wide type 1 error in genome wide association studies. BMC Genomics. 2008;9(1):516–24.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu Z, Karsi A, Li P, Cao D, Dunham R. An AFLP-based genetic linkage map of channel catfish (Ictalurus punctatus) constructed by using an interspecific hybrid resource family. Genetics. 2003;165(2):687–94.PubMed CentralPubMedGoogle Scholar
- Kucuktas H, Wang S, Li P, He C, Xu P, Sha Z, et al. Construction of genetic linkage maps and comparative genome analysis of catfish using gene-associated markers. Genetics. 2009;181(4):1649–60.PubMed CentralPubMedView ArticleGoogle Scholar
- Ninwichian P, Peatman E, Liu H, Kucuktas H, Somridhivej B, Liu S, et al. Second-generation genetic linkage map of catfish and its integration with the BAC-based physical map. G3: Genes| Genomes| Genetics. 2012;2(10):1233–41.PubMed CentralPubMedView ArticleGoogle Scholar
- Burge C, Karlin S. Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997;268(1):78–94.PubMedView ArticleGoogle Scholar
- Salamov AA, Solovyev VV. Ab initio gene finding in Drosophila genomic DNA. Genome Res. 2000;10(4):516–22.PubMed CentralPubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.PubMedView ArticleGoogle Scholar
- Muffato M, Louis A, Poisnel C-E, Crollius HR. Genomicus: a database and a browser to study gene synteny in modern and ancestral genomes. Bioinformatics. 2010;26(8):1119–21.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.