Second-generation genetic linkage map of catfish and its integration with the BAC-based physical map

Hydin as the Candidate Master Sex Determination Gene in Channel Catfish (Ictalurus punctatus) and Its Epigenetic Regulation

Sex determination is a fascinating area of research. To date, more than 20 master sex determination (SD) genes have been reported from vertebrate animals. With channel catfish (Ictalurus punctatus), much work has been conducted to determine its master SD gene, ranging from genetic linkage mapping, genome-wide association (GWA) analysis, genome sequencing, comparative genome analysis, epigenomic analysis, transcriptome analysis, and functional studies. Here in this mini review, we provide positional, expression, regulatory, and functional evidence supporting hydin (hydrocephalus-inducing protein or HYDIN axonemal central pair apparatus protein-like) as a master SD gene in channel catfish. Hydin is located within the sex determination region (SDR) within a mapped 8.9-Mb non-recombinational segment on chromosome 4 of channel catfish. It is highly expressed in genetic males, but not in genetic females. The alleles of X and Y are highly differentially methylated with the X chromosome being hypermethylated and the Y chromosome hypomethylated. The hypomethylated Y allele of hydin is expressed while the hypermethylated X allele is not expressed. Such allelic expression fits well with the XY sex determination system of channel catfish. Functional analysis using a methylation blocker, 5-aza-dC, demonstrated that demethylation, especially within the SDR, is accompanied with increased expression of hydin, which led to sex reversal of genetic females into phenotypic males. These evidences support the candidacy of hydin as a master SD gene in channel catfish. Future knockout and analysis of affected genes after hydin knockout should provide insights into how hydin functions as a master SD gene.

Chromosome-level assembly and annotation of the blue catfish Ictalurus furcatus, an aquaculture species for hybrid catfish reproduction, epigenetics, and heterosis studies

Background

The blue catfish is of great value in aquaculture and recreational fisheries. The F₁ hybrids of female channel catfish (Ictalurus punctatus) × male blue catfish (Ictalurusfurcatus) have been the primary driver of US catfish production in recent years because of superior growth, survival, and carcass yield. The channel–blue hybrid also provides an excellent model to investigate molecular mechanisms of environment-dependent heterosis. However, transcriptome and methylome studies suffered from low alignment rates to the channel catfish genome due to divergence, and the genome resources for blue catfish are not publicly available.

Results

The blue catfish genome assembly is 841.86 Mbp in length with excellent continuity (8.6 Mbp contig N50, 28.2 Mbp scaffold N50) and completeness (98.6% Eukaryota and 97.0% Actinopterygii BUSCO). A total of 30,971 protein-coding genes were predicted, of which 21,781 were supported by RNA sequencing evidence. Phylogenomic analyses revealed that it diverged from channel catfish approximately 9 million years ago with 15.7 million fixed nucleotide differences. The within-species single-nucleotide polymorphism (SNP) density is 0.32% between the most aquaculturally important blue catfish strains (D&B and Rio Grande). Gene family analysis discovered significant expansion of immune-related families in the blue catfish lineage, which may contribute to disease resistance in blue catfish.

Conclusions

We reported the first high-quality, chromosome-level assembly of the blue catfish genome, which provides the necessary genomic tool kit for transcriptome and methylome analysis, SNP discovery and marker-assisted selection, gene editing and genome engineering, and reproductive enhancement of the blue catfish and hybrid catfish.

Integration of Maps Enables a Cytogenomics Analysis of the Complete Karyotype in Solea senegalensis

The Pleuronectiformes order, which includes several commercially-important species, has undergone extensive chromosome evolution. One of these species is Solea senegalensis, a flatfish with 2n = 42 chromosomes. In this study, a cytogenomics approach and integration with previous maps was applied to characterize the karyotype of the species. Synteny analysis of S. senegalensis was carried out using two flatfish as a reference: Cynoglossus semilaevis and Scophthalmus maximus. Most S. senegalensis chromosomes (or chromosome arms for metacentrics and submetacentrics) showed a one-to-one macrosyntenic pattern with the other two species. In addition, we studied how repetitive sequences could have played a role in the evolution of S. senegalensis bi-armed (3, and 5–9) and acrocentric (11, 12 and 16) chromosomes, which showed the highest rearrangements compared with the reference species. A higher abundance of TEs (Transposable Elements) and other repeated elements was observed adjacent to telomeric regions on chromosomes 3, 7, 9 and 16. However, on chromosome 11, a greater abundance of DNA transposons was detected in interstitial BACs. This chromosome is syntenic with several chromosomes of the other two flatfish species, suggesting rearrangements during its evolution. A similar situation was also found on chromosome 16 (for microsatellites and low complexity sequences), but not for TEs (retroelements and DNA transposons). These differences in the distribution and abundance of repetitive elements in chromosomes that have undergone remodeling processes during the course of evolution also suggest a possible role for simple repeat sequences in rearranged regions.

