Genome-wide identification, phylogeny and expressional profile of the Sox gene family in channel catfish (Ictalurus punctatus)

Expression and localization of HPG axis-related genes in Carassius auratus with different ploidy

Introduction

In the Dongting water system, the Carassius auratus (Crucian carp) complex is characterized by the coexistence of diploid forms (2n=100, 2nCC) and polyploidy forms. The diploid (2nCC) and triploid C.auratus (3n=150, 3nCC) had the same fertility levels, reaching sexual maturity at one year.

Methods

The nucleotide sequence, gene expression, methylation, and immunofluorescence of the gonadotropin releasing hormone 2(Gnrh2), Gonadotropin hormone beta(Gthβ), and Gonadotropin-releasing hormone receptor(Gthr) genes pivotal genes of the hypothalamic-pituitary-gonadal (HPG) axis were analyzed.

Results

The analysis results indicated that Gnrh2, follicle-stimulating hormone receptor(Fshr), and Lethal hybrid rescue(Lhr) genes increased the copy number and distinct structural differentiation in 3nCC compared to that in 2nCC. The transcript levels of HPG axis genes in 3nCC were higher than 2nCC (P<0.05), which could promote the production and secretion of sex steroid hormones conducive to the gonadal development of 3nCC. Meanwhile, the DNA methylation levels in the promoter regions of the HPG axis genes were lower in 3nCC than in 2nCC. These results suggested that methylation of the promoter region had a potential regulatory effect on gene expression after triploidization. Immunofluorescence showed that the localization of the Fshβ, Lhβ, and Fshr genes between 3nCC and 2nCC remained unchanged, ensuring the normal expression of these genes at the corresponding sites after triploidization.

Discussion

Relevant research results provide cell and molecular biology evidence for normal reproductive activities such as gonad development and gamete maturation in triploid C. auratus, and contribute to further understanding of the genetic basis for fertility restoration in triploid C. auratus.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

Molecular Evolution and Inheritance Pattern of Sox Gene Family among Bovidae

Sox genes are an evolutionarily conserved family of transcription factors that play important roles in cellular differentiation and numerous complex developmental processes. In vertebrates, Sox proteins are required for cell fate decisions, morphogenesis, and the control of self-renewal in embryonic and adult stem cells. The Sox gene family has been well-studied in multiple species including humans but there has been scanty or no research into Bovidae. In this study, we conducted a detailed evolutionary analysis of this gene family in Bovidae, including their physicochemical properties, biological functions, and patterns of inheritance. We performed a genome-wide cataloguing procedure to explore the Sox gene family using multiple bioinformatics tools. Our analysis revealed a significant inheritance pattern including conserved motifs that are critical to the ability of Sox proteins to interact with the regulatory regions of target genes and orchestrate multiple developmental and physiological processes. Importantly, we report an important conserved motif, EFDQYL/ELDQYL, found in the SoxE and SoxF groups but not in other Sox groups. Further analysis revealed that this motif sequence accounts for the binding and transactivation potential of Sox proteins. The degree of protein–protein interaction showed significant interactions among Sox genes and related genes implicated in embryonic development and the regulation of cell differentiation. We conclude that the Sox gene family uniquely evolved in Bovidae, with a few exhibiting important motifs that drive several developmental and physiological processes.

Genome-wide identification and expression analysis of mitogen-activated protein kinase (MAPK) genes in response to salinity stress in channel catfish (Ictalurus punctatus)

The mitogen-activated protein kinase (MAPK) gene family has been systematically described in several fish species, but less so in channel catfish (Ictalurus punctatus), which is an important global aquaculture species. In this study, 16 MAPK genes were identified in the channel catfish genome and classified into three subfamilies based on phylogenetic analysis, including six extracellular signal regulated kinase (ERK) genes, six p38-MAPK genes and four C-Jun N-terminal kinase (JNK) genes. All MAPK genes were distributed unevenly across 10 chromosomes, of which three (IpMAPK8, IpMAPK12 and IpMAPK14) underwent teleost-specific whole genome duplication during evolution. Gene expression profiles in channel catfish during salinity stress were analysed using transcriptome sequencing and qRT-PCR (quantitative reverse transcription PCR). Results from reads per kilobase million (RPKM) analysis showed IpMAPK13, IpMAPK14a and IpMAPK14b genes were differentially expressed when compared with other genes between treatment and control groups. Furthermore, three of these genes were validated by qRT-PCR, of which IpMAPK14a expression levels were significantly upregulated in treatment groups (high and low salinity) when compared with the control group, with the highest expression levels in the low salinity group (P < 0.05). Therefore, IpMAPK14a may have important response roles to salinity stress in channel catfish.

