A comparison of brain gene expression levels in domesticated and wild animals

Single-cell chromatin accessibility profiling reveals regulatory mechanisms and evolution in pig brains

BACKGROUND

Pig brains serve as a valuable biomedical model for studying brain-related diseases due to their significant structural similarities to the human brain. Furthermore, the long-term domestication and artificial selection of domestic pigs have profoundly shaped their brains, making them an interesting subject for research. However, a comprehensive understanding of the regulatory mechanisms governing pig brain function and their impact on various phenotypes remains elusive due to the high degree of cellular heterogeneity present in the brain.

RESULTS

In this study, we profiled 71,798 cells from domestic pig and wild boar cerebral cortex and cerebellum, identifying nine cell types, and integrated single-cell RNA sequencing data to explore cell type-specific regulatory landscapes and oligodendrocyte developmental trajectory. Furthermore, comparative analysis of each cell type between domestic pigs and wild boars indicated that oligodendrocyte progenitor cells may potentially exhibit a faster evolutionary rate. Finally, cross-species analysis suggested that, compared to humans, the proportion of sequence-conserved and functionally conserved regulatory elements in each cell type appears to be higher in pigs than in mice. Studies on the enrichment of genetic variants associated with 15 human diseases and complex traits in conserved regulatory elements across cell types indicated that immune-related diseases were more enriched in pigs, whereas neurological diseases were somewhat more enriched in mice. However, the enrichment of Alzheimer's disease-associated variants in pigs but not in mice suggests that pigs could be a more suitable model for this condition.

CONCLUSIONS

Our research offers preliminary insights into the heterogeneity of pig brains and suggests the potential underlying regulatory mechanisms. Additionally, we explore the possible impact of nervous system differences on phenotypic changes, which could lay the groundwork for further biomedical studies.

Genetic Contributions to Aggressive Behaviour in Pigs: A Comprehensive Review

Aggressive behaviour in pigs poses significant challenges to animal welfare, production efficiency, and economic performance in the pork industry. This review explores the multifaceted causes of pig aggression, focusing on genetic, environmental, and physiological factors. Aggression in pigs is categorized into social, maternal, fear-induced, play, and redirected aggression, with early-life hierarchies and environmental stressors playing critical roles. Physiological markers, such as elevated cortisol and reduced serotonin levels, are closely linked to aggressive behaviour, while dietary interventions, including tryptophan supplementation, have shown promise in mitigating aggression. Environmental factors, such as overcrowding, noise, and heat stress, exacerbate aggressive tendencies, whereas enrichment strategies, like music and improved housing conditions, can reduce stress and aggression. Genome-wide analyses have pinpointed specific polymorphisms in neurotransmitter genes (DRD2, SLC6A4, MAOA) and stress-response loci (RYR1) as significant predictors of porcine aggression. Advances in genomic technologies, including genome-wide association studies (GWASs) and transcriptomic analyses, have further elucidated the genetic and epigenetic underpinnings of aggressive behaviour. Practical application in breeding programmes remains challenging due to aggression polygenic nature and industry hesitancy toward genomic approaches. Future research should focus on integrating genetic markers into breeding programmes, developing multitrait selection indices, and exploring epigenetic modifications to improve animal welfare and production efficiency. By addressing these challenges, the pork industry can enhance both the well-being of pigs and the sustainability of production systems.

Gene expression comparisons between captive and wild shrew brains reveal captivity effects

Compared with their free-ranging counterparts, wild animals in captivity experience different conditions with lasting physiological and behavioural effects. Although shifts in gene expression are expected to occur upstream of these phenotypes, we found no previous gene expression comparisons of captive versus free-ranging mammals. We assessed gene expression profiles of three brain regions (cortex, olfactory bulb and hippocampus) of wild shrews (Sorex araneus) compared with shrews kept in captivity for two months and undertook sample dropout to examine robustness given limited sample sizes. Consistent with captivity effects, we found hundreds of differentially expressed genes in all three brain regions, 104 overlapping across all three, that enriched pathways associated with neurodegenerative disease, oxidative phosphorylation and genes encoding ribosomal proteins. In the shrew, transcriptomic changes detected under captivity resemble responses in several human pathologies, including major depressive disorder and neurodegeneration. While interpretations of individual genes are tempered by small sample sizes, we propose captivity influences brain gene expression and function and can confound analyses of natural processes in wild individuals under captive conditions.

Modulating vertebrate physiology by genomic fine-tuning of GPCR functions

G protein-coupled receptors (GPCRs) play a crucial role as membrane receptors, facilitating the communication of eukaryotic species with their environment and regulating cellular and organ interactions. Consequently, GPCRs hold immense potential in contributing to adaptation to ecological niches and responding to environmental shifts. Comparative analyses of vertebrate genomes reveal patterns of GPCR gene loss, expansion, and signatures of selection. Integrating these genomic data with insights from functional analyses of gene variants enables the interpretation of genotype-phenotype correlations. This review underscores the involvement of GPCRs in adaptive processes, presenting numerous examples of how alterations in GPCR functionality influence vertebrate physiology or, conversely, how environmental changes impact GPCR functions. The findings demonstrate that modifications in GPCR function contribute to adapting to aquatic, arid, and nocturnal habitats, influencing camouflage strategies, and specializing in particular dietary preferences. Furthermore, the adaptability of GPCR functions provides an effective mechanism in facilitating past, recent, or ongoing adaptations in animal domestication and human evolution and should be considered in therapeutic strategies and drug development.

Candidate SNP Markers Significantly Altering the Affinity of the TATA-Binding Protein for the Promoters of Human Genes Associated with Primary Open-Angle Glaucoma

Primary open-angle glaucoma (POAG) is the most common form of glaucoma. This condition leads to optic nerve degeneration and eventually to blindness. Tobacco smoking, alcohol consumption, fast-food diets, obesity, heavy weight lifting, high-intensity physical exercises, and many other bad habits are lifestyle-related risk factors for POAG. By contrast, moderate-intensity aerobic exercise and the Mediterranean diet can alleviate POAG. In this work, we for the first time estimated the phylostratigraphic age indices (PAIs) of all 153 POAG-related human genes in the NCBI Gene Database. This allowed us to separate them into two groups: POAG-related genes that appeared before and after the phylum Chordata, that is, ophthalmologically speaking, before and after the camera-type eye evolved. Next, in the POAG-related genes' promoters, we in silico predicted all 3835 candidate SNP markers that significantly change the TATA-binding protein (TBP) affinity for these promoters and, through this molecular mechanism, the expression levels of these genes. Finally, we verified our results against five independent web services-PANTHER, DAVID, STRING, MetaScape, and GeneMANIA-as well as the ClinVar database. It was concluded that POAG is likely to be a symptom of the human self-domestication syndrome, a downside of being civilized.

Abundant genetic variation is retained in many laboratory schistosome populations