High-Density Linkage Map and QTLs for Growth in Snapper (Chrysophrys auratus)

Characterizing the genetic variation underlying phenotypic traits is a central objective in biological research. This research has been hampered in the past by the limited genomic resources available for most non-model species. However, recent advances in sequencing technologies and related genotyping methods are rapidly changing this. Here we report the use of genome-wide SNP data from the ecologically and commercially important marine fish species Chrysophrys auratus (snapper) to 1) construct the first linkage map for this species, 2) scan for growth QTL, and 3) search for putative candidate genes in the surrounding QTL regions. The newly constructed linkage map contained ∼11K SNP markers and is one of the densest maps to date in the fish family Sparidae. Comparisons with genome scaffolds of the recently assembled snapper genome indicated that marker placement was mostly consistent between the scaffolds and linkage map (R = 0.7), but that at fine scales (< 5 cM) some precision limitations occurred. Of the 24 linkage groups, which likely reflect the 24 chromosomes of this species, three were found to contain QTL with genome-wide significance for growth-related traits. A scan of 13 candidate growth genes located the growth hormone, myogenin, and parvalbumin genes within 5.3, 9.6, and 25.0 cM of these QTL, respectively. The linkage map and QTL found in this study will advance the investigation of genome structure and aquaculture breeding efforts in this and related species.

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.

Restriction site-associated DNA sequencing for SNP discovery and high-density genetic map construction in southern catfish (Silurus meridionalis)

Single-nucleotide polymorphism (SNP) markers and high-density genetic maps are important resources for marker-assisted selection, mapping of quantitative trait loci (QTLs) and genome structure analysis. Although linkage maps in certain catfish species have been obtained, high-density maps remain unavailable in the economically important southern catfish (Silurus meridionalis). Recently developed restriction site-associated DNA (RAD) markers have proved to be a promising tool for SNP detection and genetic map construction. The objective of the present study was to construct a high-density linkage map using SNPs generated by next-generation RAD sequencing in S. meridionalis for future genetic and genomic studies. An F1 population of 100 individuals was obtained by intraspecific crossing of two wild heterozygous individuals. In total, 77 634 putative high-quality bi-allelic SNPs between the parents were discovered by mapping the parents' paired-end RAD reads onto the reference contigs from both parents, of which 54.7% were transitions and 45.3% were transversions (transition/transversion ratio of 1.2). Finally, 26 714 high-quality RAD markers were grouped into 29 linkage groups by using de novo clustering methods (Stacks). Among these markers, 4514 were linked to the female genetic map, 23 718 to the male map and 6715 effective loci were linked to the integrated map spanning 5918.31 centimorgans (cM), with an average marker interval of 0.89 cM. High-resolution genetic maps are a useful tool for both marker-assisted breeding and various genome investigations in catfish, such as sequence assembly, gene localization, QTL detection and genome structure comparison. Hence, such a high-density linkage map will serve as a valuable resource for comparative genomics and fine-scale QTL mapping in catfish species.

A Genome-Wide Association Study Reveals That Genes with Functions for Bone Development Are Associated with Body Conformation in Catfish

