Expansion by whole genome duplication and evolution of the sox gene family in teleost fish. PLoS One. 2017;12:e0180936. [PMID: 28738066 PMCID: PMC5524304 DOI: 10.1371/journal.pone.0180936]

Identification and expression analysis of Sox family genes in echinoderms

The Sox gene family, a collection of transcription factors widely distributed throughout the animal kingdom, plays a crucial role in numerous developmental processes. Echinoderms occupy a pivotal position in many research fields, such as neuroscience, sex determination and differentiation, and embryonic development. However, to date, no comprehensive study has been conducted to characterize and analyze Sox genes in echinoderms. In the present study, the evolution and expression of Sox family genes across 11 echinoderms were analyzed using bioinformatics methods. The results revealed a total of 70 Sox genes, with counts ranging from 5 to 8 across different echinoderms. Phylogenetic analysis revealed that the identified Sox genes could be categorized into seven distinct classes: the SoxB1 class, SoxB2 class, SoxC class, SoxD class, SoxE class, SoxF class and SoxH class. Notably, the SoxB1, SoxB2, and SoxF genes were ubiquitously present in all the echinoderms studied, which suggests that these genes may be conserved in echinoderms. The spatiotemporal expression patterns observed for Sox genes in the three echinoderms indicated that various Sox members perform distinct functional roles. Notably, SoxB1 is likely involved in echinoderm ovary development, while SoxH may play a crucial role in testis development in starfish and sea cucumber. In general, the present investigation provides a molecular foundation for exploring the Sox gene in echinoderms, providing a valuable resource for future phylogenetic and genomic studies.

Genetic and developmental divergence in the neural crest programme between cichlid fish species

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic programme in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to inter-specific morphological differences. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared to the teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes - particularly those controlled by sox10s - might be involved in generating morphological diversification between species and lays the groundwork for further investigations into mechanisms underpinning vertebrate NC diversification.

Genome-Wide Identification, Molecular Characterization, and Expression Analysis of the HSP70 and HSP90 Gene Families in Thamnaconus septentrionalis

Heat shock proteins (HSPs) are a class of highly conserved proteins that play an important role in biological responses to various environmental stresses. The mariculture of Thamnaconus septentrionalis, a burgeoning aquaculture species in China, frequently encounters stressors such as extreme temperatures, salinity variations, and elevated ammonia levels. However, systematic identification and analysis of the HSP70 and HSP90 gene families in T. septentrionalis remain unexplored. This study conducted the first genome-wide identification of 12 HSP70 and 4 HSP90 genes in T. septentrionalis, followed by a comprehensive analysis including phylogenetics, gene structure, conserved domains, chromosomal localization, and expression profiling. Expression analysis from RNA-seq data across various tissues and developmental stages revealed predominant expression in muscle, spleen, and liver, with the highest expression found during the tailbud stage, followed by the gastrula, neurula, and juvenile stages. Under abiotic stress, most HSP70 and HSP90 genes were upregulated in response to high temperature, high salinity, and low salinity, notably hspa5 during thermal stress, hspa14 in high salinity, and hsp90ab1 under low salinity conditions. Ammonia stress led to a predominance of downregulated HSP genes in the liver, particularly hspa2, while upregulation was observed in the gills, especially for hsp90b1. Quantitative real-time PCR analysis corroborated the expression levels under environmental stresses, validating their involvement in stress responses. This investigation provides insights into the molecular mechanisms of HSP70 and HSP90 in T. septentrionalis under stress, offering valuable information for future functional studies of HSPs in teleost evolution, optimizing aquaculture techniques, and developing stress-resistant strains.

gjSOX9 Cloning, Expression, and Comparison with gjSOXs Family Members in Gekko japonicus

SOX9 plays a crucial role in the male reproductive system, brain, and kidneys. In this study, we firstly analyzed the complete cDNA sequence and expression patterns for SOX9 from Gekko japonicus SOX9 (gjSOX9), carried out bioinformatic analyses of physiochemical properties, structure, and phylogenetic evolution, and compared these with other members of the gjSOX family. The results indicate that gjSOX9 cDNA comprises 1895 bp with a 1482 bp ORF encoding 494aa. gjSOX9 was not only expressed in various adult tissues but also exhibited a special spatiotemporal expression pattern in gonad tissues. gjSOX9 was predicted to be a hydrophilic nucleoprotein with a characteristic HMG-Box harboring a newly identified unique sequence, "YKYQPRRR", only present in SOXE members. Among the 20 SOX9 orthologs, gjSOX9 shares the closest genetic relationships with Eublepharis macularius SOX9, Sphacrodactylus townsendi SOX9, and Hemicordylus capensis SOX9. gjSOX9 and gjSOX10 possessed identical physicochemical properties and subcellular locations and were tightly clustered with gjSOX8 in the SOXE group. Sixteen gjSOX family members were divided into six groups: SOXB, C, D, E, F, and H with gjSOX8, 9, and 10 in SOXE among 150 SOX homologs. Collectively, the available data in this study not only facilitate a deep exploration of the functions and molecular regulation mechanisms of the gjSOX9 and gjSOX families in G. japonicus but also contribute to basic research regarding the origin and evolution of SOX9 homologs or even sex-determination mode in reptiles.

Hsp70 Gene Family in Sebastiscus marmoratus: The Genome-Wide Identification and Transcriptome Analysis under Thermal Stress

Heat shock protein 70 kDa (Hsp70) is a highly conserved heat stress protein that is important in biotic processes and responses to abiotic stress. Hsp70 genes may be important in Sebastiscus marmoratus, for it is a kind of nearshore reef fish, and habitat temperature change is more drastic during development. However, genome-wide identification and expression analysis in the Hsp70 gene family of S. marmoratus are still lacking. Here, a total of 15 Hsp70 genes in the genome of S. marmoratus are identified, and their expression patterns were investigated using transcriptomic data from thermal stress experiments. The expansion and gene duplication events of Hsp70 genes from the Hspa4, Hspa8, and Hspa12a subfamilies in S. marmoratus are revealed by phylogenetic analysis. qRT-PCR expression patterns demonstrated that seven Hsp70 genes were significantly up-regulated and none were significantly down-regulated after heat treatment. Only the hsp70 gene was significantly up-regulated after cold treatment. The selection test further showed a purifying selection on the duplicated gene pairs, suggesting that these genes underwent subfunctionalization. Our results add novel insight to aquaculture and biological research on S. marmoratus, providing important information on how Hsp70 genes are regulated in Scorpaeniformes under thermal stress.

Genome-Wide Identification, Molecular Characterization, and Involvement in Response to Abiotic and Biotic Stresses of the HSP70 Gene Family in Turbot (Scophthalmus maximus)

Heat shock proteins 70 (HSP70s) are known to play essential roles in organisms' response mechanisms to various environmental stresses. However, no systematic identification and functional analysis has been conducted for HSP70s in the turbot (Scophthalmus maximus), a commercially important worldwide flatfish. Herein, 16 HSP70 genes unevenly distributed on nine chromosomes were identified in the turbot at the genome-wide level. Analyses of gene structure, motif composition, and phylogenetic relationships provided valuable data on the HSP70s regarding their evolution, classification, and functional diversity. Expression profiles of the HSP70 genes under five different stresses were investigated by examining multiple RNA-seq datasets. Results showed that 10, 6, 8, 10, and 9 HSP70 genes showed significantly up- or downregulated expression after heat-induced, salinity-induced, and Enteromyxum scophthalmi, Vibrio anguillarum, and Megalocytivirus infection-induced stress, respectively. Among them, hsp70 (hspa1a), hspa1b, and hspa5 showed significant responses to each kind of induced stress, and qPCR analyses further validated their involvement in comprehensive anti-stress, indicating their involvement in organisms' anti-stress mechanisms. These findings not only provide new insights into the biological function of HSP70s in turbot adapting to various environmental stresses, but also contribute to the development of molecular-based selective breeding programs for the production of stress-resistant turbot strains in the aquaculture industry.

Phylogenomic analysis of the bowfin (Amia calva) reveals unrecognized species diversity in a living fossil lineage

The Bowfin (Amia calva), as currently recognized, represents the sole living member of the family Amiidae, which dates back to approximately 150 Ma. Prior to 1896, 13 species of extant Bowfins had been described, but these were all placed into a single species with no rationale or analysis given. This situation has persisted until the present day, with little attention given to re-evaluation of those previously described nominal forms. Here, we present a phylogenomic analysis based on over 21,000 single nucleotide polymorphisms (SNPs) from 94 individuals that unambiguously demonstrates the presence of at least two independent evolutionary lineages within extant Amia populations that merit species-level standing, as well as the possibility of two more. These findings not only expand the recognizable species diversity in an iconic, ancient lineage, but also demonstrate the utility of such methods in addressing previously intractable questions of molecular systematics and phylogeography in slowly evolving groups of ancient fishes.

Improved biomarker discovery through a plot twist in transcriptomic data analysis

Background

Transcriptomic analysis is crucial for understanding the functional elements of the genome, with the classic method consisting of screening transcriptomics datasets for differentially expressed genes (DEGs). Additionally, since 2005, weighted gene co-expression network analysis (WGCNA) has emerged as a powerful method to explore relationships between genes. However, an approach combining both methods, i.e., filtering the transcriptome dataset by DEGs or other criteria, followed by WGCNA (DEGs + WGCNA), has become common. This is of concern because such approach can affect the resulting underlying architecture of the network under analysis and lead to wrong conclusions. Here, we explore a plot twist to transcriptome data analysis: applying WGCNA to exploit entire datasets without affecting the topology of the network, followed with the strength and relative simplicity of DEG analysis (WGCNA + DEGs). We tested WGCNA + DEGs against DEGs + WGCNA to publicly available transcriptomics data in one of the most transcriptomically complex tissues and delicate processes: vertebrate gonads undergoing sex differentiation. We further validate the general applicability of our approach through analysis of datasets from three distinct model systems: European sea bass, mouse, and human.

Results

In all cases, WGCNA + DEGs clearly outperformed DEGs + WGCNA. First, the network model fit and node connectivity measures and other network statistics improved. The gene lists filtered by each method were different, the number of modules associated with the trait of interest and key genes retained increased, and GO terms of biological processes provided a more nuanced representation of the biological question under consideration. Lastly, WGCNA + DEGs facilitated biomarker discovery.

Conclusions

We propose that building a co-expression network from an entire dataset, and only thereafter filtering by DEGs, should be the method to use in transcriptomic studies, regardless of biological system, species, or question being considered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01398-w.

The genome-wide identification and adaptive evolution of slc9 genes in Leuciscus waleckii under extremely alkaline conditions

The solute carrier family 9 (slc9) genes, especially slc9a isoform coding proteins contribute to electroneutral countertransport of H⁺ for Na⁺ across the plasmalemmal and organellar membranes, intracellular pH and cellular volume regulation as well as the electrolyte, acid-base, and fluid volume homeostasis at the systemic level. These functional properties determine a potential basis for organisms to challenge stressful conditions. However, these well-done researches have been reported more in mammals. Thus, in this study, a total of eleven slc9 genes were identified from the latest version genome of L. waleckii, a cyprinid fish that could tolerate extremely alkaline environments (pH 9.6). The evolutionary footprint of slc9 genes was uncovered via the analysis of copy numbers, gene structure, motif composition, chromosome location and phylogenetic relationship. More importantly, there were two SNPs located on 5' UTR and three non-synonymous mutations in the coding region of the slc9a3.2 gene by comparing freshwater with alkaline water populations attached to resequencing technology. Slc9a3.2 gene was a statistically significant low expression in gill tissue with extremely alkaline pressure. Generally, slc9 gene family in L. waleckii was highly conserved. Several important SNPs with high Fst values were identified where non-synonymous mutations occurred between freshwater and alkaline water populations, and they may play an important role in specific functional differentiation. Slc9 genes had clear tissue expression preferences and were involved in abiotic stress response, indicating their roles in physiological function and strong self-regulating capacity. Our insight into the genetic variations that take place in the individual genes under extreme conditions could provide a feasible example for studying specific molecular mechanisms based on genomic data with increasing environmental stress.

Genome-wide identification and transcriptome-based expression profiling of the Sox gene family in the spinyhead croaker (Collichthys lucidus)

Sox genes encode transcription factors with a high-mobility group (HMG) box, playing critical roles in the initiation and maintenance of a variety of developmental processes, such as sex determination and differentiation. In the present study, we identified 26 Sox genes in the genome of spinyhead croaker Collichthys lucidus (Richardson, 1844) with homology-based analysis of the HMG box. The transcriptome-based expression profiles revealed that the expression of the Sox gene in gonads began to differ between sexes when the body length was 2.74 ± 0.24 cm. At that time, three Sox genes (Sox11b, Sox8a and Sox19) were significantly upregulated, accompanied by the downregulation of 12 Sox genes in the ovary, and six Sox genes were temporarily significantly upregulated in the testis. Afterwards, the expression profile of Sox genes changed only with a small amplitude in both the ovary and testis. For adult tissues, huge differences were observed in the expression profiles of Sox genes between ovaries and testes, as well as small differences in somatic tissues between sexes. These results provide clues to further decipher the role of Sox genes in the processes of sex determination and differentiation in spinyhead croaker and other teleosts.

Gene family evolution underlies cell-type diversification in the hypothalamus of teleosts

Hundreds of cell types form the vertebrate brain but it is largely unknown how similar cellular repertoires are between or within species or how cell-type diversity evolves. To examine cell-type diversity across and within species, we performed single-cell RNA sequencing of ~130,000 hypothalamic cells from zebrafish (Danio rerio) and surface and cave morphs of Mexican tetra (Astyanax mexicanus). We found that over 75% of cell types were shared between zebrafish and Mexican tetra, which diverged from a common ancestor over 150 million years ago. Shared cell types displayed shifts in paralogue expression that were generated by subfunctionalization after genome duplication. Expression of terminal effector genes, such as neuropeptides, was more conserved than the expression of their associated transcriptional regulators. Species-specific cell types were enriched for the expression of species-specific genes and characterized by the neofunctionalization of expression patterns of members of recently expanded or contracted gene families. Comparisons between surface and cave morphs revealed differences in immune repertoires and transcriptional changes in neuropeptidergic cell types associated with genomic differences. The single-cell atlases presented here are a powerful resource to explore hypothalamic cell types and reveal how gene family evolution and shifts in paralogue expression contribute to cellular diversity.

Promoter methylation and expression of the rassf4 gene in Nile tilapia (Oreochromis niloticus)

Oreochromis niloticus RAS-association domain family 4 (rassf4) was specifically expressed in the ovaries. Immunohistochemistry results showed that Rassf4 was located in follicular cells. The methylation percentages of the promoter region (-734/-1012) of rassf4 in the ovaries and testes were 13.28% and 92.85%, respectively. Deleting the fragment of -734 to -1012 sites significantly reduced the basal activity of the rassf4 promoter to 62.33%, which indicated that this region was important for the transcription of the rassf4 gene. This study lays the foundation for further research on the function of fish rassf4.

Diverse Regulation but Conserved Function: SOX9 in Vertebrate Sex Determination

Sex determination occurs early during embryogenesis among vertebrates. It involves the differentiation of the bipotential gonad to ovaries or testes by a fascinating diversity of molecular switches. In most mammals, the switch is SRY (sex determining region Y); in other vertebrates it could be one of a variety of genes including Dmrt1 or dmy. Downstream of the switch gene, SOX9 upregulation is a central event in testes development, controlled by gonad-specific enhancers across the 2 Mb SOX9 locus. SOX9 is a ‘hub’ gene of gonadal development, regulated positively in males and negatively in females. Despite this diversity, SOX9 protein sequence and function among vertebrates remains highly conserved. This article explores the cellular, morphological, and genetic mechanisms initiated by SOX9 for male gonad differentiation.

The Evolution of Sox Gene Repertoires and Regulation of Segmentation in Arachnids

The Sox family of transcription factors regulates many processes during metazoan development, including stem cell maintenance and nervous system specification. Characterizing the repertoires and roles of these genes can therefore provide important insights into animal evolution and development. We further characterized the Sox repertoires of several arachnid species with and without an ancestral whole-genome duplication and compared their expression between the spider Parasteatoda tepidariorum and the harvestman Phalangium opilio. We found that most Sox families have been retained as ohnologs after whole-genome duplication and evidence for potential subfunctionalization and/or neofunctionalization events. Our results also suggest that Sox21b-1 likely regulated segmentation ancestrally in arachnids, playing a similar role to the closely related SoxB gene, Dichaete, in insects. We previously showed that Sox21b-1 is required for the simultaneous formation of prosomal segments and sequential addition of opisthosomal segments in P. tepidariorum. We studied the expression and function of Sox21b-1 further in this spider and found that although this gene regulates the generation of both prosomal and opisthosomal segments, it plays different roles in the formation of these tagmata reflecting their contrasting modes of segmentation and deployment of gene regulatory networks with different architectures.

Investigating the molecular guts of endoderm formation using zebrafish

The vertebrate endoderm makes major contributions to the respiratory and gastrointestinal tracts and all associated organs. Zebrafish and humans share a high degree of genetic homology and strikingly similar endodermal organ systems. Combined with a multitude of experimental advantages, zebrafish are an attractive model organism to study endoderm development and disease. Recent functional genomics studies have shed considerable light on the gene regulatory programs governing early zebrafish endoderm development, while advances in biological and technological approaches stand to further revolutionize our ability to investigate endoderm formation, function and disease. Here, we discuss the present understanding of endoderm specification in zebrafish compared to other vertebrates, how current and emerging methods will allow refined and enhanced analysis of endoderm formation, and how integration with human data will allow modeling of the link between non-coding sequence variants and human disease.

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.

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.

Analysis of structural variants in four African cichlids highlights an association with developmental and immune related genes

BACKGROUND

East African lake cichlids are one of the most impressive examples of an adaptive radiation. Independently in Lake Victoria, Tanganyika, and Malawi, several hundreds of species arose within the last 10 million to 100,000 years. Whereas most analyses in cichlids focused on nucleotide substitutions across species to investigate the genetic bases of this explosive radiation, to date, no study has investigated the contribution of structural variants (SVs) in the evolution of adaptive traits across the three Great Lakes of East Africa.

RESULTS

Here, we annotate and characterize the repertoires and evolutionary potential of different SV classes (deletion, duplication, inversion, insertions and translocations) in four cichlid species: Haplochromis burtoni, Metriaclima zebra, Neolamprologus brichardi and Pundamilia nyererei. We investigate the patterns of gain and loss evolution for each SV type, enabling the identification of lineage specific events. Both deletions and inversions show a significant overlap with SINE elements, while inversions additionally show a limited, but significant association with DNA transposons. Inverted regions are enriched for genes regulating behaviour, or involved in skeletal and visual system development. We also find that duplicated regions show enrichment for genes associated with "antigen processing and presentation" and other immune related categories. Our pipeline and results were further tested by PCR validation of selected deletions and inversions, which confirmed respectively 7 out of 10 and 6 out of 9 events.

CONCLUSIONS

Altogether, we provide the first comprehensive overview of rearrangement evolution in East African cichlids, and some important insights into their likely contribution to adaptation.

Regeneration of the central nervous system-principles from brain regeneration in adult zebrafish

