Identification of a novel Gig2 gene family specific to non-amniote vertebrates

Origin of an antifreeze protein gene in response to Cenozoic climate change

Antifreeze proteins (AFPs) inhibit ice growth within fish and protect them from freezing in icy seawater. Alanine-rich, alpha-helical AFPs (type I) have independently (convergently) evolved in four branches of fishes, one of which is a subsection of the righteye flounders. The origin of this gene family has been elucidated by sequencing two loci from a starry flounder, Platichthys stellatus, collected off Vancouver Island, British Columbia. The first locus had two alleles that demonstrated the plasticity of the AFP gene family, one encoding 33 AFPs and the other allele only four. In the closely related Pacific halibut, this locus encodes multiple Gig2 (antiviral) proteins, but in the starry flounder, the Gig2 genes were found at a second locus due to a lineage-specific duplication event. An ancestral Gig2 gave rise to a 3-kDa "skin" AFP isoform, encoding three Ala-rich 11-a.a. repeats, that is expressed in skin and other peripheral tissues. Subsequent gene duplications, followed by internal duplications of the 11 a.a. repeat and the gain of a signal sequence, gave rise to circulating AFP isoforms. One of these, the "hyperactive" 32-kDa Maxi likely underwent a contraction to a shorter 3.3-kDa "liver" isoform. Present day starry flounders found in Pacific Rim coastal waters from California to Alaska show a positive correlation between latitude and AFP gene dosage, with the shorter allele being more prevalent at lower latitudes. This study conclusively demonstrates that the flounder AFP arose from the Gig2 gene, so it is evolutionarily unrelated to the three other classes of type I AFPs from non-flounders. Additionally, this gene arose and underwent amplification coincident with the onset of ocean cooling during the Cenozoic ice ages.

Pattern of New Gene Origination in a Special Fish Lineage, the Flatfishes

Origination of new genes are of inherent interest of evolutionary geneticists for decades, but few studies have addressed the general pattern in a fish lineage. Using our recent released whole genome data of flatfishes, which evolved one of the most specialized body plans in vertebrates, we identified 1541 (6.9% of the starry flounder genes) flatfish-lineage-specific genes. The origination pattern of these flatfish new genes is largely similar to those observed in other vertebrates, as shown by the proportion of DNA-mediated duplication (1317; 85.5%), RNA-mediated duplication (retrogenes; 96; 6.2%), and de novo-origination (128; 8.3%). The emergence rate of species-specific genes is 32.1 per Mya and the whole average level rate for the flatfish-lineage-specific genes is 20.9 per Mya. A large proportion (31.4%) of these new genes have been subjected to selection, in contrast to the 4.0% in primates, while the old genes remain quite similar (66.4% vs. 65.0%). In addition, most of these new genes (70.8%) are found to be expressed, indicating their functionality. This study not only presents one example of systematic new gene identification in a teleost taxon based on comprehensive phylogenomic data, but also shows that new genes may play roles in body planning.

A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination

Fetal mammalian testes secrete Anti-Müllerian hormone (Amh), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes polycystic ovarian syndrome. Zebrafish (Danio rerio) lacks Müllerian ducts and the Amh receptor gene amhr2 but, curiously, retains amh To discover the roles of Amh in the absence of Müllerian ducts and the ancestral receptor gene, we made amh null alleles in zebrafish. Results showed that normal amh prevents female-biased sex ratios. Adult male amh mutants had enormous testes, half of which contained immature oocytes, demonstrating that Amh regulates male germ cell accumulation and inhibits oocyte development or survival. Mutant males formed sperm ducts and some produced a few offspring. Young female mutants laid a few fertile eggs, so they also had functional sex ducts. Older amh mutants accumulated nonvitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-sequencing data partitioned juveniles at 21 days postfertilization (dpf) into two groups that each contained mutant and wild-type fish. Group21-1 upregulated ovary genes compared to Group21-2, which were likely developing as males. By 35 dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly underexpressed granulosa and theca genes, and testes underexpressed Leydig cell genes. These results show that ancestral Amh functions included development of the gonadal soma in ovaries and testes and regulation of gamete proliferation and maturation. A major gap in our understanding is the identity of the gene encoding a zebrafish Amh receptor; we show here that the loss of amhr2 is associated with the breakpoint of a chromosome rearrangement shared among cyprinid fishes.

The Molecular Evolution and Functional Divergence of Lamprey Programmed Cell Death Genes

The programmed cell death (PDCD) family plays a significant role in the regulation of cell survival and apoptotic cell death. However, the evolution, distribution and role of the PDCD family in lampreys have not been revealed. Thus, we identified the PDCD gene family in the lamprey genome and classified the genes into five subfamilies based on orthologs of the genes, conserved synteny, functional domains, phylogenetic tree, and conserved motifs. The distribution of the lamprey PDCD family and the immune response of the PDCD family in lampreys stimulated by different pathogens were also demonstrated. In addition, we investigated the molecular function of lamprey PDCD2, PDCD5, and PDCD10. Our studies showed that the recombinant lamprey PDCD5 protein and transfection of the L-PDCD5 gene induced cell apoptosis, upregulated the expression of the associated X protein (BAX) and TP53 and downregulated the expression of B cell lymphoma 2 (BCL-2) independent of Caspase 3. In contrast, lamprey PDCD10 suppressed apoptosis in response to cis-diaminedichloro-platinum (II) stimuli. Our phylogenetic and functional data not only provide a better understanding of the evolution of lamprey PDCD genes but also reveal the conservation of PDCD genes in apoptosis. Overall, our results provide a novel perspective on lamprey immune regulation mediated by the PDCD family.

Non-coding RNAs Function as Immune Regulators in Teleost Fish

Non-coding RNAs (ncRNAs) are functional RNA molecules that are transcribed from DNA but not translated into proteins. ncRNAs function as key regulators of gene expression and chromatin modification. Recently, the functional role of ncRNAs in teleost fish has been extensively studied. Teleost fish are a highly diverse group among the vertebrate lineage. Fish are also important in terms of aquatic ecosystem, food production and human life, being the source of animal proteins worldwide and models of biomedical research. However, teleost fish always suffer from the invasion of infectious pathogens including viruses and bacteria, which has resulted in a tremendous economic loss to the fishing industry worldwide. Emerging evidence suggests that ncRNAs, especially miRNAs and lncRNAs, may serve as important regulators in cytokine and chemokine signaling, antigen presentation, complement and coagulation cascades, and T cell response in teleost fish. In this review, we summarize current knowledge and understanding of the roles of both miRNAs and lncRNAs in immune regulation in teleost fish. Molecular mechanism insights into the function of ncRNAs in fish immune response may contribute to the development of potential biomarkers and therapeutic targets for the prevention and treatment of fish diseases.

The role of type I interferons in innate and adaptive immunity against viruses in Atlantic salmon

Type I IFNs (IFN-I) are cytokines, which play a crucial role in innate and adaptive immunity against viruses of vertebrates. In essence, IFN-I are induced and secreted upon host cell recognition of viral nucleic acids and protect other cells against infection by inducing antiviral proteins. Atlantic salmon possesses an extraordinary repertoire of IFN-I genes encompassing at least six different classes (IFNa, IFNb, IFNc, IFNd, IFNe and IFNf) most of which are encoded by several genes. This review describes recent research on the functions of salmon IFNa, IFNb, IFNc and IFNd. As in mammals, expression of different salmon IFN-I in response to virus infection is dependent on their promoters, properties of the virus and the cell's expression of nucleic acid receptors and interferon regulatory factors (IRFs). While IFNa mainly display local antiviral activity, IFNb and IFNc show systemic antiviral activity. In addition, salmon appears to possess several IFN-I receptors, which show selectivity in binding different IFN-I. This complexity in IFN-I and receptors allows for a large variation in functions of the salmon IFN-I. Studies with intramuscular injection of IFN expression plasmids have recently provided surprising results, which may be of relevance for application of IFN-I in prophylaxis against virus infection. Firstly, injection of IFNc plasmid protected salmon presmolts against virus infection for at least 10 weeks. Secondly, IFN plasmids showed potent adjuvant activity when injected together with a DNA vaccine against infectious salmon anemia virus (ISAV).

De novo annotation of the immune-enriched transcriptome provides insights into immune system genes of Chinese sturgeon (Acipenser sinensis)

Chinese sturgeon (Acipenser sinensis), one of the oldest extant actinopterygian fishes with very high evolutionary, economical and conservation interest, is considered to be one of the critically endangered aquatic animals in China. Up to date, the immune system of this species remains largely undetermined with little sequence information publicly available. Herein, the first comprehensive transcriptome of immune tissues for Chinese sturgeon was characterized using Illumina deep sequencing. Over 67 million high-quality reads were generated and de novo assembled into the final set of 91,739 unique sequences. The annotation pipeline revealed that 25,871 unigenes were successfully annotated in the public databases, of which only 2002 had significant match to the existing sequences for the genus Acipenser. Overall 22,827 unigenes were categorized into 52 GO terms, 12,742 were classified into 26 KOG categories, and 4968 were assigned to 339 KEGG pathways. A more detailed annotation search showed the presence of a notable representation of immune-related genes, which suggests that this non-teleost actinopterygian fish harbors the same intermediates as in the well known immune pathways from mammals and teleosts, such as pattern recognition receptor (PRR) signaling pathway, JAK-STAT signaling pathway, complement and coagulation pathway, T-cell receptor (TCR) and B-cell receptor (BCR) signaling pathways. Additional genetic marker discovery led to the retrieval of 20,056 simple sequence repeats (SSRs) and 327,140 single nucleotide polymorphisms (SNPs). This immune-enriched transcriptome of Chinese sturgeon represents a rich resource that adds to the currently nascent field of chondrostean fish immunogenetics and furthers the conservation and management of this valuable fish.

Insights into the antiviral immunity against grass carp (Ctenopharyngodon idella) reovirus (GCRV) in grass carp

Global fish production from aquaculture has rapidly grown over the past decades, and grass carp shares the largest portion. However, hemorrhagic disease caused by grass carp reovirus (GCRV) results in tremendous loss of grass carp (Ctenopharyngodon idella) industry. During the past years, development of molecular biology and cellular biology technologies has promoted significant advances in the understanding of the pathogen and the immune system. Immunoprophylaxis based on stimulation of the immune system of fish has also got some achievements. In this review, authors summarize the recent progresses in basic researches on GCRV; viral nucleic acid sensors, high-mobility group box proteins (HMGBs); pattern recognition receptors (PRRs), Toll-like receptors (TLRs) and retinoic acid inducible gene I- (RIG-I-) like receptors (RLRs); antiviral immune responses induced by PRRs-mediated signaling cascades of type I interferon (IFN-I) and IFN-stimulated genes (ISGs) activation. The present review also notices the potential applications of molecule genetic markers. Additionally, authors discuss the current preventive and therapeutic strategies (vaccines, RNAi, and prevention medicine) and highlight the importance of innate immunity in long term control for grass carp hemorrhagic disease.

Promoter analysis and transcriptional regulation of a Gig2 gene in grass carp (Ctenopharyngodon idella)

Grass carp reovirus (GCRV)-induced gene 2 (Gig2) is recognized as a new antiviral factor involved in response to viral infection. However, little is known about the mechanisms behind the transcriptional regulation of Gig2 when infected by virus. In this study, the upstream promoter region of grass carp (Ctenopharyngodon idella) Gig2 gene (CiGig2) was identified by homology cloning strategy. CiGig2 promoter sequence was found to be 859 bp in length and contained three scattered IFN-stimulated response elements (ISRE). In addition, some grass carp IRFs (CiIRF1, CiIRF2 and CiIRF3) ORF sequences were subcloned into the expression plasmids pET-32a and expressed in Escherichia coli BL21, then the expressed proteins were purified by affinity chromatography with the Ni-NTA His-Bind Resin. Gel mobility shift assay was employed to screen the transcriptional regulatory factor for CiGig2. The results revealed that the recombinant polypeptides of CiIRF1, CiIRF2 and CiIRF3 bound to CiGig2 promoter with high affinity; indicating that IRF1, IRF2 and IRF3 could be the potential transcriptional regulatory factors for Gig2. Subsequently, CiGig2 promoter sequence was cloned into pGL3-Basic vector and the ORFs of CiIRF1, CiIRF2 and CiIRF3 were cloned into the expression plasmids pcDNA3.1 (+). Then, pGL3-CiGig2 promoter sequence and pcDNA3.1-CiIRFs were co-transfected into C. idella kidney (CIK) cells. The in vivo effects of CiIRFs on CiGig2 promoter were measured by dual-luciferase assays in the transfected CIK cells. Our results showed that the roles of CiIRFs were diversified in regulating CiGig2 transcription, e.g., CiIRF3 played a positive role in during this process; on the contrary CiIRF1 worked as a suppressor; however the effect of CiIRF2 on CiGig2 transcription was not obvious. For further study the roles of the three ISREs in CiGig2 transcription, we cloned three mutant CiGig2 promoters called ISRE1mut-luc (deleted ISRE1), ISRE2mut-luc (deleted ISRE2) and ISRE3mut-luc (deleted ISRE3), respectively. In vitro, gel mobility shift assays showed that all three mutant promoters also were combined with CiIRFs. CIK cells were co-transfected with CiGig2 promoter mutants (ISRE1mut-luc, ISRE2mut-luc or ISRE3mut-luc, respectively) and pcDNA3.1-IRFs. The results suggested that different ISRE played the diverse roles. ISRE2 is more important than ISRE1 and ISRE3 to the transcription of CiGig2 induced by CiIRF1. ISRE1 and ISRE3 are important to the transcription of CiGig2 induced by CiIRF2 and CiIRF3.

Zebrafish liver (ZFL) cells are able to mount an anti-viral response after stimulation with Poly (I:C)

The zebrafish (Danio rerio) is a widely used model species for biomedical research and is also starting to be a model for aquaculture research. The ZFL cell line, established from zebrafish liver, has been mostly used in toxicological and ecotoxicological studies. However, no studies have previously characterised this cell line in regard to its immunological response. The aim of this work was to study the gene expression response of the ZFL cell line after incubation with different prototypical immune stimuli, such as lipopolysaccharide (LPS), peptidoglycan (PGN), zymosan, and with a special focus on the dsRNA Poly (I:C). Using PCR, microarrays, and confocal microscopy we have explored the response of the ZFL cells against Poly (I:C). This study shows that the ZFL is able to uptake very efficiently the Poly (I:C) and mount a strong anti-viral response. We can conclude that ZFL could be used not only in toxicological studies, but also in studying anti-viral responses in zebrafish.

Genome-wide identification, phylogeny, and gonadal expression of fox genes in Nile tilapia, Oreochromis niloticus

The fox genes play important roles in various biological processes, including sexual development. In the present study, we isolated 65 fox genes, belonging to 18 subfamilies named A-R, from Nile tilapia through genome-wide screening. Twenty-four of them have two or three (foxm1) copies. Furthermore, 16, 25, 68, and 45 fox members were isolated from nematodes, protochordates, teleosts, and tetrapods, respectively. Phylogenetic analyses indicated fox gene family had undergone three expansions parallel to the three rounds of genome duplication during evolution. We also analyzed the clustered fox genes and found that apparent linkage duplication existed in teleosts, which further supported fish-specific genome duplication hypothesis. In addition, species- and lineage-specific duplication is another reason for fox gene family expansion. Based on the four pairs of XX and XY gonadal transcriptome data from four critical developmental stages, we analyzed the expression profile of all fox genes and identified sexually dimorphic fox genes at each stage. All fox genes were detected in gonads, with 15 of them at the background expression level (total read per kb per million reads, RPKM < 10), 29 at moderate expression level (10 < total RPKM < 100), and 21 at high expression level (total RPKM > 100). There are 27, 24, 28, and 9 sexually dimorphic fox genes at 5, 30, 90, and 180 days after hatching (dah), respectively. foxq1a, foxf1, foxr1, and foxr1 were identified as the most differentially expressed genes at each stage. foxl2 was characterized as XX-dominant gene, while foxd5, foxi3, foxn3, foxj1a, foxj3b, and foxo6b were characterized as XY-dominant genes. qPCR and in situ hybridization of foxh1 and foxj1a were performed to confirm the expression profiles and to validate the transcriptome data. Our results suggest that fox genes might play important roles in sex determination and gonadal development in teleosts.

Gig1 and Gig2 homologs (CiGig1 and CiGig2) from grass carp (Ctenopharyngodon idella) display good antiviral activities in an IFN-independent pathway

The virus-induced genes, Gig1 and Gig2, were identified first as IFN-stimulated genes (ISGs) from CAB cells. Previous studies suggested that Gig protein may have some potential antiviral functions. In this study, we cloned and identified the full-length cDNA sequences of Gig1 and Gig2 homologs (designated as CiGig1 and CiGig2, respectively) from grass carp (Ctenopharyngodon idella). The complete cDNA sequences of Gig1 and Gig2 contain 1231 bp and 690 bp, encoding for a 194 amino acid protein and a 158 amino acid protein, respectively. Their structure characteristics of CiGig1 and CiGig2 are highly similar to the corresponding homologues in crucian carp. The tissue-specific expressions of CiGig1 and CiGig2 in liver, spleen, kidney, intestine, gill and heart were significantly up-regulated following GCHV challenge. The results indicated that CiGig1 and CiGig2 may be involved in the antiviral immune responses in cells. To better understand the antiviral functions of CiGig1 and CiGig2 in vivo, CiGig1 or CiGig2 ORF cDNA were inserted into the plasmid pcDNA3.1, respectively. Subsequently, the recombinant plasmids were transfected into C. idellus kidney (CIK) cells. The over-expressions of CiGig1 and CiGig2 were observed in the CIK cells after treatment with GCHV. Cells with pcDNA3.1-CiGig1 or pcDNA-CiGig2 exhibited a relatively higher survival rate of (70.84% or 69.24%) than non-transfection (22.16%) and mock-vehicle controls (24.38%) following the virus infection. Our data showed that both CiGig1 and CiGig2 could exert antiviral effects effectively in vivo. Cycloheximide blocking protein synthesis demonstrated that both CiGig1 and CiGig2 mRNA expression could be induced by GCHV rather than by recombinant grass carp IFN (rCiIFN) directly, suggesting that CiGig1 and CiGig2 may not be IFN-stimulated genes since they display their antivirus activities in an IFN-independent pathway.