Skip to content

Advertisement

Virology Journal

What do you think about BMC? Take part in

Open Access

Genetic analysis of porcine circovirus type 2 from pigs affected with PMWS in Chile reveals intergenotypic recombination

  • Victor Neira1,
  • Natalia Ramos2,
  • Rodrigo Tapia1,
  • Juan Arbiza2,
  • Andrónico Neira-Carrillo3,
  • Manuel Quezada4,
  • Álvaro Ruiz4 and
  • Sergio A. Bucarey5Email author
Virology Journal201714:191

https://doi.org/10.1186/s12985-017-0850-1

Received: 6 January 2017

Accepted: 17 September 2017

Published: 4 October 2017

Abstract

Background

Porcine circovirus type 2 (PCV2) is a very small, non-enveloped and icosahedral virus, with circular single stranded DNA genome. This virus is the most ubiquitous and persistent pathogen currently affecting the swine industry worldwide. PCV2 has been implicated as the major causative agent of postweaning multisystemic wasting syndrome (PMWS), a disease which is characterized by severe immunosuppressive effects in the porcine host. Worldwide PCV2 isolates have been classified into four different genotypes, PCV2a, PCV2b, PCV2c and PCVd. The goal of this work was to conduct the first phylogenetic analysis of PCV2 in Chile.

Methods

PCV2 partial ORF2 sequences (462 nt) obtained from 29 clinical cases of PMWS in 22 Chilean intensive swine farms, covering over the 90% of the local pork-production, were analyzed.

Results

14% and 52% of sequences belonged to the genotypes PCV2a and PCV2b, respectively. Surprisingly, 34% of sequences were PCV2a/PCV2d recombinant viruses.

Conclusions

Our findings suggested that a novel cluster of Chilean sequences emerged resulting from intergenotypic recombination between PCV2a and PCV2d.

Keywords

PCV2PMWSGenetic diversityRecombination

Background

Porcine Circovirus type 2 (PCV2) is a small non-enveloped virus, member of the genus Circovirus, family Circoviridae [1]. PCV2 is primary causative agent of postweaning multisystemic wasting syndrome (PMWS) and other diseases in swine referred to as porcine circovirus associated diseases (PCVAD) [2]. PMWS was first described in Canada in 1991, and nowadays it has been described in pig population worldwide, having a serious economic impact on the global swine industry. Pigs affected with PMWS show decreased weight gain, progressive weight loss, enlarged lymph nodes, and others non-specific clinical signs such as dyspnea and diarrhea [3] .

The genome of PCV2 is a 1.76–1.77 kb single-stranded circular DNA, containing three main open reading frames (ORFs) which encode viral proteins [4]. ORF1 encodes for the replicase protein (Rep) and splice variants involved in viral replication, ORF2 encodes for the capsid (Cap), an exclusive structural protein [5], and ORF3 encodes a non-structural protein which causes apoptosis in PCV2 infected cells and has been associated with pathogenicity [6]. PCV2 isolates could be classified into four different genotypes, PCV2a, PCV2b, PCV2c and PCV2d, based on pairwise sequence comparisons (PASC) [7], and most recently based on reference sequences and/or marker positions [8]. PCV2a and PCV2b are the major groups, worldwide distributed, while PCV2c has been described only in Denmark from swine historic archived tissues [9] and recently from wild boar in Brazil [10]. PCV2d correspond to a PCV2b mutant (mPCV2b) described in some countries of Europe, Asia, North America and South America [8, 11]. Moreover, a new cluster of PCV2 has been recently described and proposed as the PCV2e genotype [12].

In South America, PMWS and PCVAD have been reported; phylogenetic studies have only been conducted in Brazil, Argentina, Uruguay and Colombia [10, 13, 14]. The objective of this study is to evaluate the genetic diversity of PCV2 using sequences in pig affected by PMWS obtained from intensive swine farms in Chile.

Methods

Samples

