Skip to main content
  • Research article
  • Open access
  • Published:

Analysis of the brain transcriptome in lines of laying hens divergently selected for feather pecking

Abstract

Background

Feather pecking (FP) in laying hens reduces animal welfare and leads to economic losses for the layer industry. FP is considered a heritable condition that is influenced by dysregulation of neurotransmitter homeostasis, the gut microbiome, and the immune system. To identify genes and biological pathways responsible for FP behavior we compared the brain transcriptomes of 48 hens divergently selected for FP. In addition, we tested if high feather peckers (HFP) and low feather peckers (LFP) respond differently to light since light has been shown to trigger FP behavior.

Results

Of approximately 48 million reads/sample an average of 98.4% were mapped to the chicken genome (GRCg6a). We found 13,070 expressed genes in the analyzed brains of which 423 showed differential expression between HFP and LFP. Genes of uncertain function and non-coding RNAs were overrepresented among those transcripts. Functional analyses revealed the involvement of cholinergic signaling, postsynaptic activity, membrane channels, and the immune system. After the light stimulus, 28 genes were found to be differentially expressed. These included an interaction cluster of core components of the circadian clock. However, differences in the response to light between HFP and LFP were not detectable.

Conclusions

Genes involved in cholinergic signaling, channel activity, synaptic transmission, and immune response were found to be involved in FP behavior. We propose a model in which the gut microbiota modulates the immune system, which in turn affects cholinergic signaling. This might have an influence on monoamine signaling with possible involvement of GABA or glutamate signaling.

Background

Feather pecking (FP) is a serious problem of the layer industry causing economic losses and massive impairments of animal welfare. The propensity to perform this damaging behavior is a complex trait influenced by numerous environmental factors as well as a genetic component. Reported heritability estimates of around 0.15 [1,2,3] indicate the possibility of genetic selection against FP. Despite intensive research during the last decades, the underlying mechanisms are still not well understood.

A longstanding theory explained FP as a redirected foraging behavior [4, 5], which was not confirmed in later studies. Instead, it has been linked to general locomotor activity [6, 7] and feather eating [7, 8]. With respect to the latter, it was reported that the inclusion of feathers in the diet affects the chicken’s gut microbiome [9] as well as the actual pecking behavior [10]. Comparisons of the cecal microbiomes between layer lines divergently selected for FP or with different actual FP behavior revealed apparent differences [11,12,13]. It is difficult to determine, whether these differences are a cause for FP or just a consequence of feather consumption. The hypothesis of a link with the gut microbiome is, however, convincing as FP is influenced by the serotonergic system [14,15,16], which comprises a central and a peripheral part [17]. Central serotonin partly controls dopamine and thus also affects reward-related behaviors. The peripheral part is mainly found in enterochromaffin cells of the gastrointestinal tract. Both systems do interact via the bloodstream or the vagus nerve [17]. Gut microbiota can affect the serotonergic system in various ways such as the activation of the immune system [17, 18]. The latter is notable as FP has also been directly linked to immune response [19, 20]. Comprehensive genomic and transcriptomic studies might help to bring ends together by identifying potentially involved genes and pathways but are still rare. Several gene expression studies linked FP to genes involved in serotonergic signaling, either directly [21, 22] or indirectly via dopaminergic [22] or GABAergic [21] signaling or via a link to obsessive-compulsive disorders [23]. Mapping studies identified loci connecting FP to dopaminergic [24], GABAergic [25] and serotonergic [2, 24, 26] signaling as well as to the immune system [26].

As comprehensively reviewed by [18], it thus seems likely that peripheral serotonin, the gut microbiome, and the immune system are interacting with central serotonin and dopamine to influence FP. Despite the huge body of evidence for the involvement of the mentioned mechanisms in FP, it is still not understood in detail how they interact. With the aim to gain further insight into the underlying mechanisms, we used RNA-sequencing (RNA-seq) to analyze the whole brain transcriptomes of 24 full-sib pairs of hens from two layer lines divergently selected for FP behavior. As light intensity is known to trigger the actual behavior given the respective motivation [27, 28], we not only compared transcriptomes between lines but also before and after light stimulation. To our knowledge, this is the first comprehensive brain transcriptome study using RNA-seq with respect to feather pecking in layers. Our key results point to an involvement of cholinergic signaling and a contribution of the immune system, which complements with previous studies and might be a step towards an integrated model for FP.

Results

Transcriptome and differential expression analyses

The average number of reads per sample was 48,502,572 (SD = 6,816,183, MIN = 36,337,624, MAX = 69,407,025) with an average mapping efficiency of 98.4% after alignment to the reference genome with TopHat (detailed summary in Additional file 1). A total of 13,070 genes were expressed in the brains of the 48 analyzed animals (average FPKM > 1, FPKM values are summarized in Additional file 2). Differential expression analysis revealed that 423 genes were differentially expressed (DE) between high feather peckers (HFP) and low feather peckers (LFP) (abs. LogFC ≥1, padj < 0.01) (Figs. 1a and 2a). Among those, we found 305 genes of uncertain function (LOC symbols), of which 156 are classified as non-coding RNAs (ncRNAs). Light stimulation led to the discovery of 28 DE genes (abs. LogFC ≥0.5, padj < 0.01) (Figs. 1b and 2b) with no differences between HFP and LFP (summary of differential expression analysis in Additional file 3). The top 40 DE genes with the lowest padj between HFP and LFP are plotted in a Heatmap (Fig. 3). A table of putative enhancer RNAs (eRNAs; Table 1) contains the closest DE protein coding genes (abs. LogFC ≥0.5, padj < 0.01) to DE ncRNAs.

Fig. 1
figure 1

Number of significantly upregulated genes between (a) high feather peckers (HFP) and low feather peckers (LFP) as well as (b) laying hens after light stimulation (light) and before light stimulation (dark)