Schistosomes are obligately sexual blood flukes that can be maintained in the laboratory using freshwater snails as intermediate and rodents as definitive hosts. The genetic composition of laboratory schistosome populations is poorly understood: whether genetic variation has been purged due to serial inbreeding or retained is unclear. We sequenced 19 - 24 parasites from each of five laboratory Schistosoma mansoni populations and compared their genomes with published exome data from four S. mansoni field populations. We found abundant genomic variation (0.897 - 1.22 million variants) within laboratory populations: these retained on average 49% (π = 3.27e-04 - 8.94e-04) of the nucleotide diversity observed in the four field parasite populations (π = 1.08e-03 - 2.2e-03). However, the pattern of variation was very different in laboratory and field populations. Tajima's D was positive in all laboratory populations except SmBRE, indicative of recent population bottlenecks, but negative in all field populations. Current effective population size estimates of laboratory populations were lower (2 - 258) compared to field populations (3,174 - infinity). The distance between markers at which linkage disequilibrium (LD) decayed to 0.5 was longer in laboratory populations (59 bp - 180 kb) compared to field populations (9 bp - 9.5 kb). SmBRE was the least variable; this parasite also shows low fitness across the lifecycle, consistent with inbreeding depression. The abundant genetic variation present in most laboratory schistosome populations has several important implications: (i) measurement of parasite phenotypes, such as drug resistance, using laboratory parasite populations will determine average values and underestimate trait variation; (ii) genome-wide association studies (GWAS) can be conducted in laboratory schistosome populations by measuring phenotypes and genotypes of individual worms; (iii) genetic drift may lead to divergence in schistosome populations maintained in different laboratories. We conclude that the abundant genetic variation retained within many laboratory schistosome populations can provide valuable, untapped opportunities for schistosome research.

The Domestication of Wild Boar Could Result in a Relaxed Selection for Maintaining Olfactory Capacity

Domesticated animals are artificially selected to exhibit desirable traits, however not all traits of domesticated animals are the result of deliberate selection. Loss of olfactory capacity in the domesticated pig (Sus scrofa domesticus) is one example. We used whole transcriptome analysis (RNA-Seq) to compare patterns of gene expression in the olfactory mucosa of the pig and two subspecies of wild boar (Sus scrofa), and investigate candidate genes that could be responsible for the loss of olfactory capacity. We identified hundreds of genes with reductions in transcript abundance in pig relative to wild boar as well as differences between the two subspecies of wild boar. These differences were detected mainly in genes involved in the formation and motility of villi, cilia and microtubules, functions associated with olfaction. In addition, differences were found in the abundances of transcripts of genes related to immune defenses, with the highest levels in continental wild boar subspecies. Overall, the loss of olfactory capacity in pigs appears to have been accompanied by reductions in the expression of candidate genes for olfaction. These changes could have resulted from unintentional selection for reduced olfactory capacity, relaxed selection for maintaining olfactory capacity, pleiotropic effects of genes under selection, or other non-selective processes. Our findings could be a cornerstone for future researches on wild boars, pigs, feral populations, and their evolutionary trajectories, aimed to provide tools to better calibrate species management as well as guidelines for breeders.

The genomic landscape of mammal domestication might be orchestrated by selected transcription factors regulating brain and craniofacial development

Domestication transforms once wild animals into tamed animals that can be then exploited by humans. The process entails modifications in the body, cognition, and behavior that are essentially driven by differences in gene expression patterns. Although genetic and epigenetic mechanisms were shown to underlie such differences, less is known about the role exerted by trans-regulatory molecules, notably transcription factors (TFs) in domestication. In this paper, we conducted extensive in silico analyses aimed to clarify the TF landscape of mammal domestication. We first searched the literature, so as to establish a large list of genes selected with domestication in mammals. From this list, we selected genes experimentally demonstrated to exhibit TF functions. We also considered TFs displaying a statistically significant number of targets among the entire list of (domestication) selected genes. This workflow allowed us to identify 5 candidate TFs (SOX2, KLF4, MITF, NR3C1, NR3C2) that were further assessed in terms of biochemical and functional properties. We found that such TFs-of-interest related to mammal domestication are all significantly involved in the development of the brain and the craniofacial region, as well as the immune response and lipid metabolism. A ranking strategy, essentially based on a survey of protein-protein interactions datasets, allowed us to identify SOX2 as the main candidate TF involved in domestication-associated evolutionary changes. These findings should help to clarify the molecular mechanics of domestication and are of interest for future studies aimed to understand the behavioral and cognitive changes associated to domestication.

RatDEGdb: a knowledge base of differentially expressed genes in the rat as a model object in biomedical research

The animal models used in biomedical research cover virtually every human disease. RatDEGdb, a knowledge base of the differentially expressed genes (DEGs) of the rat as a model object in biomedical research is a collection of published data on gene expression in rat strains simulating arterial hypertension, age-related diseases, psychopathological conditions and other human afflictions. The current release contains information on 25,101 DEGs representing 14,320 unique rat genes that change transcription levels in 21 tissues of 10 genetic rat strains used as models of 11 human diseases based on 45 original scientific papers. RatDEGdb is novel in that, unlike any other biomedical database, it offers the manually curated annotations of DEGs in model rats with the use of independent clinical data on equal changes in the expression of homologous genes revealed in people with pathologies. The rat DEGs put in RatDEGdb were annotated with equal changes in the expression of their human homologs in affected people. In its current release, RatDEGdb contains 94,873 such annotations for 321 human genes in 836 diseases based on 959 original scientific papers found in the current PubMed. RatDEGdb may be interesting first of all to human geneticists, molecular biologists, clinical physicians, genetic advisors as well as experts in biopharmaceutics, bioinformatics and personalized genomics. RatDEGdb is publicly available at https://www.sysbio.ru/RatDEGdb.

Identification and analysis of circRNAs in the prefrontal cortices of wild boar and domestic pig

Domestication caused significant differences in morphology and behavior between wild and domestic pigs. However, the regulatory role of circRNA in this event is unclear. Here, we analyzed circRNA expression patterns in the prefrontal cortices of wild boar and domestic pigs to determine the potential role of circRNAs in domestication. We identified a total of 11,375 circRNAs and found that 349 and 354 circRNAs were up-regulated in wild boar and Rongchang pig, respectively. Functional enrichment analysis showed that host genes of significantly highly-expressed circRNAs in wild boar were significantly enriched in neural synapse-related categories and the categories of 'regulation of defense response (p = 0.028)' and 'neural retina development (p = 4.32 × 10-3)'. Host genes of significantly highly-expressed circRNAs in Rongchang pig were specifically involved in 'chordate embryonic development (p = 2.38 × 10-4)'. Additionally, we constructed circRNA-miRNA-mRNA regulatory axes in wild boar and Rongchang pig and found more regulatory axes in wild boar that potentially regulate synaptic activities. We identified multiple circRNAs that may be related to domesticated characteristics, such as ssc_circ_6179 (ssc_circ_6179-ssc-miR-9847-HRH3, related to aggression) and ssc_circ_3027 (ssc_circ_3027-ssc-miR-4334-5p-HCRTR1, related to attention). This study provides a resource for further investigation of the molecular basis of pig domestication.

