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.

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.

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.

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 Drive to Driven Model of Mapping Intraspecific Interaction Networks.

Community ecology theory suggests that an individual's phenotype is determined by the phenotypes of its coexisting members to the extent at which this process can shape community evolution. Here, we develop a mapping theory to identify interaction quantitative trait loci (QTL) governing inter-individual dependence. We mathematically formulate the decision-making strategy of interacting individuals. We integrate these mathematical descriptors into a statistical procedure, enabling the joint characterization of how QTL drive the strengths of ecological interactions and how the genetic architecture of QTL is driven by ecological networks. In three fish full-sib mapping experiments, we identify a set of genome-wide QTL that control a range of societal behaviors, including mutualism, altruism, aggression, and antagonism, and find that these intraspecific interactions increase the genetic variation of body mass by about 50%. We showcase how the interaction QTL can be used as editors to reconstruct and engineer new social networks for ecological communities.

The allotetraploid origin and asymmetrical genome evolution of the common carp Cyprinus carpio.

Common carp (Cyprinus carpio) is an allotetraploid species derived from recent whole genome duplication and provides a model to study polyploid genome evolution in vertebrates. Here, we generate three chromosome-level reference genomes of C. carpio and compare to related diploid Cyprinid genomes. We identify a Barbinae lineage as potential diploid progenitor of C. carpio and then divide the allotetraploid genome into two subgenomes marked by a distinct genome similarity to the diploid progenitor. We estimate that the two diploid progenitors diverged around 23 Mya and merged around 12.4 Mya based on the divergence rates of homoeologous genes and transposable elements in two subgenomes. No extensive gene losses are observed in either subgenome. Instead, we find gene expression bias across surveyed tissues such that subgenome B is more dominant in homoeologous expression. CG methylation in promoter regions may play an important role in altering gene expression in allotetraploid C. carpio.

Patterns of Geographical and Potential Adaptive Divergence in the Genome of the Common Carp (Cyprinus carpio).

The common carp, Cyprinus carpio, is a cyprinid fish species cultured in Europe and Asia. It accounts for >70% of freshwater aquaculture production worldwide. We conducted a population genomics analysis on C. carpio using high-throughput SNP genotyping of 2,198 individuals from 14 populations worldwide to determine the genetic architecture of common carp populations and the genetic bases for environmental adaptation. Structure analyses including phylogeny and principal component analysis were also conducted, showing distinct geographical patterns in European and Asian populations. The linkage disequilibrium block average lengths of the 14 populations ranged from 3.94 kb to 36.67 kb. Genes within selective sweep regions were identified by genome scanning among the different populations, including gdf6a, bmpr1b, and opsin5. Gene Ontology and KEGG enrichment analyses revealed potential trait-related loci and genes associated with body shape, scaling patterns, and skin color. This population genomics analysis may provide valuable clues for future genome-assisted breeding of C. carpio.

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, Transcriptomic, and Epigenomic Features Differentiate Genes That Are Relevant for Muscular Polyunsaturated Fatty Acids in the Common Carp.

Polyunsaturated fatty acids (PUFAs) are a set of important nutrients that mainly include arachidonic acid (ARA4), docosahexaenoic acid (DHA), eicosapentaenoic acid (EPA), and α-linolenic acid (ALA). Recently, fish-derived PUFAs have been associated with cardiovascular health, fetal development, and improvement of brain functions. Studies have shown that fish muscular tissues are rich in PUFAs, which are influenced by various factors, including genetic variations, regulatory profiles, and methylation status of desaturase genes during fatty acid desaturation and elongation processes. However, the genetic mechanism and the pathways involved in fatty acid metabolism in fishes remain unclear. The overall aim of this study was to assess differences in gene expression responses among fishes with different fatty acid levels. To achieve this goal, we conducted genome-wide association analysis (GWAS) using a 250K SNP array in a population of 203 samples of common carp (Cyprinus carpio) and identified nine SNPs and 15 genes associated with muscular PUFA content. Then, RNA-Seq and whole genome bisulfite sequencing (WGBS) of different groups with high and low EPA, DHA, ARA4, and ALA contents in muscle, liver and brain tissues were conducted, resulting in 6,750 differentially expressed genes and 5,631 genes with differentially methylated promoters. Gene ontology and KEGG pathway enrichment analyses of RNA-Seq and WGBS results identified enriched pathways for fatty acid metabolism, which included the adipocytokine signaling pathway, ARA4 and linoleic acid metabolism pathway, and insulin signaling pathway. Integrated analysis indicated significant correlations between gene expression and methylation status among groups with high and low PUFA contents in muscular tissues. Taken together, these multi-level results uncovered candidate genes and pathways that are associated with fatty acid metabolism and paved the way for further genomic selection and carp breeding for PUFA traits.

