Discovery of Transacting Long Noncoding RNAs That Regulate Smooth Muscle Cell Phenotype.

BACKGROUND: Smooth muscle cells (SMC), the major cell type in atherosclerotic plaques, are vital in coronary artery diseases (CADs). SMC phenotypic transition, which leads to the formation of various cell types in atherosclerotic plaques, is regulated by a network of genetic and epigenetic mechanisms and governs the risk of disease. The involvement of long noncoding RNAs (lncRNAs) has been increasingly identified in cardiovascular disease. However, SMC lncRNAs have not been comprehensively characterized, and their regulatory role in SMC state transition remains unknown. METHODS: A discovery pipeline was constructed and applied to deeply strand-specific RNA sequencing from perturbed human coronary artery SMC with different disease-related stimuli, to allow for the detection of novel lncRNAs. The functional relevance of a select few novel lncRNAs were verified in vitro. RESULTS: We identified 4579 known and 13 655 de novo lncRNAs in human coronary artery SMC. Consistent with previous long noncoding RNA studies, these lncRNAs overall have fewer exons, are shorter in length than protein-coding genes (pcGenes), and have relatively low expression level. Genomic location of these long noncoding RNA is disproportionately enriched near CAD-related TFs (transcription factors), genetic loci, and gene regulators of SMC identity, suggesting the importance of their function in disease. Two de novo lncRNAs, ZIPPOR (ZEB-interacting suppressor) and TNS1-AS2 (TNS1-antisense 2), were identified by our screen. Combining transcriptional data and in silico modeling along with in vitro validation, we identified CAD gene ZEB2 as a target through which these lncRNAs exert their function in SMC phenotypic transition. CONCLUSIONS: Expression of a large and diverse set of lncRNAs in human coronary artery SMC are highly dynamic in response to CAD-related stimuli. The dynamic changes in expression of these lncRNAs correspond to alterations in transcriptional programs that are relevant to CAD, suggesting a critical role for lncRNAs in SMC phenotypic transition and human atherosclerotic disease.

ZEB2 Shapes the Epigenetic Landscape of Atherosclerosis.

BACKGROUND: Smooth muscle cells (SMCs) transition into a number of different phenotypes during atherosclerosis, including those that resemble fibroblasts and chondrocytes, and make up the majority of cells in the atherosclerotic plaque. To better understand the epigenetic and transcriptional mechanisms that mediate these cell state changes, and how they relate to risk for coronary artery disease (CAD), we have investigated the causality and function of transcription factors at genome-wide associated loci. METHODS: We used CRISPR-Cas 9 genome and epigenome editing to identify the causal gene and cells for a complex CAD genome-wide association study signal at 2q22.3. Single-cell epigenetic and transcriptomic profiling in murine models and human coronary artery smooth muscle cells were used to understand the cellular and molecular mechanism by which this CAD risk gene exerts its function. RESULTS: CRISPR-Cas 9 genome and epigenome editing showed that the complex CAD genetic signals within a genomic region at 2q22.3 lie within smooth muscle long-distance enhancers for ZEB2, a transcription factor extensively studied in the context of epithelial mesenchymal transition in development of cancer. Zeb2 regulates SMC phenotypic transition through chromatin remodeling that obviates accessibility and disrupts both Notch and transforming growth factor ß signaling, thus altering the epigenetic trajectory of SMC transitions. SMC-specific loss of Zeb2 resulted in an inability of transitioning SMCs to turn off contractile programing and take on a fibroblast-like phenotype, but accelerated the formation of chondromyocytes, mirroring features of high-risk atherosclerotic plaques in human coronary arteries. CONCLUSIONS: These studies identify ZEB2 as a new CAD genome-wide association study gene that affects features of plaque vulnerability through direct effects on the epigenome, providing a new therapeutic approach to target vascular disease.

Identification of an Epigenetically Marked Locus within the Sex Determination Region of Channel Catfish.