The samples used in this study were obtained from fixed tissues submitted between 2005 and 2008 to the Veterinary Pathology Laboratory, University of Concepcion, based on PMWS suspected clinical signs. Samples from 20 intensive pig farms located in 5 administrative regions of Chile Central, plus an extensive wild boar farm were included, and data about geographic localization was registered. Further, clinical PMWS diagnosis was confirmed by macroscopic, microscopic lesions and “in situ” hybridization (ISH) based on the technique described by Segalés and Domingo, 2002 [3]. Subsequently, a subset of 51 PCV2 positive samples were selected for PCR testing and sequencing based on its ISH intensity.

PCV2 DNA sequencing

Paraffin embedded tissues were manually deparaffinized using xylene. After, the DNA was extracted using QIAamp DNA mini kit (Qiagen, Valencia, USA) following manufacturer’s instructions. First, samples were tested using a selective PCR protocol, which is able to differentiate PCV2a and PCV2b viruses and identify mixed infections [15]. Also, samples were tested with sequencing PCR primers based on protocol described by Quintana et al., 2002 [16]. This PCR protocol amplifies a 652 bp fragment located between genomic positions 940 and 1596. PCR products were visualized in a 1.8% agarose gel and purified using QIAquick PCR Purification Kit. Sequencing was performed directly on the PCR products, which were sequenced on an ABI PRISM 3100 Avant Genetic Analyzer (Applied Biosystems) at Macrogen USA (MD, United States).

Genotyping study

Molecular characterization of PCV2 strains was performed on the basis of sequence analysis of 462 nt of the ORF2. Sequence alignment was carried out with Clustal W algorithm including reference sequences corresponding to different PCV2 genotypes [8] and sequences from neighbor countries. The degree of identity among nucleotide and amino acid sequences was determined by using BioEdit v. 7.0.5. Phylogenetic tree was reconstructed by the neighbor-joining method using MEGA v5.0 software. Bootstrap values were determined with 1000 replicates of the dataset and PCV1 sequence was included as outgroup (FJ475129).

Recombination analysis

In order to analyze the potential recombinant strains, an integrated software package, the recombination detection program (RDP v.3.44), was used [17]. The general settings were a highest acceptable p-value of 0.01 and Bonferroni correction. A dataset of selected PCV2 sequences of different genotypes (n = 18), were included in the RDP analysis excluding known recombinant sequences. Each potential recombinant strain was analyzed separately.

Additionally, SIMPLOT software v. 3.5.1 [18] was employed to further analyze the potential parental PCV2 lineages and putative recombination breakpoints previously analyzed and estimated by RDP software. Four groups of 462 nt of ORF2 sequences were included in SIMPLOT analysis: a group of potential recombinants (O1/2008, M1/2007, L1/2007, H3/2007, H4/2007, B4/2008, C1/2006, B1/2005, B2/2005), a group of a reference sequences from PCV2d genotype (AY943819; JX535296), a group of PCV2a genotype sequences (AY180397; AY146993) and an out-group (PCV1 isolate FJ475129). A similarity plot was generated comparing the percentage of identity versus genome position of the query isolate to the panel of references sequences. Analysis was performed with a window size of 100 nt and a step size of 20 nt. The analysis of potential recombinant strains was also carried out by the identification of conflicting phylogenetic trees in different regions of the alignment based on the putative recombination breakpoints observed in SIMPLOT and RDP analysis. Phylogenetic reconstructions were realized as described above. One phylogenetic tree was constructed on the basis of the minor parent region (nucleotides 100 to 389) and the other comprised the major parent region concatenated alignment of nucleotides 1 to 99 and 390 to 462.

Results

PCV2 genetic material was amplified by selective PCR in 19 out 51 samples. No mixed infections were identified. Ten samples, tested negative for selective PCR, were amplified by sequencing primers. Therefore, 29 out 51 tissues samples were sequenced, 28 corresponding to 18 intensive swine farms and 1 from an extensive wild boar farm. Genomic sequences were deposited into GenBank™ with accession numbers (KU992866 to KU992894) and each sequence included a farm code and year (Table 1).
Table 1

