Journal of Animal Science and Technology
Korean Society of Animal Sciences and Technology
RESEARCH ARTICLE

Genetic association of polymorphisms in porcine RGS16 with porcine circovirus viral load in naturally infected Yorkshire pigs

Seung-Hoon Lee1,2https://orcid.org/0000-0001-6703-7914, Kyu-Sang Lim3https://orcid.org/0000-0001-5406-266X, Ki-Chang Hong1,*https://orcid.org/0000-0001-8605-5050, Jun-Mo Kim2,*https://orcid.org/0000-0002-6934-398X
*Corresponding author: Ki-Chang Hong, Division of Biotechnology, College of Life Science, Korea University, Seoul 02841, Korea., Tel: +82-2-3290-3053, E-mail: kchong@korea.ac.kr
*Corresponding author: Jun-Mo Kim, Department of Animal Science and Technology, Chung-Ang University, Anseong 17546, Korea., Tel: +82-31-670-3263, E-mail: junmokim@cau.ac.kr

© Copyright 2021 Korean Society of Animal Science and Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Aug 19, 2021; Revised: Sep 14, 2021; Accepted: Sep 16, 2021

Published Online: Nov 30, 2021

Abstract

Regulator of G protein signaling 16 (RGS16) is known to be associated with porcine circovirus type 2 (PCV2). PCV2 associated disease (PCVAD) is a serious problem in the swine industry. The representative symptoms of PCVAD are high viral titer proliferation and decreased average daily gain. In this study, we identified single nucleotide polymorphisms (SNPs) in the RGS16 region, including the upstream region. Of the 22 identified SNPs, rs332913874, rs326071195, and rs318298586 were genotyped in 142 Yorkshire pigs. These SNPs were significantly associated with the PCV2 viral load. Moreover, the haplotype combination was also related to the PCV2 viral load. The haplotype and diplotype analysis also had a significant difference with the PCV2 viral load. Taken together, our results suggest that RGS16 SNPs considerably affect the PCV2 viral load.

Keywords: PCV2; Viral load; RGS16; Single nucleotide polymorphism; Pig

INTRODUCTION

Porcine circovirus type 2 (PCV2), a small non-enveloped DNA virus, is perceived as an important pathogen in the swine industry worldwide [1]. PCV2 has four major open reading frames (ORFs). ORF1 encodes the replicase (Rep) and splicing replicase (Rep’) proteins that participate in viral replication [2,3]. ORF2 encodes a viral capsid protein (Cap) [4,5]. The ORF3-encoded protein has been reported to be involved in the apoptosis of PCV2 infected cells [6]. The fourth ORF, named ORF4, has also been reported as a mitochondrial-targeting protein that interacts with adenine nucleotide translocase 3 via the mitochondrial pathways [7].

The proliferation of virus titers correlates with weight loss in post-weaning piglets [8,9]. The typical symptoms of PCV2 associated disease (PCVAD) are an increase in viral load and a decrease in average daily gain (ADG) [10-12]. Moreover, the appearance of PCVAD in swine herds leads to low productivity and economic losses for those farms [13-15]. In addition, PCVAD is known to be influenced by the host genetic background as to which genetically resistant or susceptible to PCV2 [16]. Therefore, the control of PCV2 viral load through host genetic effects is considered important.

Regulator of G protein signaling 16 (RGS16), a GTPase-activating protein, can regulate T-lymphocyte activation [17]. It has previously been reported that RGS16 is upregulated in inflammatory dendritic cells stimulated with lipopolysaccharide [18]. Furthermore, RGS16 is associated with the ORF3 protein of PCV2 and host response [19]. In pigs, the RGS16 gene is located in the quantitative trait loci for ADG on Sus scrofa chromosome 9 [20,21]. Accordingly, RGS16 is thought to be associated with PCVAD. The aim of this study was to identify polymorphisms in the RGS16 region and investigate the association between genotypes and the PCV2 viral load. Four-week-old unvaccinated piglets against PCV2 were used for measuring PCV2 viral load by blood collection. After 6 weeks, the viral load of PCV2 was measured again via re-blood collection.