Comparative transcriptome analysis reveals expression signatures of albino Russian sturgeon, Acipenseriformes gueldenstaedtii.

Albinism is a genetically inherited condition that is caused by a series of genetic abnormalities leading to a reduction in melanin production. Russian sturgeon is one of the most valuable freshwater fish species worldwide, and albino individuals have been found in fish farms. Due to its complicated genome and scarce genome-wide genetic resources, the underlying molecular basis of albinism in Russian sturgeon is unknown. In the present study, we first generated transcriptome profile of Acipenser gueldenstaedtii using pooled tissues, which provided reliable reference sequences for future molecular genetic studies. A total of 369,441 contigs were assembled, corresponding to 32,965 unique genes. A comparative analysis of the transcripts from the skin of albino and wildtype individuals was conducted afterwards. A total of 785 unique genes were differentially expressed, including the upregulation of 385 genes and the downregulation of 400 genes in albino individuals. The expression pattern of 16 selected differentially expressed genes was validated using qRT-PCR. Additional annotation, GO enrichment analysis and gene pathway analysis indicated that the melanogenesis pathway may be interrupted in albinism. Eight potential causative genes that were highly likely to be responsible for sturgeon albinism were identified, including Dct, Tyrp1b, Slc45a2, Ctns, Pmela, Pmelb, Cd63, and Bloc1s3, which were found to be significantly down-regulated in albino Russian sturgeon. Moreover, a sliding window analysis of the ratio of nonsynonymous to synonymous nucleotide substitution rates (Ka/Ks) ratios indicated that seven out of the eight genes underwent positive selection during evolution. Our results provide a valuable basis for understanding the molecular mechanism of albinism in fish species and will facilitate future genetic selection and breeding of sturgeon with market-favored traits in aquaculture.

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

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

Genome-wide characterization of Toll-like receptor gene family in common carp (Cyprinus carpio) and their involvement in host immune response to Aeromonas hydrophila infection.

The Toll-like receptor (TLR) gene family is a class of conserved pattern recognition receptors, which play an essential role in innate immunity providing efficient defense against invading microbial pathogens. Although TLRs have been extensively characterized in both invertebrates and vertebrates, a comprehensive analysis of TLRs in common carp is lacking. In the present study, we have conducted the first genome-wide systematic analysis of common carp (Cyprinus carpio) TLR genes. A set of 27 common carp TLR genes were identified and characterized. Sequence similarity analysis, functional domain prediction and phylogenetic analysis supported their annotation and orthologies. By examining the gene copy number of TLR genes across several vertebrates, gene duplications and losses were observed. The expression patterns of TLR genes were examined during early developmental stages and in various healthy tissues, and the results showed that TLR genes were ubiquitously expressed, indicating a likely role in maintaining homeostasis. Moreover, the differential expression of TLRs was examined after Aeromons hydrophila infection, and showed that most TLR genes were induced, with diverse patterns. TLR1, TLR4-2, TLR4-3, TLR22-2, TLR22-3 were significantly up-regulated at minimum one timepoint, whereas TLR2-1, TLR4-1, TLR7-1 and TLR7-2 were significantly down-regulated. Our results suggested that TLR genes play critical roles in the common carp immune response. Collectively, our findings provide fundamental genomic resources for future studies on fish disease management and disease-resistance selective breeding strategy development.

Draft genome of the Northern snakehead, Channa argus.