List of PCV2 sequences including farm ID, year of collection, administrative region, and genotype

Farm

Year

Region

Genotype

Sequence ID

Accession No

H

2005

VI

PCV2a

H/2005

KU992879

N

2005

VIII

PCV2b

N/2005

KU992890

B

2005

VIII

PCV2rec

B1/2005

KU992870

B

2005

VIII

PCV2rec

B2/2005

KU992871

K

2006

RM

PCV2b

K/2006

KU992885

N

2006

VIII

PCV2b

N/2006

KU992891

C

2006

VI

PCV2rec

C/2006

KU992875

J

2007

RM

PCV2a

J/2007

KU992874

A

2007

RM

PCV2b

A1/2007

KU992866

A

2007

RM

PCV2b

A2/2007

KU992867

A

2007

RM

PCV2b

A3/2007

KU992868

A

2007

RM

PCV2b

A4/2007

KU992869

B

2007

VIII

PCV2b

B/2007

KU992872

D

2007

VIII

PCV2b

D/2007

KU992876

E

2007

IX

PCV2b

E/2007

KU992877

G

2007

VI

PCV2b

G/2007

KU992892

I

2007

VII

PCV2b

I/2007

KU992884

P

2007

VIII

PCV2b

P/2007

KU992886

R

2007

VIII

PCV2b

R/2007

KU992894

H

2007

VI

PCV2rec

H1/2007

KU992880

H

2007

VI

PCV2rec

H2/2007

KU992887

H

2007

VI

PCV2rec

H3/2007

KU992881

L

2007

RM

PCV2rec

L/2007

KU992888

M

2007

VI

PCV2rec

M/2007

KU992889

F

2008

VII

PCV2b

F1/2008

KU992878

B

2008

VIII

PCV2rec

F2/2008

KU992873

O

2008

RM

PCV2rec

O/2008

KU992893

Q

2013

RM

PCV2a

Q1/2013

KU992882

Q

2013

RM

PCV2a

Q2/2013

KU992883

Phylogenetic analysis using the partial ORF2 sequences obtained (462 nt) showed that PCV2 sequences belonged to PCV2b (15 out 29), PCV2a (4 out 29), and 10 sequences segregated independent labeled as Cluster R (Fig. 1). PCV2b Chilean sequences came from 10 pig farms and a nucleotide identity of 99.9% was observed in this group. On the other hand, PCV2a Chilean sequences were obtained from 3 pig farms, sharing a nucleotide identity of 96.9% among them.
Fig. 1

Phylogenetic tree generated by neighbor-joining method using 462-nt sequence of ORF2 PCV2. Chilean strains were compared with twenty-five PCV2 reference strains. PCV1 was included as out-group. GenBank accession numbers and country of origin are specified. Genotypes are indicated with straight brackets

Ten isolates were grouped in a new cluster with a very high statistical support (bootstrap value, 100%) (Fig. 1), suggested the possible presence of recombinant strains. This cluster contained viruses from 6 different pig farms with nucleotide identity of 99.5%. These sequences were detected as recombinants using RDP3 software and the recombination events were statistically supported by BootScan, Maxchi, Chimera and 3Seq methods with average p-values of 1.9 × 10−5, 6.6 × 10−4, 1.6 × 10−2 and 4.3 × 10−3, respectively. Results showed that all recombinant sequences were originated from the same parental lineages. The parent sequences resulted to be JX535296 (major parent, PCV2d genotype) and AY180397 (minor parent, PCV2a genotype). In addition, a similarity plot analysis was carried out with putative parental sequences using SIMPLOT software to corroborate the results obtained using RDP3 software. In the Fig. 2a is shown the similarity plot generated for B2/2005 sequence as an example, indicated that the recombinant presented a nucleotide similarity with a PCV2a isolate (blue) between nucleotides 91 and 391, and with an PCV2d strain (red) between nucleotides 1 to 91 and 391 to 462. Moreover, recombination graph of the recombination event for B2/2005 is also shown using BootScan method (Fig. 2b) and nucleotides 132 and 386 were predicted by RDP software to be the potential breakpoints.
Fig. 2