Rats respond to aversive emotional arousal of human handlers with the activation of the basolateral and central amygdala

Reading danger signals may save an animal's life, and learning about threats from others allows avoiding first-hand aversive and often fatal experiences. Fear expressed by other individuals, including those belonging to other species, may indicate the presence of a threat in the environment and is an important social cue. Humans and other animals respond to conspecifics' fear with increased activity of the amygdala, the brain structure crucial for detecting threats and mounting an appropriate response to them. It is unclear, however, whether the cross-species transmission of threat information involves similar mechanisms, e.g., whether animals respond to the aversively induced emotional arousal of humans with activation of fear-processing circuits in the brain. Here, we report that when rats interact with a human caregiver who had recently undergone fear conditioning, they show risk assessment behavior and enhanced amygdala activation. The amygdala response involves its two major parts, the basolateral and central, which detect a threat and orchestrate defensive responses. Further, we show that humans who learn about a threat by observing another aversively aroused human, similar to rats, activate the basolateral and centromedial parts of the amygdala. Our results demonstrate that rats detect the emotional arousal of recently aversively stimulated caregivers and suggest that cross-species social transmission of threat information may involve similar neural circuits in the amygdala as the within-species transmission.

Enrichment of self-domestication and neural crest function loci in the heritability of neurodevelopmental disorders

Self-domestication could contribute to shaping the biology of human brain and consequently the predisposition to neurodevelopmental disorders. Leveraging genome-wide data from the Psychiatric Genomics Consortium, we tested the enrichment of self-domestication and neural crest function loci with respect to the heritability of autism spectrum disorder, schizophrenia (SCZ in East Asian and European ancestries, EAS and EUR, respectively), attention-deficit/hyperactivity disorder, obsessive-compulsive disorder, and Tourette's syndrome (TS). Considering only self-domestication and neural-crest-function annotations in the linkage disequilibrium score regression (LDSC) model, our partitioned heritability analysis revealed statistically significant enrichments across all disorders investigated. The estimates of the heritability enrichments for self-domestication loci were similar across neurodevelopmental disorders, ranging from 0.902 (EAS SCZ, p = 4.55 × 10-20) to 1.577 (TS, p = 5.85 × 10-5). Conversely, a wider spectrum of heritability enrichment estimates was present for neural crest function with the highest enrichment observed for TS (enrichment = 3.453, p = 2.88 × 10-3) and the lowest for EAS SCZ (enrichment = 1.971, p = 3.81 × 10-3). Although these estimates appear to be strong, the enrichments for self-domestication and neural crest function were null once we included additional annotations related to different genomic features. This indicates that the effect of self-domestication on the polygenic architecture of neurodevelopmental disorders is not independent of other functions of human genome.

Candidate SNP Markers Significantly Altering the Affinity of TATA-Binding Protein for the Promoters of Human Hub Genes for Atherogenesis, Atherosclerosis and Atheroprotection

Atherosclerosis is a systemic disease in which focal lesions in arteries promote the build-up of lipoproteins and cholesterol they are transporting. The development of atheroma (atherogenesis) narrows blood vessels, reduces the blood supply and leads to cardiovascular diseases. According to the World Health Organization (WHO), cardiovascular diseases are the leading cause of death, which has been especially boosted since the COVID-19 pandemic. There is a variety of contributors to atherosclerosis, including lifestyle factors and genetic predisposition. Antioxidant diets and recreational exercises act as atheroprotectors and can retard atherogenesis. The search for molecular markers of atherogenesis and atheroprotection for predictive, preventive and personalized medicine appears to be the most promising direction for the study of atherosclerosis. In this work, we have analyzed 1068 human genes associated with atherogenesis, atherosclerosis and atheroprotection. The hub genes regulating these processes have been found to be the most ancient. In silico analysis of all 5112 SNPs in their promoters has revealed 330 candidate SNP markers, which statistically significantly change the affinity of the TATA-binding protein (TBP) for these promoters. These molecular markers have made us confident that natural selection acts against underexpression of the hub genes for atherogenesis, atherosclerosis and atheroprotection. At the same time, upregulation of the one for atheroprotection promotes human health.

The Australasian dingo archetype: de novo chromosome-length genome assembly, DNA methylome, and cranial morphology

BACKGROUND

One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long-read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed.

FINDINGS

We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on chromosomes 11, 16, 25, and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and 9 previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mitochondrial DNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified 2 differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphologic data, comprising geometric morphometric assessment of cranial morphology, place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue shows she had a larger cranial capacity than a similar-sized domestic dog.

CONCLUSIONS

These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphologic characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

Differentially Expressed Genes and Molecular Susceptibility to Human Age-Related Diseases

Mainstream transcriptome profiling of susceptibility versus resistance to age-related diseases (ARDs) is focused on differentially expressed genes (DEGs) specific to gender, age, and pathogeneses. This approach fits in well with predictive, preventive, personalized, participatory medicine and helps understand how, why, when, and what ARDs one can develop depending on their genetic background. Within this mainstream paradigm, we wanted to find out whether the known ARD-linked DEGs available in PubMed can reveal a molecular marker that will serve the purpose in anyone's any tissue at any time. We sequenced the periaqueductal gray (PAG) transcriptome of tame versus aggressive rats, identified rat-behavior-related DEGs, and compared them with their known homologous animal ARD-linked DEGs. This analysis yielded statistically significant correlations between behavior-related and ARD-susceptibility-related fold changes (log₂ values) in the expression of these DEG homologs. We found principal components, PC1 and PC2, corresponding to the half-sum and the half-difference of these log₂ values, respectively. With the DEGs linked to ARD susceptibility and ARD resistance in humans used as controls, we verified these principal components. This yielded only one statistically significant common molecular marker for ARDs: an excess of Fcγ receptor IIb suppressing immune cell hyperactivation.

The Australasian dingo archetype: De novo chromosome-length genome assembly, DNA methylome, and cranial morphology

Background

One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed.

Findings

We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on Chromosomes 11, 16, 25 and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and nine previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mtDNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified two differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphological data, comprising geometric morphometric assessment of cranial morphology place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue show she had a larger cranial capacity than a similar-sized domestic dog.

Conclusions

These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphological characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

Captivity induces large and population-dependent brain transcriptomic changes in wild-caught cane toads (Rhinella marina)

