Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes)

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.

Molecular signature of excessive female aggression: study of stressed mice with genetic inactivation of neuronal serotonin synthesis

Aggression is a complex social behavior, critically involving brain serotonin (5-HT) function. The neurobiology of female aggression remains elusive, while the incidence of its manifestations has been increasing. Yet, animal models of female aggression are scarce. We previously proposed a paradigm of female aggression in the context of gene x environment interaction where mice with partial genetic inactivation of tryptophan hydroxylase-2 (Tph2+/- mice), a key enzyme of neuronal 5-HT synthesis, are subjected to predation stress resulting in pathological aggression. Using deep sequencing and the EBSeq method, we studied the transcriptomic signature of excessive aggression in the prefrontal cortex of female Tph2+/- mice subjected to rat exposure stress and food deprivation. Challenged mutants, but not other groups, displayed marked aggressive behaviors. We found 26 genes with altered expression in the opposite direction between stressed groups of both Tph2 genotypes. We identified several molecular markers, including Dgkh, Arfgef3, Kcnh7, Grin2a, Tenm1 and Epha6, implicated in neurodevelopmental deficits and psychiatric conditions featuring impaired cognition and emotional dysregulation. Moreover, while 17 regulons, including several relevant to neural plasticity and function, were significantly altered in stressed mutants, no alteration in regulons was detected in stressed wildtype mice. An interplay of the uncovered pathways likely mediates partial Tph2 inactivation in interaction with severe stress experience, thus resulting in excessive female 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.

Evolutionary endocrinology and the problem of Darwin's tangled bank

Like Darwin's tangled bank of biodiversity, the endocrine mechanisms that give rise to phenotypic diversity also exhibit nearly endless forms. This tangled bank of mechanistic diversity can prove problematic as we seek general principles on the role of endocrine mechanisms in phenotypic evolution. A key unresolved question is therefore: to what degree are specific endocrine mechanisms re-used to bring about replicated phenotypic evolution? Related areas of inquiry are booming in molecular ecology, but behavioral traits are underrepresented in this literature. Here, I leverage the rich comparative tradition in evolutionary endocrinology to evaluate whether and how certain mechanisms may be repeated hotspots of behavioral evolutionary change. At one extreme, mechanisms may be parallel, such that evolution repeatedly uses the same gene or pathway to arrive at multiple independent (or, convergent) origins of a particular behavioral trait. At the other extreme, the building blocks of behavior may be unique, such that outwardly similar phenotypes are generated via lineage-specific mechanisms. This review synthesizes existing case studies, phylogenetic analyses, and experimental evolutionary research on mechanistic parallelism in animal behavior. These examples show that the endocrine building blocks of behavior have some elements of parallelism across replicated evolutionary events. However, support for parallelism is variable among studies, at least some of which relates to the level of complexity at which we consider sameness (i.e. pathway vs. gene level). Moving forward, we need continued experimentation and better testing of neutral models to understand whether, how - and critically, why - mechanism A is used in one lineage and mechanism B is used in another. We also need continued growth of large-scale comparative analyses, especially those that can evaluate which endocrine parameters are more or less likely to undergo parallel evolution alongside specific behavioral traits. These efforts will ultimately deepen understanding of how and why hormone-mediated behaviors are constructed the way that they are.

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

Constitutive gene expression differs in three brain regions important for cognition in neophobic and non-neophobic house sparrows (Passer domesticus)

Neophobia (aversion to new objects, food, and environments) is a personality trait that affects the ability of wildlife to adapt to new challenges and opportunities. Despite the ubiquity and importance of this trait, the molecular mechanisms underlying repeatable individual differences in neophobia in wild animals are poorly understood. We evaluated wild-caught house sparrows (Passer domesticus) for neophobia in the lab using novel object tests. We then selected a subset of neophobic and non-neophobic individuals (n = 3 of each, all females) and extracted RNA from four brain regions involved in learning, memory, threat perception, and executive function: striatum, caudal dorsomedial hippocampus, medial ventral arcopallium, and caudolateral nidopallium (NCL). Our analysis of differentially expressed genes (DEGs) used 11,889 gene regions annotated in the house sparrow reference genome for which we had an average of 25.7 million mapped reads/sample. PERMANOVA identified significant effects of brain region, phenotype (neophobic vs. non-neophobic), and a brain region by phenotype interaction. Comparing neophobic and non-neophobic birds revealed constitutive differences in DEGs in three of the four brain regions examined: hippocampus (12% of the transcriptome significantly differentially expressed), striatum (4%) and NCL (3%). DEGs included important known neuroendocrine mediators of learning, memory, executive function, and anxiety behavior, including serotonin receptor 5A, dopamine receptors 1, 2 and 5 (downregulated in neophobic birds), and estrogen receptor beta (upregulated in neophobic birds). These results suggest that some of the behavioral differences between phenotypes may be due to underlying gene expression differences in the brain. The large number of DEGs in neophobic and non-neophobic birds also implies that there are major differences in neural function between the two phenotypes that could affect a wide variety of behavioral traits beyond neophobia.

Transcriptome analysis provides new insights into cold adaptation of corsac fox ( Vulpes Corsac )

Vulpesare widely distributed throughout the world and have undergone drastic physiological and phenotypic changes in response to their environment. However, little is known about the underlying genetic causes of these traits, especially Vulpes corsac. In this study, RNA‐Seq was used to obtain a comprehensive dataset for multiple pooled tissues of corsac fox, and selection analysis of orthologous genes was performed to identify the genes that may be influenced by the low‐temperature environment. More than 6.32 Gb clean reads were obtained and assembled into a total of 173,353 unigenes with an average length of 557 bp for corsac fox. Selective pressure analysis showed that 16 positively selected genes (PSGs) were identified in corsac fox, red fox, and arctic fox. Enrichment analysis of PSGs showed that the LRP11 gene was enriched in several pathways related to the low‐temperature response and might play a key role in response to environmental stimuli of foxes. In addition, several positively selected genes were related to DNA damage repair (ELP2 and CHAF1A), innate immunity (ARRDC4 and S100A12), and the respiratory chain (NDUFA5), and these positively selected genes might play a role in adaptation to harsh wild fox environments. The results of common orthologous gene analysis showed that gene flow or convergent evolution might be an important factor in promoting regional differentiation of foxes. Our study provides a valuable transcriptomic resource for the evolutionary history of the corsac fox and the adaptations to the extreme environments.

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

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.

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.

RNA-Seq analysis of the guppy immune response against Gyrodactylus bullatarudis infection

Gyrodactylids are ubiquitous ectoparasites of teleost fish, but our understanding of the host immune response against them is fragmentary. Here, we used RNA-Seq to investigate genes involved in the primary response to infection with Gyrodactylus bullatarudis on the skin of guppies, Poecilia reticulata, an important evolutionary model, but also one of the most common fish in the global ornamental trade. Analysis of differentially expressed genes identified several immune-related categories, including IL-17 signalling pathway and Th17 cell differentiation, cytokine-cytokine receptor interaction, chemokine signalling pathway, NOD-like receptor signalling pathway, natural killer cell-mediated cytotoxicity and pathways involved in antigen recognition, processing and presentation. Components of both the innate and the adaptive immune responses play a role in response to gyrodactylid infection. Genes involved in IL-17/Th17 response were particularly enriched among differentially expressed genes, suggesting a significant role for this pathway in fish responses to ectoparasites. Our results revealed a sizable list of genes potentially involved in the teleost-gyrodactylid immune response.

How research on female vertebrates contributes to an expanded challenge hypothesis

The bi-directional links between hormones and behavior have been a rich area of research for decades. Theory on the evolution of testosterone (T) was greatly advanced by the challenge hypothesis, which presented a framework for understanding interspecific, seasonal, and social variation in T levels in males, and how they are shaped by the competing demands of parental care and male-male competition. Female competition is also widespread in nature, although it is less clear whether or how the challenge hypothesis applies to females. Here, we evaluate this issue in four parts: (1) We summarize and update prior analyses of seasonal plasticity and interspecific variation in T in females. (2) We evaluate experimental links between T and female aggression on shorter timescales, asking how T manipulations affect aggression and conversely, how social manipulations affect T levels in female mammals, birds, lizards, and fishes. (3) We examine alternative mechanisms that may link aggression to the social environment independently of T levels in circulation. (4) We present a case study, including new data analyses, in an aggressive female bird (the tree swallow, Tachycineta bicolor) to explore how variation in tissue-level processing of T may bridge the gap between circulating T and variation in behavior that is visible to natural selection. We close by connecting these multivariate levels of sex steroid signaling systems alongside different temporal scales (social, seasonal, and evolutionary) to generate broadly applicable insights into how animals respond to their social environment, regardless of whether they are male or female.

Ecotype differences in aggression, neural activity and behaviorally relevant gene expression in cichlid fish

In Lake Malawi, two ecologically distinct lineages of cichlid fishes (rock- vs sand-dwelling ecotypes, each comprised of over 200 species) evolved within the last million years. The rock-dwelling species (Mbuna) are aggressively territorial year-round and males court and spawn with females over rocky substrate. In contrast, males of sand-dwelling species are not territorial and instead aggregate on seasonal breeding leks in which males construct courtship "bowers" in the sand. However, little is known about how phenotypic variation in aggression is produced by the genome. In this study, we first quantify and compare behavior in seven cichlid species, demonstrating substantial ecotype and species differences in unconditioned mirror-elicited aggression. Second, we compare neural activity in mirror-elicited aggression in two representative species, Mchenga conophoros (sand-dwelling) and Petrotilapia chitimba (rock-dwelling). Finally, we compare gene expression patterns between these two species, specifically within neurons activated during mirror aggression. We identified a large number of genes showing differential expression in mirror-elicited aggression, as well as many genes that differ between ecotypes. These genes, which may underly species differences in behavior, include several neuropeptides, genes involved in the synthesis of steroid hormones and neurotransmitter activity. This work lays the foundation for future experiments using this emerging genetic model system to investigate the genomic basis of evolved species differences in both brain and behavior.

The Australian dingo: untamed or feral?

BACKGROUND

The Australian dingo continues to cause debate amongst Aboriginal people, pastoralists, scientists and the government in Australia. A lingering controversy is whether the dingo has been tamed and has now reverted to its ancestral wild state or whether its ancestors were domesticated and it now resides on the continent as a feral dog. The goal of this article is to place the discussion onto a theoretical framework, highlight what is currently known about dingo origins and taxonomy and then make a series of experimentally testable organismal, cellular and biochemical predictions that we propose can focus future research.

DISCUSSION

We consider a canid that has been unconsciously selected as a tamed animal and the endpoint of methodical or what we now call artificial selection as a domesticated animal. We consider wild animals that were formerly tamed as untamed and those wild animals that were formerly domesticated as feralized. Untamed canids are predicted to be marked by a signature of unconscious selection whereas feral animals are hypothesized to be marked by signatures of both unconscious and artificial selection. First, we review the movement of dingo ancestors into Australia. We then discuss how differences between taming and domestication may influence the organismal traits of skull morphometrics, brain and size, seasonal breeding, and sociability. Finally, we consider cellular and molecular level traits including hypotheses concerning the phylogenetic position of dingoes, metabolic genes that appear to be under positive selection and the potential for micronutrient compensation by the gut microbiome.

CONCLUSIONS

Western Australian Government policy is currently being revised to allow the widespread killing of the Australian dingo. These policies are based on an incomplete understanding of the evolutionary history of the canid and assume the dingo is feralized. However, accumulated evidence does not definitively show that the dingo was ever domesticated and additional focused research is required. We suggest that incorporating ancient DNA data into the debate concerning dingo origins will be pivotal to understanding the evolutionary history of the canid. Further, we advocate that future morphological, behavioural and genetic studies should focus on including genetically pure Alpine and Desert dingoes and not dingo-dog hybrids. Finally, we propose that future studies critically examine genes under selection in the dingo and employ the genome from a wild canid for comparison.

Hippocampal glucocorticoid receptor and microRNA gene expression and serum cortisol concentration in foxes selected for behavior toward humans

In many cases, stress reactivity is one of the important bases of aggressive behavior. It appears as if reduced stress reactivity underlies an abrupt decrease in aggression towards man in domesticated animals. However, the mechanisms of this reduction have yet to be resolved. In this work, we used an experimental domestication model, the silver fox selected for many years for the response to humans to study cortisol stress reactivity in tame and aggressive foxes in response to immobilization in human arms. Additionally, these behavioral fox groups were explored for one of the important mechanisms of glucocorticoid negative feedback, the expression of the glucocorticoid receptor gene (NR3C1) in a portion of the dorsal hippocampus. In recent years, attention has been paid to differences in miRNA expression patterns between animals with different behavior and stress reactivity, as well as to miRNA regulation under stress. The same applies to NR3C1 mRNA as well. That is why we performed a miRNA-seq analysis on a portion of the fox dorsal hippocampus. It has been demonstrated that immobilization in human arms leads to significantly higher stressinduced cortisol levels in aggressive than tame foxes. At the same time, no differences have been found between hippocampal NR3C1 gene expression and the pattern of miRNA expression. Thus, reduced stress reactivity in foxes during selection for the absence of aggressive responses and for the presence of emotionally positive responses to humans does not seem to be associated with important mechanisms of regulation such as alterations in hippocampal NR3C1 gene expression or microRNA-mediated silencing.

Epigenetics and early domestication: differences in hypothalamic DNA methylation between red junglefowl divergently selected for high or low fear of humans

Background

Domestication of animals leads to large phenotypic alterations within a short evolutionary time-period. Such alterations are caused by genomic variations, yet the prevalence of modified traits is higher than expected if they were caused only by classical genetics and mutations. Epigenetic mechanisms may also be important in driving domesticated phenotypes such as behavior traits. Gene expression can be modulated epigenetically by mechanisms such as DNA methylation, resulting in modifications that are not only variable and susceptible to environmental stimuli, but also sometimes transgenerationally stable. To study such mechanisms in early domestication, we used as model two selected lines of red junglefowl (ancestors of modern chickens) that were bred for either high or low fear of humans over five generations, and investigated differences in hypothalamic DNA methylation between the two populations.

Results

Twenty-two 1-kb windows were differentially methylated between the two selected lines at p < 0.05 after false discovery rate correction. The annotated functions of the genes within these windows indicated epigenetic regulation of metabolic and signaling pathways, which agrees with the changes in gene expression that were previously reported for the same tissue and animals.

Conclusions

Our results show that selection for an important domestication-related behavioral trait such as tameness can cause divergent epigenetic patterns within only five generations, and that these changes could have an important role in chicken domestication.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0384-z) contains supplementary material, which is available to authorized users.

Rapid, broad-scale gene expression evolution in experimentally harvested fish populations

Gene expression changes potentially play an important role in adaptive evolution under human-induced selection pressures, but this has been challenging to demonstrate in natural populations. Fishing exhibits strong selection pressure against large body size, thus potentially inducing evolutionary changes in life history and other traits that may be slowly reversible once fishing ceases. However, there is a lack of convincing examples regarding the speed and magnitude of fisheries-induced evolution, and thus, the relevant underlying molecular-level effects remain elusive. We use wild-origin zebrafish (Danio rerio) as a model for harvest-induced evolution. We experimentally demonstrate broad-scale gene expression changes induced by just five generations of size-selective harvesting, and limited genetic convergence following the cessation of harvesting. We also demonstrate significant allele frequency changes in genes that were differentially expressed after five generations of size-selective harvesting. We further show that nine generations of captive breeding induced substantial gene expression changes in control stocks likely due to inadvertent selection in the captive environment. The large extent and rapid pace of the gene expression changes caused by both harvest-induced selection and captive breeding emphasizes the need for evolutionary enlightened management towards sustainable fisheries.

Lessons from the canine Oxtr gene: populations, variants and functional aspects

Oxytocin receptor (OXTR) acts as a key behavioral modulator of the central nervous system, affecting social behavior, stress, affiliation and cognitive functions. Variants of the Oxtr gene are known to influence behavior both in animals and humans; however, canine Oxtr polymorphisms are less characterized in terms of possible relevance to function, selection criteria in breeding and domestication. In this report, we provide a detailed characterization of common variants of the canine Oxtr gene. In particular (1) novel polymorphisms were identified by direct sequencing of wolf and dog samples, (2) allelic distributions and pairwise linkage disequilibrium patterns of several canine populations were compared, (3) neighbor joining (NJ) tree based on common single nucleotide polymorphisms (SNPs) was constructed, (4) mRNA expression features were assessed, (5) a novel splice variant was detected and (6) in vitro functional assays were performed. Results indicate marked differences regarding Oxtr variations between purebred dogs of different breeds, free-ranging dog populations, wolf subspecies and golden jackals. This, together with existence of explicitly dog-specific alleles and data obtained from the NJ tree implies that Oxtr could indeed have been a target gene during domestication and selection for human preferred aspects of temperament and social behavior. This assumption is further supported by the present observations on gene expression patterns within the brain and luciferase reporter experiments, providing a molecular level link between certain canine Oxtr polymorphisms and differences in nervous system function and behavior.

The concerted impact of domestication and transposon insertions on methylation patterns between dogs and grey wolves

The process of domestication can exert intense trait-targeted selection on genes and regulatory regions. Specifically, rapid shifts in the structure and sequence of genomic regulatory elements could provide an explanation for the extensive, and sometimes extreme, variation in phenotypic traits observed in domesticated species. Here, we explored methylation differences from >24 000 cytosines distributed across the genomes of the domesticated dog (Canis familiaris) and the grey wolf (Canis lupus). PCA and model-based cluster analyses identified two primary groups, domestic vs. wild canids. A scan for significantly differentially methylated sites (DMSs) revealed species-specific patterns at 68 sites after correcting for cell heterogeneity, with weak yet significant hypermethylation typical of purebred dogs when compared to wolves (59% and 58%, P < 0.05, respectively). Additionally, methylation patterns at eight genes significantly deviated from neutrality, with similar trends of hypermethylation in purebred dogs. The majority (>66%) of differentially methylated regions contained or were associated with repetitive elements, indicative of a genotype-mediated trend. However, DMSs were also often linked to functionally relevant genes (e.g. neurotransmitters). Finally, we utilized known genealogical relationships among Yellowstone wolves to survey transmission stability of methylation marks, from which we found a substantial fraction that demonstrated high heritability (both H(2) and h(2 ) > 0.99). These analyses provide a unique epigenetic insight into the molecular consequences of recent selection and radiation of our most ancient domesticated companion, the dog. These findings suggest selection has acted on methylation patterns, providing a new genomic perspective on phenotypic diversification in domesticated species.

Transcriptome Analysis in Domesticated Species: Challenges and Strategies

Domesticated species occupy a special place in the human world due to their economic and cultural value. In the era of genomic research, domesticated species provide unique advantages for investigation of diseases and complex phenotypes. RNA sequencing, or RNA-seq, has recently emerged as a new approach for studying transcriptional activity of the whole genome, changing the focus from individual genes to gene networks. RNA-seq analysis in domesticated species may complement genome-wide association studies of complex traits with economic importance or direct relevance to biomedical research. However, RNA-seq studies are more challenging in domesticated species than in model organisms. These challenges are at least in part associated with the lack of quality genome assemblies for some domesticated species and the absence of genome assemblies for others. In this review, we discuss strategies for analyzing RNA-seq data, focusing particularly on questions and examples relevant to domesticated species.

Natural variation in brain gene expression profiles of aggressive and nonaggressive individual sticklebacks

Within many species, some individuals are consistently more aggressive than others. We examine whether there are differences in brain gene expression between aggressive versus nonaggressive behavioural types of individuals within a natural population of male three-spined sticklebacks (Gasterosteus aculeatus). We compared gene expression profiles of aggressive male sticklebacks to nonaggressive males in four regions of the brain (brainstem, cerebellum, diencephalon and telencephalon). Relatively few genes were differentially expressed between behavioural types in telencephalon, cerebellum and diencephalon, but hundreds of genes were differentially expressed in brainstem, a brain area involved in detecting threats. Six genes that were differentially expressed in response to a territorial intrusion in a previous study were also differentially expressed between behavioural types in this study, implying primarily non-shared but some shared molecular mechanisms. Our findings offer new insights into the molecular causes and correlates of behavioural plasticity and individual variation in behaviour.

Genetic signatures of adaptation revealed from transcriptome sequencing of Arctic and red foxes

Background

The genus Vulpes (true foxes) comprises numerous species that inhabit a wide range of habitats and climatic conditions, including one species, the Arctic fox (Vulpes lagopus) which is adapted to the arctic region. A close relative to the Arctic fox, the red fox (Vulpes vulpes), occurs in subarctic to subtropical habitats. To study the genetic basis of their adaptations to different environments, transcriptome sequences from two Arctic foxes and one red fox individual were generated and analyzed for signatures of positive selection. In addition, the data allowed for a phylogenetic analysis and divergence time estimate between the two fox species.

Results

The de novo assembly of reads resulted in more than 160,000 contigs/transcripts per individual. Approximately 17,000 homologous genes were identified using human and the non-redundant databases. Positive selection analyses revealed several genes involved in various metabolic and molecular processes such as energy metabolism, cardiac gene regulation, apoptosis and blood coagulation to be under positive selection in foxes. Branch site tests identified four genes to be under positive selection in the Arctic fox transcriptome, two of which are fat metabolism genes. In the red fox transcriptome eight genes are under positive selection, including molecular process genes, notably genes involved in ATP metabolism. Analysis of the three transcriptomes and five Sanger re-sequenced genes in additional individuals identified a lower genetic variability within Arctic foxes compared to red foxes, which is consistent with distribution range differences and demographic responses to past climatic fluctuations. A phylogenomic analysis estimated that the Arctic and red fox lineages diverged about three million years ago.

Conclusions

Transcriptome data are an economic way to generate genomic resources for evolutionary studies. Despite not representing an entire genome, this transcriptome analysis identified numerous genes that are relevant to arctic adaptation in foxes. Similar to polar bears, fat metabolism seems to play a central role in adaptation of Arctic foxes to the cold climate, as has been identified in the polar bear, another arctic specialist.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1724-9) contains supplementary material, which is available to authorized users.

Genotyping-By-Sequencing (GBS) Detects Genetic Structure and Confirms Behavioral QTL in Tame and Aggressive Foxes (Vulpes vulpes)

The silver fox (Vulpes vulpes) offers a novel model for studying the genetics of social behavior and animal domestication. Selection of foxes, separately, for tame and for aggressive behavior has yielded two strains with markedly different, genetically determined, behavioral phenotypes. Tame strain foxes are eager to establish human contact while foxes from the aggressive strain are aggressive and difficult to handle. These strains have been maintained as separate outbred lines for over 40 generations but their genetic structure has not been previously investigated. We applied a genotyping-by-sequencing (GBS) approach to provide insights into the genetic composition of these fox populations. Sequence analysis of EcoT22I genomic libraries of tame and aggressive foxes identified 48,294 high quality SNPs. Population structure analysis revealed genetic divergence between the two strains and more diversity in the aggressive strain than in the tame one. Significant differences in allele frequency between the strains were identified for 68 SNPs. Three of these SNPs were located on fox chromosome 14 within an interval of a previously identified behavioral QTL, further supporting the importance of this region for behavior. The GBS SNP data confirmed that significant genetic diversity has been preserved in both fox populations despite many years of selective breeding. Analysis of SNP allele frequencies in the two populations identified several regions of genetic divergence between the tame and aggressive foxes, some of which may represent targets of selection for behavior. The GBS protocol used in this study significantly expanded genomic resources for the fox, and can be adapted for SNP discovery and genotyping in other canid species.

Microsatellite polymorphism and its association with body weight and selected morphometrics of farm red fox (Vulpes vulpes L.)

Polymorphism of 30 canine-derived microsatellites was studied in a group of 200 red foxes kept on 2 Polish farms. 22 out of 30 microsatellites were selected to study association between marker genotypes and body weight (BW), body length (BL), body circumference (BC), tail length (TL), ear height (EH), length of the right front limb (FRLL), length of the right rear limb (RRLL), length of the right front foot (FRFL) and length of the right rear foot (RRFL). A total of 112 alleles and 243 genotypes were found at 22 autosomal microsatellite loci. Three monomorphic loci deemed as uninformative were excluded from the study. The association between marker genotypes and the studied traits was analysed using general linear model (GLM) procedure and least squares means (LSM). Linkage disequilibrium (LD) was estimated to assess non-random association between microsatellite loci. Out of 19 microsatellites studied four markers showed no association with the studied traits, three markers had a significant effect on one trait, and another three markers had significant effect on two traits. Among ten microsatellites with significant effect on four economically important traits (BW, BL, BC, TL) four were associated with two characters: marker FH2613 with BW and BC, marker FH2097withBL and BC, marker ZUBECA6 with BW and BC, whereas marker REN75M10 was associated with BL and TL. The strongest LD (r² ranged from 0.15 to 0.33) was estimated between nine loci with significant effect on economically important traits (BW, BL, BC, TL).

Proximate perspectives on the evolution of female aggression: good for the gander, good for the goose?

Female-female aggression often functions in competition over reproductive or social benefits, but the proximate mechanisms of this apparently adaptive behaviour are not well understood. The sex steroid hormone testosterone (T) and its metabolites are well-established mediators of male-male aggression, and several lines of evidence suggest that T-mediated mechanisms may apply to females as well. However, a key question is whether mechanisms of female aggression primarily reflect correlated evolutionary responses to selection acting on males, or whether direct selection acting on females has made modifications to these mechanisms that are adaptive in light of female life history. Here, I examine the degree to which female aggression is mediated at the level of T production, target tissue sensitivity to T, or downstream genomic responses in order to test the hypothesis that selection favours mechanisms that facilitate female aggression while minimizing the costs of systemically elevated T. I draw heavily from avian systems, including the dark-eyed junco (Junco hyemalis), as well as other organisms in which these mechanisms have been well studied from an evolutionary/ecological perspective in both sexes. Findings reveal that the sexes share many behavioural and hormonal mechanisms, though several patterns also suggest sex-specific adaptation. I argue that greater attention to multiple levels of analysis-from hormone to receptor to gene network, including analyses of individual variation that represents the raw material of evolutionary change-will be a fruitful path for understanding mechanisms of behavioural regulation and intersexual coevolution.

Transcriptome sequencing of Zhikong scallop (Chlamys farreri) and comparative transcriptomic analysis with Yesso scallop (Patinopecten yessoensis)

Background

Bivalves play an important role in the ecosystems they inhabit and represent an important food source all over the world. So far limited genetic research has focused on this group of animals largely due to the lack of sufficient genetic or genomic resources. Here, we performed de novo transcriptome sequencing to produce the most comprehensive expressed sequence tag resource for Zhikong scallop (Chlamys farreri), and conducted the first transcriptome comparison for scallops.

Results

In a single 454 sequencing run, 1,033,636 reads were produced and then assembled into 26,165 contigs. These contigs were then clustered into 24,437 isotigs and further grouped into 20,056 isogroups. About 47% of the isogroups showed significant matches to known proteins based on sequence similarity. Transcripts putatively involved in growth, reproduction and stress/immune-response were identified through Gene ontology (GO) and KEGG pathway analyses. Transcriptome comparison with Yesso scallop (Patinopecten yessoensis) revealed similar patterns of GO representation. Moreover, 38 putative fast-evolving genes were identified through analyzing the orthologous gene pairs between the two scallop species. More than 46,000 single nucleotide polymorphisms (SNPs) and 350 simple sequence repeats (SSRs) were also detected.

Conclusion

Our study provides the most comprehensive transcriptomic resource currently available for C. farreri. Based on this resource, we performed the first large-scale transcriptome comparison between the two scallop species, C. farreri and P. yessoensis, and identified a number of putative fast-evolving genes, which may play an important role in scallop speciation and/or local adaptation. A large set of single nucleotide polymorphisms and simple sequence repeats were identified, which are ready for downstream marker development. This transcriptomic resource should lay an important foundation for future genetic or genomic studies on C. farreri.

The secret language of destiny: stress imprinting and transgenerational origins of disease

Epigenetic regulation modulates gene expression without altering the DNA sequence to facilitate rapid adjustments to dynamically changing environmental conditions. The formation of an epigenetic memory allows passing on this information to subsequent generations. Here we propose that epigenetic memories formed by adverse environmental conditions and stress represent a critical determinant of health and disease in the F3 generation and beyond. Transgenerational programming of epigenetic regulation may represent a key to understand adult-onset complex disease pathogenesis and cumulative effects of life span and familial disease etiology. Ultimately, the mechanisms of generating an epigenetic memory may become of potentially promising diagnostic and therapeutic relevance due to their reversible nature. Exploring the role of environmental factors, such as stress, in causing variations in epigenetic profiles may lead to new avenues of personalized, preventive medicine based on epigenetic signatures and interventions.

Transcriptome sequencing and characterization for the sea cucumber Apostichopus japonicus (Selenka, 1867)

Background

Sea cucumbers are a special group of marine invertebrates. They occupy a taxonomic position that is believed to be important for understanding the origin and evolution of deuterostomes. Some of them such as Apostichopus japonicus represent commercially important aquaculture species in Asian countries. Many efforts have been devoted to increasing the number of expressed sequence tags (ESTs) for A. japonicus, but a comprehensive characterization of its transcriptome remains lacking. Here, we performed the large-scale transcriptome profiling and characterization by pyrosequencing diverse cDNA libraries from A. japonicus.

Results

In total, 1,061,078 reads were obtained by 454 sequencing of eight cDNA libraries representing different developmental stages and adult tissues in A. japonicus. These reads were assembled into 29,666 isotigs, which were further clustered into 21,071 isogroups. Nearly 40% of the isogroups showed significant matches to known proteins based on sequence similarity. Gene ontology (GO) and KEGG pathway analyses recovered diverse biological functions and processes. Candidate genes that were potentially involved in aestivation were identified. Transcriptome comparison with the sea urchin Strongylocentrotus purpuratus revealed similar patterns of GO term representation. In addition, 4,882 putative orthologous genes were identified, of which 202 were not present in the non-echinoderm organisms. More than 700 simple sequence repeats (SSRs) and 54,000 single nucleotide polymorphisms (SNPs) were detected in the A. japonicus transcriptome.

Conclusion

Pyrosequencing was proven to be efficient in rapidly identifying a large set of genes for the sea cucumber A. japonicus. Through the large-scale transcriptome sequencing as well as public EST data integration, we performed a comprehensive characterization of the A. japonicus transcriptome and identified candidate aestivation-related genes. A large number of potential genetic markers were also identified from the A. japonicus transcriptome. This transcriptome resource would lay an important foundation for future genetic or genomic studies on this species.