Recombination analysis. a Similarity plot generated by SIMPLOT software based on 462 nt of ORF2 sequence using B2/2005 as query sequence and representative parent sequence of genotypes PCV2d (red) and PCV2a (blue) (JX535296 and AY180397, respectively). A PCV1 isolate was used as outgroup (green). The analysis was performed with a window size of 100 nt and a step size of 20 nt. b BootScan analysis and potential breakpoints for B2/2005 predicted by RDP3 software

On the other hand, two phylogenies were constructed based on the major and minor parental regions taking to account data obtained of SIMPLOT and RDP software (Fig. 3). In the phylogeny based on the minor parental region, Chilean potential recombinant sequences clustered within PCV2a genotype (Fig. 3a). However, in the phylogeny based on the major parental region, recombinant viruses resulted to be more closely related to a set of PCV2d sequences (Fig. 3b). The occurrence of incongruent phylogenetic trees strongly provided an additional evidence of recombination. Finally, the closest sequence publicly available presents an identity of 98% (p-distance) and was obtained from Thailand (PCV2/ Thailand/THKUF49/2007, # HQ701666), which correspond a PCV2d virus.
Fig. 3

Phylogenetic trees generated by neighbor-joining method based on a the minor parent region estimated with SIMPLOT and RDP software (nt 100–389) and b the major parent region (concatenated alignment of nt 1–99 and nt 390–462). Chilean strains were compared with twenty-five PCV2 reference strains. GenBank accession numbers and country of origin are specified. Genotypes are indicated with straight brackets. Discordant clustering patterns of potential Chilean recombinant strains with different genotypes are observed ()

Discussion

The study conducts the first phylogenetic analysis of PCV2 in Chile. PCV2 isolated from Chilean intensive swine farms, covering over the 90% of the pork production, resulting a complete national perspective of PCV2. The results confirm the high presence of PCV2b (52%) and less represented PCV2a (14%). Surprisingly, 34% of sequences were PCV2a/PCV2d recombinant viruses. The no amplification of these 10 sequences by the selective PCR was further explained due to these primers were not able to amplified PCV2d, the major parenteral of these viruses. Furthermore, we were not able to detect PCV2d variant in the Chilean PCV2 diversity even if was the major parent of the recombinant viruses, probably this event occurred years before and the recombinant remained. PCV2a and PCV2b, and the recombinant viruses were found since 2005. Nevertheless, were not observed geographical, temporal and other patterns, most likely due to there is not enough sequences per farm.

Recently, a large-scale study determinates that the PCV2 vaccination is important factor in the PCV2 ecology, driving the evolution and epidemiology of the virus. PCV2 vaccination by selective pressure generates vaccine scape-mutants, which become new genotypes waves [19]. Interestingly, the vaccination against PCV2 in Chile started in 2007 and massive in swine farms around 2008–2009. This study includes samples since 2005 to 2013, but differences before and after massive vaccination were not reported mostly due low number of isolates per farm and year.

Six farms presented PCV2 recombinant viruses; this event was observed in samples from 2005. However, even if geographical or others patterns were not identified, probably the presence of these recombinant viruses is linked between farms. There are not similar sequences reported with this particular recombination in the world; therefore, probably this recombination event occurred in Chile and then was spread to the rest of the farms.

Recombination events are common among PCV2 strains, which have been reported in several countries [20, 21]. Also, different patterns of recombination have been described, with break points in both ORF1 and ORF2, which suggest that the recombination can occur easily in PCV2 co-infections [22].

On the other hand, PCV2b was the most common genotype in Chilean farms during the PMWS outbreak, as has been found of the rest of the world, where more than 50% of the viral sequences belonged into PCV2b [11]. PCV2b sequences did not show divergence over time among samples obtained in 2005 to 2008 from several different farms and present 99.9% identity with viruses found in the rest of the world.

Conclusions

This is the first study that revealed the genetic diversity of PCV2 in Chile during a PWMS outbreak. Our findings suggested that a novel cluster of Chilean sequences emerged resulting from recombination between genotypes PCV2a and PCV2d within ORF2. However, analysis of the complete PCV2 genome and further surveillance of PCV2 is needed to confirm the real impact of these recombinant viruses.

Declarations

Acknowledgments

We thank to Michelle Lineros, Daniela Rojas and Miguel Lopez for technical support.

Funding

This work has been supported by FONDEF from Chilean Council for Science and Technology (CONICYT), program IDeA, project IT13I20021 granted to Sergio A. Bucarey.

Availability of data and materials

All DNA sequences have been submitted in genbank (http://www.ncbi.nlm.nih.gov/genbank/) under the following accession numbers: KU992879, KU992890, KU992870, KU992871, KU992885, KU992891, KU992875, KU992874, KU992866, KU992867, KU992868, KU992869, KU992872, KU992876, KU992877, KU992892, KU992884, KU992886, KU992894, KU992880, KU992887, KU992881, KU992888, KU992889, KU992878, KU992873, KU992893, KU992882, KU992883.

Authors’ contributions

VN and RT performed the tests, analyzed and interpreted sequencing analysis and contributing in writing the manuscript. NR performed molecular analysis including genotyping and recombination studies and contributing in writing the manuscript. JR interpreted molecular analysis. MQ and AR performed the sample collection, macroscopic and histological examination of PMWS samples. AN contributed in writing the manuscript. SAB performed PCR tests, analyzed and interpreted genotyping data and contributed in writing the manuscript. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Departamento de Medicina Preventiva, Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile
(2)
Sección Virología, Facultad de Ciencias, Universidad de la República
(3)
Laboratorio Polyforms, Departamento de Ciencias Biológicas Animales, Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile
(4)
Facultad de Ciencias Veterinarias, Universidad de Concepción
(5)
Centro Biotecnológico Biovetec, Departamento de Ciencias Biológicas Animales, Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile

References

  1. Tischer I, Gelderblom H, Vettermann W, Koch MA. A very small porcine virus with circular single-stranded DNA. Nature. 1982;295(5844):64–6.View ArticlePubMedGoogle Scholar
  2. Allan G, Meehan B, Todd D, Kennedy S, McNeilly F, Ellis J, Clark EG, Harding J, Espuna E, Botner A, et al. Novel porcine circoviruses from pigs with wasting disease syndromes. Vet Rec. 1998;142(17):467–8.PubMedGoogle Scholar
  3. Segales J, Domingo M. Postweaning multisystemic wasting syndrome (PMWS) in pigs. A review. Vet Q. 2002;24(3):109–24.View ArticlePubMedGoogle Scholar
  4. Zhai SL, Chen SN, Xu ZH, Tang MH, Wang FG, Li XJ, Sun BB, Deng SF, Hu J, Lv DH, et al. Porcine circovirus type 2 in China: an update on and insights to its prevalence and control. Virol J. 2014;11:88.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Cheung AK. Comparative analysis of the transcriptional patterns of pathogenic and nonpathogenic porcine circoviruses. Virology. 2003;310(1):41–9.View ArticlePubMedGoogle Scholar
  6. Choi CY, Rho SB, Kim HS, Han J, Bae J, Lee SJ, Jung WW, Chun T. The ORF3 protein of porcine circovirus type 2 promotes secretion of IL-6 and IL-8 in porcine epithelial cells by facilitating proteasomal degradation of regulator of G protein signalling 16 through physical interaction. J Gen Virol. 2015;96(Pt 5):1098–108.PubMedGoogle Scholar
  7. Cortey M, Pileri E, Sibila M, Pujols J, Balasch M, Plana J, Segales J. Genotypic shift of porcine circovirus type 2 from PCV-2a to PCV-2b in Spain from 1985 to 2008. Vet J. 2011;187(3):363–8.View ArticlePubMedGoogle Scholar
  8. Franzo G, Cortey M, Olvera A, Novosel D, Castro AM, Biagini P, Segales J, Drigo M. Revisiting the taxonomical classification of porcine Circovirus type 2 (PCV2): still a real challenge. Virol J. 2015;12:131.View ArticlePubMedPubMed CentralGoogle Scholar
  9. Dupont K, Nielsen EO, Baekbo P, Larsen LE. Genomic analysis of PCV2 isolates from Danish archives and a current PMWS case-control study supports a shift in genotypes with time. Vet Microbiology. 2008;128(1–2):56–64.View ArticleGoogle Scholar
  10. Franzo G, Cortey M, de Castro AM, Piovezan U, Szabo MP, Drigo M, Segales J, Richtzenhain LJ. Genetic characterisation of porcine circovirus type 2 (PCV2) strains from feral pigs in the Brazilian Pantanal: an opportunity to reconstruct the history of PCV2 evolution. Vet Microbiology. 2015;178(1–2):158–62.View ArticleGoogle Scholar
  11. Xiao CT, Halbur PG, Opriessnig T. Global molecular genetic analysis of porcine circovirus type 2 (PCV2) sequences confirms the presence of four main PCV2 genotypes and reveals a rapid increase of PCV2d. J Gen Virol. 2015;96(Pt 7):1830–41.View ArticlePubMedGoogle Scholar
  12. Davies B, Wang X, Dvorak CM, Marthaler D, Murtaugh MP. Diagnostic phylogenetics reveals a new porcine circovirus 2 cluster. Virus Res. 2016;217:32–7.View ArticlePubMedGoogle Scholar
  13. Pereda A, Pineyro P, Bratanich A, Quiroga MA, Bucafusco D, Craig MI, Cappuccio J, Machuca M, Rimondi A, Dibarbora M, et al. Genetic characterization of porcine circovirus type 2 from pigs with porcine circovirus associated diseases in Argentina. ISRN Vet Sci. 2011;2011:560905.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Rincon Monroy MA, Ramirez-Nieto GC, Vera VJ, Correa JJ, Mogollon-Galvis JD. Detection and molecular characterization of porcine circovirus type 2 from piglets with porcine circovirus associated diseases in Colombia. Virol J. 2014;11:143.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Hesse R, Kerrigan M, Rowland RR. Evidence for recombination between PCV2a and PCV2b in the field. Virus Res. 2008;132(1–2):201–7.View ArticlePubMedGoogle Scholar
  16. Quintana J, Balasch M, Segales J, Calsamiglia M, Rodriguez-Arrioja GM, Plana-Duran J, Domingo M. Experimental inoculation of porcine circoviruses type 1 (PCV1) and type 2 (PCV2) in rabbits and mice. Vet Res. 2002;33(3):229–37.View ArticlePubMedGoogle Scholar
  17. Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. 2010;26(19):2462–3.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73(1):152–60.PubMedPubMed CentralGoogle Scholar
  19. Franzo G, Tucciarone CM, Cecchinato M, Drigo M. Porcine circovirus type 2 (PCV2) evolution before and after the vaccination introduction: a large scale epidemiological study. Sci Rep. 2016;6:39458.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Huang Y, Shao M, Xu X, Zhang X, Du Q, Zhao X, Zhang W, Lyu Y, Tong D. Evidence for different patterns of natural inter-genotype recombination between two PCV2 parental strains in the field. Virus Res. 2013;175(1):78–86.View ArticlePubMedGoogle Scholar
  21. Ramos N, Mirazo S, Castro G, Arbiza J. Molecular analysis of porcine Circovirus type 2 strains from Uruguay: evidence for natural occurring recombination. Infect Genet Evol. 2013;19:23–31.View ArticlePubMedGoogle Scholar
  22. Ma CM, Hon CC, Lam TY, Li VY, Wong CK, de Oliveira T, Leung FC. Evidence for recombination in natural populations of porcine circovirus type 2 in Hong Kong and mainland China. J Gen Virol. 2007;88(Pt 6):1733–7.View ArticlePubMedGoogle Scholar

Copyright

© The Author(s). 2017

Advertisement