sox19 regulates ovarian steroidogenesis in common carp

In teleost, ovarian steroidogenesis governed by the neuroendocrine system is also regulated by several transcription factors of gonadal origin. Investigating the synchronized interactions between the transcriptional and the hormonal factors is vital to comprehend the mechanisms that lead to gonadal differentiation. This study signifies the role of sry-related box (sox) 19 in ovarian steroidogenesis regulation of the common carp, Cyprinus carpio. Analysis of tissue distribution displayed higher sox19 expression in brain and ovary, and gonadal ontogeny showed higher expression of sox19 at 80 days post hatch (dph). Higher sox19 mRNA expression during spawning and increase of sox19 post human chorionic gonadotropin induction substantiate gonadotropin dependency. Estradiol-17β treatment but not 17α-methyl-di-hydroxy-testosterone to 50 dph common carp for inducing mono-sex, elevated sox19 expression substantially. Sox19 protein was observed in granulosa cells of the follicular layer in common carp ovary. Higher sox19 expression was detected in isolated granulosa and theca cells, in vitro. Transient gene silencing with sox19-siRNA caused downregulation of various ovary-related genes including those specific to activator protein-1 factors, fibroblast growth factors, wnt-signaling, steroidogenic genes along with certain transcription factors. Serum 17α, 20β-dihydroxy-4-pregnen-3-one and estradiol-17β reduced significantly post sox19 silencing, in vivo. Concomitantly, a decrease in aromatase activity was detected post sox19-siRNA treatment, in vivo. This study demonstrates the impact of sox19 in the regulation of common carp ovarian growth and steroidogenesis.

Labile sex chromosomes in the Australian freshwater fish family Percichthyidae

Sex-specific ecology has management implications, but rapid sex-chromosome turnover in fishes hinders sex-marker development for monomorphic species. We used annotated genomes and reduced-representation sequencing data for two Australian percichthyids, Macquarie perch Macquaria australasica and golden perch M. ambigua, and whole genome resequencing for 50 Macquarie perch of each sex, to identify sex-linked loci and develop an affordable sexing assay. In silico pool-seq tests of 1,492,004 Macquarie perch SNPs revealed that a 275-kb scaffold was enriched for gametologous loci. Within this scaffold, 22 loci were sex-linked in a predominantly XY system, with females being homozygous for the X-linked allele at all 22, and males having the Y-linked allele at >7. Seven XY-gametologous loci (all males, but no females, are heterozygous or homozygous for the male-specific allele) were within a 146-bp region. A PCR-RFLP sexing assay targeting one Y-linked SNP, tested in 66 known-sex Macquarie perch and two of each sex of three confamilial species, plus amplicon sequencing of 400 bp encompassing the 146-bp region, revealed that the few sex-linked positions differ between species and between Macquarie perch populations. This indicates sex-chromosome lability in Percichthyidae, supported by nonhomologous scaffolds containing sex-linked loci for Macquarie- and golden perches. The present resources facilitate genomic research in Percichthyidae, including formulation of hypotheses about candidate genes of interest such as transcription factor SOX1b that occurs in the 275-kb scaffold ~38 kb downstream of the 146-bp region containing seven XY-gametologous loci. Sex-linked markers will be useful for determining genetic sex in some populations and studying sex chromosome turnover.

Genome wide identification, phylogeny, and synteny analysis of sox gene family in common carp (Cyprinus carpio)

•

27 SOX (high-mobility group HMG-box) genes were identified in the C. carp genome.

•

SOX genes ranging from 3496 (SOX6) to 924bp (SOX17b) which coded with putative protein series from 307 to 509 amino acids.

•

Gene ontology revealed SOX proteins maximum involvement is in metabolic process 49.796 %.

•

Chromosomal location and synteny analysis display all SOX gene are located on different chromosomes.

