Complete mitochondrial genome and its phylogenetic implications of Rhinogobio nasutus, an endemic species from the Yellow River.

BACKGROUND: Rhinogobio nasutus, an endemic species from the Yellow River, is listed under the second class of the National Key Protected Wildlife List in China due to its dramatically decreased population. Despite its important status, the mitochondrial genes and phylogenetic relationships of R. nasutus are unknown. METHODS AND RESULTS: The complete mitochondrial genome of R. nasutus was sequenced, assembled, and annotated for the first time. The mitochondrial genome was 16,609 bp in length, consisting of 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and 1 non-coding control region. The gene order in the mitochondrial genome of R. nasutus was identical to that of other Rhinogobios species. Analysis of synonymous and non-synonymous nucleotide substitutions showed that the Ka/Ks ratio in all tested protein-coding genes was less than 1, indicating that these genes were evolving under purifying selection. Further phylogenetic analysis showed that R. nasutus was first clustered with R. typus, then grouped with the other two Rhinogobio species, indicating the phylogenetically close relationship between R. nasutus and R. typus. CONCLUSIONS: This was the first genomic resource developed for R. nasutus, which could not only improve our understanding of its phylogenetic status, but also serve as a genomic tool for the development of genetic markers that will be used in conservation and evolutionary genetics studies.

Transcriptome signatures of the Pacific white shrimp Litopenaeus vannamei hepatopancreas in response to stress in biofloc culture systems.

Comparative transcriptome analysis via high throughput sequencing was applied to gain knowledge on the immune response in Litopenaeus vannamei reared in biofloc technology systems (BFT). Two types of carbon sources, namely, traditional carbon sources (molasses) and biodegradable polymers [hydroxybutyric acid-co-3-hydroxyvaleric acid (PHBV)] were used in BFT systems. Clear water systems without the addition of carbon sources were treated as the control. Water quality assays showed that the average concentrations of several stress factors, including nitrite, nitrate and TSS, were the highest in molasses-based BFT systems. After sequencing and comparing the transcriptome profiles of the L. vannamei hepatopancreas, 743 and 201 genes were significantly differentially expressed in molasses- and PHBV-based BFT systems, respectively. GO enrichment analysis, which was performed using the differentially expressed genes, revealed seven significantly over-represented GO terms in molasses-based BFT systems, including catabolic process, hydrolase activity, cellular localization, organic substance metabolic process, cellular metabolic process, establishment of localization and response to stress. The captured key genes were mainly involved in the pathways including cellular stress response, immune response and pathogen recognition. However, no GO terms were significantly over-represented in PHBV-based BFT systems compared with control. This study indicates that shrimp are subject to stress in BFT systems when molasses serves as the carbon source. Thus, PHBV may be a better alternative.

Genomic Basis of Adaptive Evolution: The Survival of Amur Ide (Leuciscus waleckii) in an Extremely Alkaline Environment.

The Amur ide (Leuciscus waleckii) is a cyprinid fish that is widely distributed in Northeast Asia. The Lake Dali Nur population inhabits one of the most extreme aquatic environments on Earth, with an alkalinity up to 50 mmol/L (pH 9.6), thus providing an exceptional model with which to characterize the mechanisms of genomic evolution underlying adaptation to extreme environments. Here, we developed the reference genome assembly for L. waleckii from Lake Dali Nur. Intriguingly, we identified unusual expanded long terminal repeats (LTRs) with higher nucleotide substitution rates than in many other teleosts, suggesting their more recent insertion into the L. waleckii genome. We also identified expansions in genes encoding egg coat proteins and natriuretic peptide receptors, possibly underlying the adaptation to extreme environmental stress. We further sequenced the genomes of 10 additional individuals from freshwater and 18 from Lake Dali Nur populations, and we detected a total of 7.6 million SNPs from both populations. In a genome scan and comparison of these two populations, we identified a set of genomic regions under selective sweeps that harbor genes involved in ion homoeostasis, acid-base regulation, unfolded protein response, reactive oxygen species elimination, and urea excretion. Our findings provide comprehensive insight into the genomic mechanisms of teleost fish that underlie their adaptation to extreme alkaline environments.

Inclusion Body Fusion of Human Parainfluenza Virus Type 3 Regulated by Acetylated α-Tubulin Enhances Viral Replication.