MATERIALS AND METHODS

This study was conducted in accordance with the animal care guidelines of the Institutional Animal Care and Use Committee of the National Institute of Animal Science, Korea.

Animals

A total of 142 unvaccinated Yorkshire pigs (65 castrated males and 77 females) which were born from the 46 different sows selected from a single farm were chosen to have zero viral load at 4 weeks of age and were naturally infected with PCV2. Peripheral blood was collected from each pig at 4 and 10 weeks of age to assess PCV2 viral load and genomic DNA.

Measuring porcine circovirus type 2 viral load in serum

PCV2 viral DNA was extracted from serum using the QIAGENTM virus mini elute kit (QIAGEN, MD, USA). PCV2 viral load was measured using TaqManTM probe quantitative polymerase chain reaction (qPCR) with the primers listed in Table 1. These primers were designed from the complete PCV2 genome (Korea, GenBank: FR823451.1) using Primer3 [22]. The reaction mixture contained 500 nM of each primer, 250 nM probe, 25 μL of 2× TaqManTM Universal Master Mix (Applied Biosystems, Waltham, MA, USA), and 2.5 μL template. Nuclease-free water was added to bring the final volume to 50 μL. Amplification was performed under universal cycling conditions (2 min at 50°C, 10 min at 95°C, and 45 cycles of 15 s at 95°C and 1 min at 58°C).

Table 1. Primers and probe used for analyzing PCV2 viral load in serum
Oligo Length (nt) Primer sequence Location in PCV2 genome1)
Forward 22 5′-TCG ATC TCA AGG ACA ACG GAG T-3′ 380–401
Reverse 24 5′-TTG GTC TTC CAA TCA CGC TTC TGC-3′ 530–553
Probe 27 5′-CAG AGC AGC ACC CTG TAA CGT TTG TCA-3′ 457–483

GenBank: FR823451.1.

PCV2, porcine circovirus type 2.

Download Excel Table
Identification of genomic sequence variations

To identify polymorphisms, between or within four different pig breeds (Yorkshire, Landrace, Duroc, and Berkshire; n = 20 per breed), direct sequencing analysis was performed using porcine-specific primers containing the upstream and genetic regions of porcine RGS16. The primer sets were designed using the Primer3 software [22] (Table 2). PCR was performed using a MastercyclerTM gradient (Eppendorf, Hamburg, Germany) in a total volume of 20 μL containing 100 ng porcine genomic DNA template, 10 pmol of each primer, 0.25 mM of each dNTP, 1× PCR buffer, and 1.25 U DNA polymerase (i-Max II, Intron biotechnology, Seongnam, Korea). Sequencing analysis was performed using an ABI PRISM 3730 Genetic Analyzer (Applied BiosystemsUSA). Sequence assembly and polymorphism identification were executed using SeqMan [23].