The Northern snakehead (Channa argus), a member of the Channidae family of the Perciformes, is an economically important freshwater fish native to East Asia. In North America, it has become notorious as an intentionally released invasive species. Its ability to breathe air with gills and migrate short distances over land makes it a good model for bimodal breath research. Therefore, recent research has focused on the identification of relevant candidate genes. Here, we performed whole genome sequencing of C. argus to construct its draft genome, aiming to offer useful information for further functional studies and identification of target genes related to its unusual facultative air breathing. Findings: We assembled the C. argus genome with a total of 140.3 Gb of raw reads, which were sequenced using the Illumina HiSeq2000 platform. The final draft genome assembly was approximately 615.3 Mb, with a contig N50 of 81.4 kb and scaffold N50 of 4.5 Mb. The identified repeat sequences account for 18.9% of the whole genome. The 19 877 protein-coding genes were predicted from the genome assembly, with an average of 10.5 exons per gene. Conclusion: We generated a high-quality draft genome of C. argus, which will provide a valuable genetic resource for further biomedical investigations of this economically important teleost fish.

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.

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.

Host-Pathogen Interactions in Measles Virus Replication and Anti-Viral Immunity.

The measles virus (MeV) is a contagious pathogenic RNA virus of the family Paramyxoviridae, genus Morbillivirus, that can cause serious symptoms and even fetal complications. Here, we summarize current molecular advances in MeV research, and emphasize the connection between host cells and MeV replication. Although measles has reemerged recently, the potential for its eradication is promising with significant progress in our understanding of the molecular mechanisms of its replication and host-pathogen interactions.

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 channel catfish genome sequence provides insights into the evolution of scale formation in teleosts.

Catfish represent 12% of teleost or 6.3% of all vertebrate species, and are of enormous economic value. Here we report a high-quality reference genome sequence of channel catfish (Ictalurus punctatus), the major aquaculture species in the US. The reference genome sequence was validated by genetic mapping of 54,000 SNPs, and annotated with 26,661 predicted protein-coding genes. Through comparative analysis of genomes and transcriptomes of scaled and scaleless fish and scale regeneration experiments, we address the genomic basis for the most striking physical characteristic of catfish, the evolutionary loss of scales and provide evidence that lack of secretory calcium-binding phosphoproteins accounts for the evolutionary loss of scales in catfish. The channel catfish reference genome sequence, along with two additional genome sequences and transcriptomes of scaled catfishes, provide crucial resources for evolutionary and biological studies. This work also demonstrates the power of comparative subtraction of candidate genes for traits of structural significance.

Comparative transcriptome analysis between aquatic and aerial breathing organs of Channa argus to reveal the genetic basis underlying bimodal respiration.

Aerial breathing in fish was an important adaption for successful survival in hypoxic water. All aerial breathing fish are bimodal breathers. It is intriguing that they can obtain oxygen from both air and water. However, the genetic basis underlying bimodal breathing has not been extensively studied. In this study, we performed next-generation sequencing on a bimodal breathing fish, the Northern snakehead, Channa argus, and generated a transcriptome profiling of C. argus. A total of 53,591 microsatellites and 26,378 SNPs were identified and classified. A Ka/Ks analysis of the unigenes indicated that 63 genes were under strong positive selection. Furthermore, the transcriptomes from the aquatic breathing organ (gill) and the aerial breathing organ (suprabranchial chamber) were sequenced and compared, and the results showed 1,966 genes up-regulated in the gill and 2,727 genes up-regulated in the suprabranchial chamber. A gene pathway analysis concluded that four functional categories were significant, of which angiogenesis and elastic fibre formation were up-regulated in the suprabranchial chamber, indicating that the aerial breathing organ may be more efficient for gas exchange due to its highly vascularized and elastic structure. In contrast, ion uptake and transport and acid-base balance were up-regulated in the gill, indicating that the aquatic breathing organ functions in ion homeostasis and acid-base balance, in addition to breathing. Understanding the genetic mechanism underlying bimodal breathing will shed light on the initiation and importance of aerial breathing in the evolution of vertebrates.

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.