Viral inclusion bodies (IBs), or replication factories, are unique structures generated by viral proteins together with some cellular proteins as a platform for efficient viral replication, but little is known about the mechanism underlying IB formation and fusion. Our previous study demonstrated that the interaction between the nucleoprotein (N) and phosphoprotein (P) of human parainfluenza virus type 3 (HPIV3), an enveloped virus with great medical impact, can form IBs. In this study, we found that small IBs can fuse with each other to form large IBs that enhance viral replication. Furthermore, we found that acetylated α-tubulin interacts with the N-P complex and colocalizes with IBs of HPIV3 but does not interact with the N-P complex of human respiratory syncytial virus or vesicular stomatitis virus and does not colocalize with IBs of human respiratory syncytial virus. Most importantly, enhancement of α-tubulin acetylation using the pharmacological inhibitor trichostatin A (TSA), RNA interference (RNAi) knockdown of the deacetylase enzymes histone deacetylase 6 (HDAC6) and sirtuin 2 (SIRT2), or expression of α-tubulin acetyltransferase 1 (α-TAT1) resulted in the fusion of small IBs into large IBs and effective viral replication. In contrast, suppression of acetylation of α-tubulin by overexpressing HDAC6 and SIRT2 profoundly inhibited the fusion of small IBs and viral replication. Our findings offer previously unidentified mechanistic insights into the regulation of viral IB fusion by acetylated α-tubulin, which is critical for viral replication. IMPORTANCE: Inclusion bodies (IBs) are unique structures generated by viral proteins and some cellular proteins as a platform for efficient viral replication. Human parainfluenza virus type 3 (HPIV3) is a nonsegmented single-stranded RNA virus that mainly causes lower respiratory tract disease in infants and young children. However, no vaccines or antiviral drugs for HPIV3 are available. Therefore, understanding virus-host interactions and developing new antiviral strategies are increasingly important. Acetylation on lysine (K) 40 of α-tubulin is an evolutionarily conserved modification and plays an important role in many cellular processes, but its role in viral IB dynamics has not been fully explored. To our knowledge, our findings are the first to show that acetylated α-tubulin enhances viral replication by regulating HPIV3 IB fusion.

Genome wide identification of scavenger receptors class A in common carp (Cyprinus carpio) and their expression following Aeromonas hydrophila infection.

Scavenger receptors class A (SCARAs) is a subgroup of diverse families of pattern recognition receptors that bind a range of ligands, and play important roles in innate immune processes through pathogens detection, adhesion, endocytosis, and phagocytosis. However, most studies of SCARAs have focused on mammals, and much less is known of SCARAs in fish species. In this study, we identified 7 SCARAs across the common carp genome, which were classified into four subclasses according to comparative genomic analysis including sequence similarities analysis, gene structure and functional domain prediction. Further phylogenetic and syntenic analysis supported their annotation and orthologies. Through examining gene copy number of SCARA genes across several vertebrates, SCARA2, SCARA3 and SCARA4 were found have undergone gene duplication. The expression patterns of SCARAs in common carp were examined during early developmental stages, in healthy tissues, and after Aeromonas hydrophila infection. Most SCARA genes were ubiquitously expressed during common carp early developmental stages, and presented diverse patterns in various healthy tissues, with relatively high expression levels in spleen, liver, intestine, gill and brain, indicating their critical roles likely in maintaining homeostasis and host immune response activities. After A. hydrophila infection, most SCARA genes were up-regulated at 4 h post infection in mucosal tissue intestine, while generally up-regulated at 12 h post infection in spleen, suggesting a tissue-specific pattern of regulation. Taken together, all these results suggested that SCARA genes played important roles in host immune response to A. hydrophila infection in common carp, and provided important genomic resources for future studies on fish disease management.

Transcriptome signatures in common carp spleen in response to Aeromonas hydrophila infection.

The common carp is an important aquaculture species that is worldwide distributed. Nowadays, intensive rearing in aquaculture increases the susceptibility of fish to various pathogens such as Aeromonas hydrophila, which has caused severe damage to carp production. However, systematic analysis on the host response of common carp against A. hydrophila is less studied. In order to better understand the common carp immune response process against bacteria at the global gene expression level, we examined transcriptional profiles of the common carp spleen at three timepoints following experimental infection with A. hydrophila. A total of 545 million 125-bp paired end reads were generated, and all trimmed clean reads were mapped onto the common carp whole genome sequence. Comparison of the transcriptomes between the treatment and control group fish revealed 2900 unigenes with significantly differential expression, including 732, 936, 928 genes up-regulated, and 248, 475, 700 genes down-regulated at 4 h, 12 h, 24 h post infection respectively. The captured significantly differentially expressed genes are mainly involved in the pathways including junction/adhesion, pathogen recognition, cell surface receptor signaling, and immune system process/defense response. Our study will provide fundamental information on molecular mechanism underlying the immune response of teleost against bacterial infection and might suggest strategies for selection of resistant strains of common carp in aquaculture.

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.