Gene expression levels are key molecular phenotypes at the interplay between genotype and environment. Mounting evidence suggests that short-term changes in environmental conditions, such as those encountered in captivity, can substantially affect gene expression levels. Yet, the exact magnitude of this effect, how general it is, and whether it results in parallel changes across populations are not well understood. Here, we take advantage of the well-studied cane toad, Rhinella marina, to examine the effect of short-term captivity on brain gene expression levels, and determine whether effects of captivity differ between long-colonized and vanguard populations of the cane toad's Australian invasion range. We compared the transcriptomes of wild-caught toads immediately assayed with those from toads captured from the same populations but maintained in captivity for seven months. We found large differences in gene expression levels between captive and wild-caught toads from the same population, with an over-representation of processes related to behaviour and the response to stress. Captivity had a much larger effect on both gene expression levels and gene expression variability in toads from vanguard populations compared to toads from long-colonized areas, potentially indicating an increased plasticity in toads at the leading edge of the invasion. Overall, our findings indicate that short-term captivity can induce large and population-specific transcriptomic changes, which has significant implications for studies comparing phenotypic traits of wild-caught organisms from different populations that have been held in captivity.

On taming the effect of transcript level intra-condition count variation during differential expression analysis: A story of dogs, foxes and wolves

The evolution of RNA-seq technologies has yielded datasets of scientific value that are often generated as condition associated biological replicates within expression studies. With expanding data archives opportunity arises to augment replicate numbers when conditions of interest overlap. Despite correction procedures for estimating transcript abundance, a source of ambiguity is transcript level intra-condition count variation; as indicated by disjointed results between analysis tools. We present TVscript, a tool that removes reference-based transcripts associated with intra-condition count variation above specified thresholds and we explore the effects of such variation on differential expression analysis. Initially iterative differential expression analysis involving simulated counts, where levels of intra-condition variation and sets of over represented transcripts are explicitly specified, was performed. Then counts derived from inter- and intra-study data representing brain samples of dogs, wolves and foxes (wolves vs. dogs and aggressive vs. tame foxes) were used. For simulations, the sensitivity in detecting differentially expressed transcripts increased after removing hyper-variable transcripts, although at levels of intra-condition variation above 5% detection became unreliable. For real data, prior to applying TVscript, ≈20% of the transcripts identified as being differentially expressed were associated with high levels of intra-condition variation, an over representation relative to the reference set. As transcripts harbouring such variation were removed pre-analysis, a discordance from 26 to 40% in the lists of differentially expressed transcripts is observed when compared to those obtained using the non-filtered reference. The removal of transcripts possessing intra-condition variation values within (and above) the 97^th and 95^th percentiles, for wolves vs. dogs and aggressive vs. tame foxes, maximized the sensitivity in detecting differentially expressed transcripts as a result of alterations within gene-wise dispersion estimates. Through analysis of our real data the support for seven genes with potential for being involved with selection for tameness is provided. TVscript is available at: https://sourceforge.net/projects/tvscript/.

Genome scanning of behavioral selection in a canine olfactory detection breeding cohort

Research on working dogs is growing rapidly due to increasing global demand. Here we report genome scanning of the risk of puppies being eliminated for behavioral reasons prior to entering the training phase of the US Transportation Security Administration's (TSA) canine olfactory detection breeding and training program through 2013. Elimination of dogs for behavioral rather than medical reasons was based on evaluations at three, six, nine and twelve months after birth. Throughout that period, the fostered dogs underwent standardized behavioral tests at TSA facilities, and, for a subset of tests, dogs were tested in four different environments. Using methods developed for family studies, we performed a case-control genome wide association study (GWAS) of elimination due to behavioral observation and testing results in a cohort of 528 Labrador Retrievers (2002-2013). We accounted for relatedness by including the pedigree as a covariate and maximized power by including individuals with phenotype, but not genotype, data (approximately half of this cohort). We determined genome wide significance based on Bonferroni adjustment of two quasi-likelihood score tests optimized for either small or nearly-fully penetrant effect sizes. Six loci were significant and five suggestive, with approximately equal numbers of loci for the two tests and frequencies of loci with single versus multiple mapped markers. Several loci implicate a single gene, including CHD2, NRG3 and PDE1A which have strong relevance to behavior in humans and other species. We briefly discuss how expanded studies of canine breeding programs could advance understanding of learning and performance in the mammalian life course. Although human interactions and other environmental conditions will remain critical, our findings suggest genomic breeding selection could help improve working dog populations.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

BACKGROUND

Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

RESULTS

We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

CONCLUSION

The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle

BACKGROUND

A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle.

RESULTS

Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits.

CONCLUSIONS

We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.

Stress Reactivity, Susceptibility to Hypertension, and Differential Expression of Genes in Hypertensive Compared to Normotensive Patients

Although half of hypertensive patients have hypertensive parents, known hypertension-related human loci identified by genome-wide analysis explain only 3% of hypertension heredity. Therefore, mainstream transcriptome profiling of hypertensive subjects addresses differentially expressed genes (DEGs) specific to gender, age, and comorbidities in accordance with predictive preventive personalized participatory medicine treating patients according to their symptoms, individual lifestyle, and genetic background. Within this mainstream paradigm, here, we determined whether, among the known hypertension-related DEGs that we could find, there is any genome-wide hypertension theranostic molecular marker applicable to everyone, everywhere, anytime. Therefore, we sequenced the hippocampal transcriptome of tame and aggressive rats, corresponding to low and high stress reactivity, an increase of which raises hypertensive risk; we identified stress-reactivity-related rat DEGs and compared them with their known homologous hypertension-related animal DEGs. This yielded significant correlations between stress reactivity-related and hypertension-related fold changes (log2 values) of these DEG homologs. We found principal components, PC1 and PC2, corresponding to a half-difference and half-sum of these log2 values. Using the DEGs of hypertensive versus normotensive patients (as the control), we verified the correlations and principal components. This analysis highlighted downregulation of β-protocadherins and hemoglobin as whole-genome hypertension theranostic molecular markers associated with a wide vascular inner diameter and low blood viscosity, respectively.

A bioinformatic search for correspondence between differentially expressed genes of domestic versus wild animals and orthologous human genes altering reproductive potential

One of the greatest achievements of genetics in the 20th century is D.K. Belyaev’s discovery of destabilizing selection during the domestication of animals and that this selection affects only gene expression regulation (not gene structure) and inf luences systems of neuroendocrine control of ontogenesis in a stressful environment. Among the experimental data generalized by Belyaev’s discovery, there are also f indings about accelerated extinc tion of testes’ hormonal function and disrupted seasonality of reproduction of domesticated foxes in comparison with their wild congeners. To date, Belyaev’s discovery has already been repeatedly conf irmed, for example, by independent observations during deer domestication, during the use of rats as laboratory animals, after the reintroduction of endangered species such as Przewalski’s horse, and during the creation of a Siberian reserve population of the Siberian grouse when it had reached an endangered status in natural habitats. A genome-wide comparison among humans, several domestic animals, and some of their wild congeners has given rise to the concept of self-domestication syndrome, which includes autism spectrum disorders. In our previous study, we created a bioinformatic model of human self-domestication syndrome using differentially expressed genes (DEGs; of domestic animals versus their wild congeners) orthologous to the human genes (mainly, nervous-system genes) whose changes in expression affect reproductive potential, i. e., growth of the number of humans in the absence of restrictions caused by limiting factors. Here, we applied this model to 68 human genes whose changes in expression alter the reproductive health of women and men and to 3080 DEGs of domestic versus wild animals. As a result, in domestic animals, we identif ied 16 and 4 DEGs, the expression changes of which are codirected with changes in the expression of the human orthologous genes decreasing and increasing human reproductive potential, respectively. The wild animals had 9 and 11 such DEGs, respectively. This difference between domestic and wild animals was signif icant according to Pearson’s χ2 test (p < 0.05) and Fisher’s exact test (p < 0.05). We discuss the results from the standpoint of restoration of endangered animal species whose natural habitats are subject to an anthropogenic impact.