Fig. 2
figure 2

Volcano plots of DE genes between (a) high pecker and low pecker hens and between (b) chickens before and after light stimulation

Fig. 3
figure 3

Heatmap of the top 40 DE genes between HFP and LFP hens with the lowest padj. The average log(FPKM) is shown in the heatmap

Table 1 Putative enhancer RNAs and their distance to possible targets

Detailed information on putative eRNAs are available in in Additional file 4.

Functional and pathway analyses

The results of the Gene Ontology (GO) and KEGG pathway analyses with DE genes between HFP and LFP hens with clusterProfiler [29] are shown in Table 2 (complete results in Additional file 5). Most of the GO terms assigned to the comparison of HFP and LFP lines are linked to synaptic functions and membrane channels when DE genes with an absolute LogFC < 0.5 were neglected. In particular, the DE genes CHRNA2, CHRNA9, CHRNB3, and CHRNB4 were assigned to these GO terms.

Table 2 Gene set analysis of DE genes between HFP and LFP hens with clusterProfiler. The top 5 results of each category with a q value below 0.2 are shown

Concerning the comparison of DE genes between animals before and after light stimulation (Table 3), the majority of IDs belonging to the category GO biological processes were linked to circadian rhythm. The DE genes ARNTL, NFIL3, NPAS2, and PER2 are overrepresented in these GO terms. The category GO molecular function exclusively contains IDs related to gene expression.

Table 3 Gene set analysis of DE genes before and after light stimulation with clusterProfiler. The top 5 results of each category with a q value below 0.2 are shown

Geneset analysis with STRING [30] was performed on DE genes between HFP and LFP lines, where genes with an absolute LogFC < 1 were neglected (Table 4, complete results see Additional file 6). The majority of terms found in this analysis are linked to the immune system and immune response.

Table 4 DE genes between high feather peckers (HFP) and low feather peckers (LFP) were analyzed with STRING. The top 5 results of each category are shown

Analysis of DE genes before and after light stimulation with STRING (Table 5) led to the discovery of terms related to biological/circadian rhythm and gene expression as well.

Table 5 DE genes between chickens before and after light stimulation were analyzed with STRING. The top 5 results of each category are shown

Protein interaction network analysis with STRING revealed one large and one smaller interaction cluster in the comparison of HFP and HFP DE genes (Fig. 4a). The large cluster, consisting of the DE gene FIP1L1 as well as SSU72, CPSF2, CPSF3, CPSF4, WDR33, CDC5L, PCF11, CLP1, CSTF1, CSTF2, CSTF3, and PAPOLA, exclusively contains genes, which are responsible for mRNA processing as well as mRNA binding. The slightly smaller cluster with fewer connections exclusively contains genes, which are associated with the immune system - T-cell proliferation and regulation in particular. The immune cluster consists of the DE genes CD28 and LIF as well as CD274, CTLA4, CD80, CD86, IL2B, IL12A, IL6ST, and LIFR. Analysis of DE genes before and after light stimulation led to the discovery of one interaction cluster, which contains the DE genes ARNTL, PER2, PER3, NFIL3, NPAS2, and LONRF3, all of which, except LONFRF3, are core components of the circadian clock (Fig. 4b).

Fig. 4
figure 4

Protein interaction networks of DE genes between (a) high and low feather peckers and between (b) animals before and after light stimulation

Discussion

In this study, we compared 24 full-sib pairs of laying hens selected for HFP and LFP in a transcriptomics approach. Besides the line comparison, we assessed the transcriptional changes before and after light stimulation. To our knowledge, this is the first brain transcriptomic study in FP that used RNA-seq. Comparison with the few other available transcriptome studies revealed only little overlap, which might be due to the fact that these studies are microarray-based and applied different designs. Hughes and Buitenhuis [31] used a more complex phenotyping scheme distinguishing between severe and gentle FP, Brunberg et al. [23] particularly analyzed hypothalamic RNA and Wysocki et al. [22] considered much older birds. Some genes, however, that these authors reported were also found DE in the current study. Notably, Hughes and Buitenhuis reported the expression of IRF2 (interferon regulatory factor 2) to be associated with severe FP This gene plays a major role in immune response [32] and might match our results with respect to immune related genes. In our study, however, we did not identify this gene, but found IRF5 (padj = 0.0002), IRF6 (padj = 1.17E-06), and IRF8 (padj = 0.0002) to be DE between the lines. Wysocki et al. reported DE of genes related to monoamine signaling, namely MAOA (monoamine oxidase A) and HTRR1B (5-hydroxytryptamine (serotonin) receptor 1B). We were able to confirm MAOA (padj = 7.2E-05), but not HTRR1B. Instead, we found DE of HTR1A (padj = 0.04), HTR1E (padj = 0.001) and HTR6 (padj = 2.34E-05).

Most of the GO terms found enriched in the current study in genes DE between HFP and LFP lines are linked to synaptic function when a LogFC threshold of 0.5 is used. This, in particular, is due to the DE genes CHRNA2, CHRNA9, CHRNB3, and CHRNB4 coding for subunits of nicotinic acetylcholine receptors (nAchR). Most other studies found serotonergic or dopaminergic signaling to be involved in FP [18] and we recently reported a possible involvement of GABAergic signaling based on genomic studies in the very same population as used in the current study [21, 25]. Cholinergic signaling has not been implicated in FP so far, but is linked to both monoamine as well as GABAergic signaling and is implicated in certain psychiatric disorders [33]. In particular, α7 nAChRs have a crucial role in the dysfunction of cortical parvalbumin-positive GABAergic neurons, as seen in schizophrenia [34]. The activation of the dopamine system by stimulating α7 nAChRs is crucial in drug dependence and withdrawal effects including motor function and rewarding properties of drug intake [35, 36]. In turn, it was shown that α7 nAChR blockade in mice can reduce anxiety-like behaviors as a consequence of nicotine withdrawal [37]. Nicotinic α4β2 receptors, on the other hand, have been implicated in obsessive-compulsive disorder and it has been shown that positive allosteric modulation of these receptors attenuated the obsessive behavior [38]. The actual behavioral effects of nAChRs are mediated by dopamine and glutamate and ultimately by GABAergic neurons [39, 40]. Very recently, it was shown in rats that β2 nAChRs on dopamine neurons in the ventral tegmental area mediate nicotine’s conditioned aversive effects, while those on GABA neurons mediate the conditioned rewarding effects [41]. In our study, we did not find the same subunits as described above to be differentially expressed, but all other studies so far have been conducted in humans, mice or rats and little is known in chicken. Thus, the findings might possibly point to a model of FP, in which cholinergic, monoaminergic and GABAergic signaling in the central nervous system are interacting. If this holds true, it seems likely that nAChRs might be on top of this signaling hierarchy. In a second gene set analysis, we used a more stringent set of genes by setting a LogFC threshold of 1.0. Interestingly this led to a shift in enriched terms compared to the aforementioned analysis to immunoglobulin-associated genes. We detected a protein interaction cluster consisting of 10 genes (the DE genes CD28 and LIF as well as CD274, CTLA4, CD80, CD86, IL2B, IL12A, IL6ST, and LIFR), all of which are linked to T-cell proliferation and regulation (Fig. 4a). In a recent study on post-mortem human brains, an increased number of T-cells was found in the brains of schizophrenia patients [42]. A probable mechanism, by which an increase in the number of T cells could influence the propensity to FP has been discovered in a murine cell culture system. Mashimo et al. demonstrated that T cells are able to synthesize acetylcholine (ACh) but also respond to it via nAChRs by increased cell proliferation and increased Ca2+ levels [43]. We propose a mechanism in which an overrepresentation of T cells leads to excess ACh production in the brain that has two major effects: (i) nAChR activation in neurons and (ii) further increase of T cell proliferation in a positive feedback loop. Since the immune-modulating genes show a higher LogFC compared to genes encoding nAChR subunits, we hypothesize that the dysregulation of T cell proliferation and activation is the initial cause of this neuropsychiatric phenotype. A recent extensive transcriptome study in humans suffering from autism spectrum disorder, schizophrenia, and bipolar disorder identified pathways related to the immune system and to transmembrane transporters/receptors to be differentially regulated in comparison to healthy individuals [44], which indicates that similar mechanisms are responsible for neuropsychiatric disorders between species. Interestingly these researchers found a large amount of DE ncRNAs between patients and controls, which they termed “psychiatric ncRNAs”. With 36.9%, ncRNAs make up a large part of significant DE transcripts between HFP and LFP hens. Gandal et al. already proposed these putative “psychiatric ncRNAs” to be responsible for transcriptome dysregulation by regulating local splicing events. This might be connected to the observation that HFP animals show a globally reduced variance in gene expression [31].

However, comparison of cDNA sequences of the DE ncRNAs identified in this study with the “psychiatric ncRNAs” from Gandal et al. with discontiguous megablast [45] showed sequence homology between only one pair of ncRNAs: The human ncRNA ENSG00000234773, which is a novel zinc finger protein pseudogene, is partially homologous to LOC112533587. Due to the taxonomic distance, the huge difference in genome size, and the fact that ncRNAs are less evolutionary conserved than protein coding genes [46], we did not expect to find homologous ncRNAs between Homo sapiens and Gallus gallus. By finding the closest DE expressed protein coding genes to DE ncRNAs between HFP and LFP chickens we were able to identify 21 putative eRNAs (Table 1). Among possible targets of these eRNA candidates are ARHGEF38 and SH3YL1, which have been linked to bipolar disorder and suicide [47, 48], GLRA1 and MCTP2 were found to be associated with schizophrenia [49, 50], and CHIR-B3, MHCIA7, MHCIY as well as MICA are genes involved in the immune system [51]. Thus, further experiments, like overexpression studies, will help to clarify the involvement of particular ncRNAs in FP. The second objective of this study was to clarify if HFP and LFP lines respond differently to light stimuli. Shi et al. demonstrated that red light and low light intensity was shown to reduce FP behavior during the laying period [28]. Light stimulation led to DE of five genes, which are core components of circadian clock (Fig. 4b) but there were no differences detectable in the light response between HFP and LFP.

Whereas in the current study whole brains where analyzed, it could be especially useful to restrict analyses to particular brain regions e.g. arcopallium or caudocentral nidopallium where differences in the serotonin or dopamine metabolism between LFP and HFP were found [52, 53]. Furthermore, analyses should be conducted in a larger number of birds and accounting for the actual phenotype, i.e. comparing performers against non-performer. Here, we compared divergently selected lines, which might also differ in other traits than FP [54]. Although the lines were initially derived from the same founder population and selection was solely based on the estimated breeding value for feather pecking, genetic drift may have caused allele frequency differences between the lines [25]. Here, actual behavior was highly correlated with line with on average 1.6 bouts per bird in LFP and 12.7 bouts per bird in HFP animals during the standardized observation time. A phenotype-based approach would have resulted in an almost identical group assignment. Futures studies should aim to validate the results in a purely phenotype-based approach using a single and unselected field population. This is, however, not trivial, because the occurrence of FP cannot be reliably predicted in commercial flocks and phenotyping is challenging under field conditions, which is the reason for many studies being conducted in selection lines.

Conclusions

The analysis of whole-brain transcriptomes in 24 full-sib pairs of hens from two White Leghorn layer strains divergently selected for FP behavior revealed differentially expressed genes and pathways related to cholinergic signaling and immune response. Thus, we hypothesize an involvement of these genes and pathways in the development of FP. Previous studies unambiguously identified monoamine signaling as the key component in FP behavior, but also pointed to an involvement of GABAergic signaling, immune response and the gut microbiota. Our results provide a first building block for an integrated model of FP. In this model, the gut microbiota might be involved in the interaction between the peripheral and the central serotonergic system as well as in the modulation of the immune response. Via the latter route, there could be a connection with cholinergic signaling, which in turn projects into monoamine signaling, directly or via GABA or glutamate signaling. In coming studies, the current results have to be validated using alternative approaches and appropriate functional studies need to be designed and conducted to further challenge and develop this model.

Methods

Experimental animals and sample collection

Experiments were conducted in White Leghorn layer strains divergently selected for FP behavior [2, 55]. Strains were developed from a single founder population and selection over 15 generations was solely based on estimated breeding values for feather pecking [17]. These lines were created and are maintained at the Hohenheim University and neither commercially obtained nor from a private source. Details on hatching, rearing, and husbandry have previously been described in detail [25, 56]. A total of 48 hens comprising 12 full-sib pairs from each strain (high feather peckers, HFP vs. low feather peckers, LFP) were phenotyped according to established protocols at 27 weeks of age [25]. Briefly, animals were observed on four consecutive days in sessions of 20 min by at least six different experienced observers. FP was defined as non-aggressive severe pecks or pulls directed to the plumage of conspecifics [57]. A series of pecks delivered in a short sequence without changing the behavior was recorded as a single occurrence called a bout per bird. Subsequently, birds were kept under low light conditions to prevent FP (Fig. 5). One bird from each full-sib pair was then sacrificed and brains were immediately collected for RNA isolation. Chickens were CO2-stunned and sacrificed by ventral neck cutting. The remaining birds were kept under increased light intensity (≥100 lx) for several hours until they clearly showed FP and were then sacrificed as well and brains were collected for RNA isolation. All brain samples were kept in RNA later and frozen at − 80 °C until further processing.

Fig. 5
figure 5

Experimental design of the study

RNA isolation

Whole chicken brains were pulverized with a Retsch mixer mill MM400 in liquid nitrogen at a frequency of 30 Hz for 1 min and immediately frozen at − 80 °C. For RNA isolation 900 μl Qiazol (QIAGEN), 5 μl DX reagent and 20–30 ceramic beads were added to10–30 mg of pulverized brain. The mixture was placed in a Bead Ruptor (Omin Inc.) and homogenized with the following program: 4.5 m/sec, 3 times 15 s, 10s dwell. After 5 min incubation at room temperature 100 μl gDNA Eliminator solution (QIAGEN) was added and the sample was vortexed for 15 s. After addition of 180 μl chloroform the sample was vortexed for 15 s and incubated at room temperature for 3 min. Sample was centrifuged at 12,000 g and 4 °C for 20 min and aqueous phase was transferred to a fresh reaction tube. An equal amount of 70% ethanol was added to the sample and vortexed for 15 s. RNA was purified using the RNeasy Plus Universal Mini Kit (QIAGEN) according to the manufacturer’s protocol and stored at − 80 °C.

NGS library preparation and sequencing

Input total RNA was visualized on a TapeStation 4200 (Agilent) and all but three samples showed excellent RNA quality with RIN-scores > 8 (three samples with RIN scores of 6.8, 6.9 and 7.7 respectively). Samples were quantified on a Qubit 2.0 Fluorometer (ThermoFischer) and 500 ng of total RNA was used as input for the TruSeq stranded mRNA library kit (Illumina) following the manufacturers manual. Resulting libraries showed a fragment size distribution of around 300 bp and were 2 × 75 bp paired-end. Sequencing was conducted on the HiSeq 4000 (Illumina) with 10 samples / lane.

Transcriptome analyses

Raw sequencing reads were processed with Trimmomatic version 0.36 [58] for adapter removal, trimming of low quality base calls, and removal of low-quality reads. Trimmomatic was used with the settings: PE -phred33 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36. Read pairs were discarded if one read did not survive quality control. Trimmed reads were aligned to chicken genome version GCF_000002315.5 (see Availability of data and materials) using TopHat version 2.1.0 [59] with the settings --no-novel-juncs --min-isoform-fraction 0.0 --min-anchor-length 3 -r 192 and GCF_000002315.5.gff (see Availability of data and materials) as known transcript file. Genomic features were extracted from the general feature format file and grouped with the R package GenomicFeatures (Version 1.40.0) [60]. summarizeOverlaps from the R package GenomicAlignments (Version 1.24.0) was used to count exon spanning reads [60]. Differential expression between experimental groups was analyzed with DESeq2 (Version 1.28.1) [61]. Volcano plots of differential expression analyses were created with the R package EnhancedVolcano (Kevin Blighe, Sharmila Rana and Myles Lewis (2019). EnhancedVolcano (Version 1.6.0). All R packages were obtained with Bioconductor version 3.11.

Identification of eRNA candidates

One BED file was created containing DE ncRNAs and one BED file containing all protein coding genes annotated in chicken genome version GCF_000002315.5 (see Availability of data and materials). The closest coding gene to each respective ncRNA was identified with BEDtools (Version 2.27.1) [62].

Functional analyses

Gene set analyses were conducted using the R package clusterProfiler (Version 3.16.0) [29] along with the chicken genome annotation org.Gg.eg.db (Version 3.11.4) [63] and using default settings in terms of gene set size and p-value correction. Genes differentially expressed with an absolute LogFC ≥0.5 were analyzed for enrichment in GO terms and KEGG pathways against the background of all expressed genes in our data set. For the gene set analysis with STRING (Version 11.0) [8], DE genes between HFP and LFP laying lines with an absolute LogFC ≥1 and an adjusted p-value below 0.01 were used. Protein interaction networks for these genes were also computed with STIRNG using the following settings: evidence-based coloring of network edges, inclusion of all interaction sources, medium required interaction score 0.4, max number of interactions to show 1st shell: no more than 10 interactions, 2nd shell: no more than 10 interactions. For DE genes before and after light stimulation LogFC threshold was 0.5. Settings for protein interaction networks were as follows: evidence based coloring of network edges, inclusion of all interaction sources, medium required interaction score 0.4, max number of interactions to show 1st shell: none / query proteins only, 2nd shell: none.

Availability of data and materials

RNA-seq reads were mapped to chicken genome version GCF_000002315.5 (RefSeq assembly, https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/002/315/GCF_000002315.5_GRCg6a/GCF_000002315.5_GRCg6a_genomic.fna.gz). Gene features were derived from GCF_000002315.5.gff (https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/002/315/GCF_000002315.5_GRCg6a/GCF_000002315.5_GRCg6a_genomic.gff.gz). The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. The raw sequencing data of the transcriptome experiments including phenotype and treatment information are accessible via BioProject ID PRJNA656654.

Abbreviations

ACh:

Acetylcholine

DE:

Differentially expressed

eRNAs:

Enhancer RNAs

FP:

Feather pecking

FPKM:

Fragments per kilobase of exon model per million reads mapped

GO:

Gene ontology

HFP:

High feather pecker

LFP:

Low feather pecker

LogFC:

Log fold change

nAchR:

Nicotinic acetylcholine receptors

ncRNAs:

Non-coding RNAs

padj :

Adjustet p-value

RNA-seq:

RNA-sequencing

References

  1. Rodenburg TB, Buitenhuis AJ, Ask B, Uitdehaag KA, Koene P, van der Poel JJ, Bovenhuis H. Heritability of feather pecking and open-field response of laying hens at two different ages. Poult Sci. 2003;82:861–7. https://0-doi-org.brum.beds.ac.uk/10.1093/ps/82.6.861.

    Article  CAS  PubMed  Google Scholar 

  2. Grams V, Wellmann R, Preuß S, Grashorn MA, Kjaer JB, Bessei W, Bennewitz J. Genetic parameters and signatures of selection in two divergent laying hen lines selected for feather pecking behaviour. Genet Sel Evol. 2015;47:77. https://0-doi-org.brum.beds.ac.uk/10.1186/s12711-015-0154-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Kjaer JB, Sørensen P. Feather pecking behaviour in white leghorns, a genetic study. Br Poult Sci. 1997;38:333–41. https://0-doi-org.brum.beds.ac.uk/10.1080/00071669708417999.

    Article  CAS  PubMed  Google Scholar 

  4. Blokhuis HJ. Feather-pecking in poultry: its relation with ground-pecking. Appl Anim Behav Sci. 1986;16:63–7.

    Article  Google Scholar 

  5. Wechsler B, Huber-Eicher B. The effect of foraging material and perch height on feather pecking and feather damage in laying hens. Appl Anim Behav Sci. 1998;58:131–41.

    Article  Google Scholar 

  6. Kjaer JB. Feather pecking in domestic fowl is genetically related to locomotor activity levels: implications for a hyperactivity disorder model of feather pecking. Behav Genet. 2009;39:564–70. https://0-doi-org.brum.beds.ac.uk/10.1007/s10519-009-9280-1.

    Article  PubMed  Google Scholar 

  7. Lutz V, Kjaer JB, Iffland H, Rodehutscord M, Bessei W, Bennewitz J. Quantitative genetic analysis of causal relationships among feather pecking, feather eating, and general locomotor activity in laying hens using structural equation models. Poult Sci. 2016;95:1757–63. https://0-doi-org.brum.beds.ac.uk/10.3382/ps/pew146.

    Article  CAS  PubMed  Google Scholar 

  8. Bessei W, Kjaer JB. In: Poultry Research Foundation, editor. Feather pecking in layers-state of research and implications. Sydney: Australian Poultry Science Symposium; 2005.

    Google Scholar 

  9. Meyer B, Bessei W, Bessei AW, Vahjen W, Zentek J, Harlander-Matauschek A. Dietary inclusion of feathers affects intestinal microbiota and microbial metabolites in growing Leghorn-type chickens. Poult Sci. 2012;91:1506–13. https://0-doi-org.brum.beds.ac.uk/10.3382/ps.2011-01786.

    Article  CAS  PubMed  Google Scholar 

  10. Kriegseis I, Bessei W, Meyer B, Zentek J, Würbel H, Harlander-Matauschek A. Feather-pecking response of laying hens to feather and cellulose-based rations fed during rearing. Poult Sci. 2012;91:1514–21. https://0-doi-org.brum.beds.ac.uk/10.3382/ps.2011-01865.

    Article  CAS  PubMed  Google Scholar 

  11. Birkl P, Bharwani A, Kjaer JB, Kunze W, McBride P, Forsythe P, Harlander-Matauschek A. Differences in cecal microbiome of selected high and low feather-pecking laying hens. Poult Sci. 2018;97:3009–14. https://0-doi-org.brum.beds.ac.uk/10.3382/ps/pey167.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. van der Eijk JAJ, de Vries H, Kjaer JB, Naguib M, Kemp B, Smidt H, et al. Differences in gut microbiota composition of laying hen lines divergently selected on feather pecking. Poult Sci. 2019;98:7009–21. https://0-doi-org.brum.beds.ac.uk/10.3382/ps/pez336.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Meyer B, Zentek J, Harlander-Matauschek A. Differences in intestinal microbial metabolites in laying hens with high and low levels of repetitive feather-pecking behavior. Physiol Behav. 2013;110–111:96–101. https://0-doi-org.brum.beds.ac.uk/10.1016/j.physbeh.2012.12.017.

    Article  CAS  PubMed  Google Scholar 

  14. Rodenburg TB, van Krimpen MM, de Jong IC, de Haas EN, Kops MS, Riedstra BJ, et al. The prevention and control of feather pecking in laying hens: identifying the underlying principles. Worlds Poult Sci J. 2013;69:361–74. https://0-doi-org.brum.beds.ac.uk/10.1017/S0043933913000354.

    Article  Google Scholar 

  15. Kops MS, Kjaer JB, Güntürkün O, Westphal KGC, Korte-Bouws GAH, Olivier B, et al. Serotonin release in the caudal nidopallium of adult laying hens genetically selected for high and low feather pecking behavior: an in vivo microdialysis study. Behav Brain Res. 2014;268:81–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bbr.2014.03.050.

    Article  CAS  PubMed  Google Scholar 

  16. van Hierden YM, de Boer SF, Koolhaas JM, Korte SM. The control of feather pecking by serotonin. Behav Neurosci. 2004;118:575–83. https://0-doi-org.brum.beds.ac.uk/10.1037/0735-7044.118.3.575.

    Article  CAS  PubMed  Google Scholar 

  17. Coates MD, Tekin I, Vrana KE, Mawe GM. Review article: the many potential roles of intestinal serotonin (5-hydroxytryptamine, 5-HT) signalling in inflammatory bowel disease. Aliment Pharmacol Ther. 2017;46:569–80. https://0-doi-org.brum.beds.ac.uk/10.1111/apt.14226.

    Article  CAS  PubMed  Google Scholar 

  18. de Haas EN, van der Eijk JAJ. Where in the serotonergic system does it go wrong? Unravelling the route by which the serotonergic system affects feather pecking in chickens. Neurosci Biobehav Rev. 2018;95:170–88. https://0-doi-org.brum.beds.ac.uk/10.1016/j.neubiorev.2018.07.007.

    Article  CAS  PubMed  Google Scholar 

  19. Parmentier HK, Rodenburg TB, de Vries Reilingh G, Beerda B, Kemp B. Does enhancement of specific immune responses predispose laying hens for feather pecking? Poult Sci. 2009;88:536–42. https://0-doi-org.brum.beds.ac.uk/10.3382/ps.2008-00424.

    Article  CAS  PubMed  Google Scholar 

  20. van der Eijk JAJ, Verwoolde MB, de Vries Reilingh G, Jansen CA, Rodenburg TB, Lammers A. Chicken lines divergently selected on feather pecking differ in immune characteristics. Physiol Behav. 2019;212:112680. https://0-doi-org.brum.beds.ac.uk/10.1016/j.physbeh.2019.112680.

    Article  CAS  PubMed  Google Scholar 

  21. Lutz V, Stratz P, Preuß S, Tetens J, Grashorn MA, Bessei W, Bennewitz J. A genome-wide association study in a large F2-cross of laying hens reveals novel genomic regions associated with feather pecking and aggressive pecking behavior. Genet Sel Evol. 2017;49:18. https://0-doi-org.brum.beds.ac.uk/10.1186/s12711-017-0287-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Wysocki M, Preuss S, Stratz P, Bennewitz J. Investigating gene expression differences in two chicken groups with variable propensity to feather pecking. Anim Genet. 2013;44:773–7. https://0-doi-org.brum.beds.ac.uk/10.1111/age.12050.

    Article  CAS  PubMed  Google Scholar 

  23. Brunberg E, Jensen P, Isaksson A, Keeling L. Feather pecking behavior in laying hens: hypothalamic gene expression in birds performing and receiving pecks. Poult Sci. 2011;90:1145–52. https://0-doi-org.brum.beds.ac.uk/10.3382/ps.2010-00961.

    Article  CAS  PubMed  Google Scholar 

  24. Flisikowski K, Schwarzenbacher H, Wysocki M, Weigend S, Preisinger R, Kjaer JB, Fries R. Variation in neighbouring genes of the dopaminergic and serotonergic systems affects feather pecking behaviour of laying hens. Anim Genet. 2009;40:192–9. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1365-2052.2008.01821.x.

    Article  CAS  PubMed  Google Scholar 

  25. Iffland H, Wellmann R, Schmid M, Preuß S, Tetens J, Bessei W, Bennewitz J. Genomewide Mapping of Selection Signatures and Genes for Extreme Feather Pecking in Two Divergently Selected Laying Hen Lines. Animals (Basel). 2020;10:262. https://0-doi-org.brum.beds.ac.uk/10.3390/ani10020262.

    Article  Google Scholar 

  26. Biscarini F, Bovenhuis H, van der Poel J, Rodenburg TB, Jungerius AP, van Arendonk JAM. Across-line SNP association study for direct and associative effects on feather damage in laying hens. Behav Genet. 2010;40:715–27. https://0-doi-org.brum.beds.ac.uk/10.1007/s10519-010-9370-0.

    Article  CAS  PubMed  Google Scholar 

  27. Riber AB, Guzman DA. Effects of dark brooders on behavior and fearfulness in layers. Animals (Basel). 2016. https://0-doi-org.brum.beds.ac.uk/10.3390/ani6010003.

  28. Shi H, Li B, Tong Q, Zheng W, Zeng D, Feng G. Effects of LED light color and intensity on feather pecking and fear responses of layer breeders in natural mating Colony cages. Animals (Basel). 2019. https://0-doi-org.brum.beds.ac.uk/10.3390/ani9100814.

  29. Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7. https://0-doi-org.brum.beds.ac.uk/10.1089/omi.2011.0118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–13. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gky1131.

    Article  CAS  PubMed  Google Scholar 

  31. Hughes AL, Buitenhuis AJ. Reduced variance of gene expression at numerous loci in a population of chickens selected for high feather pecking. Poult Sci. 2010;89:1858–69. https://0-doi-org.brum.beds.ac.uk/10.3382/ps.2010-00827.

    Article  CAS  PubMed  Google Scholar 

  32. Chae M, Kim K, Park S-M, Jang I-S, Seo T, Kim D-M, et al. IRF-2 regulates NF-kappaB activity by modulating the subcellular localization of NF-kappaB. Biochem Biophys Res Commun. 2008;370:519–24. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bbrc.2008.03.136.

    Article  CAS  PubMed  Google Scholar 

  33. Kim J, Yang JH, Ryu IS, Sohn S, Kim S, Choe ES. Interactions of Glutamatergic neurotransmission and brain-derived Neurotrophic factor in the regulation of behaviors after nicotine administration. Int J Mol Sci. 2019. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20122943.

  34. Lin H, Hsu F-C, Baumann BH, Coulter DA, Anderson SA, Lynch DR. Cortical parvalbumin GABAergic deficits with α7 nicotinic acetylcholine receptor deletion: implications for schizophrenia. Mol Cell Neurosci. 2014;61:163–75. https://0-doi-org.brum.beds.ac.uk/10.1016/j.mcn.2014.06.007.

    Article  CAS  PubMed  Google Scholar 

  35. Livingstone PD, Wonnacott S. Nicotinic acetylcholine receptors and the ascending dopamine pathways. Biochem Pharmacol. 2009;78:744–55. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bcp.2009.06.004.

    Article  CAS  PubMed  Google Scholar 

  36. Dani JA. Roles of dopamine signaling in nicotine addiction. Mol Psychiatry. 2003;8:255–6. https://0-doi-org.brum.beds.ac.uk/10.1038/sj.mp.4001284.

    Article  CAS  PubMed  Google Scholar 

  37. Jackson A, Papke RL, Damaj MI. Pharmacological modulation of the \alpha7 nicotinic acetylcholine receptor in a mouse model of mecamylamine-precipitated nicotine withdrawal. Psychopharmacology. 2018;235:1897–905. https://0-doi-org.brum.beds.ac.uk/10.1007/s00213-018-4879-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Mitra S, Mucha M, Khatri SN, Glenon R, Schulte MK, Bult-Ito A. Attenuation of compulsive-like behavior through positive allosteric modulation of α4β2 nicotinic acetylcholine receptors in non-induced compulsive-like mice. Front Behav Neurosci. 2016;10:244. https://0-doi-org.brum.beds.ac.uk/10.3389/fnbeh.2016.00244.

    Article  CAS  PubMed  Google Scholar 

  39. Ryu IS, Kim J, Seo SY, Yang JH, Oh JH, Lee DK, et al. Repeated Administration of Cigarette Smoke Condensate Increases Glutamate Levels and Behavioral Sensitization. Front Behav Neurosci. 2018. https://0-doi-org.brum.beds.ac.uk/10.3389/fnbeh.2018.00047.

  40. Ryu IS, Kim J, Seo SY, Yang JH, Oh JH, Lee DK, et al. Behavioral changes after nicotine challenge are associated with α7 nicotinic acetylcholine receptor-stimulated glutamate release in the rat dorsal striatum. Sci Rep. 2017;7:15009. https://0-doi-org.brum.beds.ac.uk/10.1038/s41598-017-15161-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Grieder TE, Besson M, Maal-Bared G, Pons S, Maskos U, van der Kooy D. β2* nAChRs on VTA dopamine and GABA neurons separately mediate nicotine aversion and reward. Proc Natl Acad Sci U S A. 2019;116:25968–73. https://0-doi-org.brum.beds.ac.uk/10.1073/pnas.1908724116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Sneeboer MAM, van Mierlo HC, Stotijn E, MacIntyre DJ, Smith C, Kahn RS, et al. Increased number of T-lymphocytes in post-mortem brain tissue of patients with schizophrenia. Schizophr Res. 2019. https://0-doi-org.brum.beds.ac.uk/10.1016/j.schres.2019.10.032.

  43. Mashimo M, Iwasaki Y, Inoue S, Saito S, Kawashima K, Fujii T. Acetylcholine released from T cells regulates intracellular Ca2+, IL-2 secretion and T cell proliferation through nicotinic acetylcholine receptor. Life Sci. 2017;172:13–8. https://0-doi-org.brum.beds.ac.uk/10.1016/j.lfs.2016.12.015.

    Article  CAS  PubMed  Google Scholar 

  44. Gandal MJ, Zhang P, Hadjimichael E, Walker RL, Chen C, Liu S, et al. Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder. Science. 2018. https://0-doi-org.brum.beds.ac.uk/10.1126/science.aat8127.

  45. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/25.17.3389.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Qu Z, Adelson DL. Evolutionary conservation and functional roles of ncRNA. Front Genet. 2012;3:205. https://0-doi-org.brum.beds.ac.uk/10.3389/fgene.2012.00205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Gaine ME, Seifuddin F, Sabunciyan S, Lee RS, Benke KS, Monson ET, et al. Differentially methylated regions in bipolar disorder and suicide. Am J Med Genet B Neuropsychiatr Genet. 2019;180:496–507. https://0-doi-org.brum.beds.ac.uk/10.1002/ajmg.b.32754.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Gaynor SC, Monson ET, Gaine ME, Chimenti MS, Reichman RD, Parsons M, et al. Male-specific association of the 2p25 region with suicide attempt in bipolar disorder. J Psychiatr Res. 2020;121:151–8. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jpsychires.2019.11.009.

    Article  PubMed  Google Scholar 

  49. Vora AK, Fisher AM, New AS, Hazlett EA, McNamara M, Yuan Q, et al. Dimensional traits of Schizotypy associated with Glycine receptor GLRA1 polymorphism: an exploratory candidate-Gene Association study. J Personal Disord. 2017;32:421–32. https://0-doi-org.brum.beds.ac.uk/10.1521/pedi_2017_31_303.

    Article  Google Scholar 

  50. Djurovic S, Le Hellard S, Kähler AK, Jönsson EG, Agartz I, Steen VM, et al. Association of MCTP2 gene variants with schizophrenia in three independent samples of Scandinavian origin (SCOPE). Psychiatry Res. 2009;168:256–8. https://0-doi-org.brum.beds.ac.uk/10.1016/j.psychres.2008.08.007.

    Article  CAS  PubMed  Google Scholar 

  51. Maenaka K, Jones EY. MHC superfamily structure and the immune system. Curr Opin Struct Biol. 1999;9:745–53. https://0-doi-org.brum.beds.ac.uk/10.1016/s0959-440x(99)00039-1.

    Article  CAS  PubMed  Google Scholar 

  52. Kops MS, de Haas EN, Rodenburg TB, Ellen ED, Korte-Bouws GAH, Olivier B, et al. Effects of feather pecking phenotype (severe feather peckers, victims and non-peckers) on serotonergic and dopaminergic activity in four brain areas of laying hens (Gallus gallus domesticus). Physiol Behav. 2013;120:77–82. https://0-doi-org.brum.beds.ac.uk/10.1016/j.physbeh.2013.07.007.

    Article  CAS  PubMed  Google Scholar 

  53. Kops MS, Kjaer JB, Güntürkün O, Westphal KGC, Korte-Bouws GAH, Olivier B, et al. Brain monoamine levels and behaviour of young and adult chickens genetically selected on feather pecking. Behav Brain Res. 2017;327:11–20. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bbr.2017.03.024.

    Article  CAS  PubMed  Google Scholar 

  54. Kjaer JB, Sørensen P, Su G. Divergent selection on feather pecking behaviour in laying hens (Gallus gallus domesticus). Appl Anim Behav Sci. 2001;71:229–39. https://0-doi-org.brum.beds.ac.uk/10.1016/s0168-1591(00)00184-2.

    Article  PubMed  Google Scholar 

  55. Bessei W, Bauhaus H, Bögelein S. The effect of selection for high and low feather pecking on aggression – related behaviours of laying hens. Archiv für Geflügelkunde. 2013;77:10–4.

    Google Scholar 

  56. Bennewitz J, Bögelein S, Stratz P, Rodehutscord M, Piepho HP, Kjaer JB, Bessei W. Genetic parameters for feather pecking and aggressive behavior in a large F2-cross of laying hens using generalized linear mixed models. Poult Sci. 2014;93:810–7. https://0-doi-org.brum.beds.ac.uk/10.3382/ps.2013-03638.

    Article  CAS  PubMed  Google Scholar 

  57. Grams V, Bögelein S, Grashorn MA, Bessei W, Bennewitz J. Quantitative genetic analysis of traits related to fear and feather pecking in laying hens. Behav Genet. 2015;45:228–35. https://0-doi-org.brum.beds.ac.uk/10.1007/s10519-014-9695-1.

    Article  PubMed  Google Scholar 

  58. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/btu170.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/btp120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Lawrence M, Huber W, Pagès H, Aboyoun P, Carlson M, Gentleman R, et al. Software for computing and annotating genomic ranges. PLoS Comput Biol. 2013;9:e1003118. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pcbi.1003118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. https://0-doi-org.brum.beds.ac.uk/10.1186/s13059-014-0550-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/btq033.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Carlson M. org.Gg.eg.db: Genome wide annotation for Chicken: R package; 2019.

    Google Scholar 

Download references

Acknowledgements

The DFG Competence Centre for Genome Analysis Kiel (CCGA) is greatly acknowledged for the cooperation in RNA-sequencing. We further acknowledge support by the Open Access Publication Funds of the Göttingen University. We thank Judith Beier and Hanna Iffland for participation in animal slaughtering and sample preparation. The data presented here has been previously presented at the ISAG 2019 conference.

Funding

The study was funded by the German Research Foundation (DFG) under file numbers TE622/4–2 and BE3703/8–2. The funders had no role in study design, data collection and analysis, and interpretation of data and in writing the manuscript. Publication fee was covered by the Open Access Publication Funds of the Göttingen University. Open Access funding provided by Projekt DEAL.

Author information

Authors and Affiliations

Authors

Contributions

CFG performed all bioinformatics analyses and wrote the manuscript. AM participated in animal slaughtering and brain extraction, and performed the RNA extraction. SP participated in animal slaughtering and revised the manuscript. SF performed the sequencing. WB and JB developed the project outline. JT developed the project outline, developed data analysis strategies, and wrote the manuscript. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Clemens Falker-Gieske.

Ethics declarations

Ethics approval and consent to participate

The research protocol was approved by the German Ethical Commission of Animal Welfare of the Provincial Government of Baden-Wuerttemberg, Germany (code: HOH 35/15 PG, date of approval: April 25, 2017).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Transcriptome alignment summary.

Additional file 2.

FPKM values.

Additional file 3.

Results of differential expression analyses.

Additional file 4.

Detailed information on eRNAs.

Additional file 5.

Results of clusterProfile analyses.

Additional file 6.

Results of STRING analyses.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Falker-Gieske, C., Mott, A., Preuß, S. et al. Analysis of the brain transcriptome in lines of laying hens divergently selected for feather pecking. BMC Genomics 21, 595 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-020-07002-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-020-07002-1

Keywords