A leucine residue in the C terminus of human parainfluenza virus type 3 matrix protein is essential for efficient virus-like particle and virion release.

UNLABELLED: Paramyxovirus particles, like other enveloped virus particles, are formed by budding from membranes of infected cells, and matrix (M) proteins are critical for this process. To identify the M protein important for this process, we have characterized the budding of the human parainfluenza virus type 3 (HPIV3) M protein. Our results showed that expression of the HPIV3 M protein alone is sufficient to initiate the release of virus-like particles (VLPs). Electron microscopy analysis confirmed that VLPs are morphologically similar to HPIV3 virions. We identified a leucine (L302) residue within the C terminus of the HPIV3 M protein that is critical for M protein-mediated VLP production by regulating the ubiquitination of the M protein. When L302 was mutated into A302, ubiquitination of M protein was defective, the release of VLPs was abolished, and the membrane binding and budding abilities of M protein were greatly weakened, but the ML302A mutant retained oligomerization activity and had a dominant negative effect on M protein-mediated VLP production. Furthermore, treatment with a proteasome inhibitor also inhibited M protein-mediated VLP production and viral budding. Finally, recombinant HPIV3 containing the M(L302A) mutant could not be rescued. These results suggest that L302 acts as a critical regulating signal for the ubiquitination of the HPIV3 M protein and virion release. IMPORTANCE: Human parainfluenza virus type 3 (HPIV3) is an enveloped virus with a nonsegmented negative-strand RNA genome. It can cause severe respiratory tract diseases, such as bronchiolitis, pneumonia, and croup in infants and young children. However, no valid antiviral therapy or vaccine is currently available. Thus, further elucidation of its assembly and budding will be helpful in the development of novel therapeutic approaches. Here, we show that a leucine residue (L302) located at the C terminus of the HPIV3 M protein is essential for efficient production of virus-like particles (VLPs). Furthermore, we found L302 regulated M protein-mediated VLP production via regulation of M protein ubiquitination. Recombinant HPIV3 containing the M(L302A) mutant is growth defective. These findings provide new insight into the critical role of M protein-mediated VLP production and virion release of a residue that does not belong to L domain and may advance our understanding of HPIV3 viral assembly and budding.

Transcriptional Profiling Reveals Differential Gene Expression of Amur Ide (Leuciscus waleckii) during Spawning Migration.

Amur ide (Leuciscus waleckii), an important aquaculture species, inhabits neutral freshwater but can tolerate high salinity or alkalinity. As an extreme example, the population in Dali Nor lake inhabits alkalized soda water permanently, and migrates from alkaline water to neutral freshwater to spawn. In this study, we performed comparative transcriptome profiling study on the livers of Amur ide to interrogate the expression differences between the population that permanently inhabit freshwater in Ganggeng Nor lake (FW) and the spawning population that recently migrated from alkaline water into freshwater (SM). A total of 637,234,880 reads were generated, resulting in 53,440 assembled contigs that were used as reference sequences. Comparisons of these transcriptome files revealed 444 unigenes with significant differential expression (p-value ≤ 0.01, fold-change ≥ 2), including 246 genes that were up-regulated in SM and 198 genes that were up-regulated in FW. The gene ontology (GO) enrichment analysis and KEGG pathway analysis indicated that the mTOR signaling pathway, Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling pathway, and oxidative phosphorylation were highly likely to affect physiological changes during spawning migration. Overall, this study demonstrates that transcriptome changes played a role in Amur ide spawning migration. These results provide a foundation for further analyses on the physiological and molecular mechanisms underlying Amur ide spawning migration.

Development and evaluation of the first high-throughput SNP array for common carp (Cyprinus carpio).