Channel catfish has an XY sex determination system. However, the X and Y chromosomes harbor an identical gene content of 950 genes each. In this study, we conducted comparative analyses of methylome and transcriptome of genetic males and genetic females before gonadal differentiation to provide insights into the mechanisms of sex determination. Differentially methylated CpG sites (DMCs) were predominantly identified on the sex chromosome, most notably within the sex determination region (SDR), although the overall methylation profiles across the entire genome were similar between genetic males and females. The drastic differences in methylation were located within the SDR at nucleotide position 14.0-20.3 Mb of the sex chromosome, making this region an epigenetically marked locus within the sex determination region. Most of the differentially methylated CpG sites were hypermethylated in females and hypomethylated in males, suggesting potential involvement of methylation modification in sex determination in channel catfish. Along with the differential methylation in the SDR, a number of differentially expressed genes within the SDR were also identified between genetic males and females, making them potential candidate genes for sex determination and differentiation in channel catfish.

GWAS analysis using interspecific backcross progenies reveals superior blue catfish alleles responsible for strong resistance against enteric septicemia of catfish.

Infectious diseases pose significant threats to the catfish industry. Enteric septicemia of catfish (ESC) caused by Edwardsiella ictaluri is the most devastating disease for catfish aquaculture, causing huge economic losses annually. Channel catfish and blue catfish exhibit great contrast in resistance against ESC, with channel catfish being highly susceptible and blue catfish being highly resistant. As such, the interspecific backcross progenies provide an ideal system for the identification of quantitative trait locus (QTL). We previously reported one significant QTL on linkage group (LG) 1 using the third-generation backcrosses, but the number of founders used to make the second- and third-generation backcross progenies was very small. Although the third-generation backcross progenies provided a greater power for fine mapping than the first-generation backcrosses, some major QTL for disease resistance may have been missing due to the small numbers of founders used to produce the higher generation backcrosses. In this study, we performed a genome-wide association study using first-generation backcrosses with the catfish 690 K SNP arrays to identify additional ESC disease resistance QTL, especially those at the species level. Two genomic regions on LG1 and LG23 were determined to be significantly associated with ESC resistance as revealed by a mixed linear model and family-based association test. Examination of the resistance alleles indicated their origin from blue catfish, indicating that at least two major disease resistance loci exist among blue catfish populations. Upon further validation, markers linked with major ESC disease resistance QTL should be useful for marker-assisted introgression, allowing development of highly ESC resistant breeds of catfish.

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

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

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

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

A cell and transcriptome atlas of the human arterial vasculature.

Vascular beds show different propensities for different vascular pathologies, yet mechanisms explaining these fundamental differences remain unknown. We sought to build a transcriptomic, cellular, and spatial atlas of human arterial cells across multiple different arterial segments to understand this phenomenon. We found significant cell type-specific segmental heterogeneity. Determinants of arterial identity are predominantly encoded in fibroblasts and smooth muscle cells, and their differentially expressed genes are particularly enriched for vascular disease-associated loci and genes. Adventitial fibroblast-specific heterogeneity in gene expression coincides with numerous vascular disease risk genes, suggesting a previously unrecognized role for this cell type in disease risk. Adult arterial cells from different segments cluster not by anatomical proximity but by embryonic origin, with differentially regulated genes heavily influenced by developmental master regulators. Non-coding transcriptomes across arterial cells contain extensive variation in lnc-RNAs expressed in cell type- and segment-specific patterns, rivaling heterogeneity in protein coding transcriptomes, and show enrichment for non-coding genetic signals for vascular diseases.

Molecular mechanisms of coronary artery disease risk at the PDGFD locus.

Platelet derived growth factor (PDGF) signaling has been extensively studied in the context of vascular disease, but the genetics of this pathway remain to be established. Genome wide association studies (GWAS) for coronary artery disease (CAD) have identified a risk locus at 11q22.3, and we have verified with fine mapping approaches that the regulatory variant rs2019090 and PDGFD represent the functional variant and putative functional gene. Further, FOXC1/C2 transcription factor (TF) binding at rs2019090 was found to promote PDGFD transcription through the CAD promoting allele. Employing a constitutive Pdgfd knockout allele along with SMC lineage tracing in a male atherosclerosis mouse model we mapped single cell transcriptomic, cell state, and lesion anatomical changes associated with gene loss. These studies revealed that Pdgfd promotes expansion, migration, and transition of SMC lineage cells to the chondromyocyte phenotype and vascular calcification. This is in contrast to protective CAD genes TCF21, ZEB2, and SMAD3 which we have shown to promote the fibroblast-like cell transition or perturb the pattern or extent of transition to the chondromyocyte phenotype. Further, Pdgfd expressing fibroblasts and pericytes exhibited greater expression of chemokines and leukocyte adhesion molecules, consistent with observed increased macrophage recruitment to the plaque. Despite these changes there was no effect of Pdgfd deletion on SMC contribution to the fibrous cap or overall lesion burden. These findings suggest that PDGFD mediates CAD risk through promoting SMC expansion and migration, in conjunction with deleterious phenotypic changes, and through promoting an inflammatory response that is primarily focused in the adventitia where it contributes to leukocyte trafficking to the diseased vessel wall.