Common carp (Cyprinus carpio) is a commercial fish species valuable for nutritious components and plays a vital role in human healthy nutrition. The SOX (SRY-related genes systematically characterized by a high-mobility group HMG-box) encoded important gene regulatory proteins, a family of transcription factors found in a broad range of animal taxa and extensively known for its contribution in multiple developmental processes including contribution in sex determination across phyla. In our current study, we initially accomplished a genome-wide analysis to report the SOX gene family in common carp fish based on available genomic sequences of zebrafish retrieved from gene repository databases, we focused on the global identification of the Sox gene family in Common carp among wide range of vertebrates and teleosts based on bioinformatics tools and techniques and explore the evolutionary relationships. In our results, a total of 27 SOX (high-mobility group HMG-box) domain genes were identified in the C. carp genome. The full length sequences of SOX genes ranging from 3496 (SOX6) to 924bp (SOX17b) which coded with putative proteins series from 307 to 509 amino acids and all gene having exon number expect SOX9 and SOX13. All the SOX proteins contained at least one conserved DNA-binding HMG-box domain and two (SOX7 and SOX18) were found C terminal. The Gene ontology revealed SOX proteins maximum involvement is in metabolic process 49.796 %, average in biological regulation 45.188 %, biosynthetic process (19.992 %), regulation of cellular process 39.68, 45.508 % organic substance metabolic process, multicellular organismal process 23.23 %,developmental process 21.74 %, system development 16.59 %, gene expression 16.05 % and 14.337 % of RNA metabolic process. Chromosomal location and syntanic analysis show all SOX gene are located on different chromosomes and apparently does not fallow the unique pattern. The maximum linkage of chromosome is (2) on Unplaced Scaffold region. Finally, our results provide important genomic suggestion for upcoming studies of biochemical, physiological, and phylogenetic understanding on SOX genes among teleost.

Systematic identification and expression analysis of the Sox gene family in spotted sea bass (Lateolabrax maculatus)

The Sox gene family encodes a set of transcription factors characterized by a conserved Sry-related high mobility group (HMG)-box domain, which performs a series of essential biological functions in diverse tissues and developmental processes. In this study, the Sox gene family was systematically characterized in spotted sea bass (Lateolabrax maculatus). A total of 26 Sox genes were identified and classified into eight subfamilies, namely, SoxB1, SoxB2, SoxC, SoxD, SoxE, SoxF, SoxH and SoxK. The phylogenetic relationship, exon-intron and domain structure analyses supported their annotation and classification. Comparison of gene copy numbers and chromosome locations among different species indicated that except tandem duplicated paralogs of Sox17/Sox32, duplicated Sox genes in spotted sea bass were generated from teleost-specific whole genome duplication during evolution. In addition, qRT-PCR was performed to detect the expression profiles of Sox genes during development and adulthood. The results showed that the expression of 16 out of 26 Sox genes was induced dramatically at different starting points after the multicellular stage, which is consistent with embryogenesis. At the early stage of sex differentiation, 9 Sox genes exhibited sexually dimorphic expression patterns, among which Sox3, Sox19 and Sox6b showed the most significant ovary-biased expression. Moreover, the distinct expression pattern of Sox genes was observed in different adult tissues. Our results provide a fundamental resource for further investigating the functions of Sox genes in embryonic processes, sex determination and differentiation as well as controlling the homeostasis of adult tissues in spotted sea bass.

Molecular cloning, characterization, and expression analysis of a sex-biased transcriptional factor sox9 gene of mud crab Scylla paramamosain

Sox9 gene, a crucial member of the Sox gene family, is present in various organisms and involved in many physiological processes, especially in sex determination and gonad development. In this study, we reported a sox9 gene (designated as Spsox9) from Scylla paramamosain through analyzing published gonad transcriptome data. Meanwhile, the accuracy was validated by PCR technology, and the 3' sequences were cloned with 3' RACE technology. The full-length cDNA of Spsox9 is 2843 bp, consisting of a 243 bp 5' UTR, an 1124 bp 3' UTR, and a 1476 bp ORF encoding 491 amino acids. Furthermore, to better understand its conservation among crustacean species, the sox9 gene ortholog was identified in several other crustaceans species with their published transcriptome data, respectively. All of the Sox9 proteins identified in the current study had the common feature of Sox proteins (HMG domain) and were highly conserved among analyzed crustacean species. In all examined tissues, the Spsox9 was mainly expressed in the gonad (testis and ovary), eyestalk, and cerebral ganglion. During embryo development, Spsox9 was highly expressed in 5 pairs of appendages, 7 pairs of appendages, and eye-pigment formation stage. During ovary development, the expression level of Spsox9 remained stable in the first 4 stages (O1-O4) and decreased in the tertiary vitellogenesis (O5) stage. During testis development, the expression level of Spsox9 was highest in the spermatid stage (T2) and was significantly different from that in the spermatocyte stage (T1) and mature sperm stage (T3) (p < 0.05). In addition, Spsox9 exhibited a sex-biased expression pattern in T1 and O1. These present results indicated that the Spsox9 gene might play crucial roles in the gonad and embryo development of mud crab.

Sox Gene Family Revealed Genetic Variations in Autotetraploid Carassius auratus