BACKGROUND: A large number of single nucleotide polymorphisms (SNPs) have been identified in common carp (Cyprinus carpio) but, as yet, no high-throughput genotyping platform is available for this species. C. carpio is an important aquaculture species that accounts for nearly 14% of freshwater aquaculture production worldwide. We have developed an array for C. carpio with 250,000 SNPs and evaluated its performance using samples from various strains of C. carpio. RESULTS: The SNPs used on the array were selected from two resources: the transcribed sequences from RNA-seq data of four strains of C. carpio, and the genome re-sequencing data of five strains of C. carpio. The 250,000 SNPs on the resulting array are distributed evenly across the reference C.carpio genome with an average spacing of 6.6 kb. To evaluate the SNP array, 1,072 C. carpio samples were collected and tested. Of the 250,000 SNPs on the array, 185,150 (74.06%) were found to be polymorphic sites. Genotyping accuracy was checked using genotyping data from a group of full-siblings and their parents, and over 99.8% of the qualified SNPs were found to be reliable. Analysis of the linkage disequilibrium on all samples and on three domestic C.carpio strains revealed that the latter had the longer haplotype blocks. We also evaluated our SNP array on 80 samples from eight species related to C. carpio, with from 53,526 to 71,984 polymorphic SNPs. An identity by state analysis divided all the samples into three clusters; most of the C. carpio strains formed the largest cluster. CONCLUSIONS: The Carp SNP array described here is the first high-throughput genotyping platform for C. carpio. Our evaluation of this array indicates that it will be valuable for farmed carp and for genetic and population biology studies in C. carpio and related species.

Comprehensive analysis of circRNA-miRNA-mRNA networks in the kidney of snakehead (Channa argus) response to Nocardia seriolae challenge.

Non-coding RNAs (ncRNAs) play vital roles in regulating the expression levels of genes that control essential biological functions, including immune response to bacterial infections in teleost. To dissect the roles of ncRNAs in the Channa argus (snakehead), a systematic analysis of the expression profiles of circRNA, miRNA and mRNA, as well as competing endogenous RNAs (ceRNA) regulatory networks in the kidney of snakehead following Nocardia seriolae infection were performed in the present study. A total of 111 differentially expressed circRNAs, 706 differentially expressed miRNAs, and 2548 differentially expressed mRNAs were identified in the N. seriolae infected snakehead. Based on these differently expressed RNAs, we identified 55 circRNA-mRNA pairs, 124 miRNA-mRNA pairs, and 35 circRNA-miRNA-mRNA regulatory networks, including dre-miR-103-CD302, dre-miR-27e-IGSF3, novel_circ_0005462/novel_403-IGKC, novel_circ_0001750/novel_circ_0002162-novel_477-OCLN, and novel_circ_0003847-novel_4-KCNAB3. In addition, luciferase reporter assay was employed to detect the target relationships of several circRNA-miRNA-mRNA pairs. Taken together, this study demonstrates that the genes associated with immunity and structures in the kidney of snakehead can be regulated by circRNAs and miRNAs at post-transcription levels, and provided theoretical guidance for ncRNAs studies for other teleost. However, further studies are still in great need to validate the regulatory mechanisms of ncRNAs in snakehead.

Contribution of elovl5a to Docosahexaenoic Acid (DHA) Synthesis at the Transcriptional Regulation Level in Common Carp, Cyprinus carpio.

Docosahexaenoic acid (DHA) is an essential nutrient for humans and plays a critical role in human development and health. Freshwater fish, such as the common carp (Cyprinus carpio), have a certain degree of DHA biosynthesis ability and could be a supplemental source of human DHA needs. The elongase of very-long-chain fatty acid 5 (Elovl5) is an important enzyme affecting polyunsaturated fatty acid (PUFA) biosynthesis. However, the function and regulatory mechanism of the elovl5 gene related to DHA synthesis in freshwater fish is not clear yet. Previous studies have found that there are two copies of the elovl5 gene, elovl5a and elovl5b, which have different functions. Our research group found significant DHA content differences among individuals in Yellow River carp (Cyprinus carpio var.), and four candidate genes were found to be related to DHA synthesis through screening. In this study, the expression level of elovl5a is decreased in the high-DHA group compared to the low-DHA group, which indicated the down-regulation of elovl5a in the DHA synthesis pathways of Yellow River carp. In addition, using a dual-luciferase reporter gene assay, we found that by targeting the 3'UTR region of elovl5a, miR-26a-5p could regulate DHA synthesis in common carp. After CRISPR/Cas9 disruption of elovl5a, the DHA content in the disrupted group was significantly higher than in the wildtype group; meanwhile, the expression level of elovl5a in the disrupted group was significantly reduced compared with the wildtype group. These results suggest that elovl5a may be down-regulating DHA synthesis in Yellow River carp. This study could provide useful information for future research on the genes and pathways that affect DHA synthesis.

RNA-Seq reveals expression signatures of genes involved in oxygen transport, protein synthesis, folding, and degradation in response to heat stress in catfish.