Differential expression of 10 genes in the hypothalamus of two generations of rats selected for a reaction to humans

Individual behavioral differences are due to an interaction of the genotype and the environment. Phenotypic manifestation of aggressive behavior depends on the coordinated expression of gene ensembles. Nonetheless, the identification of these genes and of combinations of their mutual influence on expression remains a difficult task. Using animal models of aggressive behavior (gray rats that were selected for a reaction to humans; tame and aggressive rat strains), we evaluated the expression of 10 genes potentially associated with aggressiveness according to the literature: Cacna1b, Cacna2d3, Drd2, Egr1, Gad2, Gria2, Mapk1, Nos1, Pomc, and Syn1. To identify the genes most important for the manifestation of aggressiveness, we analyzed the expression of these genes in two generations of rats: 88th and 90th. Assessment of gene expression levels was carried out by real-time PCR in the hypothalamus of tame and aggressive rats. This analysis confirmed that 4 out of the 10 genes differ in expression levels between aggressive rats and tame rats in both generations. Specifically, it was shown that the expression of the Cacna1b, Drd2, Egr1, and Gad2 genes does not differ between the two generations (88th vs 90th) within each strain, but significantly differs between the strains: in the tame rats of both generations, the expression levels of these genes are significantly lower as compared to those in the aggressive rats. Therefore, these genes hold promise for further studies on behavioral characteristics. Thus, we confirmed polygenic causes of phenotypic manifestation of aggressive reactions.

Brain transcriptome analysis reveals gene expression differences associated with dispersal behaviour between range-front and range-core populations of invasive cane toads in Australia

Understanding the mechanisms allowing invasive species to adapt to novel environments is a challenge in invasion biology. Many invaders demonstrate rapid evolution of behavioural traits involved in range expansion such as locomotor activity, exploration and risk‐taking. However, the molecular mechanisms that underpin these changes are poorly understood. In 86 years, invasive cane toads (Rhinella marina) in Australia have drastically expanded their geographic range westward from coastal Queensland to Western Australia. During their range expansion, toads have undergone extensive phenotypic changes, particularly in behaviours that enhance the toads’ dispersal ability. Common‐garden experiments have shown that some changes in behavioural traits related to dispersal are heritable. At the molecular level, it is currently unknown whether these changes in dispersal‐related behaviour are underlain by small or large differences in gene expression, nor is known the biological function of genes showing differential expression. Here, we used RNA‐seq to gain a better understanding of the molecular mechanisms underlying dispersal‐related behavioural changes. We compared the brain transcriptomes of toads from the Hawai'ian source population, as well as three distinct populations from across the Australian invasive range. We found markedly different gene expression profiles between the source population and Australian toads. By contrast, toads from across the Australian invasive range had very similar transcriptomic profiles. Yet, key genes with functions putatively related to dispersal behaviour showed differential expression between populations located at each end of the invasive range. These genes could play an important role in the behavioural changes characteristic of range expansion in Australian cane toads.

Sexually Dimorphic Transcriptomic Changes of Developing Fetal Brain Reveal Signaling Pathways and Marker Genes of Brain Cells in Domestic Pigs

In this study, transcriptomic changes of the developing brain of pig fetuses of both sexes were investigated on gestation days (GD) 45, 60 and 90. Pig fetal brain grows rapidly around GD60. Consequently, gene expression of the fetal brain was distinctly different on GD90 compared to that of GD45 and GD60. In addition, varying numbers of differentially expressed genes (DEGs) were identified in the male brain compared to the female brain during development. The sex of adjacent fetuses also influenced gene expression of the fetal brain. Extensive changes in gene expression at the exon-level were observed during brain development. Pathway enrichment analysis showed that the ionotropic glutamate receptor pathway and p53 pathway were enriched in the female brain, whereas specific receptor-mediated signaling pathways were enriched in the male brain. Marker genes of neurons and astrocytes were significantly differentially expressed between male and female brains during development. Furthermore, comparative analysis of gene expression patterns between fetal brain and placenta suggested that genes related to ion transportation may play a key role in the regulation of the brain-placental axis in pig. Collectively, the study suggests potential application of pig models to better understand influence of fetal sex on brain development.

Domestication Explains Two-Thirds of Differential-Gene-Expression Variance between Domestic and Wild Animals; The Remaining One-Third Reflects Intraspecific and Interspecific Variation

Belyaev's concept of destabilizing selection during domestication was a major achievement in the XX century. Its practical value has been realized in commercial colors of the domesticated fox that never occur in the wild and has been confirmed in a wide variety of pet breeds. Many human disease models involving animals allow to test drugs before human testing. Perhaps this is why investigators doing transcriptomic profiling of domestic versus wild animals have searched for breed-specific patterns. Here we sequenced hypothalamic transcriptomes of tame and aggressive rats, identified their differentially expressed genes (DEGs), and, for the first time, applied principal component analysis to compare them with all the known DEGs of domestic versus wild animals that we could find. Two principal components, PC1 and PC2, respectively explained 67% and 33% of differential-gene-expression variance (hereinafter: log2 value) between domestic and wild animals. PC1 corresponded to multiple orthologous DEGs supported by homologs; these DEGs kept the log2 value sign from species to species and from tissue to tissue (i.e., a common domestication pattern). PC2 represented stand-alone homologous DEG pairs reversing the log2 value sign from one species to another and from tissue to tissue (i.e., representing intraspecific and interspecific variation).

Brain Transcriptomics of Wild and Domestic Rabbits Suggests That Changes in Dopamine Signaling and Ciliary Function Contributed to Evolution of Tameness

Domestication has resulted in immense phenotypic changes in animals despite their relatively short evolutionary history. The European rabbit is one of the most recently domesticated animals, but exhibits distinct morphological, physiological, and behavioral differences from their wild conspecifics. A previous study revealed that sequence variants with striking allele frequency differences between wild and domestic rabbits were enriched in conserved noncoding regions, in the vicinity of genes involved in nervous system development. This suggests that a large proportion of the genetic changes targeted by selection during domestication might affect gene regulation. Here, we generated RNA-sequencing data for four brain regions (amygdala, hypothalamus, hippocampus, and parietal/temporal cortex) sampled at birth and revealed hundreds of differentially expressed genes (DEGs) between wild and domestic rabbits. DEGs in amygdala were significantly enriched for genes associated with dopaminergic function and all 12 DEGs in this category showed higher expression in domestic rabbits. DEGs in hippocampus were enriched for genes associated with ciliary function, all 21 genes in this category showed lower expression in domestic rabbits. These results indicate an important role of dopamine signaling and ciliary function in the evolution of tameness during rabbit domestication. Our study shows that gene expression in specific pathways has been profoundly altered during domestication, but that the majority of genes showing differential expression in this study have not been the direct targets of selection.