The Sox gene family encoded transcription factors that played key roles in developmental processes in vertebrates. To further understand the evolutionary fate of the Sox gene family in teleosts, the Sox genes were comprehensively characterized in fish of different ploidy levels, including blunt snout bream (2n = 48, Megalobrama amblycephala, BSB), goldfish (2n = 100, Carassius auratus red var., 2nRCC), and autotetraploid C. auratus (4n = 200, 4nRCC). The 4nRCC, which derived from the whole genome duplication (WGD) of 2nRCC, were obtained through the distant hybridization of 2nRCC (♀) × BSB (♂). Compared with the 26 Sox genes in zebrafish (2n = 50, Danio rerio), 26, 47, and 92 putative Sox genes were identified in the BSB, 2nRCC, and 4nRCC genomes, respectively, and classified into seven subfamilies (B1, B2, C, D, E, F, and K). Comparative analyses showed that 89.36% (42/47) of Sox genes were duplicated in 2nRCC compared with those in BSB, while 97.83% (90/92) of Sox genes were duplicated in 4nRCC compared with those in 2nRCC, meaning the Sox gene family had undergone an expansion in BSB, 2nRCC, and 4nRCC, respectively, following polyploidization events. In addition, potential gene loss, genetic variations, and paternal parent SNP locus insertion occurred during the polyploidization events. Our data provided new insights into the evolution of the Sox gene family in polyploid vertebrates after several rounds of WGD events.

Systematic Identification and Evolution Analysis of Sox Genes in Coturnix japonica Based on Comparative Genomics

Coturnix japonica (Japanese quail) has been extensively used as a model animal for biological studies. The Sox gene family, which was systematically characterized by a high-mobility group (HMG-box) in many animal species, encodes transcription factors that play central roles during multiple developmental processes. However, genome-wide investigations on the Sox gene family in birds are scarce. In the current study, we first performed a genome-wide study to explore the Sox gene family in galliform birds. Based on available genomic sequences retrieved from the NCBI database, we focused on the global identification of the Sox gene family in C. japonica and other species in Galliformes, and the evolutionary relationships of Sox genes. In our result, a total of 35 Sox genes in seven groups were identified in the C. japonica genome. Our results also revealed that dispersed gene duplications contributed the most to the expansion of the Sox gene family in Galliform birds. Evolutionary analyses indicated that Sox genes are an ancient gene family, and strong purifying selections played key roles in the evolution of CjSox genes of C. japonica. More interestingly, we observed that most Sox genes exhibited highly embryo-specific expression in both gonads. Our findings provided new insights into the molecular function and phylogeny of Sox gene family in birds.

A draft genome of the striped catfish, Pangasianodon hypophthalmus, for comparative analysis of genes relevant to development and a resource for aquaculture improvement

Background

The striped catfish, Pangasianodon hypophthalmus, is a freshwater and benthopelagic fish common in the Mekong River delta. Catfish constitute a valuable source of dietary protein. Therefore, they are cultured worldwide, and P. hypophthalmus is a food staple in the Mekong area. However, genetic information about the culture stock, is unavailable for breeding improvement, although genetics of the channel catfish, Ictalurus punctatus, has been reported. To acquire genome sequence data as a useful resource for marker-assisted breeding, we decoded a draft genome of P. hypophthalmus and performed comparative analyses.

Results

Using the Illumina platform, we obtained both nuclear and mitochondrial DNA sequences. Molecular phylogeny using the mitochondrial genome confirmed that P. hypophthalmus is a member of the family Pangasiidae and is nested within a clade including the families Cranoglanididae and Ictaluridae. The nuclear genome was estimated at approximately 700 Mb, assembled into 568 scaffolds with an N50 of 14.29 Mbp, and was estimated to contain ~ 28,600 protein-coding genes, comparable to those of channel catfish and zebrafish. Interestingly, zebrafish produce gadusol, but genes for biosynthesis of this sunscreen compound have been lost from catfish genomes. The differences in gene contents between these two catfishes were found in genes for vitamin D-binding protein and cytosolic phospholipase A₂, which have lost only in channel catfish. The Hox cluster in catfish genomes comprised seven paralogous groups, similar to that of zebrafish, and comparative analysis clarified catfish lineage-specific losses of A5a, B10a, and A11a. Genes for insulin-like growth factor (IGF) signaling were conserved between the two catfish genomes. In addition to identification of MHC class I and sex determination-related gene loci, the hypothetical chromosomes by comparison with the channel catfish demonstrated the usefulness of the striped catfish genome as a marker resource.

Conclusions

We developed genomic resources for the striped catfish. Possible conservation of genes for development and marker candidates were confirmed by comparing the assembled genome to that of a model fish, Danio rerio, and to channel catfish. Since the catfish genomic constituent resembles that of zebrafish, it is likely that zebrafish data for gene functions is applicable to striped catfish as well.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5079-x) contains supplementary material, which is available to authorized users.