Temperature is one of the most prominent abiotic factors affecting ectotherms. Most fish species, as ectotherms, have extraordinary ability to deal with a wide range of temperature changes. While the molecular mechanism underlying temperature adaptation has long been of interest, it is still largely unexplored with fish. Understanding of the fundamental mechanisms conferring tolerance to temperature fluctuations is a topic of increasing interest as temperature may continue to rise as a result of global climate change. Catfish have a wide natural habitat and possess great plasticity in dealing with environmental variations in temperature. However, no studies have been conducted at the transcriptomic level to determine heat stress-induced gene expression. In the present study, we conducted an RNA-Seq analysis to identify heat stress-induced genes in catfish at the transcriptome level. Expression analysis identified a total of 2,260 differentially expressed genes with a cutoff of twofold change. qRT-PCR validation suggested the high reliability of the RNA-Seq results. Gene ontology, enrichment, and pathway analyses were conducted to gain insight into physiological and gene pathways. Specifically, genes involved in oxygen transport, protein folding and degradation, and metabolic process were highly induced, while general protein synthesis was dramatically repressed in response to the lethal temperature stress. This is the first RNA-Seq-based expression study in catfish in response to heat stress. The candidate genes identified should be valuable for further targeted studies on heat tolerance, thereby assisting the development of heat-tolerant catfish lines for aquaculture.

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.

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.

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.

Gene expression changes leading extreme alkaline tolerance in Amur ide (Leuciscus waleckii) inhabiting soda lake.

BACKGROUND: Amur ide (Leuciscus waleckii) is an economically and ecologically important cyprinid species in Northern Asia. The Dali Nor population living in the soda lake Dali Nor can adapt the extremely high alkalinity, providing us a valuable material to understand the adaptation mechanism against extreme environmental stress in teleost. RESULTS: In this study, we generated high-throughput RNA-Seq data from three tissues gill, liver and kidney of L. waleckii living in the soda lake Dali Nor and the fresh water lake Ganggeng Nor, then performed parallel comparisons of three tissues. Our results showed that out of assembled 64,603 transcript contigs, 28,391 contigs had been assigned with a known function, corresponding to 20,371 unique protein accessions. We found 477, 2,761 and 3,376 differentially expressed genes (DEGs) in the gill, kidney, and liver, respectively, of Dali Nor population compared to Ganggeng Nor population with FDR ≤ 0.01 and fold-change ≥ 2. Further analysis revealed that well-known functional categories of genes and signaling pathway, which are associated with stress response and extreme environment adaptation, have been significantly enriched, including the functional categories of "response to stimulus", "transferase activity", "transporter activity" and "oxidoreductase activity", and signaling pathways of "mTOR signaling", "EIF2 signaling", "superpathway of cholesterol biosynthesis". We also identified significantly DEGs encoding important modulators on stress adaptation and tolerance, including carbonic anhydrases, heat shock proteins, superoxide dismutase, glutathione S-transferases, aminopeptidase N, and aminotransferases. CONCLUSIONS: Overall, this study demonstrated that transcriptome changes in L. waleckii played a role in adaptation to complicated environmental stress in the highly alkalized Dali Nor lake. The results set a foundation for further analyses on alkaline-responsive candidate genes, which help us understand teleost adaptation under extreme environmental stress and ultimately benefit future breeding for alkaline-tolerant fish strains.

Efficient assembly and annotation of the transcriptome of catfish by RNA-Seq analysis of a doubled haploid homozygote.

BACKGROUND: Upon the completion of whole genome sequencing, thorough genome annotation that associates genome sequences with biological meanings is essential. Genome annotation depends on the availability of transcript information as well as orthology information. In teleost fish, genome annotation is seriously hindered by genome duplication. Because of gene duplications, one cannot establish orthologies simply by homology comparisons. Rather intense phylogenetic analysis or structural analysis of orthologies is required for the identification of genes. To conduct phylogenetic analysis and orthology analysis, full-length transcripts are essential. Generation of large numbers of full-length transcripts using traditional transcript sequencing is very difficult and extremely costly. RESULTS: In this work, we took advantage of a doubled haploid catfish, which has two sets of identical chromosomes and in theory there should be no allelic variations. As such, transcript sequences generated from next-generation sequencing can be favorably assembled into full-length transcripts. Deep sequencing of the doubled haploid channel catfish transcriptome was performed using Illumina HiSeq 2000 platform, yielding over 300 million high-quality trimmed reads totaling 27 Gbp. Assembly of these reads generated 370,798 non-redundant transcript-derived contigs. Functional annotation of the assembly allowed identification of 25,144 unique protein-encoding genes. A total of 2,659 unique genes were identified as putative duplicated genes in the catfish genome because the assembly of the corresponding transcripts harbored PSVs or MSVs (in the form of pseudo-SNPs in the assembly). Of the 25,144 contigs with unique protein hits, around 20,000 contigs matched 50% length of reference proteins, and over 14,000 transcripts were identified as full-length with complete open reading frames. The characterization of consensus sequences surrounding start codon and the stop codon confirmed the correct assembly of the full-length transcripts. CONCLUSIONS: The large set of transcripts assembled in this study is the most comprehensive set of genome resources ever developed from catfish, which will provide the much needed resources for functional genome research in catfish, serving as a reference transcriptome for genome annotation, analysis of gene duplication, gene family structures, and digital gene expression analysis. The putative set of duplicated genes provide a starting point for genome scale analysis of gene duplication in the catfish genome, and should be a valuable resource for comparative genome analysis, genome evolution, and genome function studies.

Genome-wide identification of endosialin family of C-type lectins in common carp (Cyprinus carpio) and their response following Aeromonas hydrophila infection.

The endosialin family is the group XIV of C-type lectin, regulating several processes involved in innate immunity and inflammation. Endosialin family genes have been extensively studied in human and mammals, however, rarely reported in teleost. In the present study, a set of 8 endosialin family genes was identified across the entire common carp genome. Functional domain and motif prediction and phylogenetic analysis supported their annotation and orthologies. Through examining gene copy number across several vertebrates, endosialin family genes were found have undergone gene duplication. Most of the endosialin family genes were ubiquitously expressed during common carp early developmental stages, and presented tissue-specific expression patterns in various healthy tissues, with relatively high expression in intestine, liver, gill, spleen and kidney, indicating their likely essential roles in maintaining homeostasis and host immune response. After Aeromonas hydrophila infection, gene thbd-1, thbd-2 and cd93-2 were significantly up-regulated at one or more timepoints in spleen and kidney, while gene cd248a-1, cd248a-2, cd248b-1, cd248b-2, and cd93-1 were significantly down-regulated. Taken together, all these results suggested that endosialin family genes were involved in host immune response to A. hydrophila infection in common carp, and provided fundamental genomic resources for better understanding the critical roles of endosialin family on the primary innate immune processes in teleost.

A pilot study for channel catfish whole genome sequencing and de novo assembly.

BACKGROUND: Recent advances in next-generation sequencing technologies have drastically increased throughput and significantly reduced sequencing costs. However, the average read lengths in next-generation sequencing technologies are short as compared with that of traditional Sanger sequencing. The short sequence reads pose great challenges for de novo sequence assembly. As a pilot project for whole genome sequencing of the catfish genome, here we attempt to determine the proper sequence coverage, the proper software for assembly, and various parameters used for the assembly of a BAC physical map contig spanning approximately a million of base pairs. RESULTS: A combination of low sequence coverage of 454 and Illumina sequencing appeared to provide effective assembly as reflected by a high N50 value. Using 454 sequencing alone, a sequencing depth of 18 X was sufficient to obtain the good quality assembly, whereas a 70 X Illumina appeared to be sufficient for a good quality assembly. Additional sequencing coverage after 18 X of 454 or after 70 X of Illumina sequencing does not provide significant improvement of the assembly. Considering the cost of sequencing, a 2 X 454 sequencing, when coupled to 70 X Illumina sequencing, provided an assembly of reasonably good quality. With several software tested, Newbler with a seed length of 16 and ABySS with a K-value of 60 appear to be appropriate for the assembly of 454 reads alone and Illumina paired-end reads alone, respectively. Using both 454 and Illumina paired-end reads, a hybrid assembly strategy using Newbler for initial 454 sequence assembly, Velvet for initial Illumina sequence assembly, followed by a second step assembly using MIRA provided the best assembly of the physical map contig, resulting in 193 contigs with a N50 value of 13,123 bp. CONCLUSIONS: A hybrid sequencing strategy using low sequencing depth of 454 and high sequencing depth of Illumina provided the good quality assembly with high N50 value and relatively low cost. A combination of Newbler, Velvet, and MIRA can be used to assemble the 454 sequence reads and the Illumina reads effectively. The assembled sequence can serve as a resource for comparative genome analysis. Additional long reads using the third generation sequencing platforms are needed to sequence through repetitive genome regions that should further enhance the sequence assembly.