Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviours

Molecular characterization and complexity of the immunoglobulin repertoire in the silver-black fox (Vulpes vulpes)

Immunoglobulins, a class of globulins with antibody properties, play a crucial role in the body's defense against pathogens. In this study, we analyzed the gene loci structure of Silver-black fox using a comparative genomics approach. The mechanisms of expression diversity and its preferences were investigated through Next-generation sequencing (NGS). The results revealed 32 potentially functional VH genes, 9 Vκ genes, 17 Vλ genes, 9 DH genes, 3 JH genes, 6 Jκ genes, and 11 Jλ genes, located on different scaffolds. Subsequently, 5'RACE and PE300 bipartite sequencing were used to obtain the reads of the expressed antibody repertoire of Silver-black foxes gene rearrangement events. The analysis indicated a strong preference in the use of V genes, DH genes and J genes by Silver-black fox. The main ways of expression diversity were V(D)J recombination and somatic hypermutation (SHM). The hypermutated region of SHM was not only concentrated in the CDR region but also had higher mutation rate in the FR region. The main types of SHM were G > A, C > T, T > C, and A > G. The findings of this study could serve as a theoretical foundation for a deeper understanding of Silver-black fox immunoglobulins, which is significant for enriching knowledge in immunogenetics and providing theoretical support for future studies on vaccine design for the Silver-black fox.

Captive environments reshape the compositions of carbohydrate active enzymes and virulence factors in wolf gut microbiome

Species in the family Canidae occupy different spatial ecological niches, and some (e.g., wolf) can be kept in zoos. The gut microbiome may differ among various wild and captive canids. Therefore, we compared the gut microbiomes of wild canids (wolf, red fox, and corsac fox) in the Hulun Lake area, captive wolves, and domestic dogs in different regions using metagenomic data. A random forest analysis revealed significant enrichment for bacterial species producing short-chain fatty acids and the thermogenesis pathway (ko04714) in the gut microbiome of wild wolf, potentially providing sufficient energy for adaptation to a wide range of spatial ecological niches. The significantly enriched bacterial species and functional pathways in the gut microbiome of corsac foxes were related to physiological stability and adaptation to arid environments. Alpha diversity of carbohydrate-active enzymes in the gut microbiome was higher in the red fox than in the corsac fox and wild wolf, which may be related to the abundance of plant seeds (containing carbohydrates) in their diets (red foxes inhabit seed-rich willow bosk habitats). However, the influence of host genetic factors cannot be excluded, and further experimental studies are needed to verify the study results. In addition, captive environments drove similarity in carbohydrate-active enzymes (CAZymes) and virulence factors (VFs) in the gut microbiomes of captive wolf and domestic dog, and increased the diversity of CAZymes and VFs in the gut microbiome of captive wolf. Increased VFs diversity may increase the pathogenic potential of the gut microbiome in captive wolves. Therefore, it is necessary to continue monitoring the health status of captive wolves and develop appropriate management strategies.

Evolutionary consequences of domestication on the selective effects of new amino acid changing mutations in canids

The domestication of wild canids led to dogs no longer living in the wild but instead residing alongside humans. Extreme changes in behavior and diet associated with domestication may have led to the relaxation of the selective pressure on traits that may be less important in the domesticated context. Thus, here we hypothesize that strongly deleterious mutations may have become less deleterious in domesticated populations. We test this hypothesis by estimating the distribution of fitness effects (DFE) for new amino acid changing mutations using whole-genome sequence data from 24 gray wolves and 61 breed dogs. We find that the DFE is strikingly similar across canids, with 26-28% of new amino acid changing mutations being neutral/nearly neutral (|s| < 1e-5), and 41-48% under strong purifying selection (|s| > 1e-2). Our results are robust to different model assumptions suggesting that the DFE is stable across short evolutionary timescales, even in the face of putative drastic changes in the selective pressure caused by artificial selection during domestication and breed formation. On par with previous works describing DFE evolution, our data indicate that the DFE of amino acid changing mutations depends more strongly on genome structure and organismal characteristics, and less so on shifting selective pressures or environmental factors. Given the constant DFE and previous data showing that genetic variants that differentiate wolf and dog populations are enriched in regulatory elements, we speculate that domestication may have had a larger impact on regulatory variation than on amino acid changing mutations.

Novel pathways linked to the expression of temperament in Merino sheep: a genome-wide association study

Animal temperament refers to the inherent behavioural and emotional characteristics of an animal, influencing how it interacts with its environment. The selection of sheep for temperament can change the temperament traits of the selected line and improve the welfare and production (reproduction, growth, immunity) of those animals. To understand the genetics that underly variation in temperament in sheep, and how selection on temperament can affect other production traits, a genome-wide association study was carried out. Merino sheep from lines selected for traits of calm and nervous temperament, and a commercial population, on which the temperament traits had never been assessed, were used. Blood samples from the three populations were genotyped using an Illumina GGP Ovine 50 K Genotyping BeadChip. The calm and nervous populations in the selected lines presented as distinct genetic populations, and 2 729 of the 45 761 single nucleotide polymorphisms (SNPs) had significantly different proportions between the two lines. Of those 2 729 SNPs, 2 084 were also associated with temperament traits in the commercial population. A genomic annotation identified 81 candidate genes for temperament, nearly half of which are associated with disorders of social behaviour in humans. Five of those 81 candidate genes are related to production traits in sheep. Two genes were associated with personality disorders in humans and with production traits in sheep. We identified significant enrichment in genes involved in nervous system processes such as the regulation of systemic arterial blood pressure, ventricular myocyte action, multicellular organismal signalling, ion transmembrane transport, and calcium ion binding, suggesting that temperament is underpinned by variation in multiple biological systems. Our results contribute to understanding of the genetic basis of animal temperament which could be applied to the genetic evaluation of temperament in sheep and other farm animals.

Strongly deleterious mutations influence reproductive output and longevity in an endangered population

Inbreeding depression has been documented in various fitness traits in a wide range of species and taxa, however, the mutational basis is not yet well understood. We investigate how putatively deleterious variation influences fitness and is shaped by individual ancestry by re-sequencing complete genomes of 37 individuals in a natural arctic fox (Vulpes lagopus) population subjected to both inbreeding depression and genetic rescue. We find that individuals with high proportion of homozygous loss of function genotypes (LoFs), which are predicted to exert a strong effect on fitness, generally have lower lifetime reproductive success and live shorter lives compared with individuals with lower proportion of LoFs. We also find that juvenile survival is negatively associated with the proportion of homozygous missense genotypes and positively associated with genome wide heterozygosity. Our results demonstrate that homozygosity of strongly and moderately deleterious mutations can be an important cause of trait specific inbreeding depression in wild populations, and mark an important step towards making more informed decisions using applied conservation genetics.

Missing history of a modern domesticate: Historical demographics and genetic diversity in farm-bred red fox populations

The first record of captive-bred red foxes (Vulpes vulpes) dates to 1896 when a breeding enterprise emerged in the provinces of Atlantic Canada. Because its domestication happened during recent history, the red fox offers a unique opportunity to examine the genetic diversity of an emerging domesticated species in the context of documented historical and economic influences. In particular, the historical record suggests that North American and Eurasian farm-bred populations likely experienced different demographic trajectories. Here, we focus on the likely impacts of founder effects and genetic drift given historical trends in fox farming on North American and Eurasian farms. A total of 15 mitochondrial haplotypes were identified in 369 foxes from 10 farm populations that we genotyped (n = 161) or that were previously published. All haplotypes are endemic to North America. Although most haplotypes were consistent with eastern Canadian ancestry, a small number of foxes carried haplotypes typically found in Alaska and other regions of western North America. The presence of these haplotypes supports historical reports of wild foxes outside of Atlantic Canada being introduced into the breeding stock. These putative Alaskan and Western haplotypes were more frequently identified in Eurasian farms compared to North American farms, consistent with historical documentation suggesting that Eurasian economic and breeding practices were likely to maintain low-frequency haplotypes more effectively than in North America. Contextualizing inter- vs. intra-farm genetic diversity alongside the historical record is critical to understanding the origins of this emerging domesticate and the relationships between wild and farm-bred fox populations.

Allele frequency dynamics under sex-biased demography and sex-specific inheritance in a pedigreed jay population

Sex-biased demography, including sex-biased survival or migration, can alter allele frequency changes across the genome. In particular, we can expect different patterns of genetic variation on autosomes and sex chromosomes due to sex-specific differences in life histories, as well as differences in effective population size, transmission modes, and the strength and mode of selection. Here, we demonstrate the role that sex differences in life history played in shaping short-term evolutionary dynamics across the genome. We used a 25-year pedigree and genomic dataset from a long-studied population of Florida Scrub-Jays (Aphelocoma coerulescens) to directly characterize the relative roles of sex-biased demography and inheritance in shaping genome-wide allele frequency trajectories. We used gene dropping simulations to estimate individual genetic contributions to future generations and to model drift and immigration on the known pedigree. We quantified differential expected genetic contributions of males and females over time, showing the impact of sex-biased dispersal in a monogamous system. Due to female-biased dispersal, more autosomal variation is introduced by female immigrants. However, due to male-biased transmission, more Z variation is introduced by male immigrants. Finally, we partitioned the proportion of variance in allele frequency change through time due to male and female contributions. Overall, most allele frequency change is due to variance in survival and births. Males and females make similar contributions to autosomal allele frequency change, but males make higher contributions to allele frequency change on the Z chromosome. Our work shows the importance of understanding sex-specific demographic processes in characterizing genome-wide allele frequency change in wild populations.

Duplications and Retrogenes Are Numerous and Widespread in Modern Canine Genomic Assemblies

Recent years have seen a dramatic increase in the number of canine genome assemblies available. Duplications are an important source of evolutionary novelty and are also prone to misassembly. We explored the duplication content of nine canine genome assemblies using both genome self-alignment and read-depth approaches. We find that 8.58% of the genome is duplicated in the canFam4 assembly, derived from the German Shepherd Dog Mischka, including 90.15% of unplaced contigs. Highlighting the continued difficulty in properly assembling duplications, less than half of read-depth and assembly alignment duplications overlap, but the mCanLor1.2 Greenland wolf assembly shows greater concordance. Further study shows the presence of multiple segments that have alignments to four or more duplicate copies. These high-recurrence duplications correspond to gene retrocopies. We identified 3,892 candidate retrocopies from 1,316 parental genes in the canFam4 assembly and find that ∼8.82% of duplicated base pairs involve a retrocopy, confirming this mechanism as a major driver of gene duplication in canines. Similar patterns are found across eight other recent canine genome assemblies, with metrics supporting a greater quality of the PacBio HiFi mCanLor1.2 assembly. Comparison between the wolf and other canine assemblies found that 92% of retrocopy insertions are shared between assemblies. By calculating the number of generations since genome divergence, we estimate that new retrocopy insertions appear, on average, in 1 out of 3,514 births. Our analyses illustrate the impact of retrogene formation on canine genomes and highlight the variable representation of duplicated sequences among recently completed canine assemblies.

First Immunohistochemical Demonstration of the Expression of a Type-2 Vomeronasal Receptor, V2R2, in Wild Canids

The mammalian vomeronasal system enables the perception of chemical signals crucial for social communication via the receptor families V1R and V2R. These receptors are linked with the G-protein subunits, Gαi2 and Gαo, respectively. Exploring the evolutionary pathways of V1Rs and V2Rs across mammalian species remains a significant challenge, particularly when comparing genomic data with emerging immunohistochemical evidence. Recent studies have revealed the expression of Gαo in the vomeronasal neuroepithelium of wild canids, including wolves and foxes, contradicting predictions based on current genomic annotations. Our study provides detailed immunohistochemical evidence, mapping the expression of V2R receptors in the vomeronasal sensory epithelium, focusing particularly on wild canids, specifically wolves and foxes. An additional objective involves contrasting these findings with those from domestic species like dogs to highlight the evolutionary impacts of domestication on sensory systems. The employment of a specific antibody raised against the mouse V2R2, a member of the C-family of vomeronasal receptors, V2Rs, has confirmed the presence of V2R2-immunoreactivity (V2R2-ir) in the fox and wolf, but it has revealed the lack of expression in the dog. This may reflect the impact of domestication on the regression of the VNS in this species, in contrast to their wild counterparts, and it underscores the effects of artificial selection on sensory functions. Thus, these findings suggest a more refined chemical detection capability in wild species.

Maintenance of a narrow hybrid zone between native and introduced red foxes (Vulpes vulpes) despite conspecificity and high dispersal capabilities

Human-facilitated introductions of nonnative populations can lead to secondary contact between allopatric lineages, resulting in lineage homogenisation or the formation of stable hybrid zones maintained by reproductive barriers. We investigated patterns of gene flow between the native Sacramento Valley red fox (Vulpes vulpes patwin) and introduced conspecifics of captive-bred origin in California's Central Valley. Considering their recent divergence (20-70 kya), we hypothesised that any observed barriers to gene flow were primarily driven by pre-zygotic (e.g. behavioural differences) rather than post-zygotic (e.g. reduced hybrid fitness) barriers. We also explored whether nonnative genes could confer higher fitness in the human-dominated landscape resulting in selective introgression into the native population. Genetic analysis of red foxes (n = 682) at both mitochondrial (cytochrome b + D-loop) and nuclear (19,051 SNPs) loci revealed narrower cline widths than expected under a simulated model of unrestricted gene flow, consistent with the existence of reproductive barriers. We identified several loci with reduced introgression that were previously linked to behavioural divergence in captive-bred and domestic canids, supporting pre-zygotic, yet possibly hereditary, barriers as a mechanism driving the narrowness and stability of the hybrid zone. Several loci with elevated gene flow from the nonnative into the native population were linked to genes associated with domestication and adaptation to human-dominated landscapes. This study contributes to our understanding of hybridisation dynamics in vertebrates, particularly in the context of species introductions and landscape changes, underscoring the importance of considering how multiple mechanisms may be maintaining lineages at the species and subspecies level.

Lipidomic and transcriptomic profiles provide new insights into the triacylglycerol and glucose handling capacities of the Arctic fox

The Arctic fox (Vulpes lagopus) is a species indigenous to the Arctic and has developed unique lipid metabolism, but the mechanisms remain unclear. Here, the significantly increased body weight of Arctic foxes was consistent with the significantly increased serum very-low-density lipoprotein (VLDL), and the 40% crude fat diet further increased the Arctic fox body weight. The enhanced body weight gain stems primarily from increased subcutaneous adipose tissue accumulation. The adipose triacylglycerol and phosphatidylethanolamine were significantly greater in Arctic foxes. The adipose fatty-acid synthase content was significantly lower in Arctic foxes, highlighting the main role of exogenous fatty-acids in fat accumulation. Considering the same diet, liver-derived fat dominates adipose expansion in Arctic foxes. Liver transcriptome analysis revealed greater fat and VLDL synthesis in Arctic foxes, consistent with the greater VLDL. Glucose homeostasis wasn't impacted in Arctic foxes. And the free fatty-acids in adipose, which promote insulin resistance, also did not differ between groups. However, the hepatic glycogen was greater in Arctic foxes and transcriptome analysis revealed upregulated glycogen synthesis, improving glucose homeostasis. These results suggest that the superior fat accumulation capacity and distinct characteristics of hepatic and adipose lipid and glucose metabolism facilitate glucose homeostasis and massive fat accumulation in Arctic foxes.

The expression profile of genes associated with behavior, stress, and adult neurogenesis along the hippocampal dorsoventral axis in tame and aggressive foxes

The hippocampus plays the key role in stress response regulation, and stress response appears to be weakened in domesticated animals compared to their wild relatives. The hippocampus is functionally heterogeneous along its dorsoventral axis, with its ventral compartment being more closely involved in stress regulation. An earlier series of experiments was conducted with a unique breeding model of animal domestication, the farm silver fox (Vulpes vulpes), which included tame, aggressive, and unselected animals. A decrease in many indices of the hypothalamic-pituitary-adrenal activity was observed in tame animals. Also, adult hippocampal neurogenesis was more intense in tame foxes, and this fact may relate to reduced stress levels in this experimental population of foxes. Nevertheless, the molecular mechanisms responsible for the reduced stress response in tame animals remain obscure. In this study, serum cortisol levels and the mRNA levels of 13 genes in the dorsal and ventral hippocampus have been measured and compared in tame, aggressive, and unselected foxes. At the current stage of domestication, stress-induced cortisol levels in tame, aggressive, and unselected animals differ significantly from each other: tame foxes show the lowest levels, and aggressive ones, the highest. Twelve genes tested demonstrate significant gene expression differences between the dorsal and ventral hippocampi. These differences are mainly consistent with those found in rodents and humans. In tame foxes, significantly elevated mRNA levels were recorded for several genes: CYP26B1 for cytochrome P450 26B1 and ADRA1A for α1A adrenergic receptor in the dorsal hippocampus, whereas the level of NR3C2 mRNA for mineralocorticoid receptor was higher in the ventral. It is presumed that these genes constitute an important part of the mechanism reducing stress induced by contacts with humans and contribute to linking stress regulation with adult neurogenesis in tame foxes and domesticated animals in general.

Evolutionary origin of genomic structural variations in domestic yaks

Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.

Immunohistological study of the unexplored vomeronasal organ of an endangered mammal, the dama gazelle (Nanger dama)

Dama gazelle is a threatened and rarely studied species found primarily in northern Africa. Human pressure has depleted the dama gazelle population from tens of thousands to a few hundred individuals. Since 1970, a founder population consisting of the last 17 surviving individuals in Western Sahara has been maintained in captivity, reproducing naturally. In preparation for the future implementation of assisted reproductive technology, certain aspects of dama gazelle reproductive biology have been established. However, the role played by semiochemical-mediated communications in the sexual behavior of dama gazelle remains unknown due partially to a lack of a neuroanatomical or morphofunctional characterization of the dama gazelle vomeronasal organ (VNO), which is the sensory organ responsible for pheromone processing. The present study characterized the dama gazelle VNO, which appears fully equipped to perform neurosensory functions, contributing to current understanding of interspecies VNO variability among ruminants. By employing histological, lectin-histochemical, and immunohistochemical techniques, we conducted a detailed morphofunctional evaluation of the dama gazelle VNO along its entire longitudinal axis. Our findings of significant structural and neurochemical transformation along the entire VNO suggest that future studies of the VNO should take a similar approach. The present study contributes to current understanding of dama gazelle VNO, providing a basis for future studies of semiochemical-mediated communications and reproductive management in this species. RESEARCH HIGHLIGHTS: This exhaustive immunohistological study of the vomeronasal organ (VNO) of the dama gazelle provides the first evidence of notable differences in the expression of neuronal markers along the rostrocaudal axis of the VNO. This provides a morphological basis for the implementation of pheromones in captive populations of dama gazelle.

The Shepherd and the Hunter: A Genomic Comparison of Italian Dog Breeds

Shepherd and hunting dogs have undergone divergent selection for specific tasks, resulting in distinct phenotypic and behavioural differences. Italy is home to numerous recognized and unrecognized breeds of both types, providing an opportunity to compare them genomically. In this study, we analysed SNP data obtained from the CanineHD BeadChip, encompassing 116 hunting dogs (representing 6 breeds) and 158 shepherd dogs (representing 9 breeds). We explored the population structure, genomic background, and phylogenetic relationships among the breeds. To compare the two groups, we employed three complementary methods for selection signature detection: FST, XP-EHH, and ROH. Our results reveal a clear differentiation between shepherd and hunting dogs as well as between gun dogs vs. hounds and guardian vs. herding shepherd dogs. The genomic regions distinguishing these groups harbour several genes associated with domestication and behavioural traits, including gregariousness (WBSRC17) and aggressiveness (CDH12 and HTT). Additionally, genes related to morphology, such as size and coat colour (ASIP and TYRP1) and texture (RSPO2), were identified. This comparative genomic analysis sheds light on the genetic underpinnings of the phenotypic and behavioural variations observed in Italian hunting and shepherd dogs.

A missense mutation in RRM1 contributes to animal tameness

The increased tameness to reduce avoidance of human in wild animals has been long proposed as the key step of animal domestication. The tameness is a complex behavior trait and largely determined by genetic factors. However, the underlying genetic mutations remain vague and how they influence the animal behaviors is yet to be explored. Behavior tests of a wild-domestic hybrid goat population indicate the locus under strongest artificial selection during domestication may exert a huge effect on the flight distance. Within this locus, only one missense mutation RRM1I241V which was present in the early domestic goat ~6500 years ago. Genome editing of RRM1I241V in mice showed increased tameness and sociability and reduced anxiety. These behavioral changes induced by RRM1I241V were modulated by the alternation of activity of glutamatergic synapse and some other synapse-related pathways. This study established a link between RRM1I241V and tameness, demonstrating that the complex behavioral change can be achieved by mutations under strong selection during animal domestication.

Defining hippocampal area CA2 in the fox (Vulpes vulpes) brain

Since 1959, the Russian Farm-Fox study has bred foxes to be either tame or, more recently, aggressive, and scientists have used them to gain insight into the brain structures associated with these behavioral features. In mice, hippocampal area CA2 has emerged as one of the essential regulators of social aggression, and so to eventually determine whether we could identify differences in CA2 between tame and aggressive foxes, we first sought to identify CA2 in foxes (Vulpes vulpes). As no clearly defined area of CA2 has been described in species such as cats, dogs, or pigs, it was not at all clear whether CA2 could be identified in foxes. In this study, we cut sections of temporal lobes from male and female red foxes, perpendicular to the long axis of the hippocampus, and stained them with markers of CA2 pyramidal cells commonly used in tissue from rats and mice. We observed that antibodies against Purkinje cell protein 4 best stained the pyramidal cells in the area spanning the end of the mossy fibers and the beginning of the pyramidal cells lacking mossy fibers, resembling the pattern seen in rats and mice. Our findings indicate that foxes do have a "molecularly defined" CA2, and further, they suggest that other carnivores like dogs and cats might as well. With this being the case, these foxes could be useful in future studies looking at CA2 as it relates to aggression.

The brain of the silver fox (Vulpes vulpes): a neuroanatomical reference of cell-stained histological and MRI images

Although the silver fox (Vulpes vulpes) has been largely overlooked by neuroscientists, it has the potential to serve as a powerful model for the investigation of brain-behavior relationships. The silver fox is a melanistic variant of the red fox. Within this species, the long-running Russian farm-fox experiment has resulted in different strains bred to show divergent behavior. Strains bred for tameness, aggression, or without selection on behavior present an excellent opportunity to investigate neuroanatomical changes underlying behavioral characteristics. Here, we present a histological and MRI neuroanatomical reference of a fox from the conventional strain, which is bred without behavioral selection. This can provide an anatomical basis for future studies of the brains of foxes from this particular experiment, as well as contribute to an understanding of fox brains in general. In addition, this can serve as a resource for comparative neuroscience and investigations into neuroanatomical variation among the family Canidae, the order Carnivora, and mammals more broadly.

High-Quality Chromosome-Level Genome Assembly of the Corsac Fox (Vulpes corsac) Reveals Adaptation to Semiarid and Harsh Environments

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.

The olfactory limbus of the red fox (Vulpes vulpes). New insights regarding a noncanonical olfactory bulb pathway

Introduction: The olfactory system in most mammals is divided into several subsystems based on the anatomical locations of the neuroreceptor cells involved and the receptor families that are expressed. In addition to the main olfactory system and the vomeronasal system, a range of olfactory subsystems converge onto the transition zone located between the main olfactory bulb (MOB) and the accessory olfactory bulb (AOB), which has been termed the olfactory limbus (OL). The OL contains specialized glomeruli that receive noncanonical sensory afferences and which interact with the MOB and AOB. Little is known regarding the olfactory subsystems of mammals other than laboratory rodents. Methods: We have focused on characterizing the OL in the red fox by performing general and specific histological stainings on serial sections, using both single and double immunohistochemical and lectin-histochemical labeling techniques. Results: As a result, we have been able to determine that the OL of the red fox (Vulpes vulpes) displays an uncommonly high degree of development and complexity. Discussion: This makes this species a novel mammalian model, the study of which could improve our understanding of the noncanonical pathways involved in the processing of chemosensory cues.

Mitogenomic analysis of Rüppell's fox (Vulpes rueppellii) confirms phylogenetic placement within the Palaearctic clade shared with its sister species, the red fox (Vulpes vulpes)

The Rüppell's fox (Vulpes rueppellii) inhabits desert regions across North Africa, the Arabian Peninsula and southwestern Asia. Its phylogenetic relationship with other fox species, especially within the phylogeographic context of its sister species, V. vulpes, remain unclear. We here report the sequencing and de-novo assembly of the first annotated mitogenome of V. rueppellii, analysed with data from other foxes (tribe Vulpini, subfamily Caninae). We used four bioinformatic approaches to reconstruct the V. rueppellii mitogenome, obtaining identical sequences except for the incompletely assembled tandem-repeat region within the D-loop. The mitogenome displayed an identical organization, number and length of genes as V. vulpes. We found high support for clustering of both known subclades of V. rueppellii within the Palearctic clade of V. vulpes, rendering the latter species paraphyletic, consistent with previous analyses of shorter mtDNA fragments. More work is needed for a full understanding of the evolutionary drivers and consequences of hybridization in foxes.

Selection Signatures in Italian Livestock Guardian and Herding Shepherd Dogs

Livestock guardian (LGD) and herding shepherd (HSD) dogs have distinct morphological and behavioural characteristics, long selected by farmers and breeders, to accomplish different tasks. This study aimed to find the genomic regions that best differentiate and characterise Italian LGD and HSD. Genomic data of 158 dogs of four LGD and five HSD breeds, obtained with the 170K canine SNPchip, were collected. The two groups were compared using FST and XP-EHH analyses, identifying regions containing 29 genes. Moreover, 16 islands of runs of homozygosity were found in LGD, and 15 in HSD; 4 of them were partially shared. Among the genes found that better differentiated HSD and LGD, several were associated with dog domestication and behavioural aspects; particularly, MSRB3 and LLPH were linked to herding behaviour in previous studies. Others, DYSK, MAP2K5, and RYR, were related to body size and muscle development. Prick ears prevailed in sampled HSD, and drop ears in LGD; this explains the identification of WIF1 and MSRB3 genes. Unexpectedly, a number of genes were also associated with eye development and functionality. These results shed further light on the differences that human selection introduced in dogs aimed at different duties, even in a limited geographic area such as Italy.

Genomic Consequences of Fragmentation in the Endangered Fennoscandian Arctic Fox (Vulpes lagopus)

Accelerating climate change is causing severe habitat fragmentation in the Arctic, threatening the persistence of many cold-adapted species. The Scandinavian arctic fox (Vulpes lagopus) is highly fragmented, with a once continuous, circumpolar distribution, it struggled to recover from a demographic bottleneck in the late 19th century. The future persistence of the entire Scandinavian population is highly dependent on the northernmost Fennoscandian subpopulations (Scandinavia and the Kola Peninsula), to provide a link to the viable Siberian population. By analyzing 43 arctic fox genomes, we quantified genomic variation and inbreeding in these populations. Signatures of genome erosion increased from Siberia to northern Sweden indicating a stepping-stone model of connectivity. In northern Fennoscandia, runs of homozygosity (ROH) were on average ~1.47-fold longer than ROH found in Siberia, stretching almost entire scaffolds. Moreover, consistent with recent inbreeding, northern Fennoscandia harbored more homozygous deleterious mutations, whereas Siberia had more in heterozygous state. This study underlines the value of documenting genome erosion following population fragmentation to identify areas requiring conservation priority. With the increasing fragmentation and isolation of Arctic habitats due to global warming, understanding the genomic and demographic consequences is vital for maintaining evolutionary potential and preventing local extinctions.

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/.

Comparative genomics uncovers the evolutionary history, demography, and molecular adaptations of South American canids

The remarkable radiation of South American (SA) canids produced 10 extant species distributed across diverse habitats, including disparate forms such as the short-legged, hypercarnivorous bush dog and the long-legged, largely frugivorous maned wolf. Despite considerable research spanning nearly two centuries, many aspects of their evolutionary history remain unknown. Here, we analyzed 31 whole genomes encompassing all extant SA canid species to assess phylogenetic relationships, interspecific hybridization, historical demography, current genetic diversity, and the molecular bases of adaptations in the bush dog and maned wolf. We found that SA canids originated from a single ancestor that colonized South America 3.9 to 3.5 Mya, followed by diversification east of the Andes and then a single colonization event and radiation of Lycalopex species west of the Andes. We detected extensive historical gene flow between recently diverged lineages and observed distinct patterns of genomic diversity and demographic history in SA canids, likely induced by past climatic cycles compounded by human-induced population declines. Genome-wide scans of selection showed that disparate limb proportions in the bush dog and maned wolf may derive from mutations in genes regulating chondrocyte proliferation and enlargement. Further, frugivory in the maned wolf may have been enabled by variants in genes associated with energy intake from short-chain fatty acids. In contrast, unique genetic variants detected in the bush dog may underlie interdigital webbing and dental adaptations for hypercarnivory. Our analyses shed light on the evolution of a unique carnivoran radiation and how it was shaped by South American topography and climate change.

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

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.

How much biology is in the product? Role and relevance of biological evolution and function for bio-inspired design

Bio-inspired design (BID) means the concept of transferring functional principles from biology to technology. The core idea driving BID-related work is that evolution has shaped functional attributes, which are termed “adaptations” in biology, to a high functional performance by relentless selective pressure. For current methods and tools, such as data bases, it is implicitly supposed that the considered biological models are adaptations and their functions already clarified. Often, however, the identification of adaptations and their functional features is a difficult task which is not yet accomplished for numerous biological structures, as happens to be the case also for various organismic features from which successful BID developments were derived. This appears to question the relevance of the much stressed importance of evolution for BID. While it is obviously possible to derive an attractive technical principle from an observed biological effect without knowing its original functionality, this kind of BID (“analog BID”) has no further ties to biology. In contrast, a BID based on an adaptation and its function (“homolog BID”) is deeply embedded in biology. It is suggested that a serious and honest clarification of the functional background of a biological structure is an essential first step in devising a BID project, to recognize possible problems and pitfalls as well as to evaluate the need for further biological analysis.

Conservation of chromatin conformation in carnivores

We found the three-dimensional (3D) structure of chromatin at the chromosome level to be highly conserved for more than 50 million y of carnivore evolution. Intrachromosomal contacts were maintained even after chromosome rearrangements within carnivore lineages, demonstrating that the maintenance of 3D chromatin architecture is essential for conserved genome functions. These discoveries enabled the identification of orthologous chromosomal DNA segments among related species, a method we call 3D comparative scaffotyping. The method has application for putative chromosome assignment of chromosome-scale DNA sequence scaffolds produced by de novo genome sequencing. Broadly applied to biodiversity genome sequencing efforts, the approach can reduce costs associated with karyotyping and the physical mapping of DNA segments to chromosomes.

High throughput chromatin conformation capture (Hi-C) of leukocyte DNA was used to investigate the evolutionary stability of chromatin conformation at the chromosomal level in 11 species from three carnivore families: Felidae, Canidae, and Ursidae. Chromosome-scale scaffolds (C-scaffolds) of each species were initially used for whole-genome alignment to a reference genome within each family. This approach established putative orthologous relationships between C-scaffolds among the different species. Hi-C contact maps for all C-scaffolds were then visually compared and found to be distinct for a given reference chromosome or C-scaffold within a species and indistinguishable for orthologous C-scaffolds having a 1:1 relationship within a family. The visual patterns within families were strongly supported by eigenvectors from the Hi-C contact maps. Analysis of Hi-C contact maps and eigenvectors across the three carnivore families revealed that most cross-family orthologous subchromosomal fragments have a conserved three-dimensional (3D) chromatin structure and thus have been under strong evolutionary constraint for ∼54 My of carnivore evolution. The most pronounced differences in chromatin conformation were observed for the X chromosome and the red fox genome, whose chromosomes have undergone extensive rearrangements relative to other canids. We also demonstrate that Hi-C contact map pattern analysis can be used to accurately identify orthologous relationships between C-scaffolds and chromosomes, a method we termed “3D comparative scaffotyping.” This method provides a powerful means for estimating karyotypes in de novo sequenced species that have unknown karyotype and no physical mapping information.

Genomic characterization of the world's longest selection experiment in mouse reveals the complexity of polygenic traits

BACKGROUND

Long-term selection experiments are a powerful tool to understand the genetic background of complex traits. The longest of such experiments has been conducted in the Research Institute for Farm Animal Biology (FBN), generating extreme mouse lines with increased fertility, body mass, protein mass and endurance. For >140 generations, these lines have been maintained alongside an unselected control line, representing a valuable resource for understanding the genetic basis of polygenic traits. However, their history and genomes have not been reported in a comprehensive manner yet. Therefore, the aim of this study is to provide a summary of the breeding history and phenotypic traits of these lines along with their genomic characteristics. We further attempt to decipher the effects of the observed line-specific patterns of genetic variation on each of the selected traits.

RESULTS

Over the course of >140 generations, selection on the control line has given rise to two extremely fertile lines (>20 pups per litter each), two giant growth lines (one lean, one obese) and one long-distance running line. Whole genome sequencing analysis on 25 animals per line revealed line-specific patterns of genetic variation among lines, as well as high levels of homozygosity within lines. This high degree of distinctiveness results from the combined effects of long-term continuous selection, genetic drift, population bottleneck and isolation. Detection of line-specific patterns of genetic differentiation and structural variation revealed multiple candidate genes behind the improvement of the selected traits.

CONCLUSIONS

The genomes of the Dummerstorf trait-selected mouse lines display distinct patterns of genomic variation harbouring multiple trait-relevant genes. Low levels of within-line genetic diversity indicate that many of the beneficial alleles have arrived to fixation alongside with neutral alleles. This study represents the first step in deciphering the influence of selection and neutral evolutionary forces on the genomes of these extreme mouse lines and depicts the genetic complexity underlying polygenic traits.

Age Effects Aggressive Behavior: RNA-Seq Analysis in Cattle with Implications for Studying Neoteny Under Domestication

The reactive type of aggression is regulated mostly by the brain's prefrontal cortex; however, the molecular changes underlying aggressiveness in adults have not been fully characterized. We used an RNA-seq approach to investigate differential gene expression in the prefrontal cortex of bovines from the aggressive Lidia breed at different ages: young three-year old and adult four-year-old bulls. A total of 50 up and 193 down-regulated genes in the adult group were identified. Furthermore, a cross-species comparative analysis retrieved 29 genes in common with previous studies on aggressive behaviors, representing an above-chance overlap with the differentially expressed genes in adult bulls. We detected changes in the regulation of networks such as synaptogenesis, involved in maintenance and refinement of synapses, and the glutamate receptor pathway, which acts as excitatory driver in aggressive responses. The reduced reactive aggression typical of domestication has been proposed to form part of a retention of juvenile traits as adults (neoteny).

Can domestication shape Canidae brain morphology? The accessory olfactory bulb of the red fox as a case in point

BACKGROUND

The accessory olfactory bulb (AOB) is the first integrative center of the vomeronasal system (VNS), and the general macroscopic, microscopic, and neurochemical organizational patterns of the AOB differ fundamentally among species. Therefore, the low degree of differentiation observed for the dog AOB is surprising. As the artificial selection pressure exerted on domestic dogs has been suggested to play a key role in the involution of the dog VNS, a wild canid, such as the fox, represents a useful model for studying the hypothetical effects of domestication on the AOB morphology.

METHODS

A comprehensive histological, lectin-histochemical, and immunohistochemical study of the fox AOB was performed. Anti-Gαo and anti-Gαi2 antibodies were particularly useful, as they label the transduction cascade of the vomeronasal receptor types 1 (V1R) and 2 (V2R), respectively. Other employed antibodies included those against proteins such as microtubule-associated protein 2 (MAP-2), tubulin, glial fibrillary acidic protein, growth-associated protein 43 (GAP-43), olfactory marker protein (OMP), calbindin, and calretinin.

RESULTS

The cytoarchitecture of the fox AOB showed a clear lamination, with neatly differentiated layers; a highly developed glomerular layer, rich in periglomerular cells; and large inner cell and granular layers. The immunolabeling of Gαi2, OMP, and GAP-43 delineated the outer layers, whereas Gαo and MAP-2 immunolabeling defined the inner layers. MAP-2 characterized the somas of AOB principal cells and their dendritic trees. Anti-calbindin and anti-calretinin antibodies discriminated neural subpopulations in both the mitral-plexiform layer and the granular cell layer, and the lectin Ulex europeus agglutinin I (UEA-I) showed selectivity for the AOB and the vomeronasal nerves.

CONCLUSION

The fox AOB presents unique characteristics and a higher degree of morphological development compared with the dog AOB. The comparatively complex neural basis for semiochemical information processing in the fox compared with that observed in dogs suggests loss of AOB anatomical complexity during the evolutionary history of dogs and opens a new avenue of research for studying the effects of domestication on brain structures.

Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes), in Europe

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.

Genetic structure of immunologically associated candidate genes suggests arctic rabies variants exert differential selection in arctic fox populations

Patterns of local adaptation can emerge in response to the selective pressures diseases exert on host populations as reflected in increased frequencies of respective, advantageous genotypes. Elucidating patterns of local adaptation enhance our understanding of mechanisms of disease spread and the capacity for species to adapt in context of rapidly changing environments such as the Arctic. Arctic rabies is a lethal disease that largely persists in northern climates and overlaps with the distribution of its natural host, arctic fox. Arctic fox populations display little neutral genetic structure across their North American range, whereas phylogenetically unique arctic rabies variants are restricted in their geographic distributions. It remains unknown if arctic rabies variants impose differential selection upon host populations, nor what role different rabies variants play in the maintenance and spread of this disease. Using a targeted, genotyping-by-sequencing assay, we assessed correlations of arctic fox immunogenetic variation with arctic rabies variants to gain further insight into the epidemiology of this disease. Corroborating past research, we found no neutral genetic structure between sampled regions, but did find moderate immunogenetic structuring between foxes predominated by different arctic rabies variants. F_ST outliers associated with host immunogenetic structure included SNPs within interleukin and Toll-like receptor coding regions (IL12B, IL5, TLR3 and NFKB1); genes known to mediate host responses to rabies. While these data do not necessarily reflect causation, nor a direct link to arctic rabies, the contrasting genetic structure of immunologically associated candidate genes with neutral loci is suggestive of differential selection and patterns of local adaptation in this system. These data are somewhat unexpected given the long-lived nature and dispersal capacities of arctic fox; traits expected to undermine local adaptation. Overall, these data contribute to our understanding of the co-evolutionary relationships between arctic rabies and their primary host and provide data relevant to the management of this disease.

Genomic Signatures of Domestication Selection in the Australasian Snapper (Chrysophrys auratus)

Domestication of teleost fish is a recent development, and in most cases started less than 50 years ago. Shedding light on the genomic changes in key economic traits during the domestication process can provide crucial insights into the evolutionary processes involved and help inform selective breeding programmes. Here we report on the recent domestication of a native marine teleost species in New Zealand, the Australasian snapper (Chrysophrys auratus). Specifically, we use genome-wide data from a three-generation pedigree of this species to uncover genetic signatures of domestication selection for growth. Genotyping-By-Sequencing (GBS) was used to generate genome-wide SNP data from a three-generation pedigree to calculate generation-wide averages of FST between every generation pair. The level of differentiation between generations was further investigated using ADMIXTURE analysis and Principal Component Analysis (PCA). After that, genome scans using Bayescan, LFMM and XP-EHH were applied to identify SNP variants under putative selection following selection for growth. Finally, genes near candidate SNP variants were annotated to gain functional insights. Analysis showed that between generations FST values slightly increased as generational time increased. The extent of these changes was small, and both ADMIXTURE analysis and PCA were unable to form clear clusters. Genome scans revealed a number of SNP outliers, indicative of selection, of which a small number overlapped across analyses methods and populations. Genes of interest within proximity of putative selective SNPs were related to biological functions, and revealed an association with growth, immunity, neural development and behaviour, and tumour repression. Even though few genes overlapped between outlier SNP methods, gene functionalities showed greater overlap between methods. While the genetic changes observed were small in most cases, a number of outlier SNPs could be identified, of which some were found by more than one method. Multiple outlier SNPs appeared to be predominately linked to gene functionalities that modulate growth and survival. Ultimately, the results help to shed light on the genomic changes occurring during the early stages of domestication selection in teleost fish species such as snapper, and will provide useful candidates for the ongoing selective breeding in the future of this and related species.

Comparative transcriptome provides insights into the selection adaptation between wild and farmed foxes

The silver fox and blue fox are economically important fur species and were domesticated by humans from their wild counterparts, the arctic fox and red fox, respectively. Farmed foxes show obvious differences from their wild counterparts, including differences in physiology, body size, energy metabolism, and immunity. However, the molecular mechanisms underlying these differences are presently unclear. In this study, we built transcriptome libraries from multiple pooled tissues for each species of farmed fox, used RNA-seq to obtain a comprehensive dataset, and performed selection analysis and sequence-level analyses of orthologous genes to identify the genes that may be influenced by human domestication. More than 153.3, 248.0, 81.6, and 65.8 million clean reads were obtained and assembled into a total of 118,577, 401,520, 79,900, and 186,988 unigenes with an average length range from 521 to 667 bp for AF, BF, RF, and SF, respectively. Selective pressure analysis showed that 11 and 14 positively selected genes were identified, respectively, in the two groups (AF vs. BF and RF vs. SF). Several of these genes were associated with natural immunity (CFI and LRRFIP1), protein synthesis (GOLGA4, CEP19 and SLC35A2), and DNA damage repair (MDC1). Further functional enrichment analyses demonstrated that two positively selected genes (ACO1 and ACAD10) were involved in metabolic process (GO:0008152, p-value = .032), representing a significant enrichment. Sequence analysis of 117 orthologous genes shared by the two groups showed that the LEMD2, RRBP1, and IGBP1 genes might be affected by artificial selection in farmed foxes, with mutation sites located within sequences that are otherwise highly conserved across most mammals. Our results provide a valuable transcriptomic resource for future genetic studies and improvement in the assisted breeding of foxes and other farmed animals.

Neuromorphological Changes following Selection for Tameness and Aggression in the Russian Farm-Fox experiment

The Russian farm-fox experiment is an unusually long-running and well-controlled study designed to replicate wolf-to-dog domestication. As such, it offers an unprecedented window onto the neural mechanisms governing the evolution of behavior. Here we report evolved changes to gray matter morphology resulting from selection for tameness versus aggressive responses toward humans in a sample of 30 male fox brains. Contrasting with standing ideas on the effects of domestication on brain size, tame foxes did not show reduced brain volume. Rather, gray matter volume in both the tame and aggressive strains was increased relative to conventional farm foxes bred without deliberate selection on behavior. Furthermore, tame- and aggressive-enlarged regions overlapped substantially, including portions of motor, somatosensory, and prefrontal cortex, amygdala, hippocampus, and cerebellum. We also observed differential morphologic covariation across distributed gray matter networks. In one prefrontal-cerebellum network, this covariation differentiated the three populations along the tame-aggressive behavioral axis. Surprisingly, a prefrontal-hypothalamic network differentiated the tame and aggressive foxes together from the conventional strain. These findings indicate that selection for opposite behaviors can influence brain morphology in a similar way.SIGNIFICANCE STATEMENT Domestication represents one of the largest and most rapid evolutionary shifts of life on earth. However, its neural correlates are largely unknown. Here we report the neuroanatomical consequences of selective breeding for tameness or aggression in the seminal Russian farm-fox experiment. Compared with a population of conventional farm-bred control foxes, tame foxes show neuroanatomical changes in the PFC and hypothalamus, paralleling wolf-to-dog shifts. Surprisingly, though, aggressive foxes also show similar changes. Moreover, both strains show increased gray matter volume relative to controls. These results indicate that similar brain adaptations can result from selection for opposite behavior, that existing ideas of brain changes in domestication may need revision, and that significant neuroanatomical change can evolve very quickly, within the span of <100 generations.

Epigenetic mechanisms in sex determination and in the evolutionary transitions between sexual systems

The hypothesis that epigenetic mechanisms of gene expression regulation have two main roles in vertebrate sex is presented. First, and within a given generation, by contributing to the acquisition and maintenance of (i) the male or female function once during the lifetime in individuals of gonochoristic species; and (ii) the male and female function in the same individual, either at the same time in simultaneous hermaphrodites, or first as one sex and then as the other in sequential hermaphrodites. Second, if environmental conditions change, epigenetic mechanisms may have also a role across generations, by providing the necessary phenotypic plasticity to facilitate the transition: (i) from one sexual system to another, or (ii) from one sex-determining mechanism to another. Furthermore, if the environmental change lasts enough time, epimutations could facilitate assimilation into genetic changes that stabilize the new sexual system or sex-determining mechanism. Examples supporting these assertions are presented, caveats or difficulties and knowledge gaps identified, and possible ways to test this hypothesis suggested.

This article is part of the theme issue ‘Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)’.

Genomic and fitness consequences of inbreeding in an endangered carnivore

Reduced fitness through genetic drift and inbreeding is a major threat to small and isolated populations. Although previous studies have generally used genetically verified pedigrees to document effects of inbreeding and gene flow, these often fail to capture the whole inbreeding history of the species. By assembling a draft arctic fox (Vulpes lagopus) genome and resequencing complete genomes of 23 additional foxes born before and after a well-documented immigration event in Scandinavia, we here look into the genomic consequences of inbreeding and genetic rescue. We found a difference in genome-wide diversity, with 18% higher heterozygosity and 81% lower F_ROH in immigrant F1 compared to native individuals. However, more distant descendants of immigrants (F2, F3) did not show the same pattern. We also found that foxes with lower inbreeding had higher probability to survive their first year of life. Our results demonstrate the important link between genetic variation and fitness as well as the transient nature of genetic rescue. Moreover, our results have implications in conservation biology as they demonstrate that inbreeding depression can effectively be detected in the wild by a genomic approach.

Population Genetics of the Invasive Red Fox, Vulpes vulpes, in South-Eastern Australia

The use of genetic information in conservation biology has become more widespread with genetic information more readily available for non-model organisms. It has also been recognized that genetic information from invasive species can inform their management and control. The red fox poses a significant threat to Australian native fauna and the agricultural industry. Despite this, there are few recently published studies investigating the population genetics of foxes in Australia. This study investigated the population genetics of 94 foxes across the Illawarra and Shoalhaven regions of New South Wales, Australia. Diversity Array sequencing technology was used to genotype a large number of single nucleotide polymorphisms (N = 33,375). Moderate genetic diversity and relatedness were observed across the foxes sampled. Low to moderate levels of inbreeding, high-levels of identity-by-state values, as well as high identity-by-descent values were also found. There was limited evidence for population genetic structure among the foxes across the landscape sampled, supporting the presence of a single population across the study area. This indicates that there may be no barriers hindering fox dispersal across the landscape.

Chromosome-level genome assembly of the Arctic fox (Vulpes lagopus) using PacBio sequencing and Hi-C technology

The Arctic fox (Vulpes lagopus) is the only fox species occurring in the Arctic and has adapted to its extreme climatic conditions. Currently, the molecular basis of its adaptation to the extreme climate has not been characterized. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first V. lagopus genome assembly, which is assembled into chromosome fragments. The genome assembly has a total length of 2.345 Gb with a contig N50 of 31.848 Mb and a scaffold N50 of 131.537 Mb, consisting of 25 pseudochromosomal scaffolds. The V. lagopus genome had approximately 32.33% repeat sequences. In total, 21,278 protein-coding genes were predicted, of which 99.14% were functionally annotated. Compared with 12 other mammals, V. lagopus was most closely related to V. Vulpes with an estimated divergence time of ~7.1 Ma. The expanded gene families and positively selected genes potentially play roles in the adaptation of V. lagopus to Arctic extreme environment. This high-quality assembled genome will not only promote future studies of genetic diversity and evolution in foxes and other canids but also provide important resources for conservation of Arctic species.

Gene expression profiles underlying aggressive behavior in the prefrontal cortex of cattle

BACKGROUND

Aggressive behavior is an ancient and conserved trait, habitual for most animals in order to eat, protect themselves, compete for mating and defend their territories. Genetic factors have been shown to play an important role in the development of aggression both in animals and humans, displaying moderate to high heritability estimates. Although such types of behaviors have been studied in different animal models, the molecular architecture of aggressiveness remains poorly understood. This study compared gene expression profiles of 16 prefrontal cortex (PFC) samples from aggressive and non-aggressive cattle breeds: Lidia, selected for agonistic responses, and Wagyu, selected for tameness.

RESULTS

A total of 918 up-regulated and 278 down-regulated differentially expressed genes (DEG) were identified, representing above-chance overlap with genes previously identified in studies of aggression across species, as well as those implicated in recent human evolution. The functional interpretation of the up-regulated genes in the aggressive cohort revealed enrichment of pathways such as Alzheimer disease-presenilin, integrins and the ERK/MAPK signaling cascade, all implicated in the development of abnormal aggressive behaviors and neurophysiological disorders. Moreover, gonadotropins, are up-regulated as natural mechanisms enhancing aggression. Concomitantly, heterotrimeric G-protein pathways, associated with low reactivity mental states, and the GAD2 gene, a repressor of agonistic reactions associated with PFC activity, are down-regulated, promoting the development of the aggressive responses selected for in Lidia cattle. We also identified six upstream regulators, whose functional activity fits with the etiology of abnormal behavioral responses associated with aggression.

CONCLUSIONS

These transcriptional correlates of aggression, resulting, at least in part, from controlled artificial selection, can provide valuable insights into the complex architecture that underlies naturally developed agonistic behaviors. This analysis constitutes a first important step towards the identification of the genes and metabolic pathways that promote aggression in cattle and, providing a novel model species to disentangle the mechanisms underlying variability in aggressive behavior.

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).

Molecules, Mechanisms, and Disorders of Self-Domestication: Keys for Understanding Emotional and Social Communication from an Evolutionary Perspective

The neural crest hypothesis states that the phenotypic features of the domestication syndrome are due to a reduced number or disruption of neural crest cells (NCCs) migration, as these cells differentiate at their final destinations and proliferate into different tissues whose activity is reduced by domestication. Comparing the phenotypic characteristics of modern and prehistoric man, it is clear that during their recent evolutionary past, humans also went through a process of self-domestication with a simultaneous prolongation of the period of socialization. This has led to the development of social abilities and skills, especially language, as well as neoteny. Disorders of neural crest cell development and migration lead to many different conditions such as Waardenburg syndrome, Hirschsprung disease, fetal alcohol syndrome, DiGeorge and Treacher-Collins syndrome, for which the mechanisms are already relatively well-known. However, for others, such as Williams-Beuren syndrome and schizophrenia that have the characteristics of hyperdomestication, and autism spectrum disorders, and 7dupASD syndrome that have the characteristics of hypodomestication, much less is known. Thus, deciphering the biological determinants of disordered self-domestication has great potential for elucidating the normal and disturbed ontogenesis of humans, as well as for the understanding of evolution of mammals in general.