Disruptive Selection of Human Immunostimulatory and Immunosuppressive Genes Both Provokes and Prevents Rheumatoid Arthritis, Respectively, as a Self-Domestication Syndrome

Using our previously published Web service SNP_TATA_Comparator, we conducted a genome-wide study of single-nucleotide polymorphisms (SNPs) within core promoters of 68 human rheumatoid arthritis (RA)-related genes. Using 603 SNPs within 25 genes clinically associated with RA-comorbid disorders, we predicted 84 and 70 candidate SNP markers for overexpression and underexpression of these genes, respectively, among which 58 and 96 candidate SNP markers, respectively, can relieve and worsen RA as if there is a neutral drift toward susceptibility to RA. Similarly, we predicted natural selection toward susceptibility to RA for 8 immunostimulatory genes (e.g., IL9R) and 10 genes most often associated with RA (e.g., NPY). On the contrary, using 25 immunosuppressive genes, we predicted 70 and 109 candidate SNP markers aggravating and relieving RA, respectively (e.g., IL1R2 and TGFB2), suggesting that natural selection can simultaneously additionally yield resistance to RA. We concluded that disruptive natural selection of human immunostimulatory and immunosuppressive genes is concurrently elevating and reducing the risk of RA, respectively. So, we hypothesize that RA in human could be a self-domestication syndrome referring to evolution patterns in domestic animals. We tested this hypothesis by means of public RNA-Seq data on 1740 differentially expressed genes (DEGs) of pets vs. wild animals (e.g., dogs vs. wolves). The number of DEGs in the domestic animals corresponding to worsened RA condition in humans was significantly larger than that in the related wild animals (10 vs. 3). Moreover, much less DEGs in the domestic animals were accordant to relieved RA condition in humans than those in the wild animals (1 vs. 8 genes). This indicates that the anthropogenic environment, in contrast to a natural one, affects gene expression across the whole genome (e.g., immunostimulatory and immunosuppressive genes) in a manner that likely contributes to RA. The difference in gene numbers is statistically significant as confirmed by binomial distribution (p < 0.01), Pearson's χ² (p < 0.01), and Fisher's exact test (p < 0.05). This allows us to propose RA as a candidate symptom within a self-domestication syndrome. Such syndrome might be considered as a human's payment with health for the benefits received during evolution.

Transcriptome Analyses Throughout Chili Pepper Fruit Development Reveal Novel Insights into the Domestication Process

Chili pepper (Capsicum spp.) is an important crop, as well as a model for fruit development studies and domestication. Here, we performed a time-course experiment to estimate standardized gene expression profiles with respect to fruit development for six domesticated and four wild chili pepper ancestors. We sampled the transcriptomes every 10 days from flowering to fruit maturity, and found that the mean standardized expression profiles for domesticated and wild accessions significantly differed. The mean standardized expression was higher and peaked earlier for domesticated vs. wild genotypes, particularly for genes involved in the cell cycle that ultimately control fruit size. We postulate that these gene expression changes are driven by selection pressures during domestication and show a robust network of cell cycle genes with a time shift in expression, which explains some of the differences between domesticated and wild phenotypes.

Farm Animals Are Long Away from Natural Behavior: Open Questions and Operative Consequences on Animal Welfare

Simple Summary

Animal welfare is a very important issue. One of the tasks of researchers is to provide explanations and possible solutions to questions arising from non-experts. This work analyzes part of the extensive literature on relationships between selection and domestic, mainly farm, animals’ behavior and deals with some very important themes, such as the role of regulations, domestication, and selection.

Abstract

The concept of welfare applied to farm animals has undergone a remarkable evolution. The growing awareness of citizens pushes farmers to guarantee the highest possible level of welfare to their animals. New perspectives could be opened for animal welfare reasoning around the concept of domestic, especially farm, animals as partial human artifacts. Therefore, it is important to understand how much a particular behavior of a farm animal is far from the natural one of its ancestors. This paper is a contribution to better understand the role of genetics of the farm animals on their behavior. This means that the naïve approach to animal welfare regarding returning animals to their natural state should be challenged and that welfare assessment should be considered.

A Bioinformatics Model of Human Diseases on the Basis of Differentially Expressed Genes (of Domestic Versus Wild Animals) That Are Orthologs of Human Genes Associated with Reproductive-Potential Changes

Earlier, after our bioinformatic analysis of single-nucleotide polymorphisms of TATA-binding protein-binding sites within gene promoters on the human Y chromosome, we suggested that human reproductive potential diminishes during self-domestication. Here, we implemented bioinformatics models of human diseases using animal in vivo genome-wide RNA-Seq data to compare the effect of co-directed changes in the expression of orthologous genes on human reproductive potential and during the divergence of domestic and wild animals from their nearest common ancestor (NCA). For example, serotonin receptor 3A (HTR3A) deficiency contributes to sudden death in pregnancy, consistently with Htr3a underexpression in guinea pigs (Cavia porcellus) during their divergence from their NCA with cavy (C. aperea). Overall, 25 and three differentially expressed genes (hereinafter, DEGs) in domestic animals versus 11 and 17 DEGs in wild animals show the direction consistent with human orthologous gene-markers of reduced and increased reproductive potential. This indicates a reliable association between DEGs in domestic animals and human orthologous genes reducing reproductive potential (Pearson's χ2 test p < 0.001, Fisher's exact test p < 0.05, binomial distribution p < 0.0001), whereas DEGs in wild animals uniformly match human orthologous genes decreasing and increasing human reproductive potential (p > 0.1; binomial distribution), thus enforcing the norm (wild type).

Searching where the treasure is: on the emergence of human companion animal partnership (HCAP)

In our view, the discipline, often referred to as human-animal interaction (HAI), lacks a well-defined conceptual framework. It is too narrow both with respect to the animal species investigated and the nature of human-animal interactions studied. So instead, we introduce the term human-companion animal partnership (HCAP) that is not only a better descriptor for most research efforts within HAI but also helps to direct research efforts on an ethological basis. In our approach, 'companion' is a function and not a feature of some species. This means that many species had and could have a potential to form mixed social groups with humans if they evolve some capacity of social competence. This view may initiate new comparative research involving a range of species to find out how complex social engagement could be maintained in such hetero-specific social groups based on evolutionary heritage, recent selection and individual experience (socialisation). Our approach emphasises the role of human caring behaviour and social competence in the emergence of a partnership with several species, and thus could also help in setting expectations for welfare and aid in designing artificial companions for specific purposes.

DNA methylation in canine brains is related to domestication and dog-breed formation

