Next Article in Journal
Non-Target Effects of dsRNA Molecules in Hemipteran Insects
Next Article in Special Issue
Molecular Regulation of Lipogenesis, Adipogenesis and Fat Deposition in Chicken
Previous Article in Journal
KLF15 Loss-of-Function Mutation Underlying Atrial Fibrillation as well as Ventricular Arrhythmias and Cardiomyopathy
Previous Article in Special Issue
High Levels of Genetic Variation in MHC-Linked Microsatellite Markers from Native Chicken Breeds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of SNPs Associated with Stress Response Traits within High Stress and Low Stress Lines of Japanese Quail

1
Department of Poultry Science, University of Arkansas, Fayetteville, AR 72701, USA
2
Oklahoma Medical Research Foundation, Oklahoma City, OK 73104, USA
3
Department of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea
*
Author to whom correspondence should be addressed.
Genes 2021, 12(3), 405; https://doi.org/10.3390/genes12030405
Submission received: 18 February 2021 / Revised: 5 March 2021 / Accepted: 5 March 2021 / Published: 12 March 2021
(This article belongs to the Special Issue Poultry Genetics, Breeding and Biotechnology)

Abstract

:
Mitigation of stress is of great importance in poultry production, as chronic stress can affect the efficiency of production traits. Selective breeding with a focus on stress responses can be used to combat the effects of stress. To better understand the genetic mechanisms driving differences in stress responses of a selectively bred population of Japanese quail, we performed genomic resequencing on 24 birds from High Stress (HS) and Low Stress (LS) lines of Japanese quail using Illumina HiSeq 2 × 150 bp paired end read technology in order to analyze Single Nucleotide Polymorphisms (SNPs) within the genome of each line. SNPs are common mutations that can lead to genotypic and phenotypic variations in animals. Following alignment of the sequencing data to the quail genome, 6,364,907 SNPs were found across both lines of quail. 10,364 of these SNPs occurred in coding regions, from which 2886 unique, non-synonymous SNPs with a SNP% ≥ 0.90 and a read depth ≥ 10 were identified. Using Ingenuity Pathway Analysis, we identified genes affected by SNPs in pathways tied to immune responses, DNA repair, and neurological signaling. Our findings support the idea that the SNPs found within HS and LS lines of quail could direct the observed changes in phenotype.

1. Introduction

Phenotypic diversities occur due to genetic variations within and between populations. Many genotypic variations arise as the result of genetic mutations. One of these mutations is known as a Single Nucleotide Polymorphism (SNP). SNPs are a common mutation and can be considered synonymous or non-synonymous [1]. Synonymous SNPs will either manifest within a noncoding region of a chromosome or silently within a coding region. Non-synonymous SNPs occur within coding regions of chromosomes and may be either missense or nonsense mutations. Both alter the amino acid at the point of mutation, though in differing ways. Missense mutations substitute one amino acid for another, while nonsense mutations stop translation entirely. Each type of mutation alters the protein produced [2]. These protein alterations may lead to alterations in phenotype. Due to this simple mechanism, SNPs are employed as biomarkers in studies surrounding population genetics [3].
Mitigation of stress is one of the primary goals of poultry producers. Stress has been shown to suppress the immune system’s function, decrease egg production, and slow body growth in chickens [4,5]. Hormones released during stress can affect the production of cytokines and chemokines by leukocytes, which have proinflammatory functions. Immunosuppression resulting from chronic stress has also been shown to decrease the effectiveness of vaccines, which would increase the occurrence of disease with a flock. Two different approaches are used for mitigating stress in chickens: improving the production environment and genetic selection for healthier stress responses [5].
This study aimed to identify genetic variations that could influence differential stress responses using a whole genome resequencing method. Two different lines of Japanese quail bred for differential stress responses were used. SNPs resulting in amino acid changes (nonsynonymous) were identified and their potential effects on functionalities were determined.

2. Materials and Methods

2.1. Ethics Statement

This study was conducted following the recommended guidelines for the care and use of laboratory animals for the National Institutes of Health. All procedures for animal care were performed according to the animal use protocols that was reviewed and approved by the University of Arkansas Institutional Animal Care and Use Committee (IACUC Protocol #14012).

2.2. Birds and DNA Sequencing

For the purposes of this study, two different lines of Japanese quail—a high stress (HS) line and a low stress (LS) line—have been utilized. These two lines have been selected since the 1980s for divergent blood corticosterone response to restraint-induced stress in order to serve as models for stress responses in poultry [6]. The differences between the two lines include higher sociability, lower fear, higher body weight, increased egg production, and a reduced heterophil/lymphocyte ratio in the LS line compared to the HS line. Additionally, the LS line had decreased stress-induced osteoporosis, accelerated onset of puberty, heightened male sexual activity and efficiency, and a decrease in corticosterone levels of approximately one-third when compared to the HS line [6].
The development of HS and LS lines of Japanese quail up to the 12th generation is documented by Satterlee and Johnson (1988) [7]. Each line was selected for its plasma corticosterone response to immobilization [6,7]. An independent random mating condition has since been used for their maintenance [8,9,10]. These research lines were shipped to the University of Arkansas at generation 44 from Louisiana State University and maintained at the Arkansas Agricultural Experimentation Station, Fayetteville, AR [11].
Bird selection and DNA sequencing was performed in direct accordance with Khatri et al. (2019) [6]. Male HS and LS birds were chosen for this study as they have a more stable physiology as a result of having fewer hormonal changes when compared to females. Blood samples (3 ml) from 20 birds each from HS and LS lines were collected. Genomic DNA was purified from each sample using QiaAmp DNA mini kit (Qiagen, Hilden, Germany) following manufacturer’s method. DNA quality was assessed using NanoDrop 1000 (Thermo Scientific, Waltham, MA, USA) and agarose gel electrophoresis. Twelve samples showing highest quality per line were pooled to represent each line. Library preparation and Illumina sequencing for the pooled DNA samples were performed by the Research Technology Support Facility at Michigan State University (East Lansing, MI, USA) using Illumina HiSeq 2 × 150 bp paired end read technology [6].

2.3. Data Quality and Assessment

We used the FastQC program (v0.11.6) (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/; accessed on 5 January 2018) to assess the quality of raw reads obtained after sequencing in form of FASTQ files. After quality assessment, the low quality reads were trimmed out using Trimmomatic tool (v0.32) [12]. The clean reads were then mapped onto the Japanese quail reference genome obtained from NCBI (Bioproject PRJNA292031; https://www.ncbi.nlm.nih.gov/genome/113; accessed on 5 February 2019) [6].
For the reference-based genome alignment, the Seqman NGen genome sequence assembly program of the Lasergene software package (DNAStar, Madison, WI, USA) was used. Assembly parameters were as follows: File format, Binary Alignment Map (BAM); mer Size, 21; mer skip query, 2; minimum match percentage, 93; maximum gap size, 6; minimum aligned length, 35; match score, 10; mismatch penalty, 20; gap penalty, 30; SNP calculation method, diploid Bayesian; minimum SNP percentage, 5; SNP confidence threshold, 10; minimum SNP count, 2; minimum base quality score, 5. After assembly, the SeqMan Pro program of the Lasergene package (DNAStar) was used for further analyses [13,14].

2.4. SNP Detection and Analysis

Detection and analysis of SNPs followed the methods of Khatri et al. (2018) [14]. The JMP genomics (SAS Institute Inc., Cary, NC, USA) program was used for filtering unique SNPs between the HS and LS lines of Japanese quail. Single nucleotide polymorphism (SNPs) occurring in both HS and LS lines were removed, leaving behind the unique SNPs for each line. To identify highly fixed and homozygous SNPs, the SNPs were filtered based on SNP percentages (SNP%) of SNP% ≥0.90 (for example, number of SNP = 9 of read depth = 10) were chosen. SNPs that induce non-synonymous changes in CDS (coding DNA sequences; protein coding) regions were chosen and unique SNPs in either HS or LS showing ≥10 read depths were selected as reliable SNPs. To reduce false positives, reliable SNPs chosen by criteria described above were confirmed by double-checking the initial assembly results with alignment view in SeqMan Pro program of Lasergene package (DNAStar) [14].

2.5. Bioinformatic Pathway Analysis

Genes containing unique SNPs from the HS and LS lines were entered into Ingenuity Pathway Analysis (version 60467501; Ingenuity Systems; Qiagen, Redwood City, CA, USA) for variant effect analyses for functional interpretation of genes including SNPs. The top 5 canonical pathways were chosen based on p-value calculated by the Fisher’s exact test, which was provided by the IPA algorithm.

2.6. Sorting Intolerant from Tolerant (SIFT) Prediction

The SIFT prediction for amino acid changes to be deleterious was conducted following the instruction (https://sift.bii.a-star.edu.sg/; accessed on 1 March 2021). Briefly, since SIFT database for Japanese quail (Coturnix japonica) is not available, a new database for Japanese quail genome was generated using sequences in FASTA format for genome and GFF file retrieved from NCBI (https://www.ncbi.nlm.nih.gov/genome/?term=txid93934[orgn]; accessed on 1 March 2021) following instruction (https://sift.bii.a-star.edu.sg/sift4g/SIFT4G_codes.html; accessed on 1 March 2021). The newly generated database was, then, annotated with .VCF file containing SNP list causing non-synonymous mutations. The SIFT score ≤ 0.05 is considered as deleterious amino acid substitution according to the instruction (https://sift.bii.a-star.edu.sg/sift4g/AnnotateVariants.html; accessed on 3 March 2021).

3. Results

3.1. Genome Re-Sequencing and Distribution of SNPs

In accordance with Khatri et al. (2019) [6], we performed whole genome resequencing of pooled DNA samples from 12 birds each from HS and LS lines of the quail and produced ~250 and ~257 million reads of 150 bp respectively. Of those, ~85 and ~84 million reads were mapped to the reference genome (NCBI/C. japonica) and their respective depth of coverage reached to ~41× and ~42× for HS and LS (Table 1) [6].
Totally, 10,364 SNPs were uniquely found across the HS and LS line of Japanese quail with SNP% ≥ 0.90 and present within the CDS region. The HS line contained 6551 SNPs, and the LS line contained 3813 SNPs. These SNPs were further filtered to remove synonymous SNPs, leaving only those within the CDS region that would induce changes in amino acid sequences. This brought the number of SNPs to 1831 within the HS line and 1055 within the LS line (Supplementary Materials S1 and S2). This list of SNPs was further filtered to include genes that occurred uniquely within each line. There were 980 and 539 uniquely affected genes (which were used for Ingenuity Pathway Analysis) in the HS and the LS line, respectively. The SIFT prediction tool predicted 707 and 374 amino acid changes to be deleterious substitutions in HS and LS lines, respectively (Supplementary Materials S1 and S2). The 20 SNPs from each line with the highest read depths can be found in Table 2 and Table 3, while the occurrence of each type of impact of the unique SNPs can be found in Table 4. A previous study of SNPs within a chicken autoimmune vitiligo model identified similar mutations in coding regions, causing changes in amino acid sequences [13].

3.2. Bioinformatics and Pathway Analyses

IPA was used to analyze connections of genes containing SNPs to canonical pathways and changes in molecular interactions.
Glycoprotein VI (GP6) Signaling, Signaling by Rho Family GTPases, Rac Signaling, Amyotrophic Lateral Sclerosis Signaling, and Synaptogenesis Signaling were the top 5 canonical pathways found to be significantly affected within the HS line (p-value < 0.01) (Figure 1; Supplementary Material S3). In the LS line, the canonical pathways found to be significantly affected (p-value < 0.05) include GP6 Signaling, Activation of Interferon-Regulatory Factor (IRF) by Cytosolic Pattern Recognition Receptors, Th17 Activation, Role of BRCA2 (BRCA2, DNA repair associated) in DNA Damage Response, and IL-15 Production (Figure 2; Supplementary Material S4).
The Synaptogenesis signaling and Amyotrophic Lateral Sclerosis signaling pathways found in the HS line are both tied to nervous system development and signaling (Supplementary Materials S5 and S6). Affected genes within the Amyotrophic Lateral Sclerosis signaling pathway are involved in regulation of neuronal damage and degeneration, as well as motor neuron apoptosis, while affected genes within the synaptogenesis signaling pathway are tied primarily to cytoskeletal development mechanisms, such as microtubule stabilization and synaptic spine organization and development [15,16,17,18,19].
The GP6 Signaling pathway is found to contain SNPs in both lines of quail (Figure 3a,b). Each line shows SNPs within collagen, laminin, phosphoinositide 3-kinase (PI3K), and protein kinase C (PKC) gene families. In the HS line, Protein kinase C delta (PKCδ) is also affected. PKCδ is part of a signaling pathway with proinflammatory functions within the nervous system [20]. In the LS line, talin 1 (TLN1) and inositol 1,4,5-triphosphate receptor type 1 (ITPR1) are affected. TLN1 expression mediates leukocyte adhesion to platelets [21]. ITPR1, also known as IP3R1, is a regulator of B cell apoptosis [22]. GP6 is a member of the immunoglobulin superfamily. It is expressed within platelets and their precursor megakaryocytes, and it serves as the major signaling receptor for collagen, which leads to platelet activation and thrombus formation [23,24]. GP6 is activated by a few different factors, one of which being low shear stress [25,26,27]. Shear stress is one of many factors that can drive gene expression within endothelial cells of the circulatory system, and low shear stress is often an indication of improper blood flow [28].
The Th17 Activation pathway, which includes genes containing SNPs in the LS line, is an integral part of cellular immune response (Figure 4). In addition to controlling differentiation of Th17 cells, this pathway produces IL-17A, IL-17F, IL-21, and IL-22 cytokines [29,30]. Both IL-17 cytokines have pro-inflammatory functions, are important regulators of neutrophil activity, and help mount antifungal immune responses alongside IL-22 [31,32]. IL-21 induces differentiation of B cells and T cells [33]. Similarly, the activation of IRF pathway is also associated with an immune response. This pathway is primarily involved in a host’s response to viral infection, but its activation is nonetheless a result of a stimulated innate immune response (Supplementary Material S7) [34].
The canonical pathways observed within the HS line of quail show significant overlap, most of which is attributed to the presence of phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit delta (PIK3CD), phosphatidylinositol-4,5-bisphosphate 3-kinase regulatory subunit 3 (PIK3R3), and phosphatidylinositol-4,5-bisphosphate 3-kinase regulatory subunit 6 (PIK3R6), which are genes found to be active in all five pathways listed in Figure 1. These three genes code for a family of molecules known as phosphoinositide 3-kinases (PI3Ks). PI3Ks are part of the immune response. These kinases carry out phosphorylation within leukocytes and regulate the activity of Rho GTPases, which are important controllers of leukocyte motility [35].
Found within the Rac signaling, ALS signaling, signaling by Rho GTPases, and synaptogenesis signaling pathways is the gene p21-activated kinase 1 (PAK1), which encodes a p21-activated kinase (PAK) protein. PAK proteins form the bridge between Rho GTPases and the cytoskeleton, leading to changes in cell motility and growth. This bridge is necessary for successful thrombin-mediated platelet activation [36].
Mitogen-activated protein kinase 8 (MAPK8), neutrophil cytosolic factor 1 (NCF1), and nuclear factor kappa B subunit 1 (NFKB1) genes are found within both the signaling by Rho GTPases and Rac signaling pathways. MAPK8, also known as c-Jun N-terminal protein kinase (JNK1), belongs to the MAP kinase family and reacts to cellular stimuli to trigger early gene expression responses (Figure 5). In mice, this specific kinase family has been linked to T cell proliferation, differentiation, and apoptosis [37]. NCF1 is a subunit of NADPH oxidase, an enzyme complex that drives production of superoxide during an immune response [38]. NFKB1 is part of a protein complex known as nuclear factor kappa B (NF-κB) that acts as a transcriptional regulator integral in both innate and adaptive immune response. Issues with NF-κB activation can lead to inflammatory diseases [39].
Within the LS line, only one gene, Janus kinase 2 (JAK2), is found within multiple pathways. JAK2 is found within both the Th17 Activation and IL-15 production pathways (Figure 4; Supplementary Material S8). JAK2 codes for a kinase that is part of the JAK family of kinases. This kinase is integral for cell responses to interferon γ (IFNγ), a cytokine that is an important part of both innate and adaptive immune responses, thus cementing this gene within the immune response [40].
The LS line of quail contained multiple SNPs at various stages of a DNA damage response that responds to damage caused by ionizing radiation (Figure 6). Both a gene involved in the initial activation complex, Fanconi anemia complementation group F (FANCF), and genes involved in each of the end points of the pathway were found to contain SNPs. In the initial activation, FANCF was identified. BRCA2, Bloom syndrome RecQ-like helicase (BLM), SWItch/Sucrose Non-Fermentable (SWI/SNF), and E2F were genes identified at end points of the repair process. This DNA damage response mediates homologous recombination and chromatin remodeling, while also playing a role in tumor suppression [41,42,43].
Diseases and bio functions that were found to be significantly affected within the HS line (p-value < 0.05) include cancers, neurological disease, organismal injury and abnormalities, psychological disorders, and developmental disorders. In the LS line, the significantly affected (p-value < 0.05) diseases and biofunctions include cancers, and organismal injury and abnormalities, neurological disease, hematological disease, and metabolic disease (Supplementary Material S9).
The top physiological systems found to be significantly affected (p-value < 0.05) in development and function within the HS line include embryonic development, nervous system development and function, organismal development, tissue development, and cardiovascular system development and function. The top significantly affected (p-value < 0.05) systems within the LS line include nervous system development and function, tissue development, embryonic development, organismal development, and organ development (Supplementary Material S10).
As the nervous system plays an integral role in stress response signaling, further analysis of its affected developmental functions was performed. The factors involved in nervous system development and functions were listed in Supplementary Materials S11 and S12. The top 10 functions and the related genes affected by SNPs in each line are listed in Supplementary Materials S13 and S14. Significantly affected (p-value < 0.01) nervous system developmental functions in the HS line include branching, growth and outgrowth of neurites, neuritogenesis, dendritic growth and branching, proliferation of neuronal cells, polarization of hippocampal neurons, transport of synaptic vesicles, migration of cerebellar granule cells, and development of neurons (Supplementary Material S13). Significantly affected (p-value < 0.02) nervous system developmental functions in the LS line include function, maturation, and development of neurons, development of central nervous system, myelination of nerves and the sciatic nerve, formation of inhibitory synapse and the brain, long-term potentiation of mossy fibers, and quantity of dendrites (Supplementary Material S14).

4. Discussion

Our analysis found unique SNPs in pathways tied to a shear stress induced immune response, a DNA damage repair response network, and several other immune response pathways, including cytokine production and regulation and differentiation of immune cells. Additionally, pathways related to nervous system development were found to be significantly affected by SNPs.
Differing mutations within GP6 Signaling pathway could be an indication that the two lines of quail have different immune responses to changes in blood pressure [44]. It is possible these two lines may experience changes in blood pressure when exposed to stressors, but the LS line could more efficiently minimize the effects that would have affected endothelial cells. The location of SNPs in the BRCA2 DNA Damage Response is evidence that the LS line may have become better adapted at repairing cellular damage caused by stress. These two factors taken in tandem indicate that the two lines of quail primarily differ in how their bodies have adapted to control the impact of stress at the cellular level. The LS line, when exposed to the same stressors as the HS line, would not show the same negative impacts to productive traits should its physiology sufficiently compensate.
The LS line was found to have SNPs affecting multiple different pathways tied to the immune response. In particular, the presence of SNPs in the Th17 activation pathway and IL-15 production indicate a possibly direct effect on cytokine production. IL-17A, IL-17F, IL-21, and IL-22 are considered signature cytokines within the immune system with various functions each, including regulation and differentiation of immune cells. Changes in regulatory T-cell activity may protect the cells of the LS line of quail, as these cells suppress many harmful effects of immune responses [5].
The HS line also contained SNPs in genes related to immune function across multiple canonical pathways. GP6 signaling, Rho GTPase signaling, and Rac signaling pathways were found to be directly connected transitively through SNP mutations of PI3K genes and Rho GTPases. The connections of these kinases to leukocytes, as well as each pathway’s connection to inflammatory responses, once again show an association between the immune system and the SNP mutations in these birds.
Both lines of quail were selected based on their response to restraint stress, which is a psychological stressor. Psychological stressors are detected by the nervous system, which can, in turn, activate hormonal responses within the endocrine system. The hypothalamo-pituitary-adrenal (HPA) axis is one example of this type of response. The HPA axis is a driving force behind the endocrine response to stress, as it initiates a hormone cascade leading to the release of corticosterone hormone [45]. The HS and LS lines of Japanese quail show divergent blood corticosterone levels, indicating a lasting effect of genetic selection on a significant change in the response of the neuroendocrine system to stress. Our findings show that SNPs are present in genes throughout several nervous system developmental networks, though the extent to which these SNPs could alter nervous system function still remains unclear. Further studies are needed to understand the relationship of these genes containing SNPs to nervous system development associated with the differential stress responses observed in the HS and LS lines of Japanese quail. Comparison of SNPs from the brain transcriptome could also provide more insight into the effect of stress on neuronal gene expression.

5. Conclusions

Overall, these findings seem to support the idea that the higher body weight and increased egg production observed in the LS line of Japanese quail may be attributed to differing neuronal development and a stronger protection from stress induced immune responses. When paired with recent findings that have found CNVs affected genes within nervous and endocrine system development, humoral and cell-mediated immune responses, and various metabolisms, our results support the theory that genotypic variation leads to the observed phenotypic diversity between the HS and LS lines of Japanese quail [6]. Additionally, more research directed towards the immune responses and cellular activity of each line of quail could deepen our understanding of the effects chronic stress may have on differently adapted production birds.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/12/3/405/s1. S1: HS Line SNPs (Excel); S2: LS Line SNPs (Excel); S3: HS Line Canonical Pathway Genes (Excel); S4: LS Line Canonical Pathway Genes (Excel); S5: Synaptogenesis Signaling Pathway (Image); S6: Amyotrophic Lateral Sclerosis Signaling Pathway (Image); S7: Activation by IRF by Cytosolic Pattern Recognition Receptors Pathway (Image); S8: IL-15 Production Pathway (Image); S9: HS and LS Line Diseases and Bio Functions Affected by Unique SNPS (Image); S10: HS and LS Line Physiological Systems Affected by Unique SNPs (Image); S11: HS Line Nervous System Development and Function SNPs (Excel); S12: LS Line Nervous System Development and Function SNPs (Excel); S13: HS Line Top 10 Nervous System Developmental Functions Affected by Unique SNPS (Image); S14: LS Line Top 10 Nervous System Developmental Functions Affected by Unique SNPS (Image).

Author Contributions

Conceptualization, B.K. (Bhuwan Khatri), S.K., W.K. and B.K. (Byungwhi Kong); methodology, B.K. (Bhuwan Khatri), D.S., S.S. (Stephanie Shouse), S.K. and B.K. (Byungwhi Kong); software, S.S. (Steven Shumaker), D.S., S.S. (Stephanie Shouse), and B.K. (Byungwhi Kong); writing—original draft preparation, S.S. (Steven Shumaker); writing—review and editing S.S. (Steven Shumaker), B.K. (Bhuwan Khatri), S.S. (Stephanie Shouse), D.S., S.K., W.K. and B.K (Byungwhi Kong). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by USDA multistate project NC1170 Advanced Technologies for the Genetic Improvement of Poultry and Division of Agriculture, University of Arkansas System.

Institutional Review Board Statement

This study was conducted following the recommended guidelines for the care and use of laboratory animals for the National Institutes of Health. All procedures for animal care were performed according to the animal use protocols that was reviewed and approved by the University of Arkansas Institutional Animal Care and Use Committee (IACUC Protocol #14012).

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequence reads described in the manuscript are available under BioProject accession PRJNA706434 in NCBI (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA706434 accessed on 5 March 2021). Illumina sequence reads have been deposited at NCBI’s SRA under the BioProject.

Conflicts of Interest

We certify that there is no conflict of interest with any financial organization regarding the materials discussed in the manuscript.

References

  1. Shastry, B.S. SNPs: Impact on gene function and phenotype. Methods Mol. Biol. 2009, 578, 3–22. [Google Scholar] [CrossRef]
  2. Petersen, G.M.; Parmigiani, G.; Thomas, D. Missense mutations in disease genes: A bayesian approach to evaluate causality. Am. J. Hum. Genet. 1998, 62, 1516–1524. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Syvänen, A.-C. Accessing genetic variation: Genotyping single nucleotide polymorphisms. Nat. Rev. Genet. 2001, 2, 930–942. [Google Scholar] [CrossRef] [PubMed]
  4. McEwen, B.S.; Seeman, T. Protective and damaging effects of mediators of stress: Elaborating and testing the concepts of allostasis and allostatic load. Ann. N. Y. Acad. Sci. 1999, 896, 30–47. [Google Scholar] [CrossRef] [PubMed]
  5. Shini, S.; Huff, G.R.; Shini, A.; Kaiser, P. Understanding stress-induced immunosuppression: Exploration of cytokine and chemokine gene profiles in chicken peripheral leukocytes. Poult. Sci. 2010, 89, 841–851. [Google Scholar] [CrossRef]
  6. Khatri, B.; Kang, S.; Shouse, S.; Anthony, N.; Kuenzel, W.; Kong, B.C. Copy number variation study in Japanese quail associated with stress related traits using whole genome re-sequencing data. PLoS ONE 2019, 14, e0214543. [Google Scholar] [CrossRef]
  7. Satterlee, D.G.; Johnson, W.A. Selection of japanese quail for contrasting blood corticosterone response to immobilization. Poult. Sci. 1988, 67, 25–32. [Google Scholar] [CrossRef]
  8. Satterlee, D.G.; Cadd, G.G.; Jones, R.B. Developmental instability in japanese quail genetically selected for contrasting adrenocortical responsiveness. Poult. Sci. 2000, 79, 1710–1714. [Google Scholar] [CrossRef] [PubMed]
  9. Satterlee, D.G.; Marin, R.H.; Jonest, R.B. Selection of Japanese quail for reduced adrenocortical responsiveness accelerates puberty in males. Poult. Sci. 2002, 81, 1071–1076. [Google Scholar] [CrossRef] [PubMed]
  10. Suhr, C.L.; Schmidt, J.B.; Treese, S.T.; Satterlee, D.G. Short-term spatial memory responses in aged Japanese quail selected for divergent adrenocortical stress re-sponsiveness. Poult. Sci. 2010, 89, 633–642. [Google Scholar] [CrossRef] [PubMed]
  11. Huff, G.R.; Huff, W.E.; Wesley, I.V.; Anthony, N.B.; Satterlee, D.G. Response of restraint stress-selected lines of Japanese quail to heat stress and Escherichia coli challenge. Poult. Sci. 2013, 92, 603–611. [Google Scholar] [CrossRef] [Green Version]
  12. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  13. Jang, H.-M.; Erf, G.F.; Rowland, K.C.; Kong, B.-W. Genome resequencing and bioinformatic analysis of SNP containing candidate genes in the autoimmune vitiligo Smyth line chicken model. BMC Genom. 2014, 15, 707. [Google Scholar] [CrossRef] [Green Version]
  14. Khatri, B.; Hayden, A.M.; Anthony, N.B.; Kong, B.C. Whole Genome Resequencing of Arkansas Progressor and Regressor Line Chickens to Identify SNPs Associated with Tumor Regression. Genes 2018, 9, 512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Lee, M.K.; Marszalek, J.R.; Cleveland, D.W. A mutant neurofilament subunit causes massive, selective motor neuron death: Implications for the pathogenesis of human motor neuron disease. Neuron 1994, 13, 975–988. [Google Scholar] [CrossRef]
  16. Martin, L.J.; Chen, K.; Liu, Z. Adult motor neuron apoptosis is mediated by nitric oxide and fas death receptor linked by dna damage and p53 activation. J. Neurosci. 2005, 25, 6449–6459. [Google Scholar] [CrossRef] [Green Version]
  17. Chernoff, J.; Vadlamudi, R.K.; Adam, L.; Wang, R.-A.; Mandal, M.; Nguyen, D.; Sahin, A.; Hung, M.-C.; Kumar, R. Regulatable expression of p21-activated kinase-1 promotes anchorage-independent growth and abnormal organization of mitotic spindles in human epithelial breast cancer cells. J. Biol. Chem. 2000, 275, 36238–36244. [Google Scholar] [CrossRef] [Green Version]
  18. Gonzalez-Billault, C.; Avila, J.; Cáceres, A. Evidence for the role of MAP1B in axon formation. Mol. Biol. Cell 2001, 12, 2087–2098. [Google Scholar] [CrossRef] [Green Version]
  19. Beffert, U.; Dillon, G.M.; Sullivan, J.M.; Stuart, C.E.; Gilbert, J.P.; Kambouris, J.A.; Ho, A. Microtubule plus-end tracking protein CLASP2 regulates neuronal polarity and synaptic function. J. Neurosci. 2012, 32, 13906. [Google Scholar] [CrossRef] [Green Version]
  20. Pierchala, B.A.; Ahrens, R.C.; Paden, A.J.; Johnson, E.M. Nerve growth factor promotes the survival of sympathetic neurons through the cooperative function of the protein kinase C and phosphatidylinositol 3-kinase pathways. J. Biol. Chem. 2004, 279, 27986–27993. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Bromberger, T.; Klapproth, S.; Rohwedder, I.; Zhu, L.; Mittmann, L.; Reichel, C.A.; Sperandio, M.; Qin, J.; Moser, M. Direct Rap1/Talin1 interaction regulates platelet and neutrophil integrin activity in mice. Blood 2018, 132, 2754–2762. [Google Scholar] [CrossRef] [Green Version]
  22. Assefa, Z.; Bultynck, G.; Szlufcik, K.; Kasri, N.N.; Vermassen, E.; Goris, J.; Missiaen, L.; Callewaert, G.; Parys, J.B.; De Smedt, H. Caspase-3-induced truncation of type 1 inositol trisphosphate receptor accelerates apoptotic cell death and induces inositol trisphosphate-independent calcium release during apoptosis. J. Biol. Chem. 2004, 279, 43227–43236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Clemetson, J.M.; Polgar, J.; Magnenat, E.; Wells, T.N.; Clemetson, K.J. The platelet collagen receptor glycoprotein VI is a member of the immunoglobulin superfamily closely related to fcαr and the natural killer receptors. J. Biol. Chem. 1999, 274, 29019–29024. [Google Scholar] [CrossRef] [Green Version]
  24. Gibbins, J.M. Platelet adhesion signalling and the regulation of thrombus formation. J. Cell Sci. 2004, 117, 3415–3425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Nieswandt, B.; Bergmeier, W.; Eckly, A.; Schulte, V.; Ohlmann, P.; Cazenave, J.-P.; Zirngibl, H.; Offermanns, S.; Gachet, C. Evidence for cross-talk between glycoprotein VI and Gi-coupled receptors during collagen-induced platelet aggregation. Blood 2001, 97, 3829–3835. [Google Scholar] [CrossRef] [Green Version]
  26. Sakariassen, K.S.; Holme, P.A.; Ørvim, U.; Barstad, R.M.; Solum, N.O.; Brosstad, F.R. Shear-induced platelet activation and platelet microparticle formation in native human blood. Thromb. Res. 1998, 92, S33–S41. [Google Scholar] [CrossRef]
  27. Roberts, D.E.; McNicol, A.; Bose, R. Mechanism of collagen activation in human platelets. J. Biol. Chem. 2004, 279, 19421–19430. [Google Scholar] [CrossRef] [Green Version]
  28. Papaioannou, T.G.; Stefanadis, C. Vascular wall shear stress: Basic principles and methods. Hell. J. Cardiol. 2005, 46, 9–15. [Google Scholar]
  29. Wilson, N.J.; Boniface, K.; Chan, J.R.; McKenzie, B.S.; Blumenschein, W.M.; Mattson, J.D.; Basham, B.; Smith, K.J.; Chen, T.; Morel, F.; et al. Development, cytokine profile and function of human interleukin 17–producing helper T cells. Nat. Immunol. 2007, 8, 950–957. [Google Scholar] [CrossRef]
  30. Chen, Z.; Laurence, A.; O’Shea, J.J. Signal transduction pathways and transcriptional regulation in the control of Th17 differentiation. Semin. Immunol. 2007, 19, 400–408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Gaffen, S.L.; Kramer, J.M.; Yu, J.J.; Shen, F. The IL-17 cytokine family. In Vitamins and Hormones; Academic Press: Cambridge, MA, USA, 2006; pp. 255–282. [Google Scholar]
  32. Eyerich, S.; Eyerich, K.; Cavani, A.; Schmidt-Weber, C. IL-17 and IL-22: Siblings, not twins. Trends Immunol. 2010, 31, 354–361. [Google Scholar] [CrossRef] [PubMed]
  33. Leonard, W.J.; Wan, C.-K. IL-21 Signaling in Immunity. F1000Research 2016, 5, 224. [Google Scholar] [CrossRef]
  34. Weaver, B.K.; Kumar, K.P.; Reich, N.C. Interferon regulatory factor 3 and CREB-binding protein/p300 are subunits of double-stranded RNA-activated transcription factor DRAF1. Mol. Cell. Biol. 1998, 18, 1359–1368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Rickert, P.; Weiner, O.D.; Wang, F.; Bourne, H.R.; Servant, G. Leukocytes navigate by compass: Roles of PI3Kgamma and its lipid products. Trends Cell Biol. 2000, 10, 466–473. [Google Scholar] [CrossRef]
  36. Aslan, J.E.; Baker, S.M.; Loren, C.P.; Haley, K.M.; Itakura, A.; Pang, J.; Greenberg, D.L.; David, L.L.; Manser, E.; Chernoff, J.; et al. The PAK system links Rho GTPase signaling to thrombin-mediated platelet activation. Am. J. Physiol. Cell Physiol. 2013, 305, C519–C528. [Google Scholar] [CrossRef] [PubMed]
  37. Arbour, N.; Naniche, D.; Homann, D.; Davis, R.J.; Flavell, R.A.; Oldstone, M.B. c-Jun NH2-terminal kinase (JNK)1 and JNK2 signaling pathways have divergent Roles in CD8+ T Cell–mediated antiviral immunity. J. Exp. Med. 2002, 195, 801–810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Vignais, P.V. The superoxide-generating NADPH oxidase: Structural aspects and activation mechanism. Cell. Mol. Life Sci. 2002, 59, 1428–1459. [Google Scholar] [CrossRef]
  39. Liu, T.; Zhang, L.; Joo, D.; Sun, S.-C. NF-κB signaling in inflammation. Signal Transduct. Target. Ther. 2017, 2, 17023. [Google Scholar] [CrossRef] [Green Version]
  40. Ghoreschi, K.; Laurence, A.; O’Shea, J.J. Janus kinases in immune cell signaling. Immunol. Rev. 2009, 228, 273–287. [Google Scholar] [CrossRef] [Green Version]
  41. Wang, Y.; Cortez, D.; Yazdi, P.; Neff, N.; Elledge, S.J.; Qin, J. BASC, a super complex of BRCA1-associated proteins involved in the recognition and repair of aberrant DNA structures. Genes Dev. 2000, 14, 927–939. [Google Scholar] [PubMed]
  42. Bochar, D.A.; Wang, L.; Beniya, H.; Kinev, A.; Xue, Y.; Lane, W.S.; Wang, W.; Kashanchi, F.; Shiekhattar, R. BRCA1 is associated with a human SWI/SNF-related complex: Linking chromatin remodeling to breast cancer. Cell 2000, 102, 257–265. [Google Scholar] [CrossRef] [Green Version]
  43. Yu, X.; Baer, R. Nuclear localization and cell cycle-specific expression of CtIP, a protein that associates with the BRCA1 tumor suppressor. J. Biol. Chem. 2000, 275, 18541–18549. [Google Scholar] [CrossRef] [Green Version]
  44. Stegner, D.; Haining, E.J.; Nieswandt, B. Targeting glycoprotein vi and the immunoreceptor tyrosine-based activation motif signaling pathway. Arter. Thromb. Vasc. Biol. 2014, 34, 1615–1620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kuenzel, W.J.; Jurkevich, A. Molecular neuroendocrine events during stress in poultry. Poult. Sci. 2010, 89, 832–840. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Canonical pathways affected by unique SNPs in HS line of Japanese quail.
Figure 1. Canonical pathways affected by unique SNPs in HS line of Japanese quail.
Genes 12 00405 g001
Figure 2. Canonical pathways affected by unique SNPs in LS line of Japanese quail.
Figure 2. Canonical pathways affected by unique SNPs in LS line of Japanese quail.
Genes 12 00405 g002
Figure 3. (a) GP6 Signaling Pathway in HS line of Japanese quail. Molecular interactions within the given pathway are shown. Grayed symbols indicate genes containing SNP mutations while white symbols indicate genes that do not contain SNPs but are functionally associated within the pathway. Each symbol represents the given molecule’s function. (b) GP6 Signaling in LS line of Japanese quail. Molecular interactions and symbols within the pathway are the same as described in Figure 3a.
Figure 3. (a) GP6 Signaling Pathway in HS line of Japanese quail. Molecular interactions within the given pathway are shown. Grayed symbols indicate genes containing SNP mutations while white symbols indicate genes that do not contain SNPs but are functionally associated within the pathway. Each symbol represents the given molecule’s function. (b) GP6 Signaling in LS line of Japanese quail. Molecular interactions and symbols within the pathway are the same as described in Figure 3a.
Genes 12 00405 g003
Figure 4. Th17 Activation Pathway in LS line of Japanese quail. Molecular interactions and symbols within the pathway are the same as described in Figure 3.
Figure 4. Th17 Activation Pathway in LS line of Japanese quail. Molecular interactions and symbols within the pathway are the same as described in Figure 3.
Genes 12 00405 g004
Figure 5. Signaling by Rho Family GTPases and Rac Signaling Pathways in HS line of Japanese quail. Molecular interaction and symbols are the same as described in Figure 3.
Figure 5. Signaling by Rho Family GTPases and Rac Signaling Pathways in HS line of Japanese quail. Molecular interaction and symbols are the same as described in Figure 3.
Genes 12 00405 g005
Figure 6. Role of BRCA in DNA Damage Response within LS line of Japanese quail. Molecular interaction and symbols are the same as described in Figure 3.
Figure 6. Role of BRCA in DNA Damage Response within LS line of Japanese quail. Molecular interaction and symbols are the same as described in Figure 3.
Genes 12 00405 g006
Table 1. Sequencing and mapping data of HS and LS lines of Japanese quail [6].
Table 1. Sequencing and mapping data of HS and LS lines of Japanese quail [6].
LineNumber of SNPsCoverage
HS3,492,46941.45x
LS2,872,43842.59x
Table 2. Top 20 SNPs with the highest read depths in HS line.
Table 2. Top 20 SNPs with the highest read depths in HS line.
LineChrRef Pos Ref BaseCalled BaseImpactSNP %Feature Name DNA ChangeAmino Acid ChangeDepthA CntC CntG CntT CntDeletion
HS215131174ACGTNon-synonymous0.99LOC107323251c.528_529 > ACp.W177R68-06700
HS222387923GANon-synonymous1R3HCC1c.1837C > Tp.P613S56560-00
HS89560343-GCTCAAACACFrameshift1RNPC3c.945_946insp.N316fs50005000
HS534813334ATGCNon-synonymous1TTLL5c.2643_2644 > GCp.S882P50-05000
HS821795021ATNon-synonymous1ZCCHC11c.2277T > Ap.F759L49-00490
HS1214314737ACNon-synonymous0.94CNTN3c.966T > Gp.H322Q48-45000
HS11.5 × 108AGNon-synonymous1LACC1c.1099A > Gp.T367A48-04800
HS391429204TANon-synonymous0.94LOC107312240c.2044A > Tp.N682Y474430-0
HS715772589AGNon-synonymous1LOC107316692c.1156A > Gp.S386G47-04700
HS456633029CCTTNon-synonymous1PCM1c.1849_1850 > AAp.G617N470-0470
HS231086012TCNon-synonymous0.91PLEKHA8c.1304T > Cp.L435S470430-0
HS170128156CGNo-stop1EPHB6c.2447G > Cp..816Sext.?450-4500
HS212172350TCNon-synonymous1LOC107323499c.719A > Gp.Q240R440440-0
HS360212158TCNon-synonymous1REV3Lc.4339T > Cp.S1447P440440-0
HS128367256CGNon-synonymous1XPCc.310G > Cp.V104L440-4400
HS718690075TGNon-synonymous0.97C7H2orf76c.12A > Cp.L4F430042-0
HS119853879GANon-synonymous0.93CCDC79c.820G > Ap.A274T43400-00
HS343357143ACNon-synonymous1KATNA1c.711T > Gp.D237E43-43000
HS268425475AGNon-synonymous1KIF13Ac.5806T > Cp.W1936R43-04300
HS176369984CGNon-synonymous0.95INPP5Ec.744C > Gp.F248L420-4000
Quail line (High Stress), chromosome (Chr) numbers, reference position (Ref Pos), reference base (Ref Base), called (SNP) base, impact (kinds of protein mutation), SNP%, feature name (gene name), DNA change, amino acid change, Depth, and five columns for SNP counts (Cnt) are indicated.
Table 3. Top 20 SNPs with the highest read depths in LS line.
Table 3. Top 20 SNPs with the highest read depths in LS line.
LineChrRef PosRef BaseCalled BaseImpactSNP %Feature NameDNA ChangeAmino Acid ChangeDepthA CntC CntG CntT CntDeletion
LS232838673AGNon-synonymous1PP2D1c.1616T > Cp.V539A54-05400
LS98491817CTNonsense1TRIP12c.4974G > Ap.W1658.520-0520
LS79557668GANon-synonymous1C7H2orf69c.964G > Ap.V322I49490-00
LS42143170TCNon-synonymous0.91ACRCc.242A > Gp.D81G470430-0
LS1110870115TCNon-synonymous1EGFL6c.401A > Gp.K134R470470-0
LS153603348CANon-synonymous0.96PKP2c.1840G > Tp.V614L4745-000
LS164715297GCATNon-synonymous0.93LOC107317569c.401_402 > ATp.R134H45420-00
LS613674979AGNon-synonymous0.98USP54c.3248A > Gp.E1083G44-04300
LS183018090-TTGInframe insertion0.93HEXDCc.1174insCAAp.S392del43000403
LS1120607492ACNon-synonymous0.98LOC107306797c.374A > Cp.K125T43-42000
LS118576282CTNon-synonymous0.93TDRD12c.794C > Tp.S265F430-0400
LS430750477CGNon-synonymous0.95FRAS1c.1649G > Cp.R550T420-4000
LSZ44293804GANon-synonymous1PDZPH1Pc.136G > Ap.A46T42420-00
LS119982878ATNon-synonymous0.98IQUBc.1246T > Ap.S416T41-00400
LS172163526AGNon-synonymous1LOC107319872c.1789A > Gp.I597V40-04000
LS729069125ATNon-synonymous1LRP1Bc.6537T > Ap.D2179E40-00400
LS134497924AGNon-synonymous1RUFY1c.1657A > Gp.T553A40-04000
LS369683334ATNon-synonymous1SMIM8c.29T > Ap.I10N40-00400
LSLGE64150066TCNon-synonymous1LOC107325885c.47T > Cp.V16A390390-0
LS466305765CTNon-synonymous0.97RBPJc.29G > Ap.R10Q390-0380
Quail Line (Low Stress), chromosome (Chr) numbers, reference position (Ref Pos), reference base (Ref Base), called (SNP) base, impact (kinds of protein mutation), SNP%, feature name (gene name), DNA change, amino acid change, Depth, and five columns for SNP counts (Cnt) are indicated.
Table 4. Occurrence of unique SNP impacts in HS and LS lines of Japanese quail.
Table 4. Occurrence of unique SNP impacts in HS and LS lines of Japanese quail.
ImpactNumber of SNPsPercentage
HS Line:
Frameshift460.702%
Inframe deletion, conservative50.076%
Inframe deletion, disruptive50.076%
Inframe insertion, conservative20.031%
Inframe insertion, disruptive40.061%
Nonsense80.122%
Non-synonymous174626.652%
No-start110.168%
No-stop110.061%
Synonymous472072.050%
Total6551
LS Line:
Frameshift330.87%
Inframe deletion, conservative40.11%
Inframe deletion, disruptive30.08%
Inframe insertion, conservative00.00%
Inframe insertion, disruptive40.11%
Nonsense50.13%
Non-synonymous100426.56%
No-start20.05%
No-stop00.00%
Synonymous275872.96%
Total3780
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shumaker, S.; Khatri, B.; Shouse, S.; Seo, D.; Kang, S.; Kuenzel, W.; Kong, B. Identification of SNPs Associated with Stress Response Traits within High Stress and Low Stress Lines of Japanese Quail. Genes 2021, 12, 405. https://doi.org/10.3390/genes12030405

AMA Style

Shumaker S, Khatri B, Shouse S, Seo D, Kang S, Kuenzel W, Kong B. Identification of SNPs Associated with Stress Response Traits within High Stress and Low Stress Lines of Japanese Quail. Genes. 2021; 12(3):405. https://doi.org/10.3390/genes12030405

Chicago/Turabian Style

Shumaker, Steven, Bhuwan Khatri, Stephanie Shouse, Dongwon Seo, Seong Kang, Wayne Kuenzel, and Byungwhi Kong. 2021. "Identification of SNPs Associated with Stress Response Traits within High Stress and Low Stress Lines of Japanese Quail" Genes 12, no. 3: 405. https://doi.org/10.3390/genes12030405

APA Style

Shumaker, S., Khatri, B., Shouse, S., Seo, D., Kang, S., Kuenzel, W., & Kong, B. (2021). Identification of SNPs Associated with Stress Response Traits within High Stress and Low Stress Lines of Japanese Quail. Genes, 12(3), 405. https://doi.org/10.3390/genes12030405

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

Article Metrics

Back to TopTop