Table 2. Primers for polymorphism identification in porcine RGS16 gene
Primer number Primer sequence Binding region Size (bp)
P1 5′-CTT GGC AGC TGC ATA ATG A-3′ Upstream 760
P2 5′-GGA GCG TGT TTT TGT GAC TC-3′ Upstream 760
P3 5′-GGT GGA ACC TCT TCT CTG CT-3′ Upstream 828
P4 5′-TCG GAA CTA AGA GGG TGT GA-3′ Upstream 828
P5 5′-TCC AGT GAA CAG CTG ATG TG-3′ Upstream 791
P6 5′-CGT GCC CAA GAT CAG AGT TA-3′ Upstream 791
P7 5′-GTT TTA GCG AAA GGG ATT GG-3′ Exon 1 881
P8 5′-CAG CTG GAC AAG CTC CTC TT-3′ Exon 1 881
P9 5′-GAC AAC GTG CGA CTT CTC TG-3′ Exon 1 870
P10 5′-AAA AAC ACT TCC CTG CCT TC-3′ Exon 1 870
P11 5′-CTG CCT GGA GAA AGA AGT CC-3′ Intron1 704
P12 5′-AGC TGT GCC TCC ATA GGT CT-3′ Intron1 704
P13 5′-CAA TCT TTG CTG AAG CTG GA-3′ Exon 2 710
P14 5′-GCC CTG GGT AGA TTT GTC TG-3′ Exon 2 710
P15 5′-TTA GGA CAT GTC AGG GTG GA-3′ Exon 3 749
P16 5′-ACA ACA GGG TCC CAC TCA CT-3′ Exon 3 749
P17 5′-CAT TTG AGG AGT CAG CTC CA-3′ Exon 4 824
P18 5′-TTT TCA GTC TCC AGC CTC CT-3′ Exon 4 824
Download Excel Table
Genotyping

Genotyping was performed on the rs332913874, rs326071195, and rs318298586 loci. Genomic DNA was extracted from peripheral blood using a DNA isolation kit (G-DexTM IIb, Intron Biotechnology). These single nucleotide polymorphism (SNP) genotypes were identified via direct sequencing analysis. PCR was performed using the aforementioned conditions. The amplification conditions were as follows: 94°C for 10 min; 30 cycles of 94°C for 1 min, 61°C for 1 min, and 72°C for 1 min; and a final extension at 72°C for 10 min.

Statistical analysis

The significance of differences in genotype and allele frequencies within the population was determined using the Hardy–Weinberg equilibrium (HWE) calculation (chi-square test). The association of SNPs with viral load was evaluated using the MIXED model procedure in the SAS software (ver. 9.4, SAS Institute, Cary, NC, USA). The model was used as follows: yijkl = μ + Si + Bj + Dk + eijkl, where yijklmis the observation of traits, μ is the general mean, Si is the fixed effect of SNP genotypes i, Bj is the fixed effect of batch j, Dk is the random effect of dam l, and eijkl is the random error. Haplotype combinations of three SNPs were reconstructed by PHASE software [24] and analyzed with the same MIXED procedure using SAS. Results are presented as the mean of the least squares with standard error. Additive and dominance effects were calculated using the method previously described by Hasenstein et al. [25]. A value of p < 0.05 was considered statistically significant.

RESULTS

Identification of SNPs in RGS16 in pig populations

We identified 22 SNPs in the upstream and genetic regions of RGS16 inYorkshire pigs (Fig. 1). Of these, 10 SNPs, including rs332913874 and rs326071195, were located in the upstream region, whereas 12 others, including rs318298586, spread over the intron region. These 22 SNPs were registered on dbSNP and are presented in Supplementary Table S1.

jast-63-6-1223-g1
Fig. 1. Diagrammatic representation of the RGS16 gene and identified polymorphisms. Arrows indicate polymorphisms from upstream and genomic regions of RGS16; textrepresents SNP IDs thatwere identified by us and used in this study. SNP, single nucleotide polymorphism.
Download Original Figure

The allelic and genotypic frequencies of rs332913874, rs326071195, and rs318298586 areshown in Table 3. In rs332913874, the G and A allele frequencies were 0.64 and 0.36, respectively, in the Yorkshire population. Genotype frequencies for the homozygote (GG and AA) and heterozygote (GA) genotypes were 0.38, 0.10, and 0.52, respectively. The genotype and allele frequencies of rs326071195 were found to be as follows: T allele, 0.52; C allele, 0.48, TT genotype, 0.19; TC genotype, 0.67; and CC genotype, 0.14. Furthermore, the genotype and allele frequencies of another SNP (rs318298586) were as follows: A allele, 0.51; T allele, 0.49; AA genotype, 0.19; AT genotype, 0.65; and TT genotype, 0.16. Chi-square tests were used to check for deviation from the HWE. We found that the rs332913874 genotypes conformed to HWE; however, the genotypes of the other two SNPs did not follow HWE in our population. Based on these frequency data, we conducted an association analysis between genotypes and the PCV2 viral load in our population.

Table 3. Allele and genotype frequencies of polymorphisms in porcine RGS16 gene
SNP ID N Genotype Allele HWE1)p-value
rs332913874 142 GG AG AA G A 0.12
0.38 (54)2) 0.52 (74) 0.10 (14) 0.64 0.36
rs326071195 142 TT TC CC T C < 0.05
0.19 (27) 0.67 (95) 0.14 (20) 0.52 0.48
rs318298586 142 AA AT TT A T < 0.05
0.19 (27) 0.65 (92) 0.16 (23) 0.51 0.49

Hardy–Weinberg equilibrium (HWE) by the X2-test (df = 1, X20.05 = 3.841).

Number of animals.

SNP, single nucleotide polymorphism.

Download Excel Table
Haplotype analysis

Haplotype analysis of the three SNPs is presented in Table 4. We observed three major haplotypes that accounted for 97.5% of all possible haplotype combinations in the population of animals. Based on these three haplotype combinations (Haplo1, Haplo4, and Haplo6), we performed association analysis with haplotype combinations and the PCV2 viral load. Results of the analyses are presented in Table 5. The Haplo1 (GTA) and Haplo6 (ACT) combinations showed significant differences in PCV2 viral loads. Animals containing one copy of the Haplo1 combination had a significantly higher PCV2 viral load than those containing two copies (p = 0.04); However, compared to that of a single copy of Haplo1, the absence of Haplo1 combination did not significantly affect PCV2 viral load. The Haplo4 (GCT) combination did not significance. The diplotype frequencies of haplotpype combinations were showed in Supplementary Table S2. The frequency of the Haplo1/Haplo6 (GTA/ACT) combination which had the most count was 0.465. The second and the third count of diplotype combinations were Haplo1/Haplo1 (GTA/GTA) and Haplo1/Haplo4 (GTA/GCT). The frequencies of these two combination were 0.176 and 0.169 each. The forth frequency of diplotype combinations (Haplo6/Haplo6, ACT/ACT) was 0.099. The associations among the top 4 diplotype combinations against viral load were presented in Supplementary Table S3. The viral load of Haplo1/Haplo6 (GTA/ACT) combination was significantly higher than Haplo1/Haplo1 (GTA/GTA) combination (p = 0.01).

Table 4. Haplotype combination from the three SNPs in the upstream region and intron 1 of the RGS16 gene
Haplotype Allele combination Count Estimated frequency Cumulated frequency SE
Haplo1 GTA 144 0.507 0.507 0.001
Haplo2 GTT 5 0.018 0.525 0.001
Haplo3 GCA 1 0.004 0.529 0.001
Haplo4 GCT 32 0.113 0.642 0.001
Haplo5 ACA 1 0.004 0.646 0.001
Haplo6 ACT 101 0.355 1.001 0.001

SNP, single nucleotide polymorphism.

Download Excel Table
Table 5. Significant haplotype effect of polymorphisms in RGS16 on the PCV2 viral load in 142 Yorkshire pigs
Combination Significant trait Haplotype
Absent 1 copy 2 copies p-value
No. LSM ± SE No. LSM ± SE No. LSM ± SE
Haplo1 (GTA) Viral load 23 2.99 ± 0.39ab 94 3.43 ± 0.19a 25 2.36 ± 0.39b 0.04
Haplo4 (GCT) Viral load 110 3.20 ± 0.20 32 3.34 ± 0.35 0 N/A 0.71
Haplo6 (ACT) Viral load 55 2.75 ± 0.26b 73 3.64 ± 0.22a 14 2.57 ± 0.49b 0.01

Means with different superscript letters are significantly different (p < 0.05).

PCV2, porcine circovirus type 2; LSM, least squares mean.

Download Excel Table
Associations between porcine circovirus type 2 viral load and RGS16 genotypes

The animal frequency against viral load at 10-week old is showed using a histogram in Supplementary Fig. S1. The interval of the viral load was set by positive integer values. The results of the analyzed associations between the genotypes of SNPs in RGS16 and PCV2 viral loads are presented in Table 6. The rs332913874 SNP was significantly associated with PCV2 viral load (p = 0.03). Animals with the heterozygous genotype had a higher viral load than those with the homozygous GG genotype. However, the heterozygous GA genotype did not show any statistical significance with another homozygous AA genotype. The association analysis of rs326071195 and PCV2 viral load revealed that the viral load in animals with the TC genotype was higher than that in animals with the TT genotype; however, no significance was observed between the viral loads in animals with the CC and heterozygous genotypes. The results of residual normality test are showed in Supplementary Fig. S2. All of the residual was satisfied normality. In the case of rs318298586, the viral load of pigs with the AT genotype was higher than that of pigs with the AA genotype. Conversely, no significant difference against viral load was found between the heterozygous and TT genotypes.

Table 6. Significant genotype effect of polymorphisms in RGS16 on the PCV2 viral load in 142 Yorkshire pigs
Locus(Allele1/Allele2) Significant trait Genotype Additive Dominance
Allele1/Allele1 Allele1/Allele2 Allele2/Allele2 p-value
No. LSM ± SE No. LSM ± SE No. LSM ± SE
rs332913874(G/A) Viral load1) 54 2.80 ± 0.27b 74 3.59 ± 0.22a 14 2.58 ± 0.49ab 0.03 −0.11 1.81*
rs326071195 (T/C) Viral load 27 2.33 ± 0.37b 95 3.45 ± 0.19a 20 2.96 ± 0.42ab 0.02 0.31 1.61*
rs318298586 (A/T) Viral load 27 2.31 ± 0.37b 92 3.47 ± 0.19a 23 2.99 ± 0.39ab 0.02 −0.34 1.63*

The unit of viral load: Log PCV2 copies/mL.

Means with different superscript letters are significantly different (p < 0.05).

Significant levels of additive and dominance effects:

p < 0.05.

PCV2, porcine circovirus type 2; LSM, least squares mean.

Download Excel Table

DISCUSSION

In this study, we investigated the association between PCV2 viral load and SNP genotypes in porcine RGS16. PCV2, a small non-enveloped, circular ssDNA virus, leads to PCVAD with clinical signs of wasting disease, enlarged lymph nodes, and diarrhea [9,11,16,26]. PCVAD is an endemic disease that causes reduced productivity in the swine industry [1,10,27]. PCV2 is commonly known to have 11 ORFs [16,28]; however, the functions of only four of these ORF proteins have been identified. ORF1 codes for nonstructural replicase proteins known as Rep and Rep’ [5,29]. ORF2 encodes the Cap protein [30,31]. ORF3 proteins have been previously reported to induce apoptosis in PK-15 cells [32]. ORF4 is known to be a mitochondrial-targeting protein that triggers apoptosis in the host cell [7].

Various types of research on PCV2 have revealed that the viral load of experimental or naturally infected pigs is associated with PCVAD clinical signs, including postweaning multisystemic wasting syndrome (PMWS) [33-36]. To prevent PCVAD, vaccinations against PCV2 are done in the swine farms. Nevertheless, the swine industry has been consistently damaged by PCV2. Therefore, researchers have experimented on host genetic susceptibility to PCV2 using infection and genome-wide association studies [37,38]. Walker et al. reported that porcine SYNGR2 polymorphism is related to the PCV2 viral load. Zhang et al. also reported that some miRNAs are associated with PCVAD susceptibility to PCV2 infection in lung tissues. In this study, we investigated the differences in viral load based on genotype in the RGS16 polymorphism in naturally infected Yorkshire pigs. We identified SNPs in porcine RGS16 gene. A total of 22 SNPs were identified, including those in the 5′UTR and 3′UTR regions. In the statistical analysis model, we did not include sex as a fixed effect. There was no significant effect between sex in the statistical model including sex as a fixed effect. The previous study reported that the difference of sex effect against PCV2 was not founded [39]. Therefore, it is estimated that there was no difference in infection with gender in our study. Three SNPs located in the intron region were associated with the PCV2 viral load. In particular, the heterozygous genotypes of these three SNPs were related to high PCV2 viral loads. Furthermore, the reference homozygous genotypes (Allele1/Allele1) were associated with a low viral load. In addition, the association between SNP and viral load did not follow the additive effects and the complete dominance effect. This phenomenon is presumed to be the effect of hybrid weakness, but more follow-up studies seem to be needed. As previously described, PCV2 viral load is a triggering factor that causes PCVAD, including PMWS. Moreover, RGS16 has been previously reported to be involved in the activation of immune cells, with the important exception of macrophages in rats [40,41]. Another study has reported that RGS16 activated by inflammatory stimulation and the PCV2 ORF3 protein colocalize in porcine peripheral blood mononuclear cells [19]. In the present study, we did not identify differences in the expression levels of RGS16 in porcine whole blood among the genotypes. In addition, the PCV2 viral load was not measured every week in the study animals; hence, it was difficult to trace the change of aspect against viral load from the initial point to the end point.

CONCLUSION

Although the measurement of viral load in this study showed two points, our findings confirm that the polymorphisms in the RGS16 gene have considerable effects on the PCV2 viral load in naturally infected pigs. Therefore, these polymorphisms can be regarded as an effective marker for estimating the proliferation of PCV2 viral load.

SUPPLEMENTARY MATERIALS

Supplementary Figures

jast-63-6-1223-suppl1.pdf

Competing interests

No potential conflict of interest relevant to this article was reported.

Funding sources

This work was carried out with the support of ‘Cooperative Research Program for Agriculture Science & Technology Development (Project No. PJ01620403),’ Rural Development Administration, Korea. This research was supported by the Chung-Ang University Research Grants in 2020.

Acknowledgements

Not applicable.

Availability of data and material

Upon reasonable request, the datasets of this study can be available from the corresponding author.

Authors’ contributions

Conceptualization: Hong KC, Kim JM.

Data curation: Kim JM.

Formal analysis: Lim KS.

Methodology: Lim KS.

Software: Lee SH.

Validation: Lee SH.

Investigation: Lee SH.

Writing - original draft: Lee SH.

Writing - review & editing: Kim JM.

Ethics approval and consent to participate

This study was conducted in accordance with the animal care guidelines of the Institutional Animal Care and Use Committee of the National Institute of Animal Science, Korea (2015-137).

REFERENCES

1.

Segalés J, Allan GM, Domingo M. Porcine circovirus diseases. Anim Health Res Rev. 2005; 6:119-42

2.

Fort M, Sibila M, Nofrarías M, Pérez-Martín E, Olvera A, Mateu E, et al. Porcine circovirus type 2 (PCV2) Cap and Rep proteins are involved in the development of cell-mediated immunity upon PCV2 infection. Vet Immunol Immunopathol. 2010; 137:226-34

3.

Cheung AK. Transcriptional analysis of porcine circovirus type 2. Virology. 2003; 305:168-80

4.

Timmusk S, Fossum C, Berg M. Porcine circovirus type 2 replicase binds the capsid protein and an intermediate filament-like protein. J Gen Virol. 2006; 87:3215-23

5.

Cheung AK. Porcine circovirus: transcription and DNA replication. Virus Res. 2012; 164:46-53

6.

Karuppannan AK, Liu S, Jia Q, Selvaraj M, Kwang J. Porcine circovirus type 2 ORF3 protein competes with p53 in binding to Pirh2 and mediates the deregulation of p53 homeostasis. Virology. 2010; 398:1-11

7.

Lin C, Gu J, Wang H, Zhou J, Li J, Wang S, et al. Caspase-dependent apoptosis induction via viral protein ORF4 of porcine circovirus 2 binding to mitochondrial adenine nucleotide translocase 3. J Virol. 2018; 92:e00238-18

8.

Segalés J, Calsamiglia M, Olvera A, Sibila M, Badiella L, Domingo M. Quantification of porcine circovirus type 2 (PCV2) DNA in serum and tonsillar, nasal, tracheo-bronchial, urinary and faecal swabs of pigs with and without postweaning multisystemic wasting syndrome (PMWS). Vet Microbiol. 2005; 111:223-9

9.

Finsterbusch T, Mankertz A. Porcine circoviruses—small but powerful. Virus Res. 2009; 143:177-83

10.

Chae C. A review of porcine circovirus 2-associated syndromes and diseases. Vet J. 2005; 169:326-36

11.

Chae C. Porcine circovirus type 2 and its associated diseases in Korea. Virus Res. 2012; 164:107-13

12.

Opriessnig T, Fenaux M, Thomas P, Hoogland MJ, Rothschild MF, Meng XJ, et al. Evidence of breed-dependent differences in susceptibility to porcine circovirus type-2-associated disease and lesions. Vet Pathol. 2006; 43:281-93

13.

Afghah Z, Webb B, Meng XJ, Ramamoorthy S. Ten years of PCV2 vaccines and vaccination: is eradication a possibility?. Vet Microbiol. 2017; 206:21-8

14.

de Castro AMMG, Opriessnig T. PCV2 genotypes and cross-protection [Internet]. Respig. 2020 [cited 2021 Aug 4]http://www.respigdigitalmagazine.com/EN/Respig1/index.html

15.

He J, Leng C, Pan J, Li A, Zhang H, Cong F, et al. Identification of lncRNAs involved in PCV2 infection of PK-15 cells. Pathogens. 2020; 9:479

16.

Mukherjee P, Milton AAP, Karam A, Chakraborty AK, Das S, Puro K, et al. Porcine circovirus2 (PCV2): a review of disease profile and control. Indian J Comp Microbiol Immunol Infect Dis. 2018; 39:38-50

17.

Lippert E, Yowe DL, Gonzalo JA, Justice JP, Webster JM, Fedyk ER, et al. Role of regulator of G protein signaling 16 in inflammation-induced T lymphocyte migration and activation. J Immunol. 2003; 171:1542-55

18.

Shi GX, Harrison K, Han SB, Moratz C, Kehrl JH. Toll-like receptor signaling alters the expression of regulator of G protein signaling proteins in dendritic cells: implications for G protein-coupled receptor signaling. J Immunol. 2004; 172:5175-84

19.

Timmusk S, Merlot E, Lövgren T, Järvekülg L, Berg M, Fossum C. Regulator of G protein signalling 16 is a target for a porcine circovirus type 2 protein. J Gen Virol. 2009; 90:2425-36

20.

Duthie C, Simm G, Doeschl-Wilson A, Kalm E, Knap PW, Roehe R. Quantitative trait loci for chemical body composition traits in pigs and their positional associations with body tissues, growth and feed intake. Anim Genet. 2008; 39:130-40

21.

Rückert C, Bennewitz J. Joint QTL analysis of three connected F2-crosses in pigs. Genet Sel Evol. 2010; 42:40

22.

Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers.In In: Misener S, Krawetz SA, editors.editors Bioinformatics methods and protocols. Totowa, NJ: Humana Press. 2000; p p. 365-86

23.

Swindell SR, Plasterer TN. SEQMAN.In In: Swindell SR, editor.editor Sequence data analysis guidebook. Totowa, NJ: Springer. 1997; p p. 75-89

24.

Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001; 68:978-89

25.

Hasenstein JR, Hassen AT, Dekkers JCM, Lamont SJ. High density SNP analysis of genome-wide associations with post-challenge burden levels of salmonella in two advanced intercross populations of chickens.In In: Hasenstein JR, editor.editor Single, several, saturated: multilevel SNP analysis for Salmonella response in young chicks. Ames, Iowa: Iowa State University. 2007; p p. 99-119

26.

Allan GM, Ellis JA. Porcine circoviruses: a review. J Vet Diagn Invest. 2000; 12:3-14

27.

Segalés J. Porcine circovirus type 2 (PCV2) infections: clinical signs, pathology and laboratory diagnosis. Virus Res. 2012; 164:10-19

28.

Hamel AL, Lin LL, Nayar GPS. Nucleotide sequence of porcine circovirus associated with postweaning multisystemic wasting syndrome in pigs. J Virol. 1998; 72:5262-67

29.

Cheung AK. Rolling-circle replication of an animal circovirus genome in a theta-replicating bacterial plasmid in Escherichia coli. J Virol. 2006; 80:8686-94

30.

Mankertz A, Çaliskan R, Hattermann K, Hillenbrand B, Kurzendoerfer P, Mueller B, et al. Molecular biology of porcine circovirus: analyses of gene expression and viral replication. Vet Microbiol. 2004; 98:81-8

31.

Nawagitgul P, Morozov I, Bolin SR, Harms PA, Sorden SD, Paul PS. Open reading frame 2 of porcine circovirus type 2 encodes a major capsid protein. J Gen Virol. 2000; 81:2281-7

32.

Liu J, Chen I, Du Q, Chua H, Kwang J. The ORF3 protein of porcine circovirus type 2 is involved in viral pathogenesis in vivo. J Virol. 2006; 80:5065-73

33.

Carasova P, Celer V, Takacova K, Trundova M, Molinkova D, Lobova D, et al. The levels of PCV2 specific antibodies and viremia in pigs. Res Vet Sci. 2007; 83:274-8

34.

Chae C. Postweaning multisystemic wasting syndrome: a review of aetiology, diagnosis and pathology. Vet J. 2004; 168:41-9

35.

McKeown NE, Opriessnig T, Thomas P, Guenette DK, Elvinger F, Fenaux M, et al. Effects of porcine circovirus type 2 (PCV2) maternal antibodies on experimental infection of piglets with PCV2. Clin Diagn Lab Immunol. 2005; 12:1347-51

36.

McNeilly F, McNair I, Stewart G, Allan G, Green L, Waldner C, et al. Post-weaning multisystemic wasting syndrome: Studies on disease progress in relation to serum antibody levels to porcine circovirus type 2 (PCV2) in sows and piglets and PCV2 viraemia in young pigs. Pig J. 2006; 57:70-89

37.

Zhang P, Wang L, Li Y, Jiang P, Wang Y, Wang P, et al. Identification and characterization of microRNA in the lung tissue of pigs with different susceptibilities to PCV2 infection. Vet Res. 2018; 49:18

38.

Walker LR, Engle TB, Vu H, Tosky ER, Nonneman DJ, Smith TPL, et al. Synaptogyrin-2 influences replication of porcine circovirus 2. PLOS Genet. 2018; 14e1007750

39.

Vicente J, Segalés J, Höfle U, Balasch M, Plana-Durán J, Domingo M, et al. Epidemiological study on porcine circovirus type 2 (PCV2) infection in the European wild boar (Sus scrofa). Vet Res. 2004; 35:243-53

40.

Kveberg L, Ryan JC, Rolstad B, Inngjerdingen M. Expression of regulator of G protein signalling proteins in natural killer cells, and their modulation by Ly49A and Ly49D. Immunology. 2005; 115:358-65

41.

Suurväli J, Pahtma M, Saar R, Paalme V, Nutt A, Tiivel T, et al. RGS16 restricts the pro-inflammatory response of monocytes. Scand J Immunol. 2015; 81:23-30