Epigenetic factors such as DNA methylation act as mediators in the interaction between genome and environment. Variation in the epigenome can both affect phenotype and be inherited, and epigenetics has been suggested to be an important factor in the evolutionary process. During domestication, dogs have evolved an unprecedented between-breed variation in morphology and behavior in an evolutionary short period. In the present study, we explore DNA methylation differences in brain, the most relevant tissue with respect to behavior, between wolf and dog breeds. We optimized a combined method of genotype-by-sequencing (GBS) and methylated DNA immunoprecipitation (MeDIP) for its application in canines. Genomic DNA from the frontal cortex of 38 dogs of 8 breeds and three wolves was used. GBS and GBS-MeDIP libraries were prepared and sequenced on Illuma HiSeq2500 platform. The reduced sample represented 1.18 ± 0.4% of the total dog genome (2,4 billion BP), while the GBS-MeDIP covered 11,250,788 ± 4,042,106 unique base pairs. We find substantial DNA methylation differences between wolf and dog and between the dog breeds. The methylation profiles of the different groups imply that epigenetic factors may have been important in the speciation from dog to wolf, but also in the divergence of different dog breeds. Specifically, we highlight methylation differences in genes related to behavior and morphology. We hypothesize that these differences are involved in the phenotypic variation found among dogs, whereas future studies will have to find the specific mechanisms. Our results not only add an intriguing new dimension to dog breeding but are also useful to further understanding of epigenetic involvement.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

Independent domestication events in the blue-cheese fungus Penicillium roqueforti

Domestication provides an excellent framework for studying adaptive divergence. Using population genomics and phenotypic assays, we reconstructed the domestication history of the blue cheese mould Penicillium roqueforti. We showed that this fungus was domesticated twice independently. The population used in Roquefort originated from an old domestication event associated with weak bottlenecks and exhibited traits beneficial for pre‐industrial cheese production (slower growth in cheese and greater spore production on bread, the traditional multiplication medium). The other cheese population originated more recently from the selection of a single clonal lineage, was associated with all types of blue cheese worldwide except Roquefort, and displayed phenotypes more suited for industrial cheese production (high lipolytic activity, efficient cheese cavity colonization ability and salt tolerance). We detected genomic regions affected by recent positive selection and putative horizontal gene transfers. This study sheds light on the processes of rapid adaptation and raises questions about genetic resource conservation.

see also the Perspective by Brigida Gallone, Jan Steensels and Kevin J. Verstrepen.

Hypothalamic transcriptome of tame and aggressive silver foxes (Vulpes vulpes) identifies gene expression differences shared across brain regions

The underlying neurological events accompanying dog domestication remain elusive. To reconstruct the domestication process in an experimental setting, silver foxes (Vulpes vulpes) have been deliberately bred for tame vs aggressive behaviors for more than 50 generations at the Institute for Cytology and Genetics in Novosibirsk, Russia. The hypothalamus is an essential part of the hypothalamic-pituitary-adrenal axis and regulates the fight-or-flight response, and thus, we hypothesized that selective breeding for tameness/aggressiveness has shaped the hypothalamic transcriptomic profile. RNA-seq analysis identified 70 differentially expressed genes (DEGs). Seven of these genes, DKKL1, FBLN7, NPL, PRIMPOL, PTGRN, SHCBP1L and SKIV2L, showed the same direction expression differences in the hypothalamus, basal forebrain and prefrontal cortex. The genes differentially expressed across the three tissues are involved in cell division, differentiation, adhesion and carbohydrate processing, suggesting an association of these processes with selective breeding. Additionally, 159 transcripts from the hypothalamus demonstrated differences in the abundance of alternative spliced forms between the tame and aggressive foxes. Weighted gene coexpression network analyses also suggested that gene modules in hypothalamus were significantly associated with tame vs aggressive behavior. Pathways associated with these modules include signal transduction, interleukin signaling, cytokine-cytokine receptor interaction and peptide ligand-binding receptors (eg, G-protein coupled receptor [GPCR] ligand binding). Current studies show the selection for tameness vs aggressiveness in foxes is associated with unique hypothalamic gene profiles partly shared with other brain regions and highlight DEGs involved in biological processes such as development, differentiation and immunological responses. The role of these processes in fox and dog domestication remains to be determined.

A Rat Model of Human Behavior Provides Evidence of Natural Selection Against Underexpression of Aggressiveness-Related Genes in Humans

Aggressiveness is a hereditary behavioral pattern that forms a social hierarchy and affects the individual social rank and accordingly quality and duration of life. Thus, genome-wide studies of human aggressiveness are important. Nonetheless, the aggressiveness-related genome-wide studies have been conducted on animals rather than humans. Recently, in our genome-wide study, we uncovered natural selection against underexpression of human aggressiveness-related genes and proved it using F1 hybrid mice. Simultaneously, this natural selection equally supports two opposing traits in humans (dominance and subordination) as if self-domestication could have happened with its disruptive natural selection. Because there is still not enough scientific evidence that this could happen, here, we verified this natural selection pattern using quantitative PCR and two outbred rat lines (70 generations of artificial selection for aggressiveness or tameness, hereinafter: domestication). We chose seven genes—Cacna2d3, Gad2, Gria2, Mapk1, Nos1, Pomc, and Syn1—over- or underexpression of which corresponds to aggressive or domesticated behavior (in humans or mice) that has the same direction as natural selection. Comparing aggressive male rats with domesticated ones, we found that these genes are overexpressed statistically significantly in the hypothalamus (as a universal behavior regulator), not in the periaqueductal gray, where there was no aggressiveness-related expression of the genes in males. Database STRING showed statistically significant associations of the human genes homologous to these rat genes with long-term depression, circadian entrainment, Alzheimer’s disease, and the central nervous system disorders during chronic IL-6 overexpression. This finding more likely supports positive perspectives of further studies on self-domestication syndromes.

The Domestic BCO2 Allele Buffers Low-Carotenoid Diets in Chickens: Possible Fitness Increase Through Species Hybridization

Domestic animals are adapted to conditions vastly different from those of their wild ancestors, and this is particularly true for their diets. The most numerous of all domestic species, the chicken, originated from the Red Junglefowl (RJF), a native of subtropical forests in Southeast Asia. Surprisingly however, in domestic chicken breeds, a common haplotype of the β-carotene oxygenase 2 (BCO2) gene, which is involved in carotenoid metabolism, is introgressed from a related species, the Gray Junglefowl, and has been under strong selective pressure during domestication. This suggests that a hybridization event may have conferred a fitness advantage on chickens carrying the derived allele. To investigate the possible biological function of the introgressed BCO2 allele in chicken, we introgressed the ancestral BCO2 allele into domestic White Leghorn chickens. We measured gene expression as well as carotenoid accumulation in skin and eggs of chickens carrying either the ancestral or the derived BCO2 allele. The derived haplotype was associated with down-regulation of BCO2 in skin, muscle, and adipose tissue, but not in liver or duodenum, indicating that carotenoid accumulation occurred in the tissues with reduced gene expression. Most importantly, we found that hens with the derived BCO2 genotype were capable of allocating stored carotenoids to their eggs, suggesting a functional benefit through buffering any shortage in the diet during egg production. Nevertheless, it is of interest that loss of function mutations in BCO2 gene are prevalent in other domesticates including cows, rabbits, and sheep, and, given the importance of carotenoids in development, reproduction, and immunity, it is possible that derived BCO2 alleles may provide a general mechanism in multiple domestic species to deal with higher demand for carotenoids in an environment with carotenoid shortage in the diet.

