Next Article in Journal
Comparative Transcriptome Analysis of Purple and Green Non-Heading Chinese Cabbage and Function Analyses of BcTT8 Gene
Next Article in Special Issue
Genetic Determination of the Amount of White Spotting: A Case Study in Siberian Cats
Previous Article in Journal
Genome-Wide Association Study of Potential Meat Quality Trait Loci in Ducks
Previous Article in Special Issue
Genomic Diversity Profiling and Breed-Specific Evolutionary Signatures of Selection in Arunachali Yak
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fish as Model Systems to Study Epigenetic Drivers in Human Self-Domestication and Neurodevelopmental Cognitive Disorders

by
Dafni Anastasiadi
1,*,
Francesc Piferrer
2,
Maren Wellenreuther
1,3 and
Antonio Benítez Burraco
4
1
Seafood Technologies, The New Zealand Institute for Plant and Food Research, Nelson 7010, New Zealand
2
Institut de Ciències del Mar, Spanish National Research Council (CSIC), 08003 Barcelona, Spain
3
School of Biological Sciences, The University of Auckland, Auckland 1010, New Zealand
4
Department of Spanish, Linguistics, and Theory of Literature (Linguistics), Faculty of Philology, University of Seville, 41004 Seville, Spain
*
Author to whom correspondence should be addressed.
Genes 2022, 13(6), 987; https://doi.org/10.3390/genes13060987
Submission received: 25 February 2022 / Revised: 28 April 2022 / Accepted: 11 May 2022 / Published: 31 May 2022
(This article belongs to the Special Issue Animal Domestication and Breeding)

Abstract

:
Modern humans exhibit phenotypic traits and molecular events shared with other domesticates that are thought to be by-products of selection for reduced aggression. This is the human self-domestication hypothesis. As one of the first types of responses to a novel environment, epigenetic changes may have also facilitated early self-domestication in humans. Here, we argue that fish species, which have been recently domesticated, can provide model systems to study epigenetic drivers in human self-domestication. To test this, we used in silico approaches to compare genes with epigenetic changes in early domesticates of European sea bass with genes exhibiting methylation changes in anatomically modern humans (comparison 1), and neurodevelopmental cognitive disorders considered to exhibit abnormal self-domestication traits, i.e., schizophrenia, Williams syndrome, and autism spectrum disorders (comparison 2). Overlapping genes in comparison 1 were involved in processes like limb morphogenesis and phenotypes like abnormal jaw morphology and hypopigmentation. Overlapping genes in comparison 2 affected paralogue genes involved in processes such as neural crest differentiation and ectoderm differentiation. These findings pave the way for future studies using fish species as models to investigate epigenetic changes as drivers of human self-domestication and as triggers of cognitive disorders.

1. Introduction

Domestication is defined as an “evolutionary process that substantially reshapes the genetic, physiological and behavioral profile of a species to adapt to a human-made environment” [1]. Historically and contemporarily, this process has affected the evolutionary trajectories of several economically and culturally important vertebrate species. Domestication usually involves selection of less aggressive and more human-tolerant individuals, plus some other specific features of interest such as meat, wool, and milk production. Nonetheless, new additional phenotypic traits emerge repeatedly in independent vertebrate domestication events, even at the early stages of living in a human-made environment prior to deliberate selection; a phenomenon first noticed by Darwin himself and currently known as the “domestication syndrome” [2]. The domestication syndrome has been predominantly described in mammals, likely due to the large number of mammalian domesticates with a long domestication history, sometimes dating back millennia (e.g., dogs). Phenotypic traits of the domestication syndrome, which are not necessarily present in all domesticates, include a decreased size of the brain, heart and teeth, vertebrae variability, caudal vertebrae changes, shorter muzzle, more frequent estrous cycles, floppy ears, curly tail and hair, and depigmentation [3,4]. The morphological and physiological traits have all been considered to have arisen as by-products of selection for increased tameness. Since these traits are associated with the final sites of migration of neural crest cells, mild developmental deficits affecting their development, migration or differentiation have been suggested by some authors as the underlying mechanisms of the appearance of the domestication syndrome, termed the “neural crest cell hypothesis” (NCCH) [2,5].
Modern humans, compared to extant apes and extinct hominins, exhibit phenotypic traits similar to those of other domesticated vertebrates, suggesting these may have also been produced as a by-product of selection for reduced aggression and increased sociality [6,7,8]. This is called the human self-domestication hypothesis [6,7,8] (HSD). Domestication syndrome-like morphological traits in anatomically modern humans (AMH) include a decreased brain and teeth size, facial robusticity, and sexual dimorphism, as well as neoteny [6,7,9]. Behavioral traits include reduced aggression, increased sociability, prolonged playing behavior, and overall more flexible social skills [6,7,9]. To understand and evaluate the HSD hypothesis, we need to distinguish between deliberate selection for improved traits, as occurred in, e.g., agricultural animals, and non-deliberate selection for prosociality arising from adaptation to novel environments, as expected for species hypothesized to have gone through a self-domestication process, like bonobos [6]. The latter should be seen through the lens of domestication being a multi-stage process, where non-deliberate selection arises in response to the new selective environment, e.g., often involving a lack of predators and an increase in food availability [7], changed environmental conditions [10], or the colonization of new environments [11], which are all known factors promoting prosocial behavior. Empirical support for the HSD hypothesis is challenging to obtain, nevertheless, comparative genomics have provided tentative support for it [12]. Recent results of an elegant study by Zanella et al. [13] using a molecular genetics approach are consistent with both the NCCH and the process of human self-domestication, specifically with regard to changes in the skull and the face [13,14]. This paper shows how in silico analyses of paleogenetic data and neurodevelopmental studies of human-altered phenotypes (e.g., Williams syndrome) can be integrated with experimental in vitro and in vivo research for empirically validating aspects of the HSD hypothesis. Still, the HSD hypothesis is not universally accepted, since the hypothesis presents some challenges, especially in relation to testing it [7], but also to the target of selection, which has been suggested to be socially-mediated emotional control and plasticity rather than selection against aggression [15].
Neurodevelopmental disorders in humans characterized by social and cognitive impairments may present traits of the domestication syndrome in an altered manner and thus may be linked to altered self-domestication. This is consistent with the view of self-domestication as a variable phenotype in the human species, e.g., [16], with this variability depending on genetic and environmental factors. More generally, this is consistent with the view that cognitive diseases mostly result from the alteration of recently evolved brain mechanisms, as these mechanisms lack compensatory mechanisms to respond to ontogenetic damage. People with autism spectrum disorders (ASD) and schizophrenia (SZ) exhibit abnormal aggressive behavior, abnormal responses to social cues, as well as tooth, ear, and facial anomalies [17,18]. In ASD, increased head and brain size, and generalized overgrowth are also present, while in SZ, decreased brain volume and reproductive dysfunctions occur [19,20]. Accordingly, these two cognitive disorders can be regarded as “less self-domesticated” and “more self-domesticated” phenotypes, respectively [19,20]. Moreover, Williams syndrome (WS), caused by the hemideletion of 28 genes in chromosome 7, is a clear example of a “more self-domesticated” phenotype [9,13]. People with WS show hypersociability, decreased aggression, reduced head and brain size, pointy ears, small teeth and jaws, depigmentation and accelerated sexual maturity [9]. Zanella et al. [13] used cell lines derived from WS subjects to establish the molecular links of morphological and behavioral domesticated traits in humans with neural crest development and migration processes. Therefore, cognitive disorders and the gene networks associated with them may be used as models for further testing the HSD hypothesis.
Domestication is a process of adaptation to a new selective environment and has been considered likely to involve epigenetic changes [21,22,23,24,25]. Epigenetic mechanisms offer a way for novel phenotypes to emerge rapidly in response to environmental changes and to prime the offspring, when inherited, to face environments based on the parental experience [26,27,28]. In the first stages of domestication, which coincides with the emergence of domestication syndrome traits, epigenetic changes established during early development might regulate gene expression in the neural crest, and be maintained throughout adulthood and inherited by the offspring. Multigenerational epigenetic inheritance is nearly ubiquitous in diverse animal species (see [29] for review). The persistence of the domestication environment, together with the stability and small effect of epigenetic changes in mild developmental deficits of the neural crest, could be expected to accelerate the adaptation [30]. After several generations, epigenetic changes could be genetically assimilated as genetic variants [25,31,32], hardwiring these changes. Partial evidence for this process comes from studies on mammals (dogs-wolves [23]), birds (red jungle fowl-modern chickens [33]), and fish species. The same process could be hypothesized to account for the first steps of HSD, as most differences between extinct hominins and AMHs are epigenetic by nature, having impacted features that are related to the domestication syndrome, particularly the morphology of the face [34]. To date, DNA methylation, but also DNA hydroxymethylation, have been the epigenetic mechanisms of focus, even though other epigenetic mechanisms like histone modifications and variants or non-coding RNAs may also play a role in domestication. Unfortunately, the latter cannot be studied in ancient DNA remains.
Domestication of fish species has a distinct history from terrestrial vertebrates [35], although it is scientifically considered to represent a similar process [36]. Until the 20th century, the majority of seafood relied on wild animal captures, with a few exceptions like the common carp (Cyprinus cyprio) in China ~8000 years ago or Nile tilapia (Oreochromis nilocitus) in Egypt ~3500 years ago [36,37]. In the last century, the domestication of aquatic species has expanded rapidly, with an estimated number of 368 vertebrates that have been domesticated for aquaculture, including teleost fish, frogs and reptiles [38]. Nevertheless, the majority of species are at the early stages of domestication, without closed life cycles in captivity and in the absence of deliberate selection for specific traits [38]. Nonetheless, in parallel with the domestication process, phenotypic traits involving the domestication syndrome, with changes in growth, reproduction, morphology, pigmentation and behaviour, have become manifested in domesticated fish [39,40,41]. Furthermore, sequencing of fish genomes has revolutionized vertebrate comparative genomics and has greatly contributed to our understanding of selection targets, evolutionary changes and speciation. Subsequently, fish have been suggested to serve as suitable models for human biomedical research [42,43]. Recently, epigenetic patterns emerging during the first stages of domestication, in the absence of genetic differences, have been studied in salmonids [44,45], European sea bass (Dicentrarchus labrax) [39], Nile tilapia (Oreochromis niloticus) [46,47] and grass carp (Ctenopharyngodon idellus) [48]. These epigenetic patterns of domestication are present in the sperm of several species, i.e., salmonids [45,49,50,51], showing the potential of intergenerational transfer, while in the European sea bass ~20% are found in early embryos, showcasing the importance of developmental aspects during early domestication [39]. Taken together, (1) the recent domestication events in fish, (2) the high degree of parallelism between fish and human domestication, particularly, the absence of deliberate selection in both domestication events, (3) the increasing availability of fish methylomes and (4) the use of fish as animal models in biomedical research, render fish promising candidate models to identify the epigenetic mechanisms that lead to the emergence of HSD, including their abnormal manifestation in neurodevelopmental disorders.
Comparative epigenomic studies between domesticated animals and humans are expected to demonstrate parallel or contrasting processes operating in addition to traditional genetic aspects [15]. Here, we argue that fish hold great advantages as models to study epigenetic drivers in HSD. To test our argument, we use comparative approaches to epigenomic patterns, exemplified by the best-studied modification, DNA methylation, between humans and the European sea bass. The European sea bass was chosen because (1) 25 years of selective breeding resulted in selective sweeps in genes similar to those found under positive selection in all domesticates tested, e.g., glutamate receptors [52,53], (2) it presents traits of the domestication syndrome shared with those found in terrestrial vertebrates, e.g., depigmentation and cranial changes [39], and (3) epigenetic patterns of domestication have been assessed in four tissue types representative of all three embryonic layers, thus reducing bias due to tissue-specificity [39]. In the present study, we compare epigenetic patterns of domesticated sea bass with epigenetic patterns of (1) AMH as opposed to archaic hominins (Neanderthals and Denisovans), and (2) neurodevelopmental cognitive disorders with an abnormal presentation of traits parallel to the domestication syndrome (SZ, WS and ASD; Figure S1). The goal of these comparisons was to investigate whether genes or pathways were consistently altered during the steps of early domestication in European sea bass and humans, with a potential impact on our species-specific distinctive cognition and behavior.