Molecular mechanisms of coronary artery disease risk at the PDGFD locus.

Genome wide association studies for coronary artery disease (CAD) have identified a risk locus at 11q22.3. Here, we verify with mechanistic studies that rs2019090 and PDGFD represent the functional variant and gene at this locus. Further, FOXC1/C2 transcription factor binding at rs2019090 is shown to promote PDGFD transcription through the CAD promoting allele. With single cell transcriptomic and histology studies with Pdgfd knockdown in an SMC lineage tracing male atherosclerosis mouse model we find that Pdgfd promotes expansion, migration, and transition of SMC lineage cells to the chondromyocyte phenotype. Pdgfd also increases adventitial fibroblast and pericyte expression of chemokines and leukocyte adhesion molecules, which is linked to plaque macrophage recruitment. Despite these changes there is no effect of Pdgfd deletion on overall plaque burden. These findings suggest that PDGFD mediates CAD risk by promoting deleterious phenotypic changes in SMC, along with an inflammatory response that is primarily focused in the adventitia.

Genomic imprinting-like monoallelic paternal expression determines sex of channel catfish.

The X and Y chromosomes of channel catfish have the same gene contents. Here, we report allelic hypermethylation of the X chromosome within the sex determination region (SDR). Accordingly, the X-borne hydin-1 gene was silenced, whereas the Y-borne hydin-1 gene was expressed, making monoallelic expression of hydin-1 responsible for sex determination, much like genomic imprinting. Treatment with a methylation inhibitor, 5-aza-dC, erased the epigenetic marks within the SDR and caused sex reversal of genetic females into phenotypic males. After the treatment, hydin-1 and six other genes related to cell cycle control and proliferative growth were up-regulated, while three genes related to female sex differentiation were down-regulated in genetic females, providing additional support for epigenetic sex determination in catfish. This mechanism of sex determination provides insights into the plasticity of genetic sex determination in lower vertebrates and its connection with temperature sex determination where DNA methylation is broadly involved.

Allelically and Differentially Expressed Genes After Infection of Edwardsiella ictaluri in Channel Catfish as Determined by Bulk Segregant RNA-Seq.

Identification of genetic markers associated with resistance against enteric septicemia of catfish (ESC) is of great interest for genetic enhancement programs of catfish. In the present study, bulk segregant RNA-Seq analysis was applied to determine differentially expressed genes and alleles after ESC infection. Here we report three genomic regions on LG1, LG12, and LG26, containing significant single-nucleotide polymorphisms (SNPs). These genomic regions aligned well with quantitative trait loci (QTL) previously identified. Within the QTL regions, eleven genes were found to be differentially regulated between phenotypic bulks. Importantly, the QTL on linkage group 1 (LG1) were found to be expressed in the liver, whereas the QTL on LG12 and LG26 were expressed in the intestine, suggesting multiple mechanisms of ESC resistance. It is apparent that apolipoproteins may be important for ESC resistance as the QTL on LG1 included the 14-kDa apolipoprotein genes that are both allelically expressed and differentially expressed between the resistant and susceptible bulks. Traf2 and NCK-interacting protein kinase (TNIK) were found in the QTL on LG12, and it was downregulated in resistant fish, suggesting the importance of NCK downregulation in ESC resistance, as previously reported. In addition, we observed divergent gene expression patterns between the liver and intestine after infection. Immune/inflammatory-related processes were overrepresented from liver DEGs, while those DEGs identified from intestine were enriched for proteolysis and wounding processes. Taken together, the BSR-Seq analysis presented here advanced the knowledge of ESC resistance, providing information of not only positions of QTL but also genes and their differential expression between resistant and susceptible fish, making it one step closer to the identification of the causal genes for ESC resistance.

Gene Editing of the Catfish Gonadotropin-Releasing Hormone Gene and Hormone Therapy to Control the Reproduction in Channel Catfish, Ictalurus punctatus.

Transcription activator-like effector nuclease (TALEN) plasmids targeting the channel catfish gonadotropin-releasing hormone (cfGnRH) gene were delivered into fertilized eggs with double electroporation to sterilize channel catfish (Ictalurus punctatus). Targeted cfGnRH fish were sequenced and base deletion, substitution, and insertion were detected. The gene mutagenesis was achieved in 52.9% of P1 fish. P1 mutants (individuals with human-induced sequence changes at the cfGnRH locus) had lower spawning rates (20.0−50.0%) when there was no hormone therapy compared to the control pairs (66.7%) as well as having lower average egg hatch rates (2.0% versus 32.3−74.3%) except for one cfGnRH mutated female that had a 66.0% hatch rate. After low fertility was observed in 2016, application of luteinizing hormone-releasing hormone analog (LHRHa) hormone therapy resulted in good spawning and hatch rates for mutants in 2017, which were not significantly different from the controls (p > 0.05). No exogenous DNA fragments were detected in the genome of mutant P1 fish, indicating no integration of the plasmids. No obvious effects on other economically important traits were observed after the knockout of the reproductive gene in the P1 fish. Growth rates, survival, and appearance between mutant and control individuals were not different. While complete knock-out of reproductive output was not achieved, as these were mosaic P1 brood stock, gene editing of channel catfish for the reproductive confinement of gene-engineered, domestic, and invasive fish to prevent gene flow into the natural environment appears promising.

GWAS Analysis Indicated Importance of NF-κB Signaling Pathway in Host Resistance Against Motile Aeromonas Septicemia Disease in Catfish.

Motile Aeromonas septicemia (MAS) disease caused by a bacterial pathogen, Aeromonas hydrophila, is an emerging but severe disease of catfish. Genetic enhancement of disease resistance is considered to be effective to control the disease. To provide an insight into the genomic basis of MAS disease resistance, in this study, we conducted a genome-wide association study (GWAS) to identify quantitative trait loci (QTL). A total of 1820 interspecific backcross catfish of 7 families were challenged with A. hydrophila, and 382 phenotypic extremes were selected for genotyping with the catfish 690 K SNP arrays. Three QTL on linkage group (LG) 2, 26 and 29 were identified to be significantly associated with MAS resistance. Within these regions, a total of 24 genes had known functions in immunity, 10 of which were involved in NF-κB signaling pathway, suggesting the importance of NF-κB signaling pathway in MAS resistance. In addition, three suggestively significant QTL were identified on LG 11, 17, and 20. The limited numbers of QTL involved in MAS resistance suggests that marker-assisted selection may be a viable approach for catfish breeding.

The annotation of repetitive elements in the genome of channel catfish (Ictalurus punctatus).

Channel catfish (Ictalurus punctatus) is a highly adaptive species and has been used as a research model for comparative immunology, physiology, and toxicology among ectothermic vertebrates. It is also economically important for aquaculture. As such, its reference genome was generated and annotated with protein coding genes. However, the repetitive elements in the catfish genome are less well understood. In this study, over 417.8 Megabase (MB) of repetitive elements were identified and characterized in the channel catfish genome. Among them, the DNA/TcMar-Tc1 transposons are the most abundant type, making up ~20% of the total repetitive elements, followed by the microsatellites (14%). The prevalence of repetitive elements, especially the mobile elements, may have provided a driving force for the evolution of the catfish genome. A number of catfish-specific repetitive elements were identified including the previously reported Xba elements whose divergence rate was relatively low, slower than that in untranslated regions of genes but faster than the protein coding sequences, suggesting its evolutionary restrictions.

Transcriptome analysis reveals enrichment of genes associated with auditory system in swimbladder of channel catfish.

In aquatic organisms, hearing is an important sense for acoustic communications and detection of sound-emitting predators and prey. Channel catfish is a dominant aquaculture species in the United States. As channel catfish can hear sounds of relatively high frequency, it serves as a good model for study auditory mechanisms. In catfishes, Weberian ossicles connect the swimbladder to the inner ear to transfer the forced vibrations and improve hearing ability. In this study, we examined the transcriptional profiles of channel catfish swimbladder and other four tissues (gill, liver, skin, and intestine). We identified a total of 1777 genes that exhibited preferential expression pattern in swimbladder of channel catfish. Based on Gene Ontology enrichment analysis, many of swimbladder-enriched genes were categorized into sensory perception of sound, auditory behavior, response to auditory stimulus, or detection of mechanical stimulus involved in sensory perception of sound, such as coch, kcnq4, sptbn1, sptbn4, dnm1, ush2a, and col11a1. Six signaling pathways associated with hearing (Glutamatergic synapse, GABAergic synapse pathways, Axon guidance, cAMP signaling pathway, Ionotropic glutamate receptor pathway, and Metabotropic glutamate receptor group III pathway) were over-represented in KEGG and PANTHER databases. Protein interaction prediction revealed an interactive relationship among the swimbladder-enriched genes and genes involved in sensory perception of sound. This study identified a set of genes and signaling pathways associated with auditory system in the swimbladder of channel catfish and provide resources for further study on the biological and physiological roles in catfish swimbladder.

A Review of Molecular Responses of Catfish to Bacterial Diseases and Abiotic Stresses.

Catfish is one of the major aquaculture species in the United States. However, the catfish industry is threatened by several bacterial diseases such as enteric septicemia of catfish (ESC), columnaris disease and Aeromonas disease, as well as by abiotic stresses such as high temperature and low oxygen. Research has been conducted for several decades to understand the host responses to these diseases and abiotic stresses. With the development of sequencing technologies, and the application of genome-wide association studies in aquaculture species, significant progress has been made. This review article summarizes recent progress in understanding the molecular responses of catfish after bacterial infection and stress challenges, and in understanding of genomic and genetic basis for disease resistance and stress tolerance.

JAK and STAT members in channel catfish: Identification, phylogenetic analysis and expression profiling after Edwardsiella ictaluri infection.

The Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling pathway is one of the main pleiotropic cascades used to transmit information from extracellular receptors to the nucleus, which results in DNA transcription and expression of genes involved in immunity, proliferation, differentiation, migration, apoptosis, and cell survival. Members of JAK family and STAT family have been extensively studied in different mammalian species because of their important roles in innate and adaptive immune responses. However, they have not been systematically studied among teleost fish species. In this study, five JAK family members and eight STAT family members were identified and characterized from channel catfish. Phylogenetic analysis was conducted to properly annotate these genes. Syntenic analysis was also conducted to establish orthology, and confirm the results from phylogenetic analysis. Compared to mammals, more members of the JAK and STAT family were identified in channel catfish genome. Expression of JAK and STAT family members was detected in healthy catfish tissues, but was induced in gill, liver, and intestine after bacterial challenge. Notably, the significant upregulation of STAT1b gene in catfish liver, gill and intestine after Edwardsiella ictaluri infection supported the notion that high STAT1 expression are involved in defense against pathogens. Collectively, the increased expression of JAK and STAT members in tested tissues suggested their crucial function in defending the host against pathogen invasion.

Genome-Wide Association Study Reveals Multiple Novel QTL Associated with Low Oxygen Tolerance in Hybrid Catfish.

Hypoxic condition is common in aquaculture, leading to major economic losses. Genetic analysis of hypoxia tolerance, therefore, is not only scientifically significant, but also economically important. Catfish is generally regarded as being highly tolerant to low dissolved oxygen, but variations exist among various populations, strains, and species. In this study, we conducted a genome-wide association study (GWAS) using the catfish 250 K SNP array to identify quantitative trait locus (QTL) associated with tolerance to low dissolved oxygen in the channel catfish × blue catfish interspecific system. Four linkage groups (LG2, LG4, LG23, and LG29) were found to be associated with low oxygen tolerance in hybrid catfish. Multiple significant SNPs were found to be physically linked in genomic regions containing significant QTL for low oxygen tolerance on LG2 and LG23, and in those regions containing suggestively significant QTL on LG2, LG4, LG23, and LG29, suggesting that the physically linked SNPs were genuinely segregating and related with low oxygen tolerance. Analysis of genes within the associated genomic regions suggested that many of these genes were involved in VEGF, MAPK, mTOR, PI3K-Akt, P53-mediated apoptosis, and DNA damage checkpoint pathways. Comparative analysis indicated that most of the QTL at the species level, as analyzed by using the interspecific system, did not overlap with those identified from six strains of channel catfish, confirming the complexity of the genetic architecture of hypoxia tolerance in catfish.