Decrease of gene expression diversity during domestication of animals and plants

BACKGROUND

The genetic mechanisms underlying the domestication of animals and plants have been of great interest to biologists since Darwin. To date, little is known about the global pattern of gene expression changes during domestication.

RESULTS

We generated and collected transcriptome data for seven pairs of domestic animals and plants including dog, silkworm, chicken, rice, cotton, soybean and maize and their wild progenitors and compared the expression profiles between the domestic and wild species. Intriguingly, although the number of expressed genes varied little, the domestic species generally exhibited lower gene expression diversity than did the wild species, and this lower diversity was observed for both domestic plants and different kinds of domestic animals including insect, bird and mammal in the whole-genome gene set (WGGS), candidate selected gene set (CSGS) and non-CSGS, with CSGS exhibiting a higher degree of decreased expression diversity. Moreover, different from previous reports which found 2 to 4% of genes were selected by human, we identified 6892 candidate selected genes accounting for 7.57% of the whole-genome genes in rice and revealed that fewer than 8% of the whole-genome genes had been affected by domestication.

CONCLUSIONS

Our results showed that domestication affected the pattern of variation in gene expression throughout the genome and generally decreased the expression diversity across species, and this decrease may have been associated with decreased genetic diversity. This pattern might have profound effects on the phenotypic and physiological changes of domestic animals and plants and provide insights into the genetic mechanisms at the transcriptome level other than decreased genetic diversity and increased linkage disequilibrium underpinning artificial selection.

Genomic responses to selection for tame/aggressive behaviors in the silver fox (Vulpes vulpes)

Animal domestication efforts have led to a shared spectrum of striking behavioral and morphological changes. To recapitulate this process, silver foxes have been selectively bred for tame and aggressive behaviors for more than 50 generations at the Institute for Cytology and Genetics in Novosibirsk, Russia. To understand the genetic basis and molecular mechanisms underlying the phenotypic changes, we profiled gene expression levels and coding SNP allele frequencies in two brain tissue specimens from 12 aggressive foxes and 12 tame foxes. Expression analysis revealed 146 genes in the prefrontal cortex and 33 genes in the basal forebrain that were differentially expressed, with a 5% false discovery rate (FDR). These candidates include genes in key pathways known to be critical to neurologic processing, including the serotonin and glutamate receptor pathways. In addition, 295 of the 31,000 exonic SNPs show significant allele frequency differences between the tame and aggressive populations (1% FDR), including genes with a role in neural crest cell fate determination.

Rapid Evolution of Genes Involved in Learning and Energy Metabolism for Domestication of the Laboratory Rat

The laboratory rat, widely used in biomedical research, is domesticated from wild brown rat. The origin and genetic mechanism underlying domestication of the laboratory rat remain largely elusive. In the present study, large scale genomes supported a single origin for the laboratory rat, possibly from a sister group to wild rats from Europe/Africa/Middle East. Genomic and transcriptomic analyses uncovered many artificially selected genes (e.g., FOXP2, B3GAT1, and CLOCK) involved in the nervous system. These genes associate with learning ability and regulation of circadian rhythm, which likely enabled the successful domestication of the laboratory rat. Particularly, many genes, including mitochondrial genes responsible for energy metabolism, displayed a substantially increased expression in the brain of laboratory rats compared with wild rats. Our findings demystify the origin and evolution of this model animal, and provide insight into the process of its domestication.

Changes in pituitary gene expression may underlie multiple domesticated traits in chickens

Domesticated animals share a unique set of morphological and behavioral traits, jointly referred to as the domesticated phenotype. Striking similarities amongst a range of unrelated domesticated species suggest that similar regulatory mechanisms may underlie the domesticated phenotype. These include color pattern, growth, reproduction, development and stress response. Although previous studies have focused on the brain to find mechanisms underlying domestication, the potential role of the pituitary gland as a target of domestication is highly overlooked. Here, we study gene expression in the pituitary gland of the domesticated White Leghorn chicken and its wild ancestor, the Red Junglefowl. By overlapping differentially expressed genes with a previously published list of functionally important genes in the pituitary gland, we narrowed down to 34 genes. Amongst them, expression levels of genes with inhibitory function on pigmentation (ASIP), main stimulators of metabolism and sexual maturity (TSHB and DIO2), and a potential inhibitor of broodiness (PRLR), were higher in the domesticated breed. Additionally, expression of 2 key inhibitors of the stress response (NR3C1, CRHR2) was higher in the domesticated breed. We suggest that changes in the transcription of important modulatory genes in the pituitary gland can account not only for domestication of the stress response in domestic chickens, but also for changes in pigmentation, development, and reproduction. Given the pivotal role of the pituitary gland in the regulation of multiple shared domesticated traits, we suggest that similar changes in pituitary transcriptome may contribute to the domesticated phenotype in other species as well.

Neurogenomic Profiling Reveals Distinct Gene Expression Profiles Between Brain Parts That Are Consistent in Ophthalmotilapia Cichlids

The detection of external and internal cues alters gene expression in the brain which in turn may affect neural networks that underly behavioral responses. Previous studies have shown that gene expression profiles differ between major brain regions within individuals and between species with different morphologies, cognitive abilities and/or behaviors. A detailed description of gene expression in all macroanatomical brain regions and in species with similar morphologies and behaviors is however lacking. Here, we dissected the brain of two cichlid species into six macroanatomical regions. Ophthalmotilapia nasuta and O. ventralis have similar morphology and behavior and occasionally hybridize in the wild. We use 3′ mRNA sequencing and a stage-wise statistical testing procedure to identify differential gene expression between females that were kept in a social setting with other females. Our results show that gene expression differs substantially between all six brain parts within species: out of 11,577 assessed genes, 8,748 are differentially expressed (DE) in at least one brain part compared to the average expression of the other brain parts. At most 16% of these DE genes have |log₂FC| significantly higher than two. Functional differences between brain parts were consistent between species. The majority (61–79%) of genes that are DE in a particular brain part were shared between both species. Only 32 genes show significant differences in fold change across brain parts between species. These genes are mainly linked to transport, transmembrane transport, transcription (and its regulation) and signal transduction. Moreover, statistical equivalence testing reveals that within each comparison, on average 89% of the genes show an equivalent fold change between both species. The pronounced differences in gene expression between brain parts and the conserved patterns between closely related species with similar morphologies and behavior suggest that unraveling the interactions between genes and behavior will benefit from neurogenomic profiling of distinct brain regions.