2. Materials and Methods

2.1. Data Collection

Comparative epigenomic analyses were divided into two major groups including early domesticates of the European sea bass vs. (1) AMH and (2) neurodevelopmental cognitive disorders. For this, we compiled five lists of genes identified as differentially methylated in the literature (Figure S1).

2.1.1. European Sea Bass Early Domesticates

In European sea bass, we previously conducted work to generate genome-wide DNA methylation patterns (Reduced Representation Bisulfite Sequencing, RRBS) in fish captured in the wild vs. offspring of wild fish reared in hatchery [39]. DNA methylation data from the brain, muscle, liver and testis can be accessed through the NCBI Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/geo/, accessed on 30 May 2022) with accession codes GSE104366 and GSE125124. Since these data were published, the European sea bass genome has been included in the Ensembl database (https://www.ensembl.org/, accessed on 30 May 2022). The genome assembly v1.0 in Ensembl is the same used for data analysis by Anastasiadi and Piferrer [39], however, gene annotation has since been updated according to the Ensembl Gene Annotation pipelines. To facilitate comparative epigenomic analysis with humans, we converted the list of genes with differentially methylated regions (DMRs) to the Ensembl genebuild released version from April 2020. To do this, the genomic coordinates (chromosome, start, end position) of DMRs and surrounding 5000 bp regions were intersected with the genebuild Dicentrarchus_labrax.seabass_V1.0.101.gtf. Chromosome names were as in the primary assembly. A total of 1181 unique genes with DMRs were identified in early domesticates.

2.1.2. Anatomically Modern Humans (AMH)

A detailed map of the evolutionary dynamics of DNA methylation in human groups was recently published [34]. DMRs specific to the AMH-lineage as compared to other hominin lineages, i.e., Denisovan and Neanderthal, were identified using a conservative approach to minimize false positives and variability due to factors such as sex or age, as well as DNA methylation data from chimpanzee samples. AMH-lineage DMRs are a set of 873 DMRs that overlap with the gene body or the promoter up to 5000 bp upstream of 588 genes (Supplementary Dataset S2 of Gokhman et al. [34]; Figure S1). The list of genes with DMRs was supplied by [34] with UCSC (University of California Santa Cruz) identifiers (IDs) and we used the https://biotools.fr/human/ucsc_id_converter (accessed on 30 May 2022) tool to convert them to Ensembl IDs to facilitate comparative epigenomics with the European sea bass.

2.1.3. Neurodevelopmental Cognitive Disorders

WS has a clear genetic origin with the hemideletion of 28 genes at 7q11.23. Some of these genes, e.g., BAZ1B, are involved in epigenetic regulation, such as chromatin remodeling, providing a link to the impact on epigenomic patterns in WS [13]. Differential DNA methylation between patients with WS and healthy individuals as controls has been reported in at least two cases in the literature. DMRs identified using the Infinium® HumanMethylation450 BeadChip array (Illumina, San Diego, CA, USA) in the blood of 20 WS patients vs. 15 healthy controls found DMRs intersecting 551 unique genes [54]. Differentially methylated cytosines (DMCs) were detected more recently in the blood of a larger sample of 90 WS patients vs. 34 healthy controls using the same array and these intersected with 143 unique genes [55]. The two gene lists were combined for further analysis as genes differentially methylated (DM) in WS, with a total of 624 different genes.
SZ is a complex psychiatric disorder and epigenome-wide association studies (EWAS) have been carried out to explore the role of DNA methylation in SZ pathophysiology, with discordant results. Recently, a meta-analysis of five EWAS datasets was published, including samples taken from different parts of the brain (frontal cortex, cerebellum, hippocampus, and prefrontal cortex), between 3 and 47 samples per study and using either the Infinium HumanMethylation450 BeadChip or HumanMethylation27 BeadChip (Illumina, San Diego, CA, USA) [56]. A total of 513 genes were commonly DM in combinations of 4–5 EWAS and these were used here for further analysis as the DM genes in SZ.
ASD refers to a group of complex neurodevelopmental disorders with heterogeneous symptoms and underlying etiology. ASD heritability is complex and genetic variants involved are diverse with their number ranging between 1000 and 3000 genes reflecting ASD heterogeneity [57]. Other molecular aspects to better understand ASD include epigenetic variants and several studies were published in the last years. This allowed us to apply more stringent criteria for inclusion in this study, mainly a minimum number of 15 samples and identification of DMRs, which are considered more robust than DMCs only. Genes from four studies published in the last 4 years, thus, included: (a) 31 genes with DM that were at the same time differentially expressed and common in three independent studies based on blood samples [58], (b) 181 core genes with DMRs detected using all three approaches in blood cells [57], (c) 145 unique differentially expressed genes with DMRs in blood cells from three ASD subphenotypes (severe, intermediate, mild) and a group of combined cases [59], and (d) 58 genes with DMRs detected in postmortem brain samples [60]. The four datasets combined led to a list of 411 unique ASD genes.

2.2. Comparative Analyses

The BioMart data mining tool from Ensembl was used to identify orthologues of human genes from the genome assembly GRCh38.p13 in the European sea bass genome. Duplicate entries were eliminated for further analysis. Thus, we identified unique orthologues as follows: 589 for AMH, 506 for WS, 532 for SZ, and 367 for ASD (Figure S1). The BioMart tool was used to identify paralogues of the human genes involved in neurodevelopmental cognitive disorders in the human genome (GRCh38.p13), in turn used to identify orthologues in the European sea bass genome. Duplicate entries from the combined list of original orthologues and orthologues of human paralogues were eliminated and the number of homologs finally available for comparative analyses was as follows: 3460 for WS, 4000 for SZ and 2994 for ASD.
Pairwise comparisons were performed with the fish early domesticates (FED) as a reference and one human group as its pair. Thus, four pairwise comparisons occurred every time: (1) FED vs. AMH, (2) FED vs. WS, (3) FED vs. SZ, and (4) FED vs. ASD. Overlaps between gene lists were identified and visualized using the InteractiVenn tool [61]. The significance of overlap was tested using Fisher’s exact test for testing the independence of two variables represented by a contingency table. As the genomic background for gene overlap testing, the total number of 23,382 genes in the European sea bass genome (Ensembl genebuild released April 2020) was set. Furthermore, we performed Monte Carlo permutations to test whether overlaps were higher than expected by chance. Random samples of genes were drawn without replacement from the 23,882 total gene list according to the specific gene list each time, e.g., to test the overlap of orthologues FED vs. AMH, 1181 genes for FED vs. 589 genes for AMH were randomly drawn in each iteration. The process was repeated 10,000 times and each time the length of the intersection or overlap between the two genes lists was counted. The standard score of permutation was calculated as observed-mean(permuted)/sd(permuted) and the p-value as: times permuted overlap is higher than observed overlap divided by the number of permutations (10,000). Fisher’s tests and permutations were performed using R (v. 4.0.0) [62] and Rstudio (v. 1.4.1717) [63].
The Enrichr tool was used for enrichment analyses and knowledge discovery of gene sets [64,65,66]. Enrichment analyses were performed for the initial lists of genes (FED, AMH, WS, SZ and ASD). Enriched pathways from the databases BioPlanet, Wikipathway, Mammalian Phenotype and GO-terms Biological Process were kept for further comparisons, which included overlap testing as previously carried out, with background on the total number of terms found in each library on Erichr. Reduction and visualization of GO-terms were aided by REViGO [67]. IDs of pathways were entered in InteractiVenn to detect overlaps and Fisher’s exact tests were run to detect the statistical significance of the overlap. Enrichment analyses were also performed for the genes that overlapped in a pairwise manner between FED genes and homologs (combined lists of orthologues and orthologues of paralogues).

3. Results

3.1. Differentially Methylated Genes during Early Domestication in European Sea Bass and in Humans Are Shared

The early stages of domestication are expected to be associated with DNA methylation changes. To compare DNA methylation changes associated with the early stages of domestication between European sea bass and human, two gene lists were retrieved. In FED 1181 genes with DMRs were detected as compared to wild fish. For humans, based on limited availability and accessibility to early AMH domesticate samples, DNA methylation patterns of present-day AMHs compared to other hominins and primates were considered the most relevant proxy. A total of 589 genes with DMRs were detected as orthologues of AMH. We detected an overlap of 45 genes between FED and AMH and this was significant (Fisher’s test, odds ratio = 1.577, p = 0.004; Figure 1a). Furthermore, we found 1.7 times more genes in common between the two gene lists than expected by chance alone (z-score = 13.62, p = 3 × 10−4; Figure 1b).
Among the genes with DMRs in both groups (Table 1), we detected several genes that were repeatedly found to be involved in domestication in several species. For example, ADAM metallopeptidases with thrombospondin type 1 motifs, ephrin (eph) receptors, members of the integrin family (alpha or beta), or fibroblast growth factor receptors have been detected in other domesticates (see Dataset 1 from [39] Anastasiadi and Piferrer for overview and Montague et al., Pendleton et al., Zhang et al., Kukekova et al., Carneiro et al. [68,69,70,71,72] for each species). One of these genes is nuclear factor I X (NFIX in humans and nfxib in fish) which was found to be in the top 10 genes with DMRs in AMHs showing a strong correlation between methylation and expression [34]. Several lines of evidence suggest that hypermethylation of NFIX associates with its downregulation in the AMH lineage [34]. In FEDs, nfixb was hypermethylated in the testis (+29.98%) but hypomethylated in the muscle tissue (−35.88%). In other tissues, other nuclear factor 1 isoforms contained DMRs: in muscle tissue, nuclear factor 1 a-type contained 2 DMRs with opposite methylation patterns (+20.69% and −42.30%), and in brain tissue, nuclear factor 1 a-type contained 2 hypomethylated DMRs (−27.36% and −34.02%) and nuclear factor 1 b-type isoform X2 contained a hypomethylated DMR (−30.14%).
We performed enrichment analyses to get insight into the functional roles of the overlapping genes. GO Biological Process enrichment analysis highlighted processes such as limb morphogenesis (GO:0035108, p = 0.045), histone modifications (GO:0016570, p = 0.024), T cell apoptotic processes (GO:0070231, p = 0.014), or granulocyte activation (GO:0036230, p = 0.021) as common (Figure 2a; for full list Table S1). Analysis of MGI Mammalian Phenotypes showed enrichment in traits typical of the domestication syndrome, such as abnormal snout morphology (MP:0000443, p-adjusted = 0.031) or hypopigmentation (MP:0005408, p-adjusted = 0.034; Figure 2b; for full list Table S2). Enrichment of WikiPathways showed that affected pathways include endochondral ossification with skeletal dysplasia (WP4808, p = 0.008), endochondral ossification (WP474, p = 0.008), or androgen receptor signaling pathway (WP138, p = 0.015; for full list Table S3).

3.2. Early Domestication in European Sea Bass and Neurodevelopmental Cognitive Disorders Affect Paralogue Genes

Genes exhibiting DNA methylation changes in patients with neurodevelopmental cognitive disorders with traits parallel to the domestication syndrome such as SZ, WS and ASD were obtained from the literature. A total of 532 genes with DM were orthologues to SZ patients, 506 genes with DM in WS patients and 367 genes with DM in ASD patients. These gene lists of orthologues were compared to the genes of FED to evaluate whether DNA methylation in common genes was affected by these conditions. The pairwise overlaps were not significant in all cases, with 28 genes overlapping in SZ (odds ratio = 1.04, p = 0.439; Figure S2a), 31 overlapping in WS (odds ratio = 1.233, p = 0.88; Figure S2b) and 23 genes overlapping in ASD (odds ratio = 1.262, p = 0.169; Figure S2c). Permutation testing for the pairwise comparisons showed that the number of overlaps was within the range expected by chance in the case of SZ (z-score = −0.59, p = 0.216; Figure S2d) and ASD (z-score = 2.61, p = 0.062; Figure S2f), and only marginally significant in the case of WS (z-score = 3.75, p = 0.049; Figure S2e).
In an attempt to overcome the constraints of the conservative approach applied here for orthologues and since key candidate genes of domestication were present in all pairwise comparisons, e.g., protocadherins, ADAM metallopeptidases, collagens and glutamate receptors, we then focused on comparisons of functional properties. Orthologue genes were submitted for enrichment analyses and pairwise comparisons were performed at the pathway level following the reasoning that similar processes may be affected by different genes. We considered four libraries targeted by Enrichr as the most informative in our case: Bioplanet, WikiPathways, GO-terms Biological Process, and MGI Mammalian Phenotype. Terms in all four libraries were examined for enrichment according to the gene lists we provided (FED, SZ, WS, and ASD) and pairwise comparisons of terms were performed as follows: (1) FED vs. SZ, (2) FED vs. WS, and (3) FED vs. ASD (Figure S3). In 42% of the comparisons, there was no overlap of terms, while in three cases there were between 1 and 4 terms overlapping. The overlaps of terms were significant only in cases of SZ for WikiPathways (odds ratio = 4.477, p = 0.003; Figure S3d) and GO Biological Process (odds ratio = 2.442, p = 0.002; Figure S3g). WikiPathways included endochondral ossification with skeletal dysplasia (WP4808) and endochondral ossification (WP474) like in the enrichment of orthologue genes overlapping in AMH, but also neural crest differentiation (WP2064). GO Biological Process enriched included development of the renal system (GO:0072001), kidney (GO:0001822), or ureteric bud (GO:0001657), as well as regulation of immune cells such as T-helper 17 and alpha-beta T (GO:2000317, GO:0046639, or GO:2000320). Taken together these results indicate that further comparative analyses could reveal more additional similarities.
To investigate the role of gene families, we compared gene lists containing not only the orthologues but also the paralogues of genes. The FED gene list was maintained in the original format and served as the control in the pairwise comparisons completed as above. For the other three gene lists (SZ, WS, and ASD), paralogues in the human genome were obtained by Biomart, merged with the original genes and then orthologues in the European sea bass genome were identified, resulting in lists containing unique homologues (orthologues and paralogues). The gene lists contained 4000 homologues for SZ, 3460 homologues for WS, and 2994 homologues for ASD. Overlap between all pairwise comparisons was significant with 241 genes common in SZ (odds ratio = 1.258, p = 0.001; Figure 3a, Dataset S1), 236 in WS (odds ratio = 1.470, p = 4.422 × 10−7; Figure 3b, Dataset S2), and 178 overlapping in ASD (odds ratio = 1.222, p = 0.011; Figure 3c, Dataset S3). Since these gene lists contain ~8 times more genes than previously, the significance of the overlaps could be attributed to larger numbers. To test whether the number of overlaps could be expected by chance due to large number of genes, we performed Monte Carlo permutations using random sampling of genes from the whole genome as previously. We found that overlaps between gene lists were higher than expected by chance in all cases, including SZ (z-score = 49.03, p = 0; Figure 3d), WS (z-score = 33.01, p = 0; Figure 3e), and ASD (z-score = 69.46, p = 0; Figure 3f). These results confirmed that there were similarities between genes DM early during fish domestication and homologues of genes DM in neurodevelopmental cognitive disorders with domestication syndrome traits.
To evaluate the functional properties of the core overlaps between genes in FED and lists of homologous genes in cognitive disorders, we performed enrichment analysis using Enrichr as previously. Pathways affected in all pairwise comparisons included neural crest differentiation (WP2064), ectoderm differentiation (WP2858), hair follicle development: organogenesis–part 2 of 3 (WP2839), arrhytmogenic right ventricular cardiomyopathy (WP2118; Figure 4a–c, full lists in Tables S4–S6). Pathways affected in at least two pairwise comparisons included endochondral ossification with skeletal dysplasia (WP4808) like in the core overlap of FED with orthologues of AMH, or also focal adhesion (WP306) and BMP signaling in eyelid development (WP3927) among others (Figure 4a–c).
Further functional analyses included GO-terms of Biological Process. GO-terms affected in all pairwise comparisons included embryonic morphogenesis of skeletal system (GO:0048704), digestive tract (GO:0048557) and organ (GO:0048562), regulation of morphogenesis of a branching structure (GO:0060688), morphogenesis of an epithelium (GO:0002009), neuromuscular junction development (GO:0007528), odontogenesis (GO:0042476) and positive regulation of fibroblast proliferation (GO:0048146; Figure 5a–c; full lists in Tables S7–S9). In SZ and WS, the extracellular matrix organization was the most significantly enriched GO-term. In ASD, the most significantly enriched GO-term was renal system development and among the enriched GO-terms, we detected glutamatergic synaptic transmission (Figure 5c), a process involving glutamate receptors that have been recognized as affected by domestication across species [39,53].

4. Discussion

We have shown that a sizeable portion of epigenetic changes in early European sea bass domesticates occur in similar genes when compared to AMHs, and in similar gene families as in human-specific neurodevelopmental cognitive disorders. Thus, parallel epigenetic changes seem to manifest in independent processes across vertebrates involving domestication, hypothesized to be self-domestication in the case of humans. Our finding that similar genes or gene families exhibited epigenetic changes between human groups and European sea bass provides evidence for domestication as a process affecting similar functional biological properties in vertebrates. Further, it indicates that fish species can be suitable models for research on epigenomics in the context of HSD, as well as human cognitive disorders.
For the purposes of this study, we implemented in silico analyses to compare the lists of genes that exhibited epigenetic changes, measured as differences in DNA methylation. We followed a very conservative approach and included layers of statistical testing, however, some inevitable limitations associated with the nature of the study are present. Genes with epigenetic changes have been pulled from different studies which have used distinct methodologies to interrogate methylation status (e.g., arrays or sequencing) and distinct algorithms to analyze them. In the case of AMH, the data were deduced from comparisons with reconstructed methylomes using a robust methodology. However, this dataset lacks actual methylation data for early AMHs, as comparisons were performed with Neanderthal and Denisovans, which are considered as non–self-domesticated hominin species [34]. With regards to neurodevelopmental diseases, due to their often complex etiology, there may be differences in genes detected as DM by different authors who may have used different sampling strategies. Thus, we chose to include only studies which fulfilled stringent criteria. For example, studies involving a very small number of samples, such as comparisons of a pair of twins, were excluded. However, we cannot rule out the possibility that the exact gene lists with epigenetic changes may vary slightly when following consistent and unified guidelines for their detection. Furthermore, to detect homologs, the Biomart tool from the Ensembl database was used, which is one of the most transparent approaches to perform the task since versions of the genome and annotations can be traced. For the enrichment analyses, genes have to be well characterized and included in the query databases to be informative of the affected pathways. Relying on these bioinformatics resources carries the inherent risks of minor modifications in future updated versions. Nevertheless, the results of this study can be interpreted while taking these limitations into account, since they are based on conservative inclusion criteria and statistical testing and can be used as a step for further research on comparative epigenomics between phylogenetically distant vertebrates.
As noted in the Introduction, the HSD hypothesis, as well as the involvement of neural crest cells in HSD, even though attractive, remained mostly supported theoretically until recently. Genomic approaches comparing genes under positive selection between domesticated mammals and AMHs are starting to be used as supportive evidence for the HSD hypothesis [12,73]. Recently, the hypothesis was empirically validated and the role of BAZ1B, a gene within the hemideleted region in WS with an established role in neural crest induction and migration, was demonstrated [13]. The implication of this gene in morphological and behavioral phenotypes typical of the domestication syndrome via neural crest cell development was further shown using zebrafish as a model [74]. Our comparative results between AMHs and early European sea bass domesticates provide additional support for the role of specific genes in key processes with an impact on (self-)domestication features and suggest a role for epigenetic regulation of their expression. That said, one limitation of our approach is that comparative (epi)genomic approaches should be taken with more caution in comparison to conclusive mechanistic studies. However, appropriate model systems for experimental studies are impossible to obtain for the HSD, as viable cells from early AMH are unavailable, and extremely difficult to obtain for other domesticated vertebrates. Thus, in silico comparative studies, such as this one, may provide insightful information on the genes undergoing molecular changes. An additional problem arises from the fact that the number of genes may increase with the number of phenotypic traits to be considered, which is large in the case of domestication and the domestication syndrome, and which, up to a degree, is species-specific [3]. Greater phylogenetic distance, like between humans and fish, may explain more dissimilarities between traits and, thus, more genes altered overall but with lower biological significance individually to explain domestication. These limitations need to be taken into account when interpreting comparative analyses.
With all these limitations in mind, specific genes involved in key processes have been underlined by our analyses. NFIX is associated with craniofacial skeletal disease phenotypes and related to speech capabilities, and it has already been highlighted for its role in the development of the AMH face and larynx [13]. Another gene common between early domesticated European sea bass and AMH was GLI family zinc finger 3 (GLI3) which is a known transcriptional repressor involved in tissue development, including limb development, and immune system development [75]. GLI3 has a role during embryogenesis, controlling thalamic development [76], as well as calvarial suture development [77], while in ≈98% of Altaic Neanderthals and Denisovans it contains a mutation that is mildly disruptive [78]. The RUNX family transcription factor 3, RUNX3, is involved in the developing spinal cord and also has a role in the language and social regions of brain [79,80]. SMOC1, as well as SMOC2, play a role in endochondral bone formation and are regulated by another member of the RUNX family transcription factor [81]. RUNX2 encodes a master transcription factor during vertebrate development involved in the globularization of the human skull/brain. RUNX2 is also involved in the development of thalamus, which is functionally connected to many genes that are important for brain and language development, and that have experienced changes in our recent evolutionary history [74]. NCOR2 has already been identified as under selection in dogs [69], and is part of the cranial neural crest gene expression program [82]. The above-mentioned genes participate in the enriched mammalian phenotypes detected which match the domestication syndrome traits, like abnormal cranium morphology, hypopigmentation or decreased body strength, but also in human-distinctive features potentially associated with our self-domestication. Similarly, GLI3 and SMOC1 participate in the enriched GO-terms processes related to limb development, including limb morphogenesis and embryonic digit morphogenesis. The GO-term most significantly enriched according to its p-value ranking was the negative regulation of alpha-beta T cell differentiation. This is likely due to the involvement of the above-mentioned genes, i.e., RUNX3, GLI3 and SMOC1, in the immune system as well. These results together reinforce the role of epigenetics in the regulation of similar genes associated with the domestication syndrome during the early stages of domestication in the absence of deliberate selection, as is the case in both humans and fish. These results also provide support for the view that domestication constitutes an example of “developmental bias”, i.e., when perturbed by an altered environment, complex organisms pursue a limited number of developmental pathways [4].
As also highlighted in the Introduction, neurodevelopmental cognitive disorders in humans have been previously suggested as models for testing the HSD hypothesis [9,19,20]. WS has already been used to gather molecular evidence for the shaping of the human face and behavior underlying HSD [13]. Our initial analyses in search of common genes and pathways epigenetically altered in European sea bass early domesticates and cognitive disorders were unsuccessful. However, even though orthologue genes seemed to be absent, it was evident that similar gene families were affected, thus, justifying our subsequent approach in the search of paralogues. The lack of common genes could be due to the phylogenetic distance between species and to the nature of conditions tested, i.e., disease phenotypes vs. fish under farming conditions.
In SZ and WS, genes of key families were affected including, ADAM metallopeptidases, bone morphogenetic proteins, ephrins, fibroblast growth factors, homeoboxes, laminins, and members of the TBC1 domain family. ADAM metallopeptidases and laminins constitute the core members of the most significantly enriched GO-term of overlapping genes in both comparisons: extracellular structure organization. The role of DM genes of the extracellular matrix has already been highlighted in relation to early domestication in fish, and especially for DM changes established already early during development [39]. At the same time, the brain extracellular matrix is known to have multiple roles in brain development and function, and abnormal alteration of this matrix is increasingly acknowledged as a key etiological factor involved in neurological and psychiatric disorders (see Dityatev et al. [83] for review). In ASD, genes were slightly different and included bone morphogenetic proteins, glutamate receptors, laminins, protocadherins, and semaphorins. Neural crest migration depends on the interaction of receptors, e.g., ephrins and receptors for bone morphogenetic proteins, with extracellular matrix molecules, e.g., laminins and semaphorins [84]. The term “neural crest differentiation” was enriched in the overlapping groups of genes and consistently found in all three neurodevelopmental disorders, together with ectoderm differentiation, hair follicle development: organogenesis–part 2 of 3, and arrhytmogenic right ventricular cardiomyopathy. Members of this term were fgfr2, pax3, axin2, hdac10, cdh2, hes1, tfap2a, tfap2b, and tcf7l1. Disorders of the processes related to the neural crest are often regarded as underlying SZ, WS, and ASD [84]. FGF has an essential embryonic function during vertebrate development and Fgf signalling has been shown to serve as a target for selection during the domestication [85]. In ASD, paralogues of two key genes found in the AMH comparison were also identified as epigenetically altered, i.e., runx3 and gli3. This reinforces the idea that parallel processes are involved in HSD phenotype emergence, either evolutionary or pathologically, supporting the view that cognitive diseases can result from changes in genes involved in the human evolution [17,18]. Together these results show that epigenetic changes occur in similar gene families in independent models of early (self-)domestication and that several of these genes have already an established role in the neural crest and other processes recognized as affected by (self-)domestication.
Fish as animal models have long been used in basic science. Small teleost fish, like zebrafish or medaka, has been recently considered as models to study human neurological disorders including ASD [86], peripheral neuropathy [87], and behavioral neuroscience [88] since they possess several key advantages [89]. First, they consist of a phylogenetically diverse group with species that have evolved phenotypes naturally mimicking human diseases, called “evolutionary mutant models” [90,91,92]. Cross-species comparisons allow for the identification of the best models to study a specific physiological pathway [43]. Furthermore, in model species like zebrafish, genetic mutants for specific genes can be easily generated. Second, since they are vertebrates, their brain’s basic structure and function exhibit similarities to humans showing conserved neuronal circuitry [93]. Third, teleost genomes show homology with 70% of genes associated with human diseases [94,95]. Fourth, model fish species larvae are transparent, offering the opportunity for direct observation of the central nervous system during the development [96]. Thus, the use of fish models to study neurodevelopmental cognitive disorders exhibiting (self-)domestication-related features already has a sound basis in previous research. Indeed, zebrafish has been used as a model for the three disorders studied here, SCZ [97], WS [98], and ASD [99]. Further, fish species have been suggested as models to investigate evolutionary questions [100], and their potential as models for domestication has been recently recognized [101]. Our findings that homolog genes were differentially methylated in both human disorders and early European sea bass domesticates provide further evidence for the use of fish species as models to study the epigenomic regulation implicated in HSD-related phenotypes, which has proven to be a key source of the human uniqueness [6].
For research related to the HSD hypothesis, fish not only possess the above-mentioned advantages but also show a key similarity distinct from most farm animals: fish domestication and HSD took place in absence of deliberate selection. Our result that DNA methylation changes in European sea bass early domesticates and human groups manifested in overlapping genes supports the implication of epigenetic mechanisms in domestication as a process of adaptation to a human-made environment, including the environment resulting from our self-domestication. A recent study using zebrafish investigated the role of neural crest in the morphological and behavioral domesticated phenotypes in [13] HSD. They found that a loss of function of the key gene in WS and for the neural crest, baz1b, identified as important previously in humans as well [13], resulted in mild neural crest deficiencies during development and behavioral changes related to stress and sociality in adulthood [74]. Furthermore, comparative genomics using domesticated mammals have already been used to shed light to the HSD hypothesis [12]. Together these results show that fish species can be implemented in comparative (epi)genomics approaches and functional studies to further shed light on the validity of the HSD hypothesis.

5. Conclusions

We have demonstrated the occurrence of parallel epigenetic changes during independent domestication events in phylogenetically distant vertebrates. These events were driven by living in human-made environments, including the creation of the very human-specific niche through self-domestication, rather than by intentional selection. Epigenetic changes could be the first level of response to a new environment and could later be genomically integrated. An important part of these parallel epigenetic changes arises in genes associated with the neural crest, further supporting the involvement of mild deficits during neural crest development in the emergence of the domestication syndrome. Other common epigenetic changes manifest in genes with neurological or morphological functions that have been associated with the domestication phenotype, including HSD. These findings contribute to our understanding of the initial molecular changes happening during early (self-)domestication and pave the way for future studies using fish as models to investigate epigenetic changes as drivers of HSD, but also as etiological factors of human-specific cognitive diseases.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13060987/s1, Figure S1: Conceptual design of the study; Figure S2: Overlap of orthologue genes differentially methylated in fish early domesticates (FED) and cognitive disorders; Figure S3: Overlap of pathways enriched associated with orthologue genes differentially methylated in fish early domesticates (FED) and cognitive disorders; Table S1: Enrichment of Gene Ontology (GO) Biological Process (BP) terms associated with genes shared between early fish domesticates and anatomically modern humans; Table S2: Enrichment of Mammalian Phenotype (2014) terms associated with genes shared between early fish domesticates and anatomically modern humans; Table S3: Enrichment of Wikipathways associated with genes shared between early fish domesticates and anatomically modern humans; Table S4: Enrichment of Wikipathways associated with genes shared between early fish domesticates and human groups with schizophrenia; Table S5: Enrichment of Wikipathways associated with genes shared between early fish domesticates and human groups with Williams syndrome; Table S6: Enrichment of Wikipathways associated with genes shared between early fish domesticates and human groups with autism spectrum disorders; Table S7: Enrichment of Gene Ontology (GO) Biological Process (BP) terms associated with genes shared between early fish domesticates and human groups with schizophrenia; Table S8: Enrichment of Gene Ontology (GO) Biological Process (BP) terms associated with genes shared between early fish domesticates and human groups with Williams syndrome; Table S9: Enrichment of Gene Ontology (GO) Biological Process (BP) terms associated with genes shared between early fish domesticates and human groups with autism spectrum disorders; Dataset S1: homologue genes with epigenetic changes in human groups with schrizophrenia; Dataset S2: homologue genes with epigenetic changes in human groups with Williams syndrome; Dataset S3: homologue genes with epigenetic changes in human groups with autism spectrum disorders.

Author Contributions

Conceptualization, D.A., F.P., M.W., A.B.B.; methodology, D.A.; formal analysis, D.A.; investigation, D.A.; data curation, D.A.; writing—original draft preparation, D.A.; writing—review and editing, F.P., M.W., A.B.B.; visualization, D.A.; supervision, F.P., M.W., A.B.B. All authors have read and agreed to the published version of the manuscript.

Funding

D.A. was supported by a Marsden grant Te Pūtea Rangahau managed by the Royal Society of New Zealand Te Apārangi to M.W. (MFP-PAF-1801). This research was supported in part by grant PID2020-114516GB-I00 funded by MCIN/AEI/ 10.13039/501100011033 to A.B.B., and AEI grant PID2019-108888RB-I00 to F.P.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used in this study have been previously published and the details are included in the Materials and Methods section. Any new data generated from re-analysis are included as Supplementary Materials.

Acknowledgments

Author F.P. acknowledges the “Severo Ochoa Centre of Excellence” accreditation (CEX2019-000928-S).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Range, F.; Marshall-Pescini, S. Comparing Wolves and Dogs: Current Status and Implications for Human ‘Self-Domestication’. Trends Cogn. Sci. 2022, 26, 337–349. [Google Scholar] [CrossRef] [PubMed]
  2. Wilkins, A.S.; Wrangham, R.W.; Fitch, W.T. The “Domestication Syndrome” in Mammals: A Unified Explanation Based on Neural Crest Cell Behavior and Genetics. Genetics 2014, 197, 795–808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Sánchez-Villagra, M.R.; Geiger, M.; Schneider, R.A. The Taming of the Neural Crest: A Developmental Perspective on the Origins of Morphological Covariation in Domesticated Mammals. R. Soc. Open Sci. 2016, 3, 160107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Wilkins, A.S. A Striking Example of Developmental Bias in an Evolutionary Process: The “Domestication Syndrome”. Evol. Dev. 2020, 22, 143–153. [Google Scholar] [CrossRef] [Green Version]
  5. Wilkins, A.S. Revisiting two hypotheses on the “domestication syndrome” in light of genomic data. Vavilov J. Genet. Breed 2017, 21, 435–442. [Google Scholar] [CrossRef]
  6. Hare, B. Survival of the Friendliest: Homo Sapiens Evolved via Selection for Prosociality. Annu. Rev. Psychol. 2017, 68, 155–186. [Google Scholar] [CrossRef]
  7. Sánchez-Villagra, M.R.; van Schaik, C.P. Evaluating the Self-Domestication Hypothesis of Human Evolution. Evol. Anthropol. Issues News Rev. 2019, 28, 133–143. [Google Scholar] [CrossRef]
  8. Hare, B.; Wobber, V.; Wrangham, R. The Self-Domestication Hypothesis: Evolution of Bonobo Psychology Is Due to Selection against Aggression. Anim. Behav. 2012, 83, 573–585. [Google Scholar] [CrossRef]
  9. Niego, A.; Benítez-Burraco, A. Williams Syndrome, Human Self-Domestication, and Language Evolution. Front. Psychol. 2019, 10, 521. [Google Scholar] [CrossRef] [Green Version]
  10. Spikins, P.; French, J.C.; John-Wood, S.; Dytham, C. Theoretical and Methodological Approaches to Ecological Changes, Social Behaviour and Human Intergroup Tolerance 300,000 to 30,000 BP. J. Archaeol. Method Theory 2021, 28, 53–75. [Google Scholar] [CrossRef]
  11. Brooks, J.; Yamamoto, S. The Founder Sociality Hypothesis. Ecol. Evol. 2021, 11, 14392–14404. [Google Scholar] [CrossRef] [PubMed]
  12. Theofanopoulou, C.; Gastaldon, S.; O’Rourke, T.; Samuels, B.D.; Messner, A.; Martins, P.T.; Delogu, F.; Alamri, S.; Boeckx, C. Self-Domestication in Homo Sapiens: Insights from Comparative Genomics. PLoS ONE 2017, 12, e0185306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zanella, M.; Vitriolo, A.; Andirko, A.; Martins, P.T.; Sturm, S.; O’Rourke, T.; Laugsch, M.; Malerba, N.; Skaros, A.; Trattaro, S.; et al. Dosage Analysis of the 7q11.23 Williams Region Identifies BAZ1B as a Major Human Gene Patterning the Modern Human Face and Underlying Self-Domestication. Sci. Adv. 2019, 5, eaaw7908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wilkins, A.S. A Molecular Investigation of Human Self-Domestication. Trends Genet. 2020, 36, 227–228. [Google Scholar] [CrossRef]
  15. Shilton, D.; Breski, M.; Dor, D.; Jablonka, E. Human Social Evolution: Self-Domestication or Self-Control? Front. Psychol. 2020, 11, 134. [Google Scholar] [CrossRef]
  16. Gleeson, B.T.; Kushnick, G. Female Status, Food Security, and Stature Sexual Dimorphism: Testing Mate Choice as a Mechanism in Human Self-Domestication. Am. J. Phys. Anthr. 2018, 167, 458–469. [Google Scholar] [CrossRef]
  17. Gibson, G.; Lacek, K.A. Canalization and Robustness in Human Genetics and Disease. Annu. Rev. Genet. 2020, 54, 189–211. [Google Scholar] [CrossRef]
  18. Benton, M.L.; Abraham, A.; LaBella, A.L.; Abbot, P.; Rokas, A.; Capra, J.A. The Influence of Evolutionary History on Human Health and Disease. Nat. Rev. Genet. 2021, 22, 269–283. [Google Scholar] [CrossRef]
  19. Benítez-Burraco, A.; Lattanzi, W.; Murphy, E. Language Impairments in ASD Resulting from a Failed Domestication of the Human Brain. Front. Neurosci. 2016, 10, 373. [Google Scholar] [CrossRef] [Green Version]
  20. Benítez-Burraco, A.; Di Pietro, L.; Barba, M.; Lattanzi, W. Schizophrenia and Human Self-Domestication: An Evolutionary Linguistics Approach. Brain. Behav. Evol. 2017, 89, 162–184. [Google Scholar] [CrossRef]
  21. Belyaev, D.K.; Ruvinsky, A.O.; Trut, L.N. Inherited Activation-Inactivation of the Star Gene in Foxes: Its Bearing on the Problem of Domestication. J. Hered. 1981, 72, 267–274. [Google Scholar] [CrossRef] [PubMed]
  22. Trut, L.; Oskina, I.; Kharlamova, A. Animal Evolution during Domestication: The Domesticated Fox as a Model. BioEssays 2009, 31, 349–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Koch, I.J.; Clark, M.M.; Thompson, M.J.; Deere-Machemer, K.A.; Wang, J.; Duarte, L.; Gnanadesikan, G.E.; McCoy, E.L.; Rubbi, L.; Stahler, D.R.; et al. The Concerted Impact of Domestication and Transposon Insertions on Methylation Patterns between Dogs and Grey Wolves. Mol. Ecol. 2016, 25, 1838–1855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Vogt, G. Facilitation of Environmental Adaptation and Evolution by Epigenetic Phenotype Variation: Insights from Clonal, Invasive, Polyploid, and Domesticated Animals. Environ. Epigenetics 2017, 3, dvx002. [Google Scholar] [CrossRef]
  25. Vogt, G. Epigenetic Variation in Animal Populations: Sources, Extent, Phenotypic Implications, and Ecological and Evolutionary Relevance. J. Biosci. 2021, 46, 24. [Google Scholar] [CrossRef]
  26. Jablonka, E. The Evolutionary Implications of Epigenetic Inheritance. Interface Focus 2017, 7, 20160135. [Google Scholar] [CrossRef]
  27. Shea, N.; Pen, I.; Uller, T. Three Epigenetic Information Channels and Their Different Roles in Evolution: Epigenetic Mechanisms and Evolution. J. Evol. Biol. 2011, 24, 1178–1187. [Google Scholar] [CrossRef] [Green Version]
  28. Day, T.; Bonduriansky, R. A Unified Approach to the Evolutionary Consequences of Genetic and Nongenetic Inheritance. Am. Nat. 2011, 178, E18–E36. [Google Scholar] [CrossRef]
  29. Anastasiadi, D.; Venney, C.J.; Bernatchez, L.; Wellenreuther, M. Epigenetic Inheritance and Reproductive Mode in Plants and Animals. Trends Ecol. Evol. 2021, 36, 1124–1140. [Google Scholar] [CrossRef]
  30. Kronholm, I.; Collins, S. Epigenetic Mutations Can Both Help and Hinder Adaptive Evolution. Mol. Ecol. 2016, 25, 1856–1868. [Google Scholar] [CrossRef]
  31. O’Dea, R.E.; Noble, D.W.A.; Johnson, S.L.; Hesselson, D.; Nakagawa, S. The Role of Non-Genetic Inheritance in Evolutionary Rescue: Epigenetic Buffering, Heritable Bet Hedging and Epigenetic Traps. Environ. Epigenetics 2016, 2, dvv014. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Danchin, E.; Pocheville, A.; Rey, O.; Pujol, B.; Blanchet, S. Epigenetically Facilitated Mutational Assimilation: Epigenetics as a Hub within the Inclusive Evolutionary Synthesis. Biol. Rev. 2019, 94, 259–282. [Google Scholar] [CrossRef] [Green Version]
  33. Belteky, J.; Agnvall, B.; Bektic, L.; Hoglund, A.; Jensen, P.; Guerrero-Bosagna, C. Epigenetics and Early Domestication: Differences in Hypothalamic DNA Methylation between Red Junglefowl Divergently Selected for High or Low Fear of Humans. Genet. Sel. Evol. 2018, 50, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Gokhman, D.; Nissim-Rafinia, M.; Agranat-Tamir, L.; Housman, G.; García-Pérez, R.; Lizano, E.; Cheronet, O.; Mallick, S.; Nieves-Colón, M.A.; Li, H.; et al. Differential DNA Methylation of Vocal and Facial Anatomy Genes in Modern Humans. Nat. Commun. 2020, 11, 1189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Teletchea, F. Domestication and Genetics: What a Comparison Between Land and Aquatic Species Can Bring? In Evolutionary Biology: Biodiversification from Genotype to Phenotype; Pontarotti, P., Ed.; Springer International Publishing: Cham, Switzerland, 2015; pp. 389–401. ISBN 978-3-319-19932-0. [Google Scholar]
  36. Teletchea, F. Fish Domestication in Aquaculture: 10 Unanswered Questions. Anim. Front. 2021, 11, 87–91. [Google Scholar] [CrossRef] [PubMed]
  37. Teletchea, F. Fish Domestication in Aquaculture: Reassessment and Emerging Questions. Cybium 2019, 43, 7–15. [Google Scholar] [CrossRef]
  38. Houston, R.D.; Bean, T.P.; Macqueen, D.J.; Gundappa, M.K.; Jin, Y.H.; Jenkins, T.L.; Selly, S.L.C.; Martin, S.A.M.; Stevens, J.R.; Santos, E.M.; et al. Harnessing Genomics to Fast-Track Genetic Improvement in Aquaculture. Nat. Rev. Genet. 2020, 21, 389–409. [Google Scholar] [CrossRef]
  39. Anastasiadi, D.; Piferrer, F. Epimutations in Developmental Genes Underlie the Onset of Domestication in Farmed European Sea Bass. Mol. Biol. Evol. 2019, 36, 2252–2264. [Google Scholar] [CrossRef]
  40. Wringe, B.F.; Purchase, C.F.; Fleming, I.A. In Search of a “Cultured Fish Phenotype”: A Systematic Review, Meta-Analysis and Vote-Counting Analysis. Rev. Fish Biol. Fish. 2016, 26, 351–373. [Google Scholar] [CrossRef]
  41. Milla, S.; Pasquet, A.; El Mohajer, L.; Fontaine, P. How Domestication Alters Fish Phenotypes. Rev. Aquac. 2021, 13, 388–405. [Google Scholar] [CrossRef]
  42. Sarropoulou, E.; Fernandes, J.M.O. Comparative Genomics in Teleost Species: Knowledge Transfer by Linking the Genomes of Model and Non-Model Fish Species. Comp. Biochem. Physiol. Part D Genom. Proteom. 2011, 6, 92–102. [Google Scholar] [CrossRef] [PubMed]
  43. Ahmad, S.F.; Jehangir, M.; Srikulnath, K.; Martins, C. Fish Genomics and Its Impact on Fundamental and Applied Research of Vertebrate Biology. Rev. Fish Biol. Fish. 2021, 32, 357–385. [Google Scholar] [CrossRef]
  44. Le Luyer, J.; Laporte, M.; Beacham, T.D.; Kaukinen, K.H.; Withler, R.E.; Leong, J.S.; Rondeau, E.B.; Koop, B.F.; Bernatchez, L. Parallel Epigenetic Modifications Induced by Hatchery Rearing in a Pacific Salmon. Proc. Natl. Acad. Sci. USA 2017, 114, 12964–12969. [Google Scholar] [CrossRef] [Green Version]
  45. Wellband, K.; Roth, D.; Linnansaari, T.; Curry, R.A.; Bernatchez, L. Environment-Driven Reprogramming of Gamete DNA Methylation Occurs during Maturation and Is Transmitted Intergenerationally in Atlantic Salmon. G3 Genes Genomes Genet. 2021, 11, jkab353. [Google Scholar] [CrossRef] [PubMed]
  46. Konstantinidis, I.; Sætrom, P.; Mjelle, R.; Nedoluzhko, A.V.; Robledo, D.; Fernandes, J.M.O. Major Gene Expression Changes and Epigenetic Remodelling in Nile Tilapia Muscle after Just One Generation of Domestication. Epigenetics 2020, 15, 1052–1067. [Google Scholar] [CrossRef] [Green Version]
  47. Podgorniak, T.; Dhanasiri, A.; Chen, X.; Ren, X.; Kuan, P.-F.; Fernandes, J. Early Fish Domestication Affects Methylation of Key Genes Involved in the Rapid Onset of the Farmed Phenotype. Epigenetics 2022, 1–18. [Google Scholar] [CrossRef]
  48. Li, Y.; Fu, B.; Zhang, J.; Xie, J.; Wang, G.; Jiang, P.; Tian, J.; Li, H.; Xia, Y.; Yu, E. Genome-wide analysis of DNA methylation reveals selection signatures of grass carp during the domestication. bioRxiv 2021. [Google Scholar] [CrossRef]
  49. Leitwein, M.; Laporte, M.; Le Luyer, J.; Mohns, K.; Normandeau, E.; Withler, R.; Bernatchez, L. Epigenomic Modifications Induced by Hatchery Rearing Persist in Germ Line Cells of Adult Salmon after Their Oceanic Migration. Evol. Appl. 2021, 14, 2402–2413. [Google Scholar] [CrossRef]
  50. Rodriguez Barreto, D.; Garcia de Leaniz, C.; Verspoor, E.; Sobolewska, H.; Coulson, M.; Consuegra, S. DNA Methylation Changes in the Sperm of Captive-Reared Fish: A Route to Epigenetic Introgression in Wild Populations. Mol. Biol. Evol. 2019, 36, 2205–2211. [Google Scholar] [CrossRef] [Green Version]
  51. Gavery, M.R.; Nichols, K.M.; Goetz, G.W.; Middleton, M.A.; Swanson, P. Characterization of Genetic and Epigenetic Variation in Sperm and Red Blood Cells from Adult Hatchery and Natural-Origin Steelhead, Oncorhynchus Mykiss. G3 Genes Genomes Genet. 2018, 8, 3723–3736. [Google Scholar] [CrossRef] [Green Version]
  52. Bertolini, F.; Geraci, C.; Schiavo, G.; Sardina, M.T.; Chiofalo, V.; Fontanesi, L. Whole Genome Semiconductor Based Sequencing of Farmed European Sea Bass (Dicentrarchus Labrax) Mediterranean Genetic Stocks Using a DNA Pooling Approach. Mar. Genom. 2016, 28, 63–70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. O’Rourke, T.; Boeckx, C. Glutamate Receptors in Domestication and Modern Human Evolution. Neurosci. Biobehav. Rev. 2020, 108, 341–357. [Google Scholar] [CrossRef] [PubMed]
  54. Strong, E.; Butcher, D.T.; Singhania, R.; Mervis, C.B.; Morris, C.A.; De Carvalho, D.; Weksberg, R.; Osborne, L.R. Symmetrical Dose-Dependent DNA-Methylation Profiles in Children with Deletion or Duplication of 7q11.23. Am. J. Hum. Genet. 2015, 97, 216–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Kimura, R.; Lardenoije, R.; Tomiwa, K.; Funabiki, Y.; Nakata, M.; Suzuki, S.; Awaya, T.; Kato, T.; Okazaki, S.; Murai, T.; et al. Integrated DNA Methylation Analysis Reveals a Potential Role for ANKRD30B in Williams Syndrome. Neuropsychopharmacology 2020, 45, 1627–1636. [Google Scholar] [CrossRef]
  56. Forero, D.A.; González-Giraldo, Y. Integrative In Silico Analysis of Genome-Wide DNA Methylation Profiles in Schizophrenia. J. Mol. Neurosci. 2020, 70, 1887–1893. [Google Scholar] [CrossRef]
  57. Hu, V.W.; Hong, Y.; Xu, M.; Shu, H.T. Altered DNA Methylation in a Severe Subtype of Idiopathic Autism: Evidence for Sex Differences in Affected Metabolic Pathways. Autism 2021, 25, 887–910. [Google Scholar] [CrossRef]
  58. Zhang, B.; Hu, X.; Li, Y.; Ni, Y.; Xue, L. Identification of Methylation Markers for Diagnosis of Autism Spectrum Disorder. Metab. Brain Dis. 2021, 37, 219–228. [Google Scholar] [CrossRef]
  59. Lee, E.C.; Hu, V.W. Phenotypic Subtyping and Re-Analysis of Existing Methylation Data from Autistic Probands in Simplex Families Reveal ASD Subtype-Associated Differentially Methylated Genes and Biological Functions. Int. J. Mol. Sci. 2020, 21, 6877. [Google Scholar] [CrossRef]
  60. Nardone, S.; Sams, D.S.; Zito, A.; Reuveni, E.; Elliott, E. Dysregulation of Cortical Neuron DNA Methylation Profile in Autism Spectrum Disorder. Cereb. Cortex 2017, 27, 5739–5754. [Google Scholar] [CrossRef] [Green Version]
  61. Heberle, H.; Meirelles, G.V.; da Silva, F.R.; Telles, G.P.; Minghim, R. InteractiVenn: A Web-Based Tool for the Analysis of Sets through Venn Diagrams. BMC Bioinform. 2015, 16, 169. [Google Scholar] [CrossRef]
  62. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2015. [Google Scholar]
  63. R Studio Team. R Studio: Integrated Development Environment for R; R Studio, Inc.: Boston, MA, USA, 2015. [Google Scholar]
  64. Chen, E.Y.; Tan, C.M.; Kou, Y.; Duan, Q.; Wang, Z.; Meirelles, G.V.; Clark, N.R.; Ma’ayan, A. Enrichr: Interactive and Collaborative HTML5 Gene List Enrichment Analysis Tool. BMC Bioinform. 2013, 14, 128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A Comprehensive Gene Set Enrichment Analysis Web Server 2016 Update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Xie, Z.; Bailey, A.; Kuleshov, M.V.; Clarke, D.J.B.; Evangelista, J.E.; Jenkins, S.L.; Lachmann, A.; Wojciechowicz, M.L.; Kropiwnicki, E.; Jagodnik, K.M.; et al. Gene Set Knowledge Discovery with Enrichr. Curr. Protoc. 2021, 1, e90. [Google Scholar] [CrossRef] [PubMed]
  67. Supek, F.; Bošnjak, M.; Škunca, N.; Šmuc, T. REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms. PLoS ONE 2011, 6, e21800. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Montague, M.J.; Li, G.; Gandolfi, B.; Khan, R.; Aken, B.L.; Searle, S.M.J.; Minx, P.; Hillier, L.W.; Koboldt, D.C.; Davis, B.W.; et al. Comparative Analysis of the Domestic Cat Genome Reveals Genetic Signatures Underlying Feline Biology and Domestication. Proc. Natl. Acad. Sci. USA 2014, 111, 17230–17235. [Google Scholar] [CrossRef] [Green Version]
  69. Pendleton, A.L.; Shen, F.; Taravella, A.M.; Emery, S.; Veeramah, K.R.; Boyko, A.R.; Kidd, J.M. Comparison of Village Dog and Wolf Genomes Highlights the Role of the Neural Crest in Dog Domestication. BMC Biol. 2018, 16, 64. [Google Scholar] [CrossRef] [Green Version]
  70. Zhang, Z.; Jia, Y.; Almeida, P.; Mank, J.E.; van Tuinen, M.; Wang, Q.; Jiang, Z.; Chen, Y.; Zhan, K.; Hou, S. Whole-Genome Resequencing Reveals Signatures of Selection and Timing of Duck Domestication. GigaScience 2018, 7, giy027. [Google Scholar] [CrossRef]
  71. Kukekova, A.V.; Johnson, J.L.; Xiang, X.; Feng, S.; Liu, S.; Rando, H.M.; Kharlamova, A.V.; Herbeck, Y.; Serdyukova, N.A.; Xiong, Z.; et al. Red Fox Genome Assembly Identifies Genomic Regions Associated with Tame and Aggressive Behaviours. Nat. Ecol. Evol. 2018, 2, 1479–1491. [Google Scholar] [CrossRef]
  72. Carneiro, M.; Rubin, C.-J.; Di Palma, F.; Albert, F.W.; Alföldi, J.; Barrio, A.M.; Pielberg, G.; Rafati, N.; Sayyab, S.; Turner-Maier, J. Rabbit Genome Analysis Reveals a Polygenic Basis for Phenotypic Change during Domestication. Science 2014, 345, 1074–1079. [Google Scholar] [CrossRef] [Green Version]
  73. Benítez-Burraco, A.; Chekalin, E.; Bruskin, S.; Tatarinova, T.; Morozova, I. Recent Selection of Candidate Genes for Mammal Domestication in Europeans and Language Change in Europe: A Hypothesis. Ann. Hum. Biol. 2021, 48, 313–320. [Google Scholar] [CrossRef]
  74. Torres-Pérez, J.V.; Anagianni, S.; Mech, A.M.; Havelange, W.; García-González, J.; Fraser, S.E.; Vallortigara, G.; Brennan, C.H. Behavioural Analysis of Loss of Function Zebrafish Supports baz1b as Master Regulator of Domestication. bioRxiv 2022. [Google Scholar] [CrossRef]
  75. Matissek, S.J.; Elsawa, S.F. GLI3: A Mediator of Genetic Diseases, Development and Cancer. Cell Commun. Signal. CCS 2020, 18, 54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Haddad-Tóvolli, R.; Heide, M.; Zhou, X.; Blaess, S.; Alvarez-Bolado, G. Mouse Thalamic Differentiation: Gli-Dependent Pattern and Gli-Independent Prepattern. Front. Neurosci. 2012, 6, 27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Tanimoto, Y.; Veistinen, L.; Alakurtti, K.; Takatalo, M.; Rice, D.P.C. Prevention of Premature Fusion of Calvarial Suture in GLI-Kruppel Family Member 3 (Gli3)-Deficient Mice by Removing One Allele of Runt-Related Transcription Factor 2 (Runx2). J. Biol. Chem. 2012, 287, 21429–21438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Castellano, S.; Parra, G.; Sánchez-Quinto, F.A.; Racimo, F.; Kuhlwilm, M.; Kircher, M.; Sawyer, S.; Fu, Q.; Heinze, A.; Nickel, B.; et al. Patterns of Coding Variation in the Complete Exomes of Three Neandertals. Proc. Natl. Acad. Sci. USA 2014, 111, 6666–6671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Inoue, K.; Shiga, T.; Ito, Y. Runx Transcription Factors in Neuronal Development. Neural Dev. 2008, 3, 20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Boeckx, C.; Benítez-Burraco, A. Globularity and Language-Readiness: Generating New Predictions by Expanding the Set of Genes of Interest. Front. Psychol. 2014, 5, 1–22. [Google Scholar] [CrossRef] [Green Version]
  81. Takahata, Y.; Hagino, H.; Kimura, A.; Urushizaki, M.; Kobayashi, S.; Wakamori, K.; Fujiwara, C.; Nakamura, E.; Yu, K.; Kiyonari, H.; et al. Smoc1 and Smoc2 Regulate Bone Formation as Downstream Molecules of Runx2. Commun. Biol. 2021, 4, 1199. [Google Scholar] [CrossRef]
  82. Williams, S.S.; Mear, J.P.; Liang, H.-C.; Potter, S.S.; Aronow, B.J.; Colbert, M.C. Large-Scale Reprogramming of Cranial Neural Crest Gene Expression by Retinoic Acid Exposure. Physiol. Genom. 2004, 19, 184–197. [Google Scholar] [CrossRef] [Green Version]
  83. Dityatev, A.; Seidenbecher, C.; Morawski, M. Brain Extracellular Matrix: An Upcoming Target in Neurological and Psychiatric Disorders. Eur. J. Neurosci. 2021, 53, 3807–3810. [Google Scholar] [CrossRef]
  84. Šimić, G.; Vukić, V.; Kopić, J.; Krsnik, Ž.; Hof, P.R. Molecules, Mechanisms, and Disorders of Self-Domestication: Keys for Understanding Emotional and Social Communication from an Evolutionary Perspective. Biomolecules 2020, 11, 2. [Google Scholar] [CrossRef] [PubMed]
  85. Rohner, N.; Bercsényi, M.; Orbán, L.; Kolanczyk, M.E.; Linke, D.; Brand, M.; Nüsslein-Volhard, C.; Harris, M.P. Duplication of Fgfr1 Permits Fgf Signaling to Serve as a Target for Selection during Domestication. Curr. Biol. CB 2009, 19, 1642–1647. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Meshalkina, D.A.; Kizlyk, M.N.; Kysil, E.V.; Collier, A.D.; Echevarria, D.J.; Abreu, M.S.; Barcellos, L.J.G.; Song, C.; Warnick, J.E.; Kyzar, E.J.; et al. Zebrafish Models of Autism Spectrum Disorder. Exp. Neurol. 2018, 299, 207–216. [Google Scholar] [CrossRef] [PubMed]
  87. Won, S.Y.; Choi, B.-O.; Chung, K.W.; Lee, J.E. Zebrafish Is a Central Model to Dissect the Peripheral Neuropathy. Genes Genom. 2019, 41, 993–1000. [Google Scholar] [CrossRef] [PubMed]
  88. Bandmann, O.; Burton, E.A. Genetic Zebrafish Models of Neurodegenerative Diseases. Neurobiol. Dis. 2010, 40, 58–65. [Google Scholar] [CrossRef]
  89. Choi, T.-Y.; Choi, T.-I.; Lee, Y.-R.; Choe, S.-K.; Kim, C.-H. Zebrafish as an Animal Model for Biomedical Research. Exp. Mol. Med. 2021, 53, 310–317. [Google Scholar] [CrossRef]
  90. Schartl, M. Beyond the Zebrafish: Diverse Fish Species for Modeling Human Disease. Dis. Model. Mech. 2014, 7, 181–192. [Google Scholar] [CrossRef] [Green Version]
  91. Albertson, R.C.; Cresko, W.; Detrich, H.W.; Postlethwait, J.H. Evolutionary Mutant Models for Human Disease. Trends Genet. TIG 2009, 25, 74–81. [Google Scholar] [CrossRef] [Green Version]
  92. Beck, E.A.; Healey, H.M.; Small, C.M.; Currey, M.C.; Desvignes, T.; Cresko, W.A.; Postlethwait, J.H. Advancing Human Disease Research with Fish Evolutionary Mutant Models. Trends Genet. 2022, 38, 22–44. [Google Scholar] [CrossRef]
  93. Parker, M.O.; Brock, A.J.; Walton, R.T.; Brennan, C.H. The Role of Zebrafish (Danio Rerio) in Dissecting the Genetics and Neural Circuits of Executive Function. Front. Neural Circuits 2013, 7, 63. [Google Scholar] [CrossRef] [Green Version]
  94. Sullivan, C.; Kim, C.H. Zebrafish as a Model for Infectious Disease and Immune Function. Fish Shellfish. Immunol. 2008, 25, 341–350. [Google Scholar] [CrossRef] [PubMed]
  95. Howe, K.; Clark, M.D.; Torroja, C.F.; Torrance, J.; Berthelot, C.; Muffato, M.; Collins, J.E.; Humphray, S.; McLaren, K.; Matthews, L.; et al. The Zebrafish Reference Genome Sequence and Its Relationship to the Human Genome. Nature 2013, 496, 498–503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Matsui, H. The Use of Fish Models to Study Human Neurological Disorders. Neurosci. Res. 2017, 120, 1–7. [Google Scholar] [CrossRef] [PubMed]
  97. Langova, V.; Vales, K.; Horka, P.; Horacek, J. The Role of Zebrafish and Laboratory Rodents in Schizophrenia Research. Front. Psychiatry 2020, 11, 703. [Google Scholar] [CrossRef]
  98. Komoike, Y.; Fujii, K.; Nishimura, A.; Hiraki, Y.; Hayashidani, M.; Shimojima, K.; Nishizawa, T.; Higashi, K.; Yasukawa, K.; Saitsu, H.; et al. Zebrafish Gene Knockdowns Imply Roles for Human YWHAG in Infantile Spasms and Cardiomegaly. Genesis 2010, 48, 233–243. [Google Scholar] [CrossRef] [PubMed]
  99. Rea, V.; Van Raay, T.J. Using Zebrafish to Model Autism Spectrum Disorder: A Comparison of ASD Risk Genes Between Zebrafish and Their Mammalian Counterparts. Front. Mol. Neurosci. 2020, 13, 575575. [Google Scholar] [CrossRef]
  100. Eckhard Witten, P.; Leonor Cancela, M. What Aquaculture Does for Taxonomy, Evo-Devo, Palaeontology, Biomechanics and Biomedical Research. J. Appl. Ichthyol. 2018, 34, 429–430. [Google Scholar] [CrossRef]
  101. Sánchez-Villagra, M. The Process of Animal Domestication; Princeton University Press: Princeton, NJ, USA, 2022; ISBN 978-0-691-21768-0. [Google Scholar]
Figure 1. Overlap of genes with epigenetic changes in fish early domesticates (FED) and anatomically modern humans (AMH). The overlap was tested using Fisher’s exact test for count data (a) and permutations (b). The results of permutations are represented as the distribution of the number of overlaps (shaded grey area) with the mean number of permuted overlaps (black vertical line) and significance threshold set to 0.05 (red line). The observed number of overlaps is shown by the green line and the distance of observed vs. expected (random) overlaps is shown with the black arrow. The z-score and the p-value indicate the significance of the overlap.
Figure 1. Overlap of genes with epigenetic changes in fish early domesticates (FED) and anatomically modern humans (AMH). The overlap was tested using Fisher’s exact test for count data (a) and permutations (b). The results of permutations are represented as the distribution of the number of overlaps (shaded grey area) with the mean number of permuted overlaps (black vertical line) and significance threshold set to 0.05 (red line). The observed number of overlaps is shown by the green line and the distance of observed vs. expected (random) overlaps is shown with the black arrow. The z-score and the p-value indicate the significance of the overlap.
Genes 13 00987 g001
Figure 2. Enrichment analysis of overlapping genes with epigenetic changes in fish early domesticates and anatomically modern humans. (a) GO Biological Process terms enrichment where for each GO-term the color indicates the log10-transformed p-value of enrichment. The semantic space x (y-axis) and the semantic space y (x-axis) are the result of multidimensional scaling done by REViGO and represent semantic similarities between GO-terms. (b) Enrichment of mammalian phenotypes (MGI Mammalian Phenotype 2014). Each phenotype term (y-axis) is attributed significance values after enrichment tests (performed by Enrichr) which include the −log10-transformed p-value of enrichment (x-axis and length of grey lines) and the combined score estimated by Enrichr (color of the bubble and legend). To facilitate visualization according to the significance of enrichment based on p-values, terms are ranked in descending order from top to bottom.
Figure 2. Enrichment analysis of overlapping genes with epigenetic changes in fish early domesticates and anatomically modern humans. (a) GO Biological Process terms enrichment where for each GO-term the color indicates the log10-transformed p-value of enrichment. The semantic space x (y-axis) and the semantic space y (x-axis) are the result of multidimensional scaling done by REViGO and represent semantic similarities between GO-terms. (b) Enrichment of mammalian phenotypes (MGI Mammalian Phenotype 2014). Each phenotype term (y-axis) is attributed significance values after enrichment tests (performed by Enrichr) which include the −log10-transformed p-value of enrichment (x-axis and length of grey lines) and the combined score estimated by Enrichr (color of the bubble and legend). To facilitate visualization according to the significance of enrichment based on p-values, terms are ranked in descending order from top to bottom.
Genes 13 00987 g002
Figure 3. Overlap of homologous genes with epigenetic changes in fish early domesticates (FED) and cognitive disorders. Pairwise comparisons are shown for FED vs. schizophrenia (SZ; (a,d)), Williams syndrome (WS; (b,e)) and autism spectrum disorders (ASD; (c,f)). Significance of overlaps were tested using Fisher’s exact test for count data (ac) and permutations (df). The results of permutations are represented as the distribution of number of overlaps (shaded grey areas) with mean number of permuted overlaps (black vertical lines) and significance threshold set to 0.05 (red lines). Observed number of overlaps is indicated by the green lines and the distance of observed vs. expected (random) overlaps are shown with the black arrow. The z-scores and the p-values indicate the significance of the overlaps.
Figure 3. Overlap of homologous genes with epigenetic changes in fish early domesticates (FED) and cognitive disorders. Pairwise comparisons are shown for FED vs. schizophrenia (SZ; (a,d)), Williams syndrome (WS; (b,e)) and autism spectrum disorders (ASD; (c,f)). Significance of overlaps were tested using Fisher’s exact test for count data (ac) and permutations (df). The results of permutations are represented as the distribution of number of overlaps (shaded grey areas) with mean number of permuted overlaps (black vertical lines) and significance threshold set to 0.05 (red lines). Observed number of overlaps is indicated by the green lines and the distance of observed vs. expected (random) overlaps are shown with the black arrow. The z-scores and the p-values indicate the significance of the overlaps.
Genes 13 00987 g003
Figure 4. Pathway enrichment of genes with epigenetic changes in fish early domesticates (FED) and homologues of neurodevelopmental cognitive disorders. Pathways of the library Wikipathways enriched in schizophrenia (SZ; (a)), Williams syndrome (WS; (b)) and autism spectrum disorders (ASD; (c)). Each pathway (y-axis) is attributed significance values after enrichment tests (performed by Enrichr) which include the −log10-transformed p-value of enrichment (x-axis and length of grey lines) and the combined score estimated by Enrichr (color of the bubble and legend). To facilitate visualization according to significance of enrichment based on p-values, terms are ranked in descending order from top to bottom.
Figure 4. Pathway enrichment of genes with epigenetic changes in fish early domesticates (FED) and homologues of neurodevelopmental cognitive disorders. Pathways of the library Wikipathways enriched in schizophrenia (SZ; (a)), Williams syndrome (WS; (b)) and autism spectrum disorders (ASD; (c)). Each pathway (y-axis) is attributed significance values after enrichment tests (performed by Enrichr) which include the −log10-transformed p-value of enrichment (x-axis and length of grey lines) and the combined score estimated by Enrichr (color of the bubble and legend). To facilitate visualization according to significance of enrichment based on p-values, terms are ranked in descending order from top to bottom.
Genes 13 00987 g004
Figure 5. Enrichment of the Gene Ontology (GO) terms of genes with epigenetic changes in fish early domesticates (FED) and homologs of neurodevelopmental cognitive disorders. GO Biological Process terms enrichment in schizophrenia (SZ; (a)), Williams syndrome (WS; (b)) and autism spectrum disorders (ASD; (c)). For each GO-term the color indicates the log10-transformed p-value of enrichment which is also represented by the x-axis. The semantic space x (y-axis) is the result of multidimensional scaling done by REViGO and represent semantic similarities between GO-terms.
Figure 5. Enrichment of the Gene Ontology (GO) terms of genes with epigenetic changes in fish early domesticates (FED) and homologs of neurodevelopmental cognitive disorders. GO Biological Process terms enrichment in schizophrenia (SZ; (a)), Williams syndrome (WS; (b)) and autism spectrum disorders (ASD; (c)). For each GO-term the color indicates the log10-transformed p-value of enrichment which is also represented by the x-axis. The semantic space x (y-axis) is the result of multidimensional scaling done by REViGO and represent semantic similarities between GO-terms.
Genes 13 00987 g005
Table 1. Common genes differentially methylated in European sea bass early domesticates and anatomically modern humans.
Table 1. Common genes differentially methylated in European sea bass early domesticates and anatomically modern humans.
Gene NameGene DescriptionEnsembl Gene Stable ID
adamts17
agap1
atp7b
ADAM metallopeptidase with thrombospondin type 1 motif, 17ENSDLAG00005007818
ArfGAP with GTPase domain, ankyrin repeat and PH domain 1ENSDLAG00005018378
ATPase copper transporting betaENSDLAG00005026064
bcr
carm1
BCR activator of RhoGEF and GTPaseENSDLAG00005004082
coactivator-associated arginine methyltransferase 1ENSDLAG00005025319
celsr1a
cemip
coro7
dab2ipb
cadherin EGF LAG seven-pass G-type receptor 1aENSDLAG00005009488
cell migration inducing hyaluronidase 1ENSDLAG00005002105
coronin 7ENSDLAG00005014078
DAB2 interacting protein bENSDLAG00005020932
DIP2C
ephb3a
disco-interacting protein 2 homolog CaENSDLAG00005023732
eph receptor B3aENSDLAG00005000091
eps8l2EPS8 like 2ENSDLAG00005011013
EYA2EYA transcriptional coactivator and phosphatase 2ENSDLAG00005013401
fbrsl1fibrosin-like 1ENSDLAG00005019385
fgfrl1afibroblast growth factor receptor like 1aENSDLAG00005002545
galnt18aUDP-N-acetylalphaDgalactosamine:polypeptideN-acetylgalactosaminyltransferase 18aENSDLAG00005020537
gli3GLI family zinc finger 3ENSDLAG00005018034
itga11bintegrin, alpha 11bENSDLAG00005013142
kaznbkazrin, periplakin interacting protein bENSDLAG00005001674
lasp1LIM and SH3 protein 1ENSDLAG00005018795
lhppphospholysine phosphohistidine inorganic pyrophosphate phosphataseENSDLAG00005006011
lmx1bbLIM homeobox transcription factor 1, beta bENSDLAG00005025877
magi1bmembrane associated guanylate kinase, WW and PDZ domain containing 1bENSDLAG00005022108
mast2microtubule associated serine/threonine kinase 2ENSDLAG00005007444
meis2aMeis homeobox 2aENSDLAG00005007335
msmo1methylsterol monooxygenase 1ENSDLAG00005023171
ncor2nuclear receptor corepressor 2ENSDLAG00005024501
neurl1aaneuralized E3 ubiquitin protein ligase 1AaENSDLAG00005018019
nfixbnuclear factor I XENSDLAG00005016844
pacs2phosphofurin acidic cluster sorting protein 2ENSDLAG00005000298
parvbparvin, betaENSDLAG00005021030
phactr3bphosphatase and actin regulator 3bENSDLAG00005012177
prex1phosphatidylinositol-3,4,5-trisphosphate-dependent Rac exchange factor 1ENSDLAG00005024474
rab3il1RAB3A interacting protein (rabin3)-like 1ENSDLAG00005012050
runx3RUNX family transcription factor 3ENSDLAG00005000657
sh3pxd2aaSH3 and PX domains 2AaENSDLAG00005018046
si:ch211–243o19.4si:ch211–243o19.4ENSDLAG00005011826
smoc1SPARC related modular calcium binding 1ENSDLAG00005010838
sorcs2sortilin-related VPS10 domain containing receptor 2ENSDLAG00005020908
tbc1d22aTBC1 domain family, member 22aENSDLAG00005010430
tgfbr2btransforming growth factor beta receptor 2bENSDLAG00005010792
ZNF423zinc finger protein 423ENSDLAG00005008914
ENSDLAG00005005197
ENSDLAG00005012304
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Anastasiadi, D.; Piferrer, F.; Wellenreuther, M.; Benítez Burraco, A. Fish as Model Systems to Study Epigenetic Drivers in Human Self-Domestication and Neurodevelopmental Cognitive Disorders. Genes 2022, 13, 987. https://doi.org/10.3390/genes13060987

AMA Style

Anastasiadi D, Piferrer F, Wellenreuther M, Benítez Burraco A. Fish as Model Systems to Study Epigenetic Drivers in Human Self-Domestication and Neurodevelopmental Cognitive Disorders. Genes. 2022; 13(6):987. https://doi.org/10.3390/genes13060987

Chicago/Turabian Style

Anastasiadi, Dafni, Francesc Piferrer, Maren Wellenreuther, and Antonio Benítez Burraco. 2022. "Fish as Model Systems to Study Epigenetic Drivers in Human Self-Domestication and Neurodevelopmental Cognitive Disorders" Genes 13, no. 6: 987. https://doi.org/10.3390/genes13060987

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop