Novel Polymorphisms and Genetic Characteristics of the Prion Protein Gene in Pheasants

Transmissible spongiform encephalopathies (TSEs) also known as prion diseases, are fatal neurodegenerative diseases. Prion diseases are caused by abnormal prion protein (PrP^Sc) derived from normal prion protein (PrP^C), which is encoded by the prion protein gene (PRNP). Prion diseases have been reported in several mammals. Notably, chickens, one species of bird, have not been reported to develop prion diseases and showed resistance to bovine spongiform encephalopathy (BSE) infection. However, genetic polymorphisms of the PRNP gene and protein structure of the prion protein (PrP) related to vulnerability to prion diseases have not been investigated in pheasants, another species of bird. We performed amplicon sequencing of the pheasant PRNP gene to identify genetic polymorphisms in 148 pheasants. We analyzed the genotype, allele and haplotype frequencies of the pheasant PRNP polymorphisms. In addition, we evaluated the effect of genetic polymorphisms of the pheasant PRNP gene on pheasant PrP by the AMYCO, PROVEAN, PolyPhen-2 and PANTHER softwares. Furthermore, we compared the amino acid sequences of tandem repeat domains and secondary and tertiary structures of prion proteins (PrPs) among several animals. Finally, we investigated the impact of non-synonymous single nucleotide polymorphisms (SNPs) on hydrogen bonds and tertiary structures of pheasant PrP by Swiss PDB viewer software. We identified 34 novel genetic polymorphisms of the pheasant PRNP gene including 8 non-synonymous SNPs and 6 insertion/deletion polymorphisms. Among the non-synonymous SNPs, the L23F, G33C and R177Q SNPs showed that they could have a deleterious effect on pheasant PrP. In addition, the R177Q SNP was predicted to show an increase in amyloid propensity and a reduction in hydrogen bonds of pheasant PrP. Among the insertion/deletion polymorphisms, c.163_180delAACCCGGGGTATCCCCAC showed that it could have a detrimental effect on pheasant PrP. Furthermore, secondary and tertiary structures of pheasant PrP were predicted to have structures similar to those of chicken PrP. To the best of our knowledge, this is the first study on genetic polymorphisms of the pheasant PRNP gene.

Introduction

Transmissible spongiform encephalopathies (TSEs), also known as prion diseases, are incurable and malignant neurodegenerative disorders. TSEs are caused by misfolded abnormal prion protein (PrP^Sc) derived from normal prion protein (PrP^C) encoded by the prion protein gene (PRNP) (1). PrP^Sc is characterized by a higher amyloid propensity and a higher proportion of β-sheet structure compared to PrP^C. Prion diseases have been classified as various types including Creutzfeldt-Jakob disease (CJD), Gerstmann-Sträussler-Scheinker syndrome (GSS) and fatal familial insomnia (FFI) in humans, bovine spongiform encephalopathy (BSE) in cattle, feline spongiform encephalopathy (FSE) in cats and cheetahs, transmissible mink encephalopathy in minks, scrapie in sheep and goats and chronic wasting disease (CWD) in elk and deer (2–7).

Several single nucleotide polymorphisms (SNPs) of the PRNP gene play a pivotal role in the conversion of PrP^C to PrP^Sc and have been related to susceptibility to prion diseases. In humans, PRNP codons 129 and 219 SNPs are associated with susceptibility to CJD. In particular, M129V and E219K heterozygotes confer resistance to CJD (2–4). In sheep, the susceptibility to scrapie is associated with haplotypes of PRNP polymorphisms at codons 136, 154 and 171. A₁₃₆R₁₅₄R₁₇₁/A₁₃₆R₁₅₄R₁₇₁ diplotype is strongly associated with resistance to scrapie. However, V₁₃₆R₁₅₄Q₁₇₁/V₁₃₆R₁₅₄Q₁₇₁ diplotype is related to susceptibility to scrapie (5, 8). In goats, 9 polymorphisms including codons T110P, G127S, M142I, G145D, N146D,S, R154H, R211Q and Q222K are associated with scrapie (5, 9–11). In addition, the M132L SNP of the PRNP gene in elk plays a pivotal role in the conversion of PrP^Sc in CWD (12–14). Conversely, representative prion disease-resistant animals, including dogs and horses have specific amino acids related to resistance in the conversion of PrP^C to PrP^Sc. In dog PrP, D163 located on the α1-β2 loop provides stability to the PrP structure and prolongs the incubation period in prion-infected transgenic mice carrying this allele (15, 16). In horse PrP, S167 located on the β2-α2 loop confers structural stability of equine PrP (17–19). Thus, the genetic characteristics of the PRNP gene are important in the susceptibility or resistance to prion diseases.

In birds, chickens showed resistance to experimental infection of BSE (20). Amino acid sequences of the PrP were predicted to have low amyloid propensity. In addition, there were no SNPs in the open reading frame (ORF) of the PRNP gene in 4 breeds including Dekalb White, Ross, Ogolgye and Korean native chickens (21, 22). However, Pekin duck showed a higher amyloid propensity of PrP than chickens. Furthermore, the specific amino acids of Pekin duck PrP showed an additional two β-sheet structures compared to chicken PrP and an increase in amyloid propensity (21). Since poultry are closely related to the food industry and, prion disease is infectious, investigation of prion disease-related genetic properties is very important in several avian species.

Pheasant (Phasianus colchicus) is a representative domesticated bird used for eggs and meats within the family Phasianidae. In Korea, more than 320,000 pheasants (Phasianus colchicus karpowi) are raised on over 90 farms. Pheasants and chickens belong to the order Galliformes and showed close evolutionary relationships in the phylogenetic tree. However, ducks also belong to the order Anseriformes and showed a close evolutionary relationship with Galliformes (23). Thus, there is a question of whether the pheasant has chicken-like characteristics or Pekin duck-like characteristics, and it is necessary to investigate the prion disease-related characteristics of the pheasant PRNP gene.

In this study, we performed amplicon sequencing in the pheasant PRNP gene to identify genetic polymorphisms in 148 pheasants and investigated the genotype, allele and haplotype frequencies of the pheasant PRNP polymorphisms. In addition, we analyzed the impact of genetic polymorphisms of the pheasant PRNP gene on pheasant PrP using the AMYCO, PROVEAN, PolyPhen-2 and PANTHER programs. Furthermore, we compared the amino acid sequences of the tandem repeat domains of PrPs among several species. Finally, we investigated the effect of nonsynonymous SNPs on tertiary structures of pheasant PrP and compared secondary and tertiary structures of pheasant PrP with avian PrPs.

Materials and Methods

Ethical Statements

A total of 148 pheasants (Phasianus colchicus) were obtained from a slaughterhouse provided from adjacent pheasant farms in Korea. To minimize internal relationships, we randomly selected samples 3 times. All experimental protocols were approved by the Institutional Animal Care and Use Committee (IACUC) of Jeonbuk National University (JBNU 2020-209). All efforts were made to minimize the number of animals used.

Genomic DNA Extraction

Genomic DNA was extracted from 20 mg cerebral cortex of 148 pheasants using a Bead Genomic DNA Prep kit (Biofact, Daejeon, Korea) following the manufacturer's instructions.

Genetic Analysis

Polymerase chain reaction (PCR) was carried out to amplify the pheasant PRNP gene with gene-specific primers, including PRNP-F (ATAAAGGAGGTGGGGATGGG) and PRNP-R (CGTGGACACGATGTCATCTC). These primers were designed based on the pheasant PRNP gene (GenBank ID; 116238382). The PCR was performed using Taq DNA polymerase kit (Biofact, Daejeon, Korea) contained 0.2 μl of Taq DNA polymerase, 5 μl of 5x Band Helper, 2.5 μl of 10x Taq DNA polymerase buffer, 0.5 μl of 10 mM dNTP mixture, 1 μl (10 pmol) of each primer and 1 μl (50-70 ng/μl) pheasant genomic DNA and sterile deionized water in a total volume of 25 μl. The PCR conditions were followed by the manufacturer's instructions. All PCR products (917 bp) were purified using a FavorPrep gel/PCR Purification Mini Kit (FAVORGEN, Pingtung County, Taiwan). Purified PCR products were directly sequenced in both directions with an ABI 3730xl (ABI, Foster City, CA, USA). Genotyping was performed using Finch TV software (Geospiza Inc., Seattle, WA, USA).

Statistical Analysis

Genotype, allele and haplotype frequencies were calculated by SAS 9.4 software (SAS Institute Inc., Cary, NC, USA). Linkage disequilibrium (LD), Hardy-Weinberg equilibrium (HWE) and haplotype analyses were performed by Haploview version 4.2 (Broad Institute, Cambridge, MA, USA).

In silico Analysis

PROVEAN (http://provean.jcvi.org/index.php), AMYCO (http://bioinf.uab.es/amycov04/index_CompSeq.html), PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/) and PANTHER (http://www.pantherdb.org/tools/csnpScore.do) were used to analyze the effect of pheasant PRNP polymorphisms on pheasant PrP. The PROVEAN score was calculated as the pairwise alignment scores between the query protein sequences and the single locus variations of the other proteins (24). If the PROVEAN score is equal to or below −2.5, it is called a “deleterious” effect, and if it is higher than −2.5, it is called a “neutral” effect. AMYCO using the PAPA score and the pWALTZ score for calculation indicates the protein amyloid propensity (25). An AMYCO score below 0.45 indicates low amyloid propensity, and an AMYCO score above 0.78 indicates high amyloid propensity. PolyPhen-2 was used to predict the possible impact of an amino acid substitution on the structure and function of a protein using direct physical and comparative considerations. The PolyPhen-2 score was calculated from the position-specific independent count (PSIC) score of the wild-type and mutant amino acids (26). The score ranged from 0.0 to 1.0, and the results were divided into “benign,” “probably damaging” and “possibly damaging.” PANTHER measured the preservation time for PSEP (position-specific evolutionary preservation). PSEP was estimated to preserve a position in the current pheasant protein by tracing back to its reconstructed direct ancestors. A longer preservation time indicated that it had a deleterious effect (27). The score for preservation time was “Probably damaging [time>450 millions of years (my)],” “possibly damaging (450>time>200 my)” and “probably benign (time <200 my).”

3D Structure Modeling of Avian PrP

The 3D structure of avian PrPs was predicted by the SWISS-model program (https://swissmodel.expasy.org/) and IntFOLD program (https://www.Reading.ac.uk/bioinf/IntFOLD/). Homology-based modeling was performed by the SWISS-model program based on BLAST and HHblits from the SWISS-model template library (SMTL). The 3D structure of chicken PrP based on nuclear magnetic resonance spectroscopy (NMR) was obtained from a protein data bank (PDB ID:1U3M). The 3D structure of codons 95–230 of Pekin duck PrP was predicted by the SWISS-model program based on the NMR structure of human PrP (PDB ID:2lft.1.A). The 3D structure of codons 134–248 of wild type of pheasant PrP was predicted by the SWISS-model program based on the NMR structure of chicken PrP (PDB ID:1U3M). 3D modeling was performed using three criteria global model quality estimate (GMQE), qualitative model energy analysis (QMEAN) and sequence identity. IntFOLD is an integrated web resource for the prediction of the 3D structure of proteins. The prediction was performed by the ModFOLD6 server with intuitive local and global quality scores of 3D modeling. The 3D structure of codons 1–133 and 248–272 of pheasant PrP was predicted by the IntFOLD program. The predicted 3D structure of pheasant PrP was visualized by Swiss PDB viewer 4.1 software (Swiss Institute of Bioinformatics, Lausanne, Switzerland). The secondary structure and hydrogen bond of pheasant PrP were analyzed by PDB viewer 4.1 software. Hydrogen bonds were predicted if a hydrogen was in the range from 2.195 to 3.3 Å of an accessible donor atom.

Results

Identification of Novel Genetic Polymorphisms of the Pheasant PRNP Gene

The pheasant PRNP gene is comprised of three exons and the ORF is located on exon 3. We amplified the pheasant PRNP gene using gene-specific primers, and the sequencing result was identical to the registered pheasant PRNP gene on GenBank (Gene ID: 116238382). To investigate polymorphisms of the pheasant PRNP gene, PCR and amplicon sequencing were performed in 148 pheasants. We found a total of 34 novel polymorphisms including 28 SNPs and 6 insertion/deletion polymorphisms (Figure 1, Supplementary Figures 1, 2). Of the 28 SNPs, 8 SNPs, c.61G>T (V21F), c.61G>C (V21L), c.67C>T (L23F), c.97G>A (G33C), c.530G>A (R177Q), c.750C>G (I250M), c.766G>A (D256N) and c.781G>A (V261I), were non-synonymous SNPs. We investigated the genotype and allele frequencies and HWE of the 34 polymorphisms of the pheasant PRNP gene (Table 1). Except for c.750C>G and c.766G>A, all genetic polymorphisms were in HWE. The genotype and allele frequencies of the 34 polymorphisms are described in Table 1. In addition, we analyzed LD among all polymorphisms of the pheasant PRNP gene with their D' and r² value. Notably, all genetic polymorphisms have linkage disequilibrium (D' > 0). In addition, a total of 25 strong LDs (r² > 0.333) were found (Table 2). Furthermore, we performed haplotype analysis of all polymorphisms of the pheasant PRNP gene and identified 14 major haplotypes (Table 3). Among the 14 haplotypes, the ACGCGTCWtGGInsWtATWtGCCWtCCGGCGCTWtACGA (13.1%) haplotype had the highest frequency, followed by the GCGCGTCWtGGWtWtATWtGCTWtCCGGTGTCWtGGAG (10.1%) and ACGCGTTWtGGWtWtATInsGCCInsCCGGCGCTWtGCGA (8.4%) haplotypes.

FIGURE 1

Figure 1. Novel genetic polymorphisms of the pheasant prion protein gene (PRNP) found in this study. The open reading frame (ORF) is indicated by a shaded block and the 5' and 3' untranslated regions (UTRs) are indicated by white blocks. The edged horizontal bar indicates the region sequenced. Asterisks denote nonsynonymous single nucleotide polymorphisms (SNPs). Upper lines indicate the novel SNPs. The lower lines indicate the novel insertion/deletion polymorphisms.

TABLE 1

Table 1. Genotype and allele frequencies of pheasant prion protein gene (PRNP) polymorphisms in 148 pheasants.

TABLE 2

Table 2. Linkage disequilibrium (LD) among genetic polymorphisms of the pheasant PRNP gene in pheasants.

TABLE 3

Table 3. Haplotype frequency of 34 PRNP polymorphisms in pheasants.

In silico Analysis of the Effect of Polymorphisms in the Pheasant PRNP Gene

We evaluated the functional and structural effects of PRNP polymorphisms on pheasant PrP using PolyPhen-2, PANTHER and PROVEAN (Table 4). Among the 8 non-synonymous SNPs, 5 non-synonymous SNPs including V21F, V21L, I250M, D256N and V261I were predicted to have benign effects on pheasant PrP in the PANTHER and PROVEAN programs. The L23F SNP was predicted to be “Neutral” and “Probably Damaging” by PROVEAN and PANTHER, respectively. G33C was predicted to be “Deleterious” and “Probably benign” by PROVEAN and PANTHER, respectively. In addition, R177Q was predicted to be “Neutral,” “Probably benign” and “Possibly Damaging” by PROVEAN, PANTHER and PolyPhen-2, respectively. Notably, 4 insertion/deletion polymorphisms, c.180_181insAACCCGGGGTATCCCCAC, c.180_181insAACCCGGGGTATCCCCACAACCCGGGGTATCCCCAC, c.198_199insAACCCAGGATATCCCCAC, c.216_217insAACCCCGGCTATCCCCACAACCCCGGCTATCCCCAC and c.624_626delGAA were predicted to have benign effects on pheasant PrP by PROVEAN. Notably, c.163_180delAACCCGGGGTATCCCCAC was predicted to be “Deleterious” by PROVEAN. We evaluated the amyloid propensity of pheasant PrP according to the alleles of the PRNP polymorphisms by AMYCO (Table 4). All polymorphisms were predicted to show a low amyloid propensity. However, R177Q showed an increase in AMYCO score (0.23) compared to wild type of pheasant PrP (0).

TABLE 4

Table 4. In silico evaluation of the effect of PRNP polymorphisms of pheasant prion protein (PrP).

Next, we analyzed the effect of nonsynonymous SNPs on the 3D structure of pheasant PrP (Figure 2). The V21, L21 and F21 alleles showed two hydrogen bonds with S24 (3.08 Å and 3.25 Å) (Figure 2A). The L23 and F23 alleles showed two hydrogen bonds, with K25 (3.12 Å) and with D20 (3.12 Å) (Figure 2B). The G33 and C33 alleles did not have hydrogen bonds (Figure 2C). The R177 allele showed two hydrogen bonds with Y141 (2.33 Å) and D178 (2.83 Å) (Figure 2D). Notably, the Q177 allele showed only one hydrogen bond with D178 (2.83 Å) (Figure 2D). The I250 and M250 alleles did not have hydrogen bonds (Figure 2E). The D256 allele did not have a hydrogen bond. However, the N256 allele had one hydrogen bond with S248 (3.23 Å) (Figure 2F). The V261 and L261 alleles showed 6 identical hydrogen bonds with N79 (3.07 Å), T257 (2.60 Å), W258 (2.84 Å and 3.19 Å), L259 (3.11 Å) and L265 (3.13 Å) (Figure 2G).

FIGURE 2

Figure 2. Prediction of the tertiary structure and hydrogen bonds of the pheasant prion protein (PrP) according to the alleles of the non-synonymous SNPs of the pheasant PRNP gene. (A) 3D structure of pheasant PrP with the V21, L21, and F21 alleles. (B) 3D structure of the pheasant PrP with the L23 and F23 alleles. (C) 3D structure of the pheasant PrP with the G33 and C33 alleles. (D) 3D structure of the pheasant PrP with the R177 and Q177 alleles. (E) 3D structure of the pheasant PrP with the I250 and M250 alleles. (F) 3D structure of the pheasant PrP with the D256 and N256 alleles. (G) 3D structure of the pheasant PrP with the V261 and L261 alleles. The green dotted line indicates hydrogen bonds. The green numbers indicate the distance of the hydrogen bonds.

Comparison of Tandem Repeat Domains of Pheasants and Several Species

We compared amino acid sequences of tandem repeat domains among human, sheep, goat, cattle, dog, Pekin duck, chicken and pheasant PrPs (Figure 3). Mammalian and avian PrPs showed octapeptide repeat domains (PHGGGWGQ) and hexapeptide repeat domains (NPGYPH), respectively. Pheasant PrP has 9 units of hexapeptide repeat domains, which are identical in length to those of chickens. Except for the second tandem repeat unit (chicken: QPGYPH, pheasant: QPSYPH), the amino acid sequences of the tandem repeat domains of pheasant PrP are identical to those of chicken. However, Pekin duck PrP showed shorter tandem repeat domains than pheasant and chicken PrPs (Pekin duck: 36 aa, pheasant: 54 aa, chicken: 54 aa).

FIGURE 3

Figure 3. Comparison of tandem repeat sequences in humans, sheep, goats, cattle, dogs, Pekin ducks, chickens and pheasants. The amino acid sequences of the tandem repeat region were obtained from GenBank at the National Center for Biotechnology Information (NCBI), including human (BAG32277.1), sheep (NP_001009481.1), goat (NP_001301176.1), cattle (NP_001258555.1), dog (AGA63676.1), Pekin duck (21), chicken (NP_990796.2) and pheasant (In this study). The arrows indicate the insertion/deletion polymorphisms found in this study. U 1-9: hexapeptide repeat 1-9.

Comparison of Secondary and Tertiary Structures of Pheasants and Avian PrPs

We compared the secondary and tertiary structures of avian PrPs by Swiss-model and Swiss PDB viewer. Chicken PrP was predicted to have two β-sheet structures (codons 136–138, 168–170) and three α-helix structures (codons 157–167, 185–200 and 219–246) (Figure 4A). Pekin duck PrP was predicted to have four β-sheet structures (codons 125–127, 157–159, 195–197 and 200–202) and five α-helix structures (140–147, 151, 168, 172–182 and 206–228) (Figure 4B). Pheasant PrP was predicted to have two β-sheet structures (codons 142–144 and 174–176) and three α-helix structures (codons 157–167, 185–200 and 219–246) (Figure 4C). The tertiary structures of chicken and pheasant PrPs showed three α-helices and two antiparallel β-sheet structures (Figures 4D,F) However, the tertiary structure of Pekin duck PrP showed five α-helices and four antiparallel β-sheet structures (β1-β2 and β3-β4) (Figure 4E).

FIGURE 4

Figure 4. The secondary and tertiary structures of avian PrPs. (A) The secondary structure of the chicken PrP. (B) The secondary structure of the Pekin duck PrP. (C) The secondary structure of the pheasant PrP. (D) The tertiary structure of the chicken PrP. (E) The tertiary structure of the Pekin duck PrP. (F) The tertiary structure of the pheasant PrP. The colors indicate α-helices (red), β-sheets (yellow) and coils (white). The secondary and tertiary structures of avian PrPs were analyzed by SWISS-MODEL and Swiss PDB Viewer.

Discussion

In this study, we found a total of 34 novel polymorphisms of the pheasant PRNP gene including 8 non-synonymous SNPs and 6 insertion/deletion polymorphisms (Figure 1, Table 1). Except for c.750C>G (I250M) and c.766G>A (D256N), all genetic polymorphisms were in HWE (Table 1). Since c.750C>G (I250M) and c.766G>A (D256N) are non-synonymous SNPs, functional alteration caused by the non-synonymous SNPs may be related to genetic selection and induced the HWE violation. Among the 8 non-synonymous SNPs, 2 SNPs including G33C and R177Q were predicted to have deleterious effects on pheasant PrP by in silico analyses (Table 4). In addition, L23F was predicted to have damaging effect on pheasant PrP by PANTHER. Since the L23F and G33C SNPs, which are located in the interspecies-conserved signal peptide of PrP, were predicted to have a deleterious effect on pheasant PrP, further analysis of the signal peptide of pheasant PrP according to the allele of the pheasant PRNP gene is highly desirable. The R177Q SNP, which is located in the structured region of pheasant PrP, confers an increase in amyloid propensity (Table 4). Although R177Q has a higher AMYCO score than other polymorphisms, it does not exceed the threshold with a low amyloid propensity. It indicates that genetic polymorphisms of the pheasant PRNP gene may be induced low amyloid propensity and pheasant PrP has a relatively stable protein structure. In addition, pheasant PrP with Q177 showed a reduction in the number of hydrogen bonds compared to pheasant PrP with R177. Since hydrogen bonds play a pivotal role in the stability of PrP and the stability of PrP confers resistance to the conformational changes in PrP (28), the R177Q SNP was thought to be a potential risk factor for prion diseases. Further investigation of the in vitro and/or in vivo evaluation of the amyloid propensity of pheasant PrP with the R177 allele and Q177 allele is needed in the future. However, although several genetic polymorphisms have strong genetic linkages, the effect of each SNP and insertion/deletion were only evaluated individually due to the internal setting of the software. Further evaluation of sequence variations of the pheasant PRNP gene according to haplotypes is highly desirable in the future.

We identified a highly polymorphic tandem repeat region of the pheasant PRNP gene (Figures 1, 3). The tandem repeat region of PrP is related to the activity of copper-related enzymes via the binding of copper ions in mammals (histidine) and birds (histidine, tyrosine) (29, 30). Since pheasant has additional copper binding sites, the insertion/deletion polymorphisms of the pheasant PRNP gene would not cause critical dysfunction of the copper-related activity of pheasant PrP and genetic polymorphisms of the pheasant tandem repeat region seem to be more frequently observed than in mammals. Previous studies have reported that insertion/deletion polymorphisms in the octapeptide repeat region of the PRNP gene are associated with vulnerability to CJD (31). Since insertion/deletion polymorphisms are a potent risk factor for prion disease, further analysis on relationship between insertion/deletion polymorphisms of the pheasant PRNP gene and susceptibility to prion disease is needed in the future.

We also found that the secondary and tertiary structures of the pheasant PrP have two β-sheets and three α-helices, similar to chicken PrP. Conversely, the pheasant PrP showed different secondary and tertiary structures from those of Pekin duck, which has 4 β-sheets and 5 α-helices. Since chickens showed resistance to BSE infection and pheasant PrP has a PrP structure similar to that of chicken PrP, pheasant PrP is predicted to have a relatively prion-resistant structure. Further validation is needed in the future using transgenic mice carrying pheasant PrP.

Conclusion

In the present study, we identified 34 novel genetic polymorphisms of the pheasant PRNP gene including 8 non-synonymous SNPs and 6 insertion/deletion polymorphisms in 148 pheasants. Among the 8 non-synonymous SNPs, the L23F, G33C and R177Q SNPs were predicted to have a deleterious effect on pheasant PrP. In addition, the R177Q SNP induced an increase in amyloid propensity and a reduction in hydrogen bonds. Among the 6 insertion/deletion polymorphisms, the c.163_180delAACCCGGGGTATCCCCAC polymorphism was predicted to have a deleterious effect on pheasant PrP. Furthermore, the secondary and tertiary structures of the pheasant PrP are very similar to those of chicken PrP. To the best of our knowledge, this is the first report on genetic polymorphisms of the pheasant PRNP gene.

Data Availability Statement

Data are available on reasonable request. Requests may be made to bhjeong@jbnu.ac.kr.

Ethics Statement

The animal study was reviewed and approved by Institutional Animal Care and Use Committee (IACUC) of Jeonbuk National University.

Author Contributions

KK, Y-CK, and B-HJ conceived, designed the experiment, analyzed the data, and wrote the paper. KK and Y-CK performed the experiments. All authors read and approved the final manuscript.

Funding

KK was supported by the BK21 Plus Program in the Department of Bioactive Material Sciences. This research was supported by the Basic Science Research Program through the National Research Foundation (NRF) of Korea funded by the Ministry of Education (2017R1A6A1A03015876, 2021R1A6A3A010864). This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (2021R1A2C1013213, 2022R1C1C2004792) (NRF-2019-Fostering Core Leaders of the Future Basic Science Program/Global Ph.D. Fellowship Program).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.935476/full#supplementary-material

Supplementary Figure 1. Electropherograms of 28 novel single nucleotide polymorphisms (SNPs) of the PRNP gene in pheasants. Four colors indicate the individual bases of the DNA sequence using an ABI3730xl automatic sequencer (blue: cytosine, red: thymine, black: guanine, green: adenine). The upper panels indicate major homozygotes, the middle panels indicate heterozygotes, and the lower panels indicate minor homozygotes.

Supplementary Figure 2. Electropherograms of the insertion/deletion polymorphisms of the pheasant PRNP gene. (A) c.163_180delAACCCGGGGTATCCCCAC (p.55_60delNPGYPH) (B) c.180_181insAACCCGGGGTATCCCCAC (p.60_61insNPGYPH) (C) c.180_181insAACCCGGGGTATCCCCACAACCCGGGGTATCCCCAC (p.60_61insNPGYPHNPGYPH) (D) c.198_199insAACCCAGGATATCCCCAC (p.66_67insNPGYPH) (E) c.216_217insAACCCCGGCTATCCCCACAACCCCGGCTATCCCCAC (p.72_73insNPGYPHNPGYPH) (F) c.624_626delGAA (p.208delK). Four colors indicate individual bases of the DNA sequence (blue: cytosine, red: thymine, black: guanine, green: adenine). WT: wild type of the pheasant prion protein (PRNP) gene.

References

1. Prusiner SB. Prions. Sci Am. (1984) 251:50–9 doi: 10.1038/scientificamerican1084-50

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Jeong BH, Kim YS. Genetic studies in human prion diseases. J Korean Med Sci. (2014) 29:623–32. doi: 10.3346/jkms.2014.29.5.623

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Kim YC, Park KJ, Hwang JY, Park HC, Kang HE, Sohn HJ, et al. In-depth examination of PrP(Sc) in Holstein cattle carrying the E211K somatic mutation of the bovine prion protein gene (PRNP). Transbound Emerg Dis. (2021). doi: 10.1111/tbed.14309. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Kim YC, Won SY, Jeong MJ, Jeong BH. Absence of proteinase K-resistant PrP in Korean Holstein cattle carrying potential bovine spongiform encephalopathy-related E211K somatic mutation. Transbound Emerg Dis. (2021) 69:805–12. doi: 10.1111/tbed.14053

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Baylis M, Goldmann W. The genetics of scrapie in sheep and goats. Curr Mol Med. (2004) 4:385–96. doi: 10.2174/1566524043360672

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Murdoch BM, Murdoch GK. Genetics of prion disease in cattle. Bioinform Biol Insights. (2015) 9:1–10. doi: 10.4137/BBI.S29678

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Williams ES, Young S. Chronic wasting disease of captive mule deer: a spongiform encephalopathy. J Wildl Dis. (1980) 16:89–98. doi: 10.7589/0090-3558-16.1.89

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Bossers A, de Vries R, Smits MA. Susceptibility of sheep for scrapie as assessed by in vitro conversion of nine naturally occurring variants of PrP. J Virol. (2000) 74:1407–14. doi: 10.1128/JVI.74.3.1407-1414.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Hazards E, Ricci A, Allende A, Bolton D, Chemaly M, Davies R, et al. Genetic resistance to transmissible spongiform encephalopathies (TSE) in goats. EFSA J. (2017) 15:e04962. doi: 10.2903/j.efsa.2017.4962

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Kim SK, Kim YC, Won SY, Jeong BH. Potential scrapie-associated polymorphisms of the prion protein gene (PRNP) in Korean native black goats. Sci Rep. (2019) 9:15293. doi: 10.1038/s41598-019-51621-y

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kim YC, Kim SK, Jeong BH. Scrapie susceptibility-associated indel polymorphism of shadow of prion protein gene (SPRN) in Korean native black goats. Sci Rep. (2019) 9:15261. doi: 10.1038/s41598-019-51625-8

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Green KM, Browning SR, Seward TS, Jewell JE, Ross DL, Green MA, et al. The elk PRNP codon 132 polymorphism controls cervid and scrapie prion propagation. J Gen Virol. (2008) 89:598–608. doi: 10.1099/vir.0.83168-0

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Moore SJ, Vrentas CE, Hwang S, West Greenlee MH, Nicholson EM, Greenlee JJ. Pathologic and biochemical characterization of PrP(Sc) from elk with PRNP polymorphisms at codon 132 after experimental infection with the chronic wasting disease agent. BMC Vet Res. (2018)14:80. doi: 10.1186/s12917-018-1400-9

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Roh IS, Kim YC, Kim HJ, Won SY, Jeong MJ, Hwang JY, et al. Polymorphisms of the prion-related protein gene are strongly associated with cervids' susceptibility to chronic wasting disease. Vet Rec. (2021) 190:e940. doi: 10.1002/vetr.940

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Fernandez-Borges N, Parra B, Vidal E, Erana H, Sanchez-Martin MA, de Castro J, et al. Unraveling the key to the resistance of canids to prion diseases. PLoS Pathog. (2017) 13:e1006716. doi: 10.1371/journal.ppat.1006716

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kim D-J, Kim Y-C, Kim A-D, Jeong B-H. Novel polymorphisms and genetic characteristics of the prion protein gene (PRNP) in dogs—a resistant animal of prion disease. Int J Mol Sci. (2020) 21:4160. doi: 10.3390/ijms21114160

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Perez DR, Damberger FF, Wuthrich K. Horse prion protein NMR structure and comparisons with related variants of the mouse prion protein. J Mol Biol. (2010) 400:121–8. doi: 10.1016/j.jmb.2010.04.066

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Kim Y-C, Jeong B-H. The first report of polymorphisms and genetic characteristics of the prion protein gene (PRNP) in horses. Prion. (2018) 12:245–52 doi: 10.1080/19336896.2018.1513316

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Kim Y-C, Won S-Y, Do K, Jeong B-H. Identification of the novel polymorphisms and potential genetic features of the prion protein gene (PRNP) in horses, a prion disease-resistant animal. Sci Rep. (2020) 10:8926. doi: 10.1038/s41598-020-65731-5

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Moore J, Hawkins SA, Austin AR, Konold T, Green RB, Blamire IW, et al. Studies of the transmissibility of the agent of bovine spongiform encephalopathy to the domestic chicken. BMC Res Notes. (2011) 4:501. doi: 10.1186/1756-0500-4-501

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Jeong MJ, Kim YC, Jeong BH. The first report of the prion protein gene (PRNP) sequence in pekin ducks (Anas platyrhynchos domestica): the potential prion disease susceptibility in ducks. Genes. (2021) 12:193. doi: 10.3390/genes12020193

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Kim YC, Won SY, Jeong BH. Absence of single nucleotide polymorphisms (SNPs) in the open reading frame (ORF) of the prion protein gene (PRNP) in a large sampling of various chicken breeds. BMC Genomics. (2019) 20:922. doi: 10.1186/s12864-019-6315-8

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Prum RO, Berv JS, Dornburg A, Field DJ, Townsend JP, Lemmon EM, et al. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature. (2015) 526:569–73. doi: 10.1038/nature15697

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Choi Y, Chan AP. PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. (2015) 31:2745–7. doi: 10.1093/bioinformatics/btv195

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Iglesias V, Conchillo-Sole O, Batlle C, Ventura S. AMYCO evaluation of mutational impact on prion-like proteins aggregation propensity. BMC Bioinformatics. (2019) 20:24. doi: 10.1186/s12859-019-2601-3

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. (2013) 7:20. doi: 10.1002/0471142905.hg0720s76

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Tang H, Thomas PD. PANTHER-PSEP predicting disease-causing genetic variants using position-specific evolutionary preservation. Bioinformatics. (2016) 32:2230–2. doi: 10.1093/bioinformatics/btw222

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Myers JK, Pace CN. Hydrogen bonding stabilizes globular proteins. Biophys J. (1996) 71:2033–9. doi: 10.1016/S0006-3495(96)79401-8

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Miura T, Hori-i A, Mototani H, Takeuchi H. Raman spectroscopic study on the copper(II) binding mode of prion octapeptide and its pH dependence. Biochemistry. (1999) 38:11560–9. doi: 10.1021/bi9909389

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Stanczak P, Luczkowski M, Juszczyk P, Grzonka Z, Kozlowski H. Interactions of Cu2+ ions with chicken prion tandem repeats. Dalton Trans. (2004) 2102–7. doi: 10.1039/B405753H

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Moore RA, Herzog C, Errett J, Kocisko DA, Arnold KM, Hayes SF, et al. Octapeptide repeat insertions increase the rate of protease-resistant prion protein formation. Protein Sci. (2006) 15:609–19. doi: 10.1110/ps.051822606

PubMed Abstract | CrossRef Full Text | Google Scholar