Body conformation is of great scientific and commercial interest for aquaculture fish species because it affects biological adaptation of the organism to environments, and is of economic importance to the aquaculture industry considering its direct effect on fillet yield. Catfish is the primary aquaculture species in the USA. Two major species used in the aquaculture industry, channel catfish and blue catfish, differ in body shape and therefore the backcross progenies serve as a good model for quantitative trait locus (QTL) analysis. Here, a genome-wide association study (GWAS) with hybrid catfish was conducted to identify the QTL for body conformation, including deheaded body length (DBL), body length (BL), body depth (BD), and body breadth (BB), which were all standardized by cubic root of body weight. Overall, the results indicate that the traits are polygenic. For DBL, linkage group (LG) 2 and LG 24 contain significant QTL, and LG 13 and LG 26 contain suggestively associated QTL (-log₁₀(P value) > 4.5). Compared with DBL, additional SNPs were identified to be associated with body length on LG 2, LG 7, and LG 18. Although no significant QTL for body depth was found, three suggestively associated QTLs were identified on LG 5, LG 13, and LG 14. No SNP for body breadth reached the threshold for suggestive association. Genes close to the associated SNPs were determined, many of which are known to be involved in bone development. This work therefore provides the basis for future identification of causal genes for the control of body conformation.

A High-Density Genetic Linkage Map and QTL Fine Mapping for Body Weight in Crucian Carp (Carassius auratus) Using 2b-RAD Sequencing

A high-resolution genetic linkage map is essential for a wide range of genetics and genomics studies such as comparative genomics analysis and QTL fine mapping. Crucian carp (Carassius auratus) is widely distributed in Eurasia, and is an important aquaculture fish worldwide. In this study, a high-density genetic linkage map was constructed for crucian carp using 2b-RAD technology. The consensus map contains 8487 SNP markers, assigning to 50 linkage groups (LGs) and spanning 3762.88 cM, with an average marker interval of 0.44 cM and genome coverage of 98.8%. The female map had 4410 SNPs, and spanned 3500.42 cM (0.79 cM/marker), while the male map had 4625 SNPs and spanned 3346.33 cM (0.72 cM/marker). The average recombination ratio of female to male was 2.13:1, and significant male-biased recombination suppressions were observed in LG47 and LG49. Comparative genomics analysis revealed a clear 2:1 syntenic relationship between crucian carp LGs and chromosomes of zebrafish and grass carp, and a 1:1 correspondence, but extensive chromosomal rearrangement, between crucian carp and common carp, providing evidence that crucian carp has experienced a fourth round of whole genome duplication (4R-WGD). Eight chromosome-wide QTL for body weight at 2 months after hatch were detected on five LGs, explaining 10.1-13.2% of the phenotypic variations. Potential candidate growth-related genes, such as an EGF-like domain and TGF-β, were identified within the QTL intervals. This high-density genetic map and QTL analysis supplies a basis for genome evolutionary studies in cyprinid fishes, genome assembly, and QTL fine mapping for complex traits in crucian carp.

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.

Development of a 690 K SNP array in catfish and its application for genetic mapping and validation of the reference genome sequence

Single nucleotide polymorphisms (SNPs) are capable of providing the highest level of genome coverage for genomic and genetic analysis because of their abundance and relatively even distribution in the genome. Such a capacity, however, cannot be achieved without an efficient genotyping platform such as SNP arrays. In this work, we developed a high-density SNP array with 690,662 unique SNPs (herein 690 K array) that were relatively evenly distributed across the entire genome, and covered 98.6% of the reference genome sequence. Here we also report linkage mapping using the 690 K array, which allowed mapping of over 250,000 SNPs on the linkage map, the highest marker density among all the constructed linkage maps. These markers were mapped to 29 linkage groups (LGs) with 30,591 unique marker positions. This linkage map anchored 1,602 scaffolds of the reference genome sequence to LGs, accounting for over 97% of the total genome assembly. A total of 1,007 previously unmapped scaffolds were placed to LGs, allowing validation and in few instances correction of the reference genome sequence assembly. This linkage map should serve as a valuable resource for various genetic and genomic analyses, especially for GWAS and QTL mapping for genes associated with economically important traits.

A Genome-Wide Association Study Identifies Multiple Regions Associated with Head Size in Catfish

Skull morphology is fundamental to evolution and the biological adaptation of species to their environments. With aquaculture fish species, head size is also important for economic reasons because it has a direct impact on fillet yield. However, little is known about the underlying genetic basis of head size. Catfish is the primary aquaculture species in the United States. In this study, we performed a genome-wide association study using the catfish 250K SNP array with backcross hybrid catfish to map the QTL for head size (head length, head width, and head depth). One significantly associated region on linkage group (LG) 7 was identified for head length. In addition, LGs 7, 9, and 16 contain suggestively associated regions for head length. For head width, significantly associated regions were found on LG9, and additional suggestively associated regions were identified on LGs 5 and 7. No region was found associated with head depth. Head size genetic loci were mapped in catfish to genomic regions with candidate genes involved in bone development. Comparative analysis indicated that homologs of several candidate genes are also involved in skull morphology in various other species ranging from amphibian to mammalian species, suggesting possible evolutionary conservation of those genes in the control of skull morphologies.

High-density interspecific genetic linkage mapping provides insights into genomic incompatibility between channel catfish and blue catfish

Catfish is the leading aquaculture species in the United States. The interspecific hybrid catfish produced by mating female channel catfish with male blue catfish outperform both of their parent species in a number of traits. However, mass production of the hybrids has been difficult because of reproductive isolation. Investigations of genome structure and organization of the hybrids provide insights into the genetic basis for maintenance of species divergence in the face of gene flow, thereby helping develop strategies for introgression and efficient production of the hybrids for aquaculture. In this study, we constructed a high-density genetic linkage map using the hybrid catfish system with the catfish 250K SNP array. A total of 26,238 SNPs were mapped to 29 linkage groups, with 12,776 unique marker positions. The linkage map spans approximately 3240 cM with an average intermarker distance of 0.25 cM. A fraction of markers (986 of 12,776) exhibited significant deviation from the expected Mendelian ratio of segregation, and they were clustered in major genomic blocks across 15 LGs, most notably LG9 and LG15. The distorted markers exhibited significant bias for maternal alleles among the backcross progenies, suggesting strong selection against the blue catfish alleles. The clustering of distorted markers within genomic blocks should lend insights into speciation as marked by incompatibilities between the two species. Such findings should also have profound implications for understanding the genomic evolution of closely related species as well as the introgression of hybrid production programs in aquaculture.

Construction of the High-Density Genetic Linkage Map and Chromosome Map of Large Yellow Croaker (Larimichthys crocea)

High-density genetic maps are essential for genome assembly, comparative genomic analysis and fine mapping of complex traits. In this study, 31,191 single nucleotide polymorphisms (SNPs) evenly distributed across the large yellow croaker (Larimichthys crocea) genome were identified using restriction-site associated DNA sequencing (RAD-seq). Among them, 10,150 high-confidence SNPs were assigned to 24 consensus linkage groups (LGs). The total length of the genetic linkage map was 5451.3 cM with an average distance of 0.54 cM between loci. This represents the densest genetic map currently reported for large yellow croaker. Using 2889 SNPs to target specific scaffolds, we assigned 533 scaffolds, comprising 421.44 Mb (62.04%) of the large yellow croaker assembled sequence, to the 24 linkage groups. The mapped assembly scaffolds in large yellow croaker were used for genome synteny analyses against the stickleback (Gasterosteus aculeatus) and medaka (Oryzias latipes). Greater synteny was observed between large yellow croaker and stickleback. This supports the hypothesis that large yellow croaker is more closely related to stickleback than to medaka. Moreover, 1274 immunity-related genes and 195 hypoxia-related genes were mapped to the 24 chromosomes of large yellow croaker. The integration of the high-resolution genetic map and the assembled sequence provides a valuable resource for fine mapping and positional cloning of quantitative trait loci associated with economically important traits in large yellow croaker.

Phenoscape: Identifying Candidate Genes for Evolutionary Phenotypes

Phenotypes resulting from mutations in genetic model organisms can help reveal candidate genes for evolutionarily important phenotypic changes in related taxa. Although testing candidate gene hypotheses experimentally in nonmodel organisms is typically difficult, ontology-driven information systems can help generate testable hypotheses about developmental processes in experimentally tractable organisms. Here, we tested candidate gene hypotheses suggested by expert use of the Phenoscape Knowledgebase, specifically looking for genes that are candidates responsible for evolutionarily interesting phenotypes in the ostariophysan fishes that bear resemblance to mutant phenotypes in zebrafish. For this, we searched ZFIN for genetic perturbations that result in either loss of basihyal element or loss of scales phenotypes, because these are the ancestral phenotypes observed in catfishes (Siluriformes). We tested the identified candidate genes by examining their endogenous expression patterns in the channel catfish, Ictalurus punctatus. The experimental results were consistent with the hypotheses that these features evolved through disruption in developmental pathways at, or upstream of, brpf1 and eda/edar for the ancestral losses of basihyal element and scales, respectively. These results demonstrate that ontological annotations of the phenotypic effects of genetic alterations in model organisms, when aggregated within a knowledgebase, can be used effectively to generate testable, and useful, hypotheses about evolutionary changes in morphology.

Second generation physical and linkage maps of yellowtail (Seriola quinqueradiata) and comparison of synteny with four model fish

Background

Physical and linkage maps are important aids for the assembly of genome sequences, comparative analyses of synteny, and to search for candidate genes by quantitative trait locus analysis. Yellowtail, Seriola quinqueradiata, is an economically important species in Japanese aquaculture, and genetic information will be useful for DNA-assisted breeding. We report the construction of a second generation radiation hybrid map, its synteny analysis, and a second generation linkage map containing SNPs (single nucleotide polymorphisms) in yellowtail.

Results

Approximately 1.4 million reads were obtained from transcriptome sequence analysis derived from 11 tissues of one individual. To identify SNPs, cDNA libraries were generated from a pool of 500 whole juveniles, and the gills and kidneys of 100 adults. 9,356 putative SNPs were detected in 6,025 contigs, with a minor allele frequency ≥25%. The linkage and radiation hybrid maps were constructed based on these contig sequences. 2,081 markers, including 601 SNPs markers, were mapped onto the linkage map, and 1,532 markers were mapped in the radiation hybrid map.

Conclusions

The second generation linkage and physical maps were constructed using 6,025 contigs having SNP markers. These maps will aid the de novo assembly of sequencing reads, linkage studies and the identification of candidate genes related to important traits. The comparison of marker contigs in the radiation hybrid map indicated that yellowtail is evolutionarily closer to medaka than to green-spotted pufferfish, three-spined stickleback or zebrafish. The synteny analysis may aid studies of chromosomal evolution in yellowtail compared with model fish.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1600-7) contains supplementary material, which is available to authorized users.

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.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

A high-density linkage map for Astyanax mexicanus using genotyping-by-sequencing technology

The Mexican tetra, Astyanax mexicanus, is a unique model system consisting of cave-adapted and surface-dwelling morphotypes that diverged >1 million years (My) ago. This remarkable natural experiment has enabled powerful genetic analyses of cave adaptation. Here, we describe the application of next-generation sequencing technology to the creation of a high-density linkage map. Our map comprises more than 2200 markers populating 25 linkage groups constructed from genotypic data generated from a single genotyping-by-sequencing project. We leveraged emergent genomic and transcriptomic resources to anchor hundreds of anonymous Astyanax markers to the genome of the zebrafish (Danio rerio), the most closely related model organism to our study species. This facilitated the identification of 784 distinct connections between our linkage map and the Danio rerio genome, highlighting several regions of conserved genomic architecture between the two species despite ~150 My of divergence. Using a Mendelian cave-associated trait as a proof-of-principle, we successfully recovered the genomic position of the albinism locus near the gene Oca2. Further, our map successfully informed the positions of unplaced Astyanax genomic scaffolds within particular linkage groups. This ability to identify the relative location, orientation, and linear order of unaligned genomic scaffolds will facilitate ongoing efforts to improve on the current early draft and assemble future versions of the Astyanax physical genome. Moreover, this improved linkage map will enable higher-resolution genetic analyses and catalyze the discovery of the genetic basis for cave-associated phenotypes.

Construction of a high-density, high-resolution genetic map and its integration with BAC-based physical map in channel catfish

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.

Identification and analysis of genome-wide SNPs provide insight into signatures of selection and domestication in channel catfish (Ictalurus punctatus)

Domestication and selection for important performance traits can impact the genome, which is most often reflected by reduced heterozygosity in and surrounding genes related to traits affected by selection. In this study, analysis of the genomic impact caused by domestication and artificial selection was conducted by investigating the signatures of selection using single nucleotide polymorphisms (SNPs) in channel catfish (Ictalurus punctatus). A total of 8.4 million candidate SNPs were identified by using next generation sequencing. On average, the channel catfish genome harbors one SNP per 116 bp. Approximately 6.6 million, 5.3 million, 4.9 million, 7.1 million and 6.7 million SNPs were detected in the Marion, Thompson, USDA103, Hatchery strain, and wild population, respectively. The allele frequencies of 407,861 SNPs differed significantly between the domestic and wild populations. With these SNPs, 23 genomic regions with putative selective sweeps were identified that included 11 genes. Although the function for the majority of the genes remain unknown in catfish, several genes with known function related to aquaculture performance traits were included in the regions with selective sweeps. These included hypoxia-inducible factor 1β· HIFιβ ¨ and the transporter gene ATP-binding cassette sub-family B member 5 (ABCB5). HIF1β· is important for response to hypoxia and tolerance to low oxygen levels is a critical aquaculture trait. The large numbers of SNPs identified from this study are valuable for the development of high-density SNP arrays for genetic and genomic studies of performance traits in catfish.

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

BACKGROUND

Cytochrome P450s (CYPs) encode one of the most diverse enzyme superfamily in nature. They catalyze oxidative reactions of endogenous molecules and exogenous chemicals.

METHODS

We identified CYPs genes through in silico analysis using EST, RNA-Seq and genome databases of channel catfish. Phylogenetic analyses and conserved syntenic analyses were conducted to determine their identities and orthologies. Meta-analysis of RNA-Seq databases was conducted to analyze expression profile of CYP genes following bacterial infection.

RESULTS

A full set of 61 CYP genes was identified and characterized in channel catfish. Phylogenetic tree and conserved synteny provided strong evidence of their identities and orthorlogy. Lineage-specific gene duplication was evident in a number of clans in channel catfish. CYP46A1 is missing in the catfish genome as observed with syntenic analysis and RT-PCR analysis. Thirty CYPs were found up- or down-regulated in liver, while seven and eight CYPs were observed regulated in intestine and gill following bacterial infection.

CONCLUSION

We systematically identified and characterized a full set of 61 CYP genes in channel catfish and studied their expression profiles after bacterial infection. While bacterial challenge altered the expression of large numbers of CYP genes, the mechanisms and significance of these changes are not known.

GENERAL SIGNIFICANCE

This work provides an example to systematically study CYP genes in non-model species. Moreover, it provides a basis for further toxicological and physiological studies in channel catfish.

Development of the catfish 250K SNP array for genome-wide association studies

BACKGROUND

Quantitative traits, such as disease resistance, are most often controlled by a set of genes involving a complex array of regulation. The dissection of genetic basis of quantitative traits requires large numbers of genetic markers with good genome coverage. The application of next-generation sequencing technologies has allowed discovery of over eight million SNPs in catfish, but the challenge remains as to how to efficiently and economically use such SNP resources for genetic analysis.

RESULTS

In this work, we developed a catfish 250K SNP array using Affymetrix Axiom genotyping technology. The SNPs were obtained from multiple sources including gene-associated SNPs, anonymous genomic SNPs, and inter-specific SNPs. A set of 640K high-quality SNPs obtained following specific requirements of array design were submitted. A panel of 250,113 SNPs was finalized for inclusion on the array. The performance evaluated by genotyping individuals from wild populations and backcross families suggested the good utility of the catfish 250K SNP array.

CONCLUSIONS

This is the first high-density SNP array for catfish. The array should be a valuable resource for genome-wide association studies (GWAS), fine QTL mapping, high-density linkage map construction, haplotype analysis, and whole genome-based selection.

Consolidation of the genetic and cytogenetic maps of turbot (Scophthalmus maximus) using FISH with BAC clones