Poor recovery of neuronal functions is one of the most common healthcare challenges for patients with different types of brain injuries and/or neurodegenerative diseases. Therapeutic interventions face two major challenges: (1) How to generate neurons de novo to replenish the neuronal loss caused by injuries or neurodegeneration (restorative neurogenesis) and (2) How to prevent or limit the secondary tissue damage caused by long-term accumulation of glial cells, including microglia, at injury site (glial scar). In contrast to mammals, zebrafish have extensive regenerative capacity in numerous vital organs, including the brain, thus making them a valuable model to improve the existing therapeutic approaches for human brain repair. In response to injuries to the central nervous system (CNS), zebrafish have developed specific mechanisms to promote the recovery of the lost tissue architecture and functionality of the damaged CNS. These mechanisms include the activation of a restorative neurogenic program in a specific set of glial cells (ependymoglia) and the resolution of both the glial scar and inflammation, thus enabling proper neuronal specification and survival. In this review, we discuss the cellular and molecular mechanisms underlying the regenerative ability in the adult zebrafish brain and conclude with the potential applicability of these mechanisms in repair of the mammalian CNS.

Identification of genomic regions regulating sex determination in Atlantic salmon using high density SNP data

Background

A complete understanding of the genetic basis for sexual determination and differentiation is necessary in order to implement efficient breeding schemes at early stages of development. Atlantic salmon belongs to the family Salmonidae of fishes and represents a species of great commercial value. Although the species is assumed to be male heterogametic with XY sex determination, the precise genetic basis of sexual development remains unclear. The complexity is likely associated to the relatively recent salmonid specific whole genome duplication that may be responsible for certain genome instability. This instability together with the capacity of the sex-determining gene to move across the genome as reported by previous studies, may explain that sexual development genes are not circumscribed to the same chromosomes in all members of the species. In this study, we have used a 220 K SNP panel developed for Atlantic salmon to identify the chromosomes explaining the highest proportion of the genetic variance for sex as well as candidate regions and genes associated to sexual development in this species.

Results

Results from regional heritability analysis showed that the chromosomes explaining the highest proportion of variance in these populations were Ssa02 (heritability = 0.42, SE = 0.12) and Ssa21 (heritability = 0.26, SE = 0.11). After pruning by linkage disequilibrium, genome-wide association analyses revealed 114 SNPs that were significantly associated with sex, being Ssa02 the chromosome containing a greatest number of regions. Close examination of the candidate regions evidenced important genes related to sex in other species of Class Actinopterygii, including SDY, genes from family SOX, RSPO1, ESR1, U2AF2A, LMO7, GNRH-R, DND and FIGLA.

Conclusions

The combined results from regional heritability analysis and genome-wide association have provided new advances in the knowledge of the genetic regulation of sex determination in Atlantic salmon, supporting that Ssa02 is the candidate chromosome for sex in this species and suggesting an alternative population lineage in Spanish wild populations according to the results from Ssa21.

Using zebrafish to study skeletal genomics

While genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of skeletal diseases, animal models are required to identify causal mechanisms and to translate underlying biology into new therapies. Despite large-scale knockout mouse phenotyping efforts, the skeletal functions of most genes residing at GWAS-identified loci remain unknown, highlighting a need for complementary model systems to accelerate gene discovery. Over the past several decades, zebrafish (Danio rerio) has emerged as a powerful system for modeling the genetics of human diseases. In this review, our goal is to outline evidence supporting the utility of zebrafish for accelerating our understanding of human skeletal genomics, as well as gaps in knowledge that need to be filled for this purpose. We do this by providing a basic foundation of the zebrafish skeletal morphophysiology and phenotypes, and surveying evidence of skeletal gene homology and the use of zebrafish for post-GWAS analysis in other tissues and organs. We also outline challenges in translating zebrafish mutant phenotypes. Finally, we conclude with recommendations of future directions and how to leverage the large body of tools and knowledge of skeletal genetics in zebrafish for the needs of human skeletal genomic exploration. Due to their amenability to rapid genetic approaches, as well as the large number of conserved genetic and phenotypic features, there is a strong rationale supporting the use of zebrafish for human skeletal genomic studies.

The p53 gene family in vertebrates: Evolutionary considerations

The origin of the p53 gene family predates multicellular life since TP53 members of this gene family have been found in unicellular eukaryotes. In invertebrates one or two genes attributable to a TP53-like or TP63/73-like gene are present. The radiation into three genes, TP53, TP63, and TP73, has been reported as a vertebrate invention. TP53 is considered the "guardian of the genome" given its role in protecting cells against the DNA damage and cellular stressors. TP63 and TP73 play a role in epithelial development and neurogenesis, respectively. The evolution of the p53 gene family has been the subject of considerable analyses even if several questions remain still open. In this study we addressed the evolutionary history of the p53 gene family in vertebrates performing an extended microsyntenic investigation coupled with a phylogenetic analysis, together with protein domain organization and structure assessment. On the basis of our results we discussed a possible evolutionary scenario according to which a TP53/63/73 ancestor form gave rise to the current TP53 and a TP63/73 form, which in turn independently duplicated into two genes in agnathe and gnathostome lineages.

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.

The Gene Toolkit Implicated in Functional Sex in Sparidae Hermaphrodites: Inferences From Comparative Transcriptomics

Sex-biased gene expression is the mode through which sex dimorphism arises from a nearly identical genome, especially in organisms without genetic sex determination. Teleost fishes show great variations in the way the sex phenotype forms. Among them, Sparidae, that might be considered as a model family displays a remarkable diversity of reproductive modes. In this study, we sequenced and analyzed the sex-biased transcriptome in gonads and brain (the tissues with the most profound role in sexual development and reproduction) of two sparids with different reproductive modes: the gonochoristic common dentex, Dentex dentex, and the protandrous hermaphrodite gilthead seabream, Sparus aurata. Through comparative analysis with other protogynous and rudimentary protandrous sparid transcriptomes already available, we put forward common male and female-specific genes and pathways that are probably implicated in sex-maintenance in this fish family. Our results contribute to the understanding of the complex processes behind the establishment of the functional sex, especially in hermaphrodite species and set the groundwork for future experiments by providing a gene toolkit that can improve efforts to control phenotypic sex in finfish in the ever-increasingly important field of aquaculture.

Cluster expansion of apolipoprotein D (ApoD) genes in teleost fishes

BACKGROUND

Gene and genome duplication play important roles in the evolution of gene function. Compared to individual duplicated genes, gene clusters attract particular attention considering their frequent associations with innovation and adaptation. Here, we report for the first time the expansion of the apolipoprotein D (ApoD) ligand-transporter genes in a cluster manner specific to teleost fishes.

RESULTS

Based on comparative genomic and transcriptomic analyses, protein 3D structure comparison, positive selection detection and breakpoints detection, the single ApoD gene in the ancestor expanded into two clusters following a dynamic evolutionary pattern in teleost fishes. Orthologous genes show conserved expression patterns, whereas lineage-specific duplicated genes show tissue-specific expression patterns and even evolve new gene expression profiles. Positive selection occurred in branches before and after gene duplication, especially for lineage-specific duplicated genes. Cluster analyses based on protein 3D structure comparisons, especially comparisons of the four loops at the opening side, show gene duplication-segregating patterns. Duplicated ApoD genes are predicted to be associated with forkhead transcription factors and MAPK genes. ApoD clusters are located next to the breakpoints of genome rearrangements.

CONCLUSIONS

Here, we report the expansion of ApoD genes specific to teleost fishes in a cluster manner for the first time. Neofunctionalization and subfunctionalization were observed at both the protein and expression levels after duplication. Evidence from different aspects-i.e., abnormal expression-induced disease in humans, fish-specific expansion, predicted associations with forkhead transcription factors and MAPK genes, specific expression patterns in tissues related to sexual selection and adaptation, duplicated genes under positive selection and their location next to the breakpoints of genome rearrangements-suggests the potentially advantageous roles of ApoD genes in teleost fishes. The cluster expansion of ApoD genes specific to teleost fishes provides thus an ideal evo-devo model for studying gene duplication, cluster maintenance and new gene function emergence.

Duplication and expression of Sox genes in spiders

BACKGROUND

The Sox family of transcription factors is an important part of the genetic 'toolbox' of all metazoans examined to date and is known to play important developmental roles in vertebrates and insects. However, outside the commonly studied Drosophila model little is known about the repertoire of Sox family transcription factors in other arthropod species. Here we characterise the Sox family in two chelicerate species, the spiders Parasteatoda tepidariorum and Stegodyphus mimosarum, which have experienced a whole genome duplication (WGD) in their evolutionary history.

RESULTS

We find that virtually all of the duplicate Sox genes have been retained in these spiders after the WGD. Analysis of the expression of Sox genes in P. tepidariorum embryos suggests that it is likely that some of these genes have neofunctionalised after duplication. Our expression analysis also strengthens the view that an orthologue of vertebrate Group B1 genes, SoxNeuro, is implicated in the earliest events of CNS specification in both vertebrates and invertebrates. In addition, a gene in the Dichaete/Sox21b class is dynamically expressed in the spider segment addition zone, suggestive of an ancient regulatory mechanism controlling arthropod segmentation as recently suggested for flies and beetles. Together with the recent analysis of Sox gene expression in the embryos of other arthropods, our findings support the idea of conserved functions for some of these genes, including a potential role for SoxC and SoxD genes in CNS development and SoxF in limb development.

CONCLUSIONS

Our study provides a new chelicerate perspective to understanding the evolution and function of Sox genes and how the retention of duplicates of such important tool-box genes after WGD has contributed to different aspects of spider embryogenesis. Future characterisation of the function of these genes in spiders will help us to better understand the evolution of the regulation of important developmental processes in arthropods and other metazoans including neurogenesis and segmentation.

The Role of Transposable Elements in Speciation

Understanding the phenotypic and molecular mechanisms that contribute to genetic diversity between and within species is fundamental in studying the evolution of species. In particular, identifying the interspecific differences that lead to the reduction or even cessation of gene flow between nascent species is one of the main goals of speciation genetic research. Transposable elements (TEs) are DNA sequences with the ability to move within genomes. TEs are ubiquitous throughout eukaryotic genomes and have been shown to alter regulatory networks, gene expression, and to rearrange genomes as a result of their transposition. However, no systematic effort has evaluated the role of TEs in speciation. We compiled the evidence for TEs as potential causes of reproductive isolation across a diversity of taxa. We find that TEs are often associated with hybrid defects that might preclude the fusion between species, but that the involvement of TEs in other barriers to gene flow different from postzygotic isolation is still relatively unknown. Finally, we list a series of guides and research avenues to disentangle the effects of TEs on the origin of new species.

Teleost Fish-Specific Preferential Retention of Pigmentation Gene-Containing Families After Whole Genome Duplications in Vertebrates

Vertebrate pigmentation is a highly diverse trait mainly determined by neural crest cell derivatives. It has been suggested that two rounds (1R/2R) of whole-genome duplications (WGDs) at the basis of vertebrates allowed changes in gene regulation associated with neural crest evolution. Subsequently, the teleost fish lineage experienced other WGDs, including the teleost-specific Ts3R before teleost radiation and the more recent Ss4R at the basis of salmonids. As the teleost lineage harbors the highest number of pigment cell types and pigmentation diversity in vertebrates, WGDs might have contributed to the evolution and diversification of the pigmentation gene repertoire in teleosts. We have compared the impact of the basal vertebrate 1R/2R duplications with that of the teleost-specific Ts3R and salmonid-specific Ss4R WGDs on 181 gene families containing genes involved in pigmentation. We show that pigmentation genes (PGs) have been globally more frequently retained as duplicates than other genes after Ts3R and Ss4R but not after the early 1R/2R. This is also true for non-pigmentary paralogs of PGs, suggesting that the function in pigmentation is not the sole key driver of gene retention after WGDs. On the long-term, specific categories of PGs have been repeatedly preferentially retained after ancient 1R/2R and Ts3R WGDs, possibly linked to the molecular nature of their proteins (e.g., DNA binding transcriptional regulators) and their central position in protein-protein interaction networks. Taken together, our results support a major role of WGDs in the diversification of the pigmentation gene repertoire in the teleost lineage, with a possible link with the diversity of pigment cell lineages observed in these animals compared to other vertebrates.