- Research article
- Open Access
The Y chromosome sequence of the channel catfish suggests novel sex determination mechanisms in teleost fish
- Lisui Bao1,
- Changxu Tian1,
- Shikai Liu1,
- Yu Zhang1,
- Ahmed Elaswad1,
- Zihao Yuan1,
- Karim Khalil1,
- Fanyue Sun1,
- Yujia Yang1,
- Tao Zhou1,
- Ning Li1,
- Suxu Tan1,
- Qifan Zeng1,
- Yang Liu1,
- Yueru Li1,
- Yun Li1,
- Dongya Gao1,
- Rex Dunham1,
- Kenneth Davis2,
- Geoffrey Waldbieser2 and
- Zhanjiang Liu3Email authorView ORCID ID profile
© The Author(s). 2019
- Received: 6 November 2018
- Accepted: 9 January 2019
- Published: 25 January 2019
Sex determination mechanisms in teleost fish broadly differ from mammals and birds, with sex chromosomes that are far less differentiated and recombination often occurring along the length of the X and Y chromosomes, posing major challenges for the identification of specific sex determination genes. Here, we take an innovative approach of comparative genome analysis of the genomic sequences of the X chromosome and newly sequenced Y chromosome in the channel catfish.
Using a YY channel catfish as the sequencing template, we generated, assembled, and annotated the Y genome sequence of channel catfish. The genome sequence assembly had a contig N50 size of 2.7 Mb and a scaffold N50 size of 26.7 Mb. Genetic linkage and GWAS analyses placed the sex determination locus within a genetic distance less than 0.5 cM and physical distance of 8.9 Mb. However, comparison of the channel catfish X and Y chromosome sequences showed no sex-specific genes. Instead, comparative RNA-Seq analysis between females and males revealed exclusive sex-specific expression of an isoform of the breast cancer anti-resistance 1 (BCAR1) gene in the male during early sex differentiation. Experimental knockout of BCAR1 gene converted genetic males (XY) to phenotypic females, suggesting BCAR1 as a putative sex determination gene.
We present the first Y chromosome sequence among teleost fish, and one of the few whole Y chromosome sequences among vertebrate species. Comparative analyses suggest that sex-specific isoform expression through alternative splicing may underlie sex determination processes in the channel catfish, and we identify BCAR1 as a potential sex determination gene.
- Y chromosome
- Sex determination
The processes of sex determination are tremendously diverse in teleosts and include hermaphroditism and environmental sex determination as well as genetic sex determination . In most cases, the mechanisms of genetic sex determination of teleost fish are quite different from those of tetrapod, even though they originated from the same lineage about 450 million years ago, and share approximately 70% genome sequence similarity . Unlike in mammals and birds, where distinguishable sex chromosomes and common master sex-determining genes are present [3, 4], heterogametic sex chromosomes have only observed in approximately 270 species (less than 1%) of teleost fish; of those, approximately 70% are male heterogametic (XX females and XY males) and 30% are female heterogametic (ZZ males and ZW females) [5–7].
A wide variety of genetic sex-determining mechanisms have been identified in teleost species, with various genes serving as the “master sex-determining genes” such as DMRT1 in medaka species Oryzias latipes and O. curvinotus [8, 9], GSDF in the medaka species Oryzias luzonensis , SDY in rainbow trout (Oncorhynchus mykiss) , and AMHY in the Patagonian pejerrey (Odontesthes hatcheri) . In addition to specific sex-determining genes, a single nucleotide polymorphism (SNP) was reported to be responsible for sex determination in fugu (Takifugu rubripes) . Analysis of sex determination patterns in half-smooth tongue sole (Cynoglossus semilaevis) revealed that this species has a ZW type sex determination system, and that the dmrt1 gene exhibited characteristics of sex-determining genes such as sex chromosome linkage, male-specific expression, and essentiality for testis development , and that its knockout in ZZ fish (male karyotype) led to the development of female phenotypes .
In contrast to highly differentiated sex chromosomes in mammals, where the Y chromosome has lost most of its genes compared to the X chromosome, sex chromosomes of fish species are generally less differentiated. The mammalian Y chromosome has been genetically isolated without recombination with the X chromosome beyond the pseudoautosomal region (PAR) [16, 17], but recombination occurs along the length of the catfish Y and X chromosomes. Significant differences between the W and Z chromosomes have been reported in the female heterogametic half-smooth tongue sole where the W chromosome is more than 8 Mb larger than the Z chromosome. However, sex chromosome karyotypes, and presumably sex chromosome lengths, are highly similar in size in most fish species studied to date . Therefore, even though genetic studies can clearly map the sex determination gene to a chromosomal region, identification of sex determination gene is still a daunting task with teleost fish.
In this study, we adopted an innovative approach for the identification of the sex determination gene in channel catfish by comparative genome analysis of the genomic sequences of the X chromosome and the Y chromosome. Through sex reversal, we were able to produce XY phenotypic females, and mating of XY phenotypic females with the normal XY males allowed the production of offsprings carrying XX, XY, or YY sex chromosome sets. The YY fish offered a unique template for the sequencing and assembly of the sex chromosome sequences. Here, we report the generation of a high-quality whole genome assembly of a YY male fish, comparative genome and transcriptome analyses, and identification of the breast cancer anti-resistance 1 (BCAR1) gene as a candidate for the sex-determining locus in channel catfish.
Strategy for the identification of the sex determination gene
Resources for the channel catfish genome included a dense genetic map with > 253,000 SNP markers  and the reference genome sequence assembly . On average, 1 cM on the genetic map was equivalent to 230 Kb in the reference genome sequence assembly. Because the reference genome was generated using a gynogenetic female, only the X chromosome was represented; therefore, we used hormonal sex reversal- and marker-assisted mating to produce males with two Y chromosomes. A large number of individuals were analyzed to reduce the area of the sex determination locus to a small genomic region which was then located on the X and Y chromosome-specific sequence assemblies.
Mapping of the sex determination locus
Realizing the potential for long haplotype blocks in the F2 offspring, we performed GWAS on unrelated catfish to take advantage of historical recombination which is accumulated over multiple generations within populations to associate traits with individual sequence-level polymorphisms.
Sequencing and assembly of the YY genome
Sequencing and assembly of the genome of YY channel catfish
A. Sequences generated for the assembly of YY channel catfish genome
Number of reads
Average read length
Raw PacBio reads
Contig assembly and assembly polishing
Corrected PacBio reads
B. Assembly statistics
Total number of contigs
Total number of scaffold
Contig N50 (bp)
Scaffold N50 (bp)
Total assembled genome size in scaffolds (bp)
C. Assembly assessment
Total mapped BUSCOs
A total of 96.5% of the channel catfish genes matched the BUSCO set, of which 94.4% were completely matched, and 2.1% partially matched (Table 1c); the YY sequence assembly had a much greater contiguity than the reference genome sequence, as expected for a long sequence assembly, and the gene content also suggested a complete coverage of the YY genome.
Comparison of sequences and gene contents of X and Y chromosomes
Annotation of the channel catfish Y chromosome genes and their comparison with those on the X chromosome
X chromosome (Liu et al. 2016)
Number of scaffolds
Size without gaps (bp)
Size with repeats masked(bp)
SNPs on genetic map
A total of 950 genes were annotated on the assembled X chromosome, 950 genes were annotated on the Y chromosome, and gene identities were completely shared between the X and Y chromosomes. These results pointed to variation in genic or non-genic sequences in the sex determination region or differential gene expression as a trigger of sex determination.
Sexually differentially expressed genes during early stages of sex determination
Summary of de novo assembly of catfish transcriptome during sex determination stages using Trinity
Contigs (≥ 200 bp)
Large contigs (≥ 1000 bp)
Maximum length (bp)
Average length (bp)
Reads mapped to final reference (%)
Differentially expressed genes between male and female catfish 10–19 days post fertilization (dpf). Fish were sexed by using a microsatellite marker . Expression levels were determined using RNA-Seq datasets
High/low in male 10–14 dpf
High/low in male 15–19 dpf
Differentially expressed genes across whole genome
Differentially expressed genes on LG4
BCAR1 knockout of male genotypes led to female phenotypes
A great variety of sex determination mechanisms exist in fish among individuals within a population, contrary to the apparent conservation of the gene differentiation network in vertebrates . Generally, the sex determination mechanisms in fish include genetic and environmental regulation. In genetic systems, certain components become dominant in influencing the direction of sex determination . In environmental regulation, the environment in which fish live crucially affects the embryo development. In most cases, this influence is associated with the change of temperature . For example, normally, the sex of channel catfish is determined genetically by the XY system, but female-skewed sex ratios could be generated by applying high temperature during the critical period for sex determination . It is worthwhile to mention that high temperature induce expression of aromatase CYP19A1, a cytochrome that converts testosterone to estradiol, in turn skewing sex ratios toward females. If the same mechanism controls both the sex determination and sex reversal processes in catfish, BCAR1 may interact with estrogen receptor alpha and modulate an estrogen signaling pathway . This observation suggests that the sex of channel catfish is a combination of both genetic sex determination (GSD) and environmental sex determination (ESD), and similar phenomenon has been seen more widespread in fishes than previously believed. In teleost fish, both heterogametic male (XY) and heterogametic female systems exist, with XY system being the major sex determination system. For example, medaka , channel catfish , and rainbow trout  have an XY system. However, many species have female heterogamety (ZW). Fish such as the turbot (Scophthalmus maximus) , half-smooth tongue sole , have a ZW sex determination system.
Much effort has been devoted to analysis of sex determination genes [33–38]. In the genomic era, analysis of sex determination can be accelerated by understanding of sex chromosomes, their structure and organization, and functions [39–41]. However, because of technical difficulties, the whole Y chromosome sequences were generated for only a handful of organisms, all mammals up to date, including human , chimpanzee , rhesus macaque , mouse , swine , and gorilla .
Even though the X and Y chromosomes in mammals are well differentiated, sequencing the Y chromosome in mammalian species is complicated by high levels of repetitive sequence. Teleost X and Y chromosomes are less differentiated and highly similar, and in some cases, sex is determined by a single nucleotide polymorphism , making determination of sex chromosome sequences complicated. Sex reversal and production of YY males provided pure sequencing templates of the Y chromosome without homologous X chromosome sequence that could confound assembly . We previously generated a well-assembled XX reference genome from XX gynogen that is a doubled haploid without Y chromosome. The reference genome sequence assembled from this gynogen fish provided the X chromosome sequence. In the current project, we generated whole genome sequence assembly using the YY male fish as the sequencing template. As a result, the sex chromosome sequence assembled from the YY fish are from the Y chromosome without any confusion of the X chromosome sequences. These materials provided a unique system for comparative genome analysis of the X and Y sequences. This innovative approach should be applicable to most, if not all, fish species.
The assembled lengths of sequences from X and Y chromosomes differed by 3.9 Mb (34.59 Mb for X chromosome assembly and 38.49 Mb for the Y chromosome assembly). This difference in length between the X and Y chromosomes as indicated in Table 2 are caused by differences of the sequence assemblies, not by the different lengths of the chromosomes. It must be noted that the X and Y chromosome sequence assemblies were generated by application of different sequencing platforms. The X chromosome sequence was assembled by sequencing a gynogen fish (doubled haploid, no Y chromosome) using mostly the second-generation sequencing platform, Illumina sequencing, supplemented with a low level of PacBio sequencing, while the Y chromosome sequence was assembled by sequencing a YY male fish (no X chromosome) using the third-generation sequencing technology, the PacBio Sequencing. With short Illumina reads, repetitive sequences were not assembled into the scaffolds, leading to the shorter total sequences assembled from X chromosome. With PacBio long reads, the repetitive sequences were mostly read through, and thus, such repetitive sequences were assembled into the final genome assembly, resulted in longer assembled sequences of the Y chromosome. It is the unassembled short repetitive sequences in the sex determination region that made a difference of 3.9 Mb in the final assembly.
Although comparative analysis of sequences generated from different sequencing platforms can be less conclusive, we have the firsthand knowledge of the sequence assemblies of the X and Y chromosomes. In an effort to validate the whole genome sequence assembly, we have developed the 690 K SNP arrays that allowed us to determine the marker positions on the genetic linkage map and on the reference genome sequence assemblies . The collinearity of markers on the linkage map and on the reference genome sequence, and the collinearity of markers on the X and Y chromosomes provided significant confidence of the quality of the assemblies. In addition, to detect sex-specific gene sequences, we took two approaches: First, we conducted a reciprocal mapping of X and Y sequences by alignment of sequences generated from the gynogen and the YY fish after repeat masking. The results (Fig. 3) showed an even coverage across the whole chromosome, with no gaps or leftover gene sequences. Second, we annotated the genes in the assembled X and Y chromosomes, and the gene contents were identical. The unassembled short repetitive sequence pools were also thoroughly analyzed, without detection of any sex-specific sequences. Additionally, the RNA-Seq analysis also did not detect any sex-specific genes.
In this work, we produced a whole genome reference sequence using a YY male as sequencing template, allowing the generation of the first fully assembled Y chromosome sequence outside of mammals, as the very first Y chromosome sequence in teleost fish, and the seventh of all organisms. In addition, channel catfish is the leading aquaculture species in the USA, understanding of its sex determination not only important for research, but also have practical implications for aquaculture. In addition, the Y chromosome sequence should be a valuable resource for studies of sex chromosome biology and evolution.
With a well-assembled female reference genome sequence generated from a XX gynogen , and the YY male genome sequence produced in the current research, we were able to combine genome comparative analysis, genetic mapping, and GWAS to narrow the sex determination region to 8.9 Mb. Coupling the positional candidate genes with expression candidate genes during ovary differentiation stage revealed that only a specific transcript of one gene, BCAR1, demonstrated male-specific expression during the phenocritical period (10–19 dpf). Knockout of the BCAR1 gene showed evidence of genetic males with a female phenotype, which provided functional support for BCAR1 as a candidate sex determination gene in channel catfish. What is worth to mention, another transcript, Golgin subfamily B member 1 gene, were found to be expressed highly in female during 10–14 dpf, which means it highly likely participated in ovary differentiation. As it is involved in forming intercisternal cross-bridges of the Golgi complex, the activity of Golgi complex should increase in female during ovary differentiation, which acts to process and package the macromolecules such as proteins and lipids that are synthesized by the cell. Future studies are needed to annotate the function of this gene.
Six master sex-determining genes have been isolated in fish to date: DM-domain gene on the Y chromosome (dmY), the major testis-determining factor in the Japanese medaka (XX-XY) , anti-Müllerian hormone (amhY) in the Patagonian pejerrey  and anti-Müllerian hormone receptor, type II (amhr2) in fugu , gonadal somatic cell-derived factor (gsdf) in the Luzon ricefish (O. luzonensis) , sexually dimorphic on the Y chromosome gene (sdY) in the rainbow trout , and finally SRY-related HMG-Box gene 3 (sox3) in the medaka (Oryzias dancena) . Until 2011, all four vertebrate master SD genes were known to code for transcription factors which could have been construed as evidence that gonadal sex determination in vertebrates is always triggered by transcription factors. However, the three novel candidates for the master SD genes in the Patagonian pejerrey, Oryzias luzonensis, and fugu code for growth factors or one of their receptors. Thus, these findings suggest alternative mechanisms of genotypic sex determination, in which the main trigger is not constrained to be a transcription factor [37, 51].
In many non-mammalian vertebrate species of fish, newts, frogs, and turtles, exposure to estrogen during early gonadal development can override genetic sex and cause male-to-female sex reversal [52–57]. It is proposed that estrogen acts in the native pathway by which the bipotential gonad forms an ovary in these species . These results suggest that endogenous estrogen feminizes the medulla of the bipotential turtle gonad by inhibiting SOX9 expression .
The two amino acid changes observed between the X and Y haplotypes of BCAR1 are conservative, one contains the acidic glutamic acid or aspartic acid, the other contains the hydroxylic threonine or serine. The Thr735Ser is included in the Crk-associated substrate domain of BCAR1which is crucial for protein interaction of estrogen receptor. One hypothesis is that this substitution alters the interaction of BCAR1 with the estrogen receptor. It is reported that in breast cancer cells, overexpression of BCAR1 increases c-Src kinase activity and modulates non-genomic estrogen signaling . Interestingly, BCAR1 was identified in germ cell-Sertoli cell junction signaling pathway using pathway analysis with differentially expressed genes during ovary differentiation (Additional file 5: Figure S2). Sertoli cells are supporting cells within the developing testis, and the germ cell-Sertoli cell junction signaling is regulated by testosterone and TGF-β signaling pathway. Testosterone could be converted into estrogen by P450 aromatase; in channel catfish, such a conversion was observed to be almost always the case because both estrogens and testosterone led to sex reversal into females. TGF-β family members such as amh and gsdf are also of great importance in sex determination and gonadal somatic cell development. It is well known that sex determination gene initiates the differentiation of the supporting cells in the indifferent gonad to Sertoli cells, which subsequently form tight and adherent junctions to establish the testis cord . Taken all together, it is possible that the activation of BCAR1 in XY cells interacts with estrogen receptor-α and suppresses estrogen-based ovarian developmental signaling and drives the bipotential genital ridge toward testis (Fig. 7). Future research is warranted to test this hypothesis and to investigate the mechanisms of how BCAR1 may function as a sex determination gene.
Alternately, sex determination in catfish can be regulated at the RNA level, and the additional 102 bp in the male-specific transcript that is spliced in females contains 3 potential miRNA targeting sites. Future studies are warranted to determine the mechanisms of alternative splicing and its involvement in sex determination in catfish. The importance of miRNA target sites suggested that the alternatively spliced RNA may have a different fate upon interactions with miRNA, but these works are beyond the scope of this paper.
We sequenced and assembled the Y chromosome of channel catfish with a unique strategy of creating a male with two Y chromosomes (normal males are XY), thus avoiding any X chromosome sequences in the assembly. This is the first Y chromosome sequence outside of only a handful of mammalian Y chromosome sequences. We identified a candidate gene for sex determination in catfish through positional and functional evidence. Because BCAR1 was both a positional and expression candidate, and experimental knockout of the BCAR1 gene converted genetic males (XY) to phenotypic females, we posit BCAR1 as a candidate locus for sex determination in channel catfish.
Fish sources and sampling for whole genome sequencing
All procedures involving the handling and treatment of used fish during this study was approved by the Auburn University Institutional Animal Care and Use Committee (AU-IACUC) prior to initiation of the project. 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 USA. The wild channel catfish used in this project were obtained from 12 populations of 7 major watersheds in Alabama . The YY catfish genome donor was derived from the USDA103 catfish strain .
Generation of the YY male channel catfish genome donor
Experimental animals and tissue collection for RNA-Seq
Three-year-old channel catfish that displayed physical characteristics of sexual maturity were harvested from ponds (E.W. Shell Fisheries Research Center, Auburn University, Alabama, USA) and held in tanks. The general artificial spawning, fertilization, and husbandry procedures were performed as previously described [64, 65]. In brief, females were injected with luteinizing hormone-releasing hormone analog (LHRHa) at 100 μg/kg body weight, and then, eggs were hand-stripped and fertilized artificially. After hatching, fry were reared in a trough (304 cm × 60 cm × 30 cm) with daily water quality monitoring to maintain a constant water flow of 8–10 L/min, temperature at 28 °C, and dissolved oxygen at 6 ppm. Fry were fed daily to satiation with the Aquamax Fry Starter100 (Cat#: 000-5553, Aquamax, MO, USA). Approximately 50 fry were collected each day starting from the 10th day post fertilization (dpf) to 19th dpf. Fry were euthanized with tricaine methanesulfonate (MS 222) at 300 mg/L before sample collection. Each fry was cut into two portions of head and body. The samples were flash-frozen and stored at − 80 °C until DNA extraction and RNA extraction.
Thirty catfish fry samples were selected each day for sex-linked genotyping using the microsatellite marker AUEST0678 . Genomic DNA was extracted from the head sample of each catfish fry using the Gentra Puregene Tissue Kit (QIAGEN, USA), following manufacturer’s instructions. A touchdown PCR profile was performed as previously described , and amplification products were analyzed on a 7% polyacrylamide gel using a LICOR 4300 DNA Analyzer (LICOR Biosciences, Lincoln, NE). After gel electrophoresis, the sex of a catfish fry was determined based on the specific banding pattern that is one shared band (~ 212 bp) was observed in both males and females, while a 205-bp product was present only in male channel catfish .
DNA extraction, library preparation, and PacBio sequencing
Peripheral blood was collected from a YY male catfish and stored at − 80 °C until DNA extraction. To obtain high molecular weight genomic DNA, 15-μl whole blood was digested in a 500-μl solution containing 10-mM Tris-Cl [pH 8.0], 25-mM EDTA, 100-mM NaCl, 0.5% SDS, and 100-μg/ml proteinase K overnight at room temperature. Proteins were precipitated by addition of ½ volume 7.5 M ammonium acetate, gentle mixing for 10 min, and centrifugation for 10 min at 18,000 × g, 4 °C. The supernatant was decanted into a 2-ml microcentrifuge tube, and nucleic acids were precipitated by addition of 0.6 volumes isopropanol. The nucleic acids were pelleted by brief centrifugation, the supernatant was removed, and the pellet was dissolved in 400-μl TE buffer. Carbohydrates were removed by addition of 0.3 volumes ethanol, gentle mixing for 10 min, and then centrifugation for 10 min at 18,000 × g. The supernatant was decanted to a new microcentrifuge tube, and nucleic acids were re-precipitated with 1.5 volumes of ethanol, washed in 70% ethanol, and resuspended in TE buffer. The purified genomic DNA (A260/280 > 1.8, A260/230 > 2.0) was verified as high molecular weight via agarose gel electrophoresis. Long read sequencing was performed by Pacific Biosciences (Menlo Park, CA, USA) with two Sequel SMRT Cell 1 m v2 Trays on the Sequel platform and by the USDA, ARS, Genomics and Bioinformatics Research Unit (Stoneville, MS) with 20 SMRT cells on the Pacific Biosciences RS II platform.
De novo assembly of the YY genome
The bam files from the Sequel system were first converted to FASTQ format and assembled with sequences generated from PacBio RS II system. Both Canu  and Miniasm  were used in parallel to generate the assembly of YY genome of channel catfish. The best assembly was evaluated according to the number of contigs generated, N50 size, average contig size, and maximum contig length. Racon  was used to further omit the sequencing errors and call high-quality consensus for the Miniasm assembly. Quickmerge was used to merge the Canu and Miniasm assemblies to further improve contiguity and completeness. Finally, Arrow (Pacific Biosciences), an upgrade from Quiver , was used to polish the consensus sequences and further improve the accuracy of the assembly.
The YY genome scaffolds were produced using Chromosomer  to align and scaffold YY contigs against the XX reference genome . Where appropriate, estimated gap distances were transferred from the XX reference genome as ambiguous nucleotides (Ns). Chromosome level scaffolding and further correction was conducted by integrating the YY genome assembly with the high-density SNP-based genetic map . The final genetic linkage map was comprised of 54,342 markers (at 29,081 unique positions) distributed across a genetic distance of 3505.4 cM. The 70 bp of sequence flanking each SNP locus on the genetic map was mapped to the assembly scaffolds by BLASTN alignment (E value maximal threshold 1e−10, > 95% sequence similarity, > 65-bp alignment length). Scaffolds that were adjacent on a chromosome were manually joined with a string of 100 “Ns” to represent the gaps between the two adjacent scaffolds based on the high-resolution genetic map.
The assembly quality was assessed by comparison with the benchmarking universal single-copy orthologs (BUSCO; ) with the “vertebrate” data set and with the channel catfish reference genome.
QTL mapping and GWAS of the sex determination region
For QTL mapping of the sex determination (SD) locus, we first constructed a genetic map by genotyping 187 fish from a single family using the catfish 250 K SNP array . The sex phenotype of each fish was determined by examination of reproductive tissues after anesthetization and recorded as a binary trait. Next, 3418 informative SNPs on the genetic map were selected based on the pedigree information and quality value. Linkage between each marker and the SD locus was estimated using the “scanone” function with EM algorithm and binary model integrated in R/qtl software . Significant LOD threshold was obtained by permutation tests with 1000 replicates. The 95% Bayes confidence interval was determined using the “bayesint” function.
For GWAS analysis of the SD locus, 199 wild catfish were genotyped using the catfish 250 K SNP array and their genders were determined using sex-linked markers. The analysis was performed with a linear regression model using SNP & Variation Suite 7 (Golden Helix, Inc., Bozeman, MT). Samples passing the quality control (dish value > 0.85), and SNP call rate threshold (> 95%) were retained for analysis. SNPs with missing genotypes > 5% and minor allele frequency < 5% were removed. The Manhattan plot of significant SNPs was generated by Circos .
Annotation of the YY genome assembly
A library of catfish-specific repetitive sequences was first produced using RepeatModeler v1.0.8 (http://www.repeatmasker.org/RepeatModeler/). Next, the repeat-annotated library was used to mask the genome assembly using RepeatMasker with the option of -nolow and -s. Finally, Augustus  and Fgenesh (http://www.softberry.com) were used to predict genes within the repeat-masked assembly. Gene model parameters for Augustus were trained from conserved genes from vertebrate species using CEGMA . The amino acids predicted from the same genomic positon shared by both algorithms were retained and integrated for the downstream homolog-based annotation. Then, the amino acid sequences of the predicted genes were retrieved and queried against the Uniprot and NCBI non-redundant (NR) databases using BLASTP to identify the homologous genes with an E value maximum threshold of 1e−5. The names of predicted channel catfish genes were assigned based on their homologous proteins. All of the genes with names of retrotransposable element, transposon, reverse transcriptase, RNA-directed DNA polymerase from transposon, transposase, and transposable elements were treated as potential repetitive elements and removed from the gene annotation. All of the genes with names of uncharacterized, unnamed, unknown, hypothetical, and predicted protein were treated as unnamed genes and searched against Uniprot and NR database as above to return the top 10 matches. Any matches with names of potential repetitive elements were excluded from the annotation. The chromosomal position of each gene was assigned based on the chromosome level assembly, and genes with adjacent locations and identical gene names were labeled. Then, the tandem duplicated genes among the labeled genes were identified using MCScanX , and the remaining genes were treated as fragments of a same gene and subgrouped separately based on their position and names. Only the longest protein of each subgroup was kept representing the corresponding gene.
Comparative analysis of X and Y chromosome sequences
In order to identify the Y-specific fragments, the X chromosome from the reference assembly and the Y chromosome from the present assembly were aligned using MUMmer . Also, the male and female genome assemblies were compared using BLASTN (E value < 1e−10, and identity > 95%. MUMmer was also used to perform the genome-wide alignment.
High throughput RNA sequencing
A total of 10 female and 10 male fry with conclusive determination of sexes were randomly selected from each day for RNA-Seq sequencing. The tissues of 10 male catfish fry collected in the same day were pooled and homogenized with mortar and pestle in the presence of liquid nitrogen. Similarly, the tissues of 10 female catfish fry that were collected in the same day were also pooled and homogenized. Total RNA was extracted using the RNeasy Kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s instructions. RNA concentration and integrity were measured on an Agilent 2100 Bioanalyzer. Equal amounts of RNA from each 5-day period were pooled from 10 to 14 dpf and 15 to 19 dpf for male and female, respectively. These four pooled samples were outsourced for RNA-Seq at the HudsonAlpha Genomic Services Lab (Huntsville, AL, USA). The RNA-Seq libraries were prepared following the standard TruSeq protocols and were sequenced on an Illumina HiSeq 2500 instrument for 100 bp paired-end reads.
De novo transcriptome assembly and annotation
The de novo assembly of short reads was performed using Trinity (2012-10-05 release) with RNA-Seq reads sequenced from all eight samples. Before assembly, raw reads were trimmed by removing adaptor sequences and ambiguous nucleotides. Reads with quality scores < 20 and length < 30 bp were filtered, and the resulting high-quality sequences were used in the subsequent assembly. The short reads were assembled into unique transcripts using “Inchworm” with a greedy k-mer extension (k-mer = 25). After mapping reads to Inchworm-assembled contigs, “Chrysalis” incorporated reads into deBruijn graphs and the “Butterfly” module processed the individual graphs to generate full-length transcripts. In order to reduce redundancy, the assembly was passed to the CD-Hit (version 4.5.4) for multiple alignments . The threshold of sequence identity was set at 1.0, with the alignment covering > 90% of the shorter sequence. The assembled transcript contigs were annotated by searching against the NCBI zebrafish RefSeq protein database, the UniProtKB/SwissProt database, and the non-redundant (nr) protein database using the BLASTX program with E value ≤ 1e−5.
Identification of differentially expressed genes
High-quality reads from each sample were mapped onto the assembled transcript sequences using Bowtie software  with default parameters. The RSEM program  was then used to estimate the expression abundance of the transcripts. The total number of mapped reads for each transcript was determined and then normalized to determine FPKM (fragments per kilobase of transcript per million mapped reads). Differential expression analysis was performed using the “edgeR” package . All the samples were normalized by the trimmed mean of M values (TMM)  which equated overall gene expression levels between samples under the assumption that the majority were not differentially expressed. Transcripts with absolute fold change values of greater than 1.5 and a false discovery rate-adjusted P value < 0.05 were included in the analysis as differentially expressed genes. Ingenuity Pathway Analysis (IPA) approaches were used to analyze the differential expression data.
Identification of miRNA in 3′-UTR of BCAR1
The sequences of miRNA of channel catfish were retrieved from previous study . Potential miRNA targeting sites on 3′-UTR of BCAR1 were identified using RNA22 v2 microRNA target detection  with minimum number of 15 paired-up bases in heteroduplex and maximum folding energy as − 15 for heteroduplex (Kcal/mol).
Single guide RNA (sgRNA) sequences targeting the BCAR1 gene were searched against the whole genome and any other hits with identity larger than 95% and alignment length larger than 16 bp were labeled as potential off-targets. SgRNA with 0 off-targets were selected, designed and obtained from Transposagen (Lexington, KY), and microinjected together with the Cas9 nuclease mRNA (target sequence1 5′-GGAGAGGGACACGCAAGGCC-3′; target sequence2 5′-GGGCGCCAGGGCATTGTACC-3′). Channel catfish eggs were manually stripped from female broodstock, fertilized with normal sperm, and incubated in fresh water for 30 min. One hundred eggs were placed on a 10-ml Petri dish and 100 pg sgRNA, and 1 ng Cas9 mRNA were injected into the blastodisc using an Eppendorf Microinjector 5242 system (Hamburg, Germany). The control group was microinjected with the same solution without mRNAs. Then, embryos were moved into 10-L tubs filled with Holfreter’s solution containing 10 ppm doxycycline and incubated until hatch. Dead embryos were removed, and water was changed daily. Channel catfish fry were then transferred into a recirculating system.
TA clone and sequencing
To identify specific modifications of the BCAR1 gene after the knockout experiment, primers were designed flanking the target sites and PCR products were amplified from each individual using the IsoPure DNA Purification Kit (Denville Scientific, Holliston, MA, USA). Purified PCR products were cloned into the TOPO TA Cloning Kit for Sequencing (Invitrogen, Grand Island, NY), and ten transformed recombinant colonies per individual fish were selected for Sanger sequencing (forward primer: 5′-TGGCCAAGGCTCTGTATGAC-3′; reverse primer 5′-CACTTGTGGGTGGTACCTGG-3′) (Eurofins Genomics, Louisville, KY).
We thank Milla Kaltenboeck, Ruijia Wang, Luyang Sun and all the other members of The Fish Molecular Genetics and Biotechnology Laboratory for fruitful discussions.
This project was supported by a USDA grant (grant no. 2017-67015-26295) from the Agriculture and Food Research Initiative (AFRI) under the Animal Breeding, Genetics and Genomics Program and USDA, ARS, CRIS 6402-31000-008-00D. Lisui Bao was supported by a scholarship from the China Scholarship Council (CSC), and a fellowship from the Alabama EPSCoR Graduate Research Scholars Program.
Availability of data and materials
Raw data sets and genome assembly from this study have been submitted to the NCBI BioProject (http://www.ncbi.nlm.nih.gov/bioproject) under accession number PRJNA427343.
This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession PKUT00000000. The version described in this paper is version PKUT01000000.
All authors contributed to the project design. LB, CT, SL, AE, ZY, KK, YJ, TZ, NL, ST, QZ, YL, YL, and RD performed the CRISPR experiments and collected the data. YZ and FY performed the RNA-Seq analysis. KD and GW generated the YY catfish. LB assembled the genome and performed the association analysis and comparison analysis. LB and ZL drafted the full manuscript. All authors critically revised and approved the final manuscript.
Ethics approval and consent to participate
All procedures involving the handling and treatment of used fish during this study was approved by the Auburn University Institutional Animal Care and Use Committee (AU-IACUC) prior to initiation of the project.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Charlesworth B. Sex determination: primitive Y chromosomes in fish. Curr Biol. 2004;14(18):R745–7.View ArticleGoogle Scholar
- Le Page Y, Diotel N, Vaillant C, Pellegrini E, Anglade I, Merot Y, Kah O. Aromatase, brain sexualization and plasticity: the fish paradigm. Eur J Neurosci. 2010;32(12):2105–15.View ArticleGoogle Scholar
- Sinclair AH, Berta P, Palmer MS, Hawkins JR, Griffiths BL, Smith MJ, Foster JW, Frischauf AM, Lovellbadge R, Goodfellow PN. A gene from the human sex-determining region encodes a protein with homology to a conserved DNA-binding motif. Nature. 1990;346(6281):240–4.View ArticleGoogle Scholar
- Smith CA, Roeszler KN, Ohnesorg T, Cummins DM, Farlie PG, Doran TJ, Sinclair AH. The avian Z-linked gene DMRT1 is required for male sex determination in the chicken. Nature. 2009;461(7261):267–71.View ArticleGoogle Scholar
- Mank JE, Avise JC: Evolutionary diversity and turn-over of sex determination in teleost fishes. Sex Dev 2009, 3(2-3):60-67Google Scholar
- Kottler VA, Schartl M. The colorful sex chromosomes of teleost fish. Genes (Basel). 2018;9(5):233.Google Scholar
- Mank JE, Promislow DEL, Avise JC. Evolution of alternative sex-determining mechanisms in teleost fishes. Biol J Linn Soc. 2006;87(1):83–93.View ArticleGoogle Scholar
- Matsuda M, Nagahama Y, Shinomiya A, Sato T, Matsuda C, Kobayashi T, Morrey CE, Shibata N, Asakawa S, Shimizu N, et al. DMY is a Y-specific DM-domain gene required for male development in the medaka fish. Nature. 2002;417(6888):559–63.Google Scholar
- Matsuda M, Nagahama Y, Kobayashi T, Matsuda C, Hamaguchi S, Sakaizumi M. The sex determining gene of medaka: a Y-specific DM domain gene (DMY) is required for male development. Fish Physiol Biochem. 2003;28(1–4):135–9.View ArticleGoogle Scholar
- Shibata Y, Paul-Prasanth B, Suzuki A, Usami T, Nakamoto M, Matsuda M, Nagahama Y. Expression of gonadal soma derived factor (GSDF) is spatially and temporally correlated with early testicular differentiation in medaka. Gene Expr Patterns. 2010;10(6):283–9.Google Scholar
- Yano A, Guyomard R, Nicol B, Jouanno E, Quillet E, Klopp C, Cabau C, Bouchez O, Fostier A, Guiguen Y. An immune-related gene evolved into the master sex-determining gene in rainbow trout, Oncorhynchus mykiss. Curr Biol. 2012;22(15):1423–8.View ArticleGoogle Scholar
- Hattori RS, Murai Y, Oura M, Masuda S, Majhi SK, Sakamoto T, Fernandino JI, Somoza GM, Yokota M, Strussmann CA. A Y-linked anti-Mullerian hormone duplication takes over a critical role in sex determination. Proc Natl Acad Sci U S A. 2012;109(8):2955–9.View ArticleGoogle Scholar
- Kamiya T, Kai W, Tasumi S, Oka A, Matsunaga T, Mizuno N, Fujita M, Suetake H, Suzuki S, Hosoya S, et al. A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (Fugu). PLoS Genet. 2012;8(7): e1002798.Google Scholar
- Chen SL, Zhang GJ, Shao CW, Huang QF, Liu G, Zhang P, Song WT, An N, Chalopin D, Volff JN, et al. Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle. Nat Genet. 2014;46(3):253.View ArticleGoogle Scholar
- Cui Z, Liu Y, Wang W, Wang Q, Zhang N, Lin F, Wang N, Shao C, Dong Z, Li Y, et al. Genome editing reveals dmrt1 as an essential male sex-determining gene in Chinese tongue sole (Cynoglossus semilaevis). Sci Rep. 2017;7:42213.View ArticleGoogle Scholar
- Graves JA, Wakefield MJ, Toder R. The origin and evolution of the pseudoautosomal regions of human sex chromosomes. Hum Mol Genet. 1998;7(13):1991–6.View ArticleGoogle Scholar
- Bachtrog D. Y-chromosome evolution: emerging insights into processes of Y-chromosome degeneration. Nat Rev Genet. 2013;14(2):113–24.View ArticleGoogle Scholar
- Zeng Q, Fu Q, Li Y, Waldbieser G, Bosworth B, Liu S, Yang Y, Bao L, Yuan Z, Li N, et al. Development of a 690 K SNP array in catfish and its application for genetic mapping and validation of the reference genome sequence. Sci Rep. 2017;7:40347.View ArticleGoogle Scholar
- Liu Z, Liu S, Yao J, Bao L, Zhang J, Li Y, Jiang C, Sun L, Wang R, Zhang Y, et al. The channel catfish genome sequence provides insights into the evolution of scale formation in teleosts. Nat Commun. 2016;7:11757.Google Scholar
- Li Y, Liu S, Qin Z, Waldbieser G, Wang R, Sun L, Bao L, Danzmann RG, Dunham R, Liu Z. 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.View ArticleGoogle Scholar
- Simmons M, Mickett K, Kucuktas H, Li P, Dunham R, Liu ZJ. Comparison of domestic and wild channel catfish (Ictalurus punctatus) populations provides no evidence for genetic impact. Aquaculture. 2006;252(2–4):133–46.View ArticleGoogle Scholar
- Liu S, Sun L, Li Y, Sun F, Jiang Y, Zhang Y, Zhang J, Feng J, Kaltenboeck L, Kucuktas H, et al. Development of the catfish 250K SNP array for genome-wide association studies. BMC Res Notes. 2014;7:135.View ArticleGoogle Scholar
- Tiersch TR, Simco BA, Davis KB, Chandler RW, Wachtel SS, Carmichael GJ. Stability of genome size among stocks of the channel catfish. Aquaculture. 1990;87(1):15–22.View ArticleGoogle Scholar
- Sun FY, Liu SK, Gao XY, Jiang YL, Perera D, Wang XL, Li C, Sun LY, Zhang JR, Kaltenboeck L, et al. Male-biased genes in catfish as revealed by RNA-Seq analysis of the testis transcriptome. PLoS One. 2013;8(7):e68452.Google Scholar
- Patino R, Davis KB, Schoore JE, Uguz C, Strussmann CA, Parker NC, Simco BA, Goudie CA. Sex differentiation of channel catfish gonads: normal development and effects of temperature. J Exp Zool. 1996;276(3):209–18.View ArticleGoogle Scholar
- Ninwichian P, Peatman E, Perera D, Liu S, Kucuktas H, Dunham R, Liu Z. Identification of a sex-linked marker for channel catfish. Anim Genet. 2012;43(4):476–7.View ArticleGoogle Scholar
- Graves JAM, Peichel CL. Are homologies in vertebrate sex determination due to shared ancestry or to limited options? Genome Biol. 2010;11(4):205.Google Scholar
- Devlin RH, Nagahama Y. Sex determination and sex differentiation in fish: an overview of genetic, physiological, and environmental influences. Aquaculture. 2002;208(3–4):191–364.View ArticleGoogle Scholar
- Diaz N, Piferrer F. Lasting effects of early exposure to temperature on the gonadal transcriptome at the time of sex differentiation in the European sea bass, a fish with mixed genetic and environmental sex determination. BMC Genomics. 2015;16:679.Google Scholar
- Cabodi S, Moro L, Baj G, Smeriglio M, Di Stefano P, Gippone S, Surico N, Silengo L, Turco E, Tarone G, et al. p130Cas interacts with estrogen receptor alpha and modulates non-genomic estrogen signaling in breast cancer cells. J Cell Sci. 2004;117(Pt 8):1603–11.View ArticleGoogle Scholar
- Martinez P, Bouza C, Hermida M, Fernandez J, Toro MA, Vera M, Pardo B, Millan A, Fernandez C, Vilas R, et al. Identification of the major sex-determining region of turbot (Scophthalmus maximus). Genetics. 2009;183(4):1443–52.View ArticleGoogle Scholar
- Liao X, Xu G, Chen SL: Molecular method for sex identification of half-smooth tongue sole (Cynoglossus semilaevis) using a novel sex-linked microsatellite marker. Int J Mol Sci 2014, 15(7):12952-12958.Google Scholar
- Foster JW, Brennan FE, Hampikian GK, Goodfellow PN, Sinclair AH, Lovell-Badge R, Selwood L, Renfree MB, Cooper DW, Graves JA. Evolution of sex determination and the Y chromosome: SRY-related sequences in marsupials. Nature. 1992;359(6395):531–3.View ArticleGoogle Scholar
- Crews D, Bergeron JM, McLachlan JA. The role of estrogen in turtle sex determination and the effect of PCBs. Environ Health Perspect. 1995;103 Suppl 7:73–7.PubMedGoogle Scholar
- Morais da Silva S, Hacker A, Harley V, Goodfellow P, Swain A, Lovell-Badge R. Sox9 expression during gonadal development implies a conserved role for the gene in testis differentiation in mammals and birds. Nat Genet. 1996;14(1):62–8.View ArticleGoogle Scholar
- Kondo M, Hornung U, Nanda I, Imai S, Sasaki T, Shimizu A, Asakawa S, Hori H, Schmid M, Shimizu N, et al. Genomic organization of the sex-determining and adjacent regions of the sex chromosomes of medaka. Genome Res. 2006;16(7):815–26.View ArticleGoogle Scholar
- Kikuchi K, Hamaguchi S. Novel sex-determining genes in fish and sex chromosome evolution. Dev Dyn. 2013;242(4):339–53.View ArticleGoogle Scholar
- Shao C, Li Q, Chen S, Zhang P, Lian J, Hu Q, Sun B, Jin L, Liu S, Wang Z, et al. Epigenetic modification and inheritance in sexual reversal of fish. Genome Res. 2014;24(4):604–15.View ArticleGoogle Scholar
- Karmin M, Saag L, Vicente M, Wilson Sayres MA, Jarve M, Talas UG, Rootsi S, Ilumae AM, Magi R, Mitt M, et al. A recent bottleneck of Y chromosome diversity coincides with a global change in culture. Genome Res. 2015;25(4):459–66.View ArticleGoogle Scholar
- Small CM, Bassham S, Catchen J, Amores A, Fuiten AM, Brown RS, Jones AG, Cresko WA. The genome of the Gulf pipefish enables understanding of evolutionary innovations. Genome Biol. 2016;17(1):258.View ArticleGoogle Scholar
- Davidson WS, Koop BF, Jones SJ, Iturra P, Vidal R, Maass A, Jonassen I, Lien S, Omholt SW. Sequencing the genome of the Atlantic salmon (Salmo salar). Genome Biol. 2010;11(9):403.PubMedPubMed CentralGoogle Scholar
- Skaletsky H, Kuroda-Kawaguchi T, Minx PJ, Cordum HS, Hillier L, Brown LG, Repping S, Pyntikova T, Ali J, Bieri T, et al. The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes. Nature. 2003;423(6942):825–37.View ArticleGoogle Scholar
- Hughes JF, Skaletsky H, Pyntikova T, Graves TA, van Daalen SK, Minx PJ, Fulton RS, McGrath SD, Locke DP, Friedman C, et al. Chimpanzee and human Y chromosomes are remarkably divergent in structure and gene content. Nature. 2010;463(7280):536–9.View ArticleGoogle Scholar
- Hughes JF, Skaletsky H, Page DC. Sequencing of rhesus macaque Y chromosome clarifies origins and evolution of the DAZ (deleted in AZoospermia) genes. Bioessays. 2012;34(12):1035–44.View ArticleGoogle Scholar
- Soh YQ, Alfoldi J, Pyntikova T, Brown LG, Graves T, Minx PJ, Fulton RS, Kremitzki C, Koutseva N, Mueller JL, et al. Sequencing the mouse Y chromosome reveals convergent gene acquisition and amplification on both sex chromosomes. Cell. 2014;159(4):800–13.View ArticleGoogle Scholar
- Skinner BM, Sargent CA, Churcher C, Hunt T, Herrero J, Loveland JE, Dunn M, Louzada S, Fu B, Chow W, et al. The pig X and Y chromosomes: structure, sequence, and evolution. Genome Res. 2016;26(1):130–9.View ArticleGoogle Scholar
- Tomaszkiewicz M, Rangavittal S, Cechova M, Campos Sanchez R, Fescemyer HW, Harris R, Ye D, O'Brien PC, Chikhi R, Ryder OA, et al. A time- and cost-effective strategy to sequence mammalian Y chromosomes: an application to the de novo assembly of gorilla Y. Genome Res. 2016;26(4):530–40.View ArticleGoogle Scholar
- Nanda I, Kondo M, Hornung U, Asakawa S, Winkler C, Shimizu A, Shan ZH, Haaf T, Shimizu N, Shima A, et al. A duplicated copy of DMRT1 in the sex-determining region of the Y chromosome of the medaka, Oryzias latipes. Proc Natl Acad Sci U S A. 2002;99(18):11778–83.View ArticleGoogle Scholar
- Myosho T, Otake H, Masuyama H, Matsuda M, Kuroki Y, Fujiyama A, Naruse K, Hamaguchi S, Sakaizumi M. Tracing the emergence of a novel sex-determining gene in medaka, Oryzias luzonensis. Genetics. 2012;191(1):163.View ArticleGoogle Scholar
- Takehana Y, Matsuda M, Myosho T, Suster ML, Kawakami K, Shin IT, Kohara Y, Kuroki Y, Toyoda A, Fujiyama A, et al. Co-option of Sox3 as the male-determining factor on the Y chromosome in the fish Oryzias dancena. Nat Commun. 2014;5:4157.View ArticleGoogle Scholar
- Graves JA. How to evolve new vertebrate sex determining genes. Dev Dyn. 2013;242(4):354–9.View ArticleGoogle Scholar
- Edwards TM, Moore BC, Guillette LJ Jr. Reproductive dysgenesis in wildlife: a comparative view. Int J Androl. 2006;29(1):109–21.View ArticleGoogle Scholar
- Matthiessen P, Sumpter JP. Effects of estrogenic substances in the aquatic environment. Exs. 1998;86:319–35.PubMedGoogle Scholar
- Crews D, Bergeron JM. Role of reductase and aromatase in sex determination in the red-eared slider (Trachemys scripta), a turtle with temperature-dependent sex determination. J Endocrinol. 1994;143(2):279–89.View ArticleGoogle Scholar
- Wibbels T, Crews D. Putative aromatase inhibitor induces male sex determination in a female unisexual lizard and in a turtle with temperature-dependent sex determination. J Endocrinol. 1994;141(2):295–9.View ArticleGoogle Scholar
- Chardard D, Dournon C. Sex reversal by aromatase inhibitor treatment in the newt Pleurodeles waltl. J Exp Zool. 1999;283(1):43–50.View ArticleGoogle Scholar
- Olmstead AW, Kosian PA, Korte JJ, Holcombe GW, Woodis KK, Degitz SJ. Sex reversal of the amphibian, Xenopus tropicalis, following larval exposure to an aromatase inhibitor. Aquat Toxicol. 2009;91(2):143–50.View ArticleGoogle Scholar
- Pieau C, Dorizzi M. Oestrogens and temperature-dependent sex determination in reptiles: all is in the gonads. J Endocrinol. 2004;181(3):367–77.View ArticleGoogle Scholar
- Barske LA, Capel B. Estrogen represses SOX9 during sex determination in the red-eared slider turtle Trachemys scripta. Dev Biol. 2010;341(1):305–14.View ArticleGoogle Scholar
- Shupnik MA. Crosstalk between steroid receptors and the c-Src-receptor tyrosine kinase pathways: implications for cell proliferation. Oncogene. 2004;23:7979.View ArticleGoogle Scholar
- Waldbieser GC, Wolters WR. SHORT COMMUNICATION: definition of the USDA103 strain of channel catfish (Ictalurus punctatus). Anim Genet. 2007;38(2):180–3.View ArticleGoogle Scholar
- Cheryl A, Goudie BDR, Simco BA, Davis KB. Feminization of channel catfish by oral administration of steroid sex hormones. Trans Am Fish Soc. 1983;112(5):3.Google Scholar
- Waldbieser GC, Bosworth BG. A standardized microsatellite marker panel for parentage and kinship analyses in channel catfish, Ictalurus punctatus. Anim Genet. 2013;44(4):476–9.View ArticleGoogle Scholar
- Dunham RA, Lambert DM, Argue BJ, Ligeon C, Yant DR, Liu ZJ. Comparison of manual stripping and pen spawning for production of channel catfish × blue catfish hybrids and aquarium spawning of channel catfish. N Am J Aquac. 2000;62(4):260–5.View ArticleGoogle Scholar
- Su BF, Perera DA, Zohar Y, Abraham E, Stubblefield J, Fobes M, Beam R, Argue B, Ligeon C, Padi J, et al. Relative effectiveness of carp pituitary extract, luteininzing hormone releasing hormone analog (LHRHa) injections and LHRHa implants for producing hybrid catfish fry. Aquaculture. 2013;372:133–6.View ArticleGoogle Scholar
- Koren S, Walenz BP, Berlin K, Miller JR, Bergman NH, Phillippy AM. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017;27(5):722–36.View ArticleGoogle Scholar
- Li H. Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences. Bioinformatics. 2016;32(14):2103–10.View ArticleGoogle Scholar
- Vaser R, Sovic I, Nagarajan N, Sikic M. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 2017;27(5):737–46.View ArticleGoogle Scholar
- Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, Clum A, Copeland A, Huddleston J, Eichler EE, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10(6):563.View ArticleGoogle Scholar
- Tamazian G, Dobrynin P, Krasheninnikova K, Komissarov A, Koepfli KP, O'Brien SJ. Chromosomer: a reference-based genome arrangement tool for producing draft chromosome sequences. GigaScience. 2016;5(1):38.View ArticleGoogle Scholar
- Simao FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3210–2.View ArticleGoogle Scholar
- Broman KW, Wu H, Sen S, Churchill GA. R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003;19(7):889–90.View ArticleGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.Google Scholar
- Stanke M, Steinkamp R, Waack S, Morgenstern B. AUGUSTUS: a web server for gene finding in eukaryotes. Nucleic Acids Res. 2004;32(Web Server issue):W309–12.View ArticleGoogle Scholar
- Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23(9):1061–7.View ArticleGoogle Scholar
- Wang Y, Li J, Paterson AH. MCScanX-transposed: detecting transposed gene duplications based on multiple colinearity scans. Bioinformatics. 2013;29(11):1458–60.View ArticleGoogle Scholar
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL. Versatile and open software for comparing large genomes. Genome Biol. 2004;5(2):R12.View ArticleGoogle Scholar
- Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.View ArticleGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.View ArticleGoogle Scholar
- Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.View ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.View ArticleGoogle Scholar
- Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11(3):R25.View ArticleGoogle Scholar
- Xu Z, Chen J, Li X, Ge J, Pan J, Xu X. Identification and characterization of microRNAs in channel catfish (Ictalurus punctatus) by using Solexa sequencing technology. PLoS One. 2013;8(1):e54174.View ArticleGoogle Scholar
- Miranda KC, Huynh T, Tay Y, Ang YS, Tam WL, Thomson AM, Lim B, Rigoutsos I. A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell. 2006;126(6):1203–17.Google Scholar