Bacterial artificial chromosomes (BAC) have been widely used for fluorescence in situ hybridization (FISH) mapping of chromosome landmarks in different organisms, including a few in teleosts. In this study, we used BAC-FISH to consolidate the previous genetic and cytogenetic maps of the turbot (Scophthalmus maximus), a commercially important pleuronectiform. The maps consisted of 24 linkage groups (LGs) but only 22 chromosomes. All turbot LGs were assigned to specific chromosomes using BAC probes obtained from a turbot 5× genomic BAC library. It consisted of 46,080 clones with inserts of at least 100 kb and <5 % empty vectors. These BAC probes contained gene-derived or anonymous markers, most of them linked to quantitative trait loci (QTL) related to productive traits. BAC clones were mapped by FISH to unique marker-specific chromosomal positions, which showed a notable concordance with previous genetic mapping data. The two metacentric pairs were cytogenetically assigned to LG2 and LG16, and the nucleolar organizer region (NOR)-bearing pair was assigned to LG15. Double-color FISH assays enabled the consolidation of the turbot genetic map into 22 linkage groups by merging LG8 with LG18 and LG21 with LG24. In this work, a first-generation probe panel of BAC clones anchored to the turbot linkage and cytogenetical map was developed. It is a useful tool for chromosome traceability in turbot, but also relevant in the context of pleuronectiform karyotypes, which often show small hardly identifiable chromosomes. This panel will also be valuable for further integrative genomics of turbot within Pleuronectiformes and teleosts, especially for fine QTL mapping for aquaculture traits, comparative genomics, and whole-genome assembly.

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.

Construction of high-density genetic linkage maps for orange-spotted grouper Epinephelus coioides using multiplexed shotgun genotyping

Background

Orange-spotted grouper, Epinephelus coioides, is one of the most valuable fish species in China. Commercial production of orange-spotted grouper could be increased by developing higher growth rates and improving commercially important traits. Information on genetic markers associated with quantitative trait loci (QTL) can be used in breeding programs to identify and select individuals carrying desired traits. A high-density genetic linkage map is the basis for QTL study, and multiplexed shotgun genotyping (MSG) facilitates the development of single nucleotide polymorphisms (SNPs) and genotyping. In this study, the first high-density genetic linkage maps for groupers were generated on the basis of the MSG method.

Results

The sex-averaged map contained a total of 4,608 SNPs, which spanned 1581.7 cM, with a mean distance between SNPs of 0.34 cM. The 4,608 SNPs were located in 2,849 unique locations on the linkage map, with an average inter-location space at 0.56 cM. There were 2,516 SNPs on the female map, and the number of unique locus was 1,902. However, the male map contained more numbers of SNP (2,939) and unique locations (2,005). The total length of the female and male maps was 1,370.9 and 1,335.5 cM, respectively.

Conclusions

The high-resolution genetic linkage maps will be very useful for QTL analyses and marker-assisted selection (MAS) for economically important traits in molecular breeding of the orange-spotted grouper.

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.

Generation of physical map contig-specific sequences useful for whole genome sequence scaffolding

Along with the rapid advances of the nextgen sequencing technologies, more and more species are added to the list of organisms whose whole genomes are sequenced. However, the assembled draft genome of many organisms consists of numerous small contigs, due to the short length of the reads generated by nextgen sequencing platforms. In order to improve the assembly and bring the genome contigs together, more genome resources are needed. In this study, we developed a strategy to generate a valuable genome resource, physical map contig-specific sequences, which are randomly distributed genome sequences in each physical contig. Two-dimensional tagging method was used to create specific tags for 1,824 physical contigs, in which the cost was dramatically reduced. A total of 94,111,841 100-bp reads and 315,277 assembled contigs are identified containing physical map contig-specific tags. The physical map contig-specific sequences along with the currently available BAC end sequences were then used to anchor the catfish draft genome contigs. A total of 156,457 genome contigs (~79% of whole genome sequencing assembly) were anchored and grouped into 1,824 pools, in which 16,680 unique genes were annotated. The physical map contig-specific sequences are valuable resources to link physical map, genetic linkage map and draft whole genome sequences, consequently have the capability to improve the whole genome sequences assembly and scaffolding, and improve the genome-wide comparative analysis as well. The strategy developed in this study could also be adopted in other species whose whole genome assembly is still facing a challenge.

Male-biased genes in catfish as revealed by RNA-Seq analysis of the testis transcriptome

Background

Catfish has a male-heterogametic (XY) sex determination system, but genes involved in gonadogenesis, spermatogenesis, testicular determination, and sex determination are poorly understood. As a first step of understanding the transcriptome of the testis, here, we conducted RNA-Seq analysis using high throughput Illumina sequencing.

Methodology/Principal Findings

A total of 269.6 million high quality reads were assembled into 193,462 contigs with a N50 length of 806 bp. Of these contigs, 67,923 contigs had hits to a set of 25,307 unigenes, including 167 unique genes that had not been previously identified in catfish. A meta-analysis of expressed genes in the testis and in the gynogen (double haploid female) allowed the identification of 5,450 genes that are preferentially expressed in the testis, providing a pool of putative male-biased genes. Gene ontology and annotation analysis suggested that many of these male-biased genes were involved in gonadogenesis, spermatogenesis, testicular determination, gametogenesis, gonad differentiation, and possibly sex determination.

Conclusion/Significance

We provide the first transcriptome-level analysis of the catfish testis. Our analysis would lay the basis for sequential follow-up studies of genes involved in sex determination and differentiation in catfish.

Comparative genomic analysis of catfish linkage group 8 reveals two homologous chromosomes in zebrafish and other teleosts with extensive inter-chromosomal rearrangements

Background

Comparative genomics is a powerful tool to transfer genomic information from model species to related non-model species. Channel catfish (Ictalurus punctatus) is the primary aquaculture species in the United States. Its existing genome resources such as genomic sequences generated from next generation sequencing, BAC end sequences (BES), physical maps, linkage maps, and integrated linkage and physical maps using BES-associated markers provide a platform for comparative genomic analysis between catfish and other model teleost fish species. This study aimed to gain understanding of genome organizations and similarities among catfish and several sequenced teleost genomes using linkage group 8 (LG8) as a pilot study.

Results

With existing genome resources, 287 unique genes were identified in LG8. Comparative genome analysis indicated that most of these 287 genes on catfish LG8 are located on two homologous chromosomes of zebrafish, medaka, stickleback, and three chromosomes of green-spotted pufferfish. Large numbers of conserved syntenies were identified. Detailed analysis of the conserved syntenies in relation to chromosome level similarities revealed extensive inter-chromosomal and intra-chromosomal rearrangements during evolution. Of the 287 genes, 35 genes were found to be duplicated in the catfish genome, with the vast majority of the duplications being interchromosomal.

Conclusions

Comparative genome analysis is a powerful tool even in the absence of a well-assembled whole genome sequence. In spite of sequence stacking due to low resolution of the linkage and physical maps, conserved syntenies can be identified although the exact gene order and orientation are unknown at present. Through chromosome-level comparative analysis, homologous chromosomes among teleosts can be identified. Syntenic analysis should facilitate annotation of the catfish genome, which in turn, should facilitate functional inference of genes based on their orthology.

A dense genetic linkage map for common carp and its integration with a BAC-based physical map

Background

Common carp (Cyprinus carpio) is one of the most important aquaculture species with an annual global production of 3.4 million metric tons. It is also an important ornamental species as well as an important model species for aquaculture research. To improve the economically important traits of this fish, a number of genomic resources and genetic tools have been developed, including several genetic maps and a bacterial artificial chromosome (BAC)-based physical map. However, integrated genetic and physical maps are not available to study quantitative trait loci (QTL) and assist with fine mapping, positional cloning and whole genome sequencing and assembly. The objective of this study was to integrate the currently available BAC-based physical and genetic maps.

Results

The genetic map was updated with 592 novel markers, including 312 BAC-anchored microsatellites and 130 SNP markers, and contained 1,209 genetic markers on 50 linkage groups, spanning 3,565.9 cM in the common carp genome. An integrated genetic and physical map of the common carp genome was then constructed, which was composed of 463 physical map contigs and 88 single BACs. Combined lengths of the contigs and single BACs covered a physical length of 498.75 Mb, or around 30% of the common carp genome. Comparative analysis between common carp and zebrafish genomes was performed based on the integrated map, providing more insights into the common carp specific whole genome duplication and segmental rearrangements in the genome.

Conclusion

We integrated a BAC-based physical map to a genetic linkage map of common carp by anchoring BAC-associated genetic markers. The density of the genetic linkage map was significantly increased. The integrated map provides a tool for both genetic and genomic studies of common carp, which will help us to understand the genomic architecture of common carp and facilitate fine mapping and positional cloning of economically important traits for genetic improvement and modification.