Next Article in Journal
African Swine Fever Virus (ASFV): Biology, Genomics and Genotypes Circulating in Sub-Saharan Africa
Previous Article in Journal
Heterogeneity in the Response of Different Subtypes of Drosophila melanogaster Midgut Cells to Viral Infections
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 11
10.3390/v13112282
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
Comment published on 5 May 2022, see Viruses 2022, 14(5), 962.
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Epidemiology and Evolution of Emerging Porcine Circovirus-like Viruses in Pigs with Hemorrhagic Dysentery and Diarrhea Symptoms in Central China from 2018 to 2021

by
Kankan Yang
1,2,†,
Menghuan Zhang
1,2,†,
Qi Liu
1,2,†,
Yingli Cao
1,2,
Wuyin Zhang
1,2,
Yueqiao Liang
1,2,
Xiangjun Song
1,2,
Kaiyuan Ji
1,2,
Ying Shao
1,2,
Kezong Qi
1,2,* and
Jian Tu
1,2,*
1
Anhui Province Engineering Laboratory for Animal Food Quality and Bio-Safety, College of Animal Science and Technology, Anhui Agricultural University, Hefei 230036, China
2
Anhui Province Key Laboratory of Veterinary Pathobiology and Disease Control, Anhui Agricultural University, Hefei 230036, China
*
Authors to whom correspondence should be addressed.
†
Menghuan Zhang, Qi Liu and Kankan Yang contributed equally to this work.
Viruses 2021, 13(11), 2282; https://doi.org/10.3390/v13112282
Submission received: 3 October 2021 / Revised: 8 November 2021 / Accepted: 10 November 2021 / Published: 15 November 2021
(This article belongs to the Topic Veterinary Infectious Diseases)
Browse Figures
Versions Notes

Abstract

:
Porcine circovirus-like virus (PCLV) is a type of circular Rep-encoding single-stranded DNA virus and may be associated with the development of diarrheal symptoms in pigs. In this study, we retrospectively analyzed three years of past cases in Anhui, China, and reported a case of hemorrhagic enteritis and death in a pregnant sow possibly caused by PCLV. In addition, we analyzed the evolutionary characteristics of PCLV and found that mutation, recombination and selective pressure all played an important role in the evolution of PCLV. We identified N15D and T17S as well as L56T, T58R, K59Q, M62R, L75I and R190K mutations in two different branches, and we noted recombination events in the Rep of a group of Chinese strains. Analysis of selection pressure revealed that PCLV gained more positive selection, indicating that the virus is in a continuous evolutionary state. The PR2 plot, ENC-plot and neutrality analysis showed a greater role of natural selection than that of mutational pressure in the formation of codon usage patterns. This study is the first to identify PCLV in sows with hemorrhagic dysentery and death, and it provides new epidemiological information on PCLV infection in pigs in China.

1. Introduction

Circular Rep-encoding single-stranded DNA (CRESS DNA) viruses form a highly diverse group of small viruses that have been found worldwide in prokaryotic, eukaryotic and even environmental samples [1,2]. The first eukaryotic CRESS DNA viruses were not discovered until the 1970s, although symptoms consistent with CRESS DNA virus infections were described more than a thousand years ago in plants [3]. Currently, CRESS DNA viruses are divided into seven families: Bacilladnaviridae, Nanoviridae, Smacoviridae, Geminiviridae, Genomoviridae, Redondoviridae and Circoviridae [4]. Members of the Circoviridae, Smacoviridae and Redondoviridae families are known to infect animals, including mosquitoes, rats, bats, ducks, cattle, pigs, dogs, humans, turkeys, and forest musk deer [5,6,7,8]. In addition, a new family, Kirkoviridae, was proposed by The International Committee on Taxonomy of Viruses (ICTV). The genome structure of members of this family contains only one open reading frame (ORF), encoding only the Rep, rather than containing two ORFs as in the case of members of other families [2,9,10].
Porcine circovirus-like virus (PCLV) is a closed circular, non-enveloped, single-stranded DNA virus, similar to porcine circovirus (PCV), and it was identified as a member of the family Kirkoviridae. As mentioned above, PCLV, containing only one ORF, is significantly different from PCV in gene structure [11]. The virus was discovered in fecal samples in the United States in 2011 [12]. In 2020 and 2021, the virus was found in Guang-xi and Guangdong provinces, located in the southern part of China [13,14]. At present, most studies show that the virus is found in piglets and is less likely to be reported in adult pigs; it causes hemorrhagic enteritis in piglets, seriously compromising their health [12,14]. In addition to causing clinical signs in piglets, the virus can also cause co-infections with other viruses, such as porcine circovirus 2 (PCV-2) and porcine epidemic diarrhea virus (PEDV), making prevention and control of the virus infection challenging.
Based on this, 140 diarrheal and non-diarrheal clinical samples (including lung, intestine, spleen, fecal, etc.) were collected from 13 pig farms in five cities in Anhui Province, China, spanning from June 2018 to August 2021. The samples were retrospectively investigated in this study to evaluate and enrich knowledge of the epidemiological characteristics and genetic diversity of PCLV. An unexplained case of sow diarrhea and death on a pig farm in Anhui Province on 15 July 2021 was also analyzed. Polymerase chain reaction (PCR) was used to monitor the prevalence of PCLV in central China, to amplify the full genome length of three strains, and to analyze the homology and evolutionary characteristics between different PCLV genomes. In this study, basic information on PCLV is provided to facilitate further understanding of its characteristics; however, its clinical significance, epidemiology and disease prevention deserve further in-depth study.

2. Materials and Methods

2.1. Sample Collection

To investigate the molecular genetic diversity and epidemiological characteristics of PCLV circulating in Anhui Province, China, and whether they correlate with signs of diarrhea, a total of 140 clinical samples from diarrheic and non-diarrheic adult pigs and piglets were collected from thirteen pig farms in Anhui Province from June 2018 to August 2021. All samples were transported packed on ice and stored at −80 °C.

2.2. Viral Nucleic Acid Extraction and PCR Detection

Fecal material (2 g) was dissolved in an Eppendorf tube containing 10% phosphate-buffered saline (PBS) and clarified by centrifugation for 3 min at 12,000× g. Tissue samples were added to PBS at a ratio of 1:3 and homogenized in a mortar, followed by centrifugation for 5 min at 13,000× g. Following this, 200 µL of supernatant was transferred to sterilized 2 mL centrifuge tubes, and viral DNA was extracted using the TIANamp Virus DNA/RNA Kit (Tiangen, Beijing, China), according to the manufacturer’s instructions. The known genomic sequence of PCLV was retrieved from Genbank and analyzed using the Lasergene package (DNAStar Inc., Madison, WI, USA). A pair of detection primers was designed for the conserved region (Table S1), and PCR was performed under the following conditions: an initial denaturation step at 94 °C for 5 min; 35 cycles at 94 °C for 30 s, at 49 °C for 30 s, and at 72 °C for 30 s; and finally extension at 72 °C for 5 min. The PCR products were analyzed by agarose gel electrophoresis, and the target bands were cut off for DNA extraction using a gel extraction kit and sent to TsingKe Co., Ltd. (Nanjing, China), China for sequencing before blast search verification.

2.3. PCLV Sequencing of the Full-Length Genome

Full-length primers were designed using the primers in previously published studies and for the reference strain (accession number: NC_025682.1) to amplify the genomes of samples that tested positive for PCLV. The primers are shown in Table S1. The geographical distribution of positive samples was marked on a map (Figure 1). All amplification products were visualized by agarose gel electrophoresis and then purified using a DNA purification kit (TIANGEN, Beijing, China) according to the manufacturer’s instructions. The PCR products were cloned into the pMD19-T vector (TaKaRa, Kusatsu, Japan), and the recombinant plasmids were transformed into Trelief™ 5α chemically competent cells (TsingKe, Nanjing, China), extracted using a QIAGEN Plasmid Mini Kit, and sent to TsingKe Co., Ltd. (Nanjing, China) for bidirectional sequencing. The reference strains were obtained from the NCBI database (https://www.ncbi.nlm.nih.gov/nucleotide/; accessed on 2 September 2021 (Table S3). Multiple sequence alignments were constructed using Clustal W.

2.4. Sequence Alignment, Phylogenetic and Recombination Analysis

The reference strains were sourced from the NCBI database, as shown in Table S3. Sequence identity was analyzed using Megalign (DNASTAR, Madison, WI, USA). Sequence alignment was performed based on nucleotide sequences via MAFFT (https://mafft.cbrc.jp/alignment/software/, accessed on 29 August 2021) and ESPript Software (https://espript.ibcp.fr/ESPript/ESPript/; accessed on 2 September 2021). Phylogenetic analysis was performed using PhyloSuite software [15]. Maximum likelihood phylogenies were inferred using IQ-TREE for 5000 ultrafast bootstraps and the Shimodaira–Hasegawa-like approximate likelihood-ratio test [16,17]. ModelFinder was used to select the best substitution model [18]. Phylogenetic trees were visualized using FigTree software.
Recombination events were identified using seven methods in the Recombination Detection Program v.4.39 (RDP 4.39), namely, RDP, GENECONV, Chimaera, MaxChi, BootScan, SiScan and 3Seq. SimPlot software v.3.5.1 was used to validate and visualize reorganization events.

2.5. Selection Pressure, B-Cell Epitope Prediction and Codon Usage Bias Analyses

Analysis of the selection pressure on Rep was performed using four methods, namely, Single-Likelihood Ancestor Counting (SLAC), Fixed Effects Likelihood (FEL), Fast Unconstrained Bayesian AppRoximation for inferring selection (FUBAR), and Mixed Effects Model of Evolution (MEME) on DATAMONKEY (http://www.datamonkey.org/; accessed on 20 August 2021) [19,20,21]. Significant recombinant strains were excluded. The B-cell epitope was predicted via SVMTriP (http://sysbio.unl.edu/SVMTriP/ (accessed on 10 August 2021) [22]. Protein structures were visualized using PyMol v4.60.
Analysis of codon usage bias was performed online using the CAIcal SERVER (https://ppuigbo.me/programs/CAIcal/; accessed on 15 August 2021). The host’s (Sus scrofa) codon usage pattern was obtained from the codon usage database (http://www.kazusa.or.jp/codon/; accessed on 26 August 2021). A Parity Rule 2 (PR2) plot was used to analyze whether the effects of mutation pressure and natural selection on the use of the Rep codon were consistent. In a PR2 plot, [A3/(A3 + T3)] is the ordinate and [G3/(G3 + C3)] is the abscissa. The center of the plot shows A = T and G = C, indicating that there is no deviation between mutation pressure and the influence of natural selection [23,24,25].
The value of ENC ranged from 20 to 60, negatively correlating with codon usage bias. ENC ≤ 35 indicates high codon usage bias, and ENC > 50 indicates that the codon usage is low [26]. The ENC-GC3s plot (ENC-plot) was analyzed with the ENC values plotted against the GC3S values to determine the main factors affecting codon use bias. The expected ENC values for each GC3S were estimated as follows [26]:
y = 2 + x + 29 x 2 + 1 − x 2
where x represents the value of G + C at the third codon position (GC3S).
Neutrality analysis was used to evaluate the effects of mutation pressure and natural selection on codon usage patterns. With GC12 as ordinate and GC3 as abscissa, the neutrality plot was generated. If the slope of the regression line was close to 1 (the point on the diagonal), there was no or weak selection pressure. In addition, the deviation from the diagonal of the regression line indicated the magnitude of the influence of natural selection [27]. Data visualization was completed using Origin 2018.

3. Results

3.1. Prevalence of PCLV

A total of 140 samples collected from different regions (Suzhou, Huaibei, Bozhou, Bengbu and Hefei) of Anhui province were tested, and three positive samples, collected in the years 2018, 2019 and 2021 from three different pig farms, A, B and C, were found (Figure 1), with a positive rate of 2.14% (3/140). As mixed infections of pathogens are common in clinically affected pigs, we investigated the mixed infections in three samples. Sample 1 (AH-23, accession no. MZ773067) from farm A was a fecal sample from a weaned piglet that showed signs of diarrhea, and PCLV was detected mixed with porcine astrovirus and PEDV. The age of the pig from which sample 2 was collected (AH-25, accession number: MZ773068) on farm B could not be determined, and the fecal sample was co-infected with PCV-2 and porcine parvovirus (PPV). Notably, sample 3 (AH-HB-2021, accession no. MZ960935) from farm C was collected from a pregnant sow with hemorrhagic enteritis, which eventually led to the death of the sow (Figure S1), and eight sows from the batch showed signs of diarrhea, with four cases resulting in death. Sample 3 was negative for other viruses routinely infecting pigs (African swine fever virus (ASFV), pseudorabies virus (PRV), porcine circovirus 2 (PCV-2), porcine circovirus 3 (PCV-3), porcine reproductive and respiratory syndrome virus (PRRSV), porcine epidemic diarrhea virus (PEDV), transmissible gastroenteritis virus (TGEV), porcine parvovirus (PPV), porcine rotavirus A (PRV-A) and classical swine fever virus (CSFV), among others) using commercially available kits, and primer sequences for pathogens (porcine astrovirus (PAstV), porcine kobuvirus (PKV), porcine deltacoronavirus (PDCoV), porcine bufavirus (PBuVs) and porcine sapovirus (PSaV)) without commercially available kits are shown in Table S2.

3.2. Full-Sequence Alignment of PCLV

We successfully amplified the full length of three PCLV strains, which had nucleotide lengths of 3955, 3946 and 3832 bp, respectively. In AH-23 and AH-HB-2021, five potential ORFs were identified, while only four potential ORFs existed in AH-25 (Figure 2). Sequence alignment showed that the similarity between the three strains obtained in this study was 86.8~91.1%, and the similarity between the three strains and reference strains was 77.7~94.6%. The nucleotide similarity and amino acid similarity of Rep among all strains were 83.4~97.5% and 84.7~99.6%, respectively. Analysis of amino acid mutations revealed that most amino acids in Rep were relatively consistent, with a small number of non-synonymous substitutions, while there were relatively more mutations in the C-terminus of Rep. Notably, N15D and T17S mutations were found in five strains. In addition, a wider range of consistent mutations was found in four other strains. The mutation sites were L56T, T58R, K59Q, M62R, L75I and R190K (Figure 3). Interestingly, from the phylogenetic tree for the Rep (Figure 4A), we found that Group 1 had mutations of N15D and T17S as mentioned above, while Group 2 had mutations of L56T, T58R, K59Q, M62R, L75I and R190K.

3.3. Phylogenetic Analysis and Genetic Divergence

Phylogenetic analysis based on complete sequences showed that all strains were divided into three groups, tentatively named Groups A, B and C. Among them, strains AH-23 and AH-HB-2021 were in a branch with the Po-Circo-like virus 22 strain found in the United States in 2010, while AH-25 was closer to the Po-Circo-like virus isolate CZQ11 isolated in China (Figure 4B). We compared the average between-group genetic distance for each genotype of PCV-2, and the results showed that the genetic distance between the three groups of PCLV was similar to that between PCV-2 types 2e and 2d, less than for 2a–2b, 2a–2c, 2a–2d, 2a-2e, 2a-2f and 2a-2g, and larger than that between other genotypes (Figure 5).

3.4. Recombination, Selection Pressure and B-Cell Antigenic Epitope Analysis

In the recombination analysis, extensive recombination events in Rep were found in a group of Chinese strains (Recombination Group, Figure 4A and Figure 6), and these were detected by four different methods (MaxChi, Chimaera, SiScan, 3Seq). The represented major parent strain and the minor parent strain were strain PCLV-AH-25 (accession number: MZ773068/China) and strain 21 (accession number: JF713716.1/USA), respectively. Analysis of selection pressure showed that PCLV was subject to positive selection (Table 1), which may lead to the tendency of the virus to mutate and the continuous evolution of the virus, possibly prompting the emergence of new varieties [29,30]. In order to explore whether the selection would affect the immunogenicity of Rep, the B-cell epitopes were predicted (Table S4), and the results showed that two positive selection sites overlapped with immune epitopes in AH-23 (codons 142 and 211), six in AH-25 (codons 29, 142, 81, 211, 290 and 72) and two in AH-HB-2021(codons 47 and 211) (Figure 7).

3.5. Codon Usage Bias Analysis

The ENC average of PCLV was 45.296 (41.2~50.3), indicating that PCLV has a low codon bias. PR2 plot analysis found that no strains clustered in the center of the coordinate axis (Figure 8A), indicating that the influence of mutation pressure and selection pressure on codon bias was inconsistent. To explore which factor was more biased towards codons, ENC-plot analysis was performed (Figure 8B). The results showed that all strains were below the expected curve, indicating that selection pressure plays a more important role in Rep during evolution. Neutrality analysis also suggested that natural selection had a greater effect on codon bias than mutation pressure (Figure 8C); the contribution of natural selection was about 73.9%.

4. Discussion

Circular Rep-encoding (replication-associated protein encoding) single-stranded DNA (CRESS DNA) viruses are a major component of the earth virome and can be detected and show high diversity and abundance in prokaryotic and eukaryotic organisms as well as in environmental samples worldwide. Currently, a new CRESS DNA virus similar to PCVs, with a circular genome but without the typical capsid protein (Cap), has been confirmed in diarrheic pigs. The epidemiology of PCLV in China needs further investigation. In this study, we performed an in-depth analysis and study of the genomic characteristics and genetic diversity of PCLV in Anhui Province, central China. PCLV infection in piglets has been confirmed in a previous report [13]. However, the results of this study suggest that PCLV may not only infect piglets, causing diarrhea, but may also be able to infect sows, causing severe diarrhea and mortality.
Rep is a relatively conserved gene in ssDNA viruses [2], which was verified in this study. Potential functional analysis of the ORFs revealed that a helicase domain (IPR027417, SSF52540, G3DSA:3.40.50.300) was found in ORF1 (Rep), whereas no domain was found in other ORFs. These domains are closely related to virus replication, and similar domains were also found in PCV [31], again confirming the function of ORF1 in encoding replicase. In addition, we found multiple amino acid mutations by complete genome sequence alignment. These mutations may be the underlying cause of these strains forming separate clades and may be key mutations that promote genotypic differentiation. Therefore, more serious consequences should be anticipated. In PCV, mutation in the Rep affect the virus’s ability to replicate, which can increase or decrease depending on the mutation site [32,33,34]. In some viruses, mutations also cause unpredictable effects, such as increased drug resistance, vaccine failure and changing abilities of replication and transmission [34,35,36,37]. Common mutational events in different branches of PCLV may indicate that the hazard of the virus to pigs may be more serious and unpredictable.
Currently, there is no standard and specific classification method for PCLV. In circoviruses, the Cap is considered to be very important for classification based on the phylogenetic tree and genetic distance [6,38,39]. However, Cap cannot be identified in PCLV [27], so the phylogenetic tree structure and genetic distances were investigated in this study. However, whether these three groups can be considered as three genotypes is currently lacking strong evidence, including a suitable genotyping basis and a large number of reference strains. At present, genotypic differentiation is known to be common in members of the genus Circovirus, such as PCV-2, PCV-3, duck circovirus and canine circovirus [40]. Given the close relationship between PCLV and PCV-2, more data are needed to prove whether the differentiation of three groups in PCLV also represents different pathogenicity and epidemic potential.
The topological differences between the tree based on Rep and the complete genome mentioned above indicate the occurrence of recombination events [41]. Rep is considered to be essential for viral replication and plays a crucial role in cellular immunity [6,42]. Recombination events are also widely found in other viruses, such as porcine circovirus, canine circovirus, pigeon circovirus, duck circovirus and goose circovirus [41,43,44,45,46]. The hotspot of recombination events is Rep, but it is still relatively conserved in ssDNA viruses, indicating that recombination may have little influence on the diversification of Rep. However, given recombination is a driving force in the evolution of viruses, we need to be concerned about the unpredictable effects of recombination on Rep.
Selection pressure can come from a wide range of sources, such as the effects of drugs, the host’s own immune system, or “competitors”. In this study, selection pressure on the virus promotes evolution. Unfortunately, we were unable to determine what brings about the selective effect of the virus, which is often not a single factor. The results of selection pressure and antigenic epitopes suggest that continuous changes in the immunogenicity of the virus may result in the diversity of PCLV and reinfection of domestic pigs in which antibodies have been acquired, causing further epidemics. Although Rep is not the main gene determining the immunogenicity of the virus, it is related to humoral immunity [42]. The same phenomenon was also found in PCV, especially PCV-2, which led to the worldwide epidemic and had a negative impact on the swine industry [47,48]. Codon usage bias is a pattern of codon use formed by viruses during evolution and influenced by mutational pressure, natural selection and other factors [25,49]. Low codon bias may promote efficient replication by reducing competition between virus and host during protein synthesis and may increase virulence. The virulence of PCLV may be underestimated currently or may increase in the future due to this low codon bias. Therefore, we should be more vigilant about the harm PCLV causes to the swine industry. In addition, the codon bias is always influenced by external factors to varying degrees. Natural selection has been shown to play a major role in PCV-2 and canine circovirus [50,51], but mutation is more important in duck circovirus codon bias [52]. Although these viruses are closely related, different factors have influenced their evolution. Therefore, as both natural selection and mutation pressure appear to influence the evolution of PCLV, further attention should be paid to these factors.
In conclusion, PCLV is capable of causing diarrhea not only in piglets but also in pregnant sows and can also cause death, as elucidated in this study. The pathogenicity of PCLV may be enhanced and revealed gradually with further evolution of amino acid mutations. However, further experiments on the construction of PCLV-based infectious clones are needed to elucidate the biological properties and pathogenic mechanisms of the mutation of this virus locus in vivo and in vitro.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13112282/s1, Figure S1: Pregnant sows developed diarrhea and hemorrhagic enteritis. Table S1: Primer sequences for detection and amplification of the whole genome of PCLV. Table S2: The primers for porcine virus detection without commercial kits. Table S3: The reference strain information used in this study. Table S4: The results of B-cell epitope prediction.

Author Contributions

K.Y., K.Q. and J.T. conceived and designed the study. M.Z., Q.L. and W.Z. conducted experiments. M.Z., W.Z. and Y.C. prepared experimental reagent materials. Y.S., K.J. and X.S. analyzed data. K.Y., Y.L. and Y.C. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported financially by the National Natural Science Foundation of China (no. 31972642) and the University Synergy Innovation Program of Anhui Province (No. GXXT-2019-035).

Institutional Review Board Statement

The study was approved by the Committee on the Ethics of Animal Care and Use at Anhui Agriculture University (Anhui, China) (No. AHAU2021008).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data set supporting the conclusions of this article is available in GenBank.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kazlauskas, D.; Varsani, A.; Koonin, E.V.; Krupovic, M. Multiple origins of prokaryotic and eukaryotic single-stranded DNA viruses from bacterial and archaeal plasmids. Nat. Commun. 2019, 10, 3425. [Google Scholar] [CrossRef] [Green Version]
  2. Zhao, L.; Rosario, K.; Breitbart, M.; Duffy, S. Eukaryotic circular rep-encoding single-stranded DNA (CRESS DNA) viruses: Ubiquitous viruses with small genomes and a diverse host range. Adv. Virus Res. 2018, 103, 71–133. [Google Scholar] [CrossRef] [PubMed]
  3. Saunders, K.; Bedford, I.D.; Yahara, T.; Stanley, J. The earliest recorded plant virus disease. Nature 2003, 422, 831. [Google Scholar] [CrossRef] [PubMed]
  4. Krupovic, M.; Varsani, A.; Kazlauskas, D.; Breitbart, M.; Delwart, E.; Rosario, K.; Yutin, N.; Wolf, Y.I.; Harrach, B.; Zerbini, F.M.; et al. Cressdnaviricota: A virus phylum unifying seven families of rep-encoding viruses with single-stranded, circular DNA genomes. J. Virol. 2020, 94, 12. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, Q.; Wang, H.; Ling, Y.; Yang, S.-X.; Wang, X.-C.; Zhou, R.; Xiao, Y.-Q.; Chen, X.; Yang, J.; Fu, W.-G.; et al. Viral metagenomics revealed diverse CRESS-DNA virus genomes in faeces of forest musk deer. Virol. J. 2020, 17, 61. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Rosario, K.; Breitbart, M.; Harrach, B.; Segalés, J.; Delwart, E.; Biagini, P.; Varsani, A. Revisiting the taxonomy of the family Circoviridae: Establishment of the genus Cyclovirus and removal of the genus Gyrovirus. Arch. Virol. 2017, 162, 1447–1463. [Google Scholar] [CrossRef] [Green Version]
  7. Varsani, A.; Krupovic, M. Correction to: Smacoviridae: A new family of animal-associated single-stranded DNA viruses. Arch. Virol. 2018, 163, 3213–3214. [Google Scholar] [CrossRef] [Green Version]
  8. Varsani, A.; Krupovic, M. Smacoviridae: A new family of animal-associated single-stranded DNA viruses. Arch. Virol. 2018, 163, 2005–2015. [Google Scholar] [CrossRef]
  9. Li, L.; Giannitti, F.; Ullmann, L.S.U.; Deng, X.; Pesavento, P.A.; Delwart, E.; Pusterla, N.; Keyes, C.; Low, J.; Aleman, M. Exploring the virome of diseased horses. J. Gen. Virol. 2015, 96, 2721–2733. [Google Scholar] [CrossRef] [Green Version]
  10. Rybicki, E. A phylogenetic and evolutionary justification for three genera of Geminiviridae. Arch. Virol. 1994, 139, 49–77. [Google Scholar] [CrossRef]
  11. Guo, Z.; He, Q.; Tang, C.; Zhang, B.; Yue, H. Identification and genomic characterization of a novel CRESS DNA virus from a calf with severe hemorrhagic enteritis in China. Virus Res. 2018, 255, 141–146. [Google Scholar] [CrossRef]
  12. Shan, T.; Li, L.; Simmonds, P.; Wang, C.; Moeser, A.; Delwart, E. The fecal virome of pigs on a high-density farm. J. Virol. 2011, 85, 11697–11708. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, X.; Zhang, X.; Xu, G.; Wang, Z.; Shen, H.; Lian, K.; Lin, Y.; Zheng, J.; Liang, P.; Zhang, L.; et al. Emergence of porcine circovirus-like viruses associated with porcine diarrheal disease in China. Transbound. Emerg. Dis. 2021. [Google Scholar] [CrossRef] [PubMed]
  14. Sun, W.; Wang, W.; Cao, L.; Zheng, M.; Zhuang, X.; Zhang, H.; Yu, N.; Tian, M.; Lu, H.; Jin, N. Genetic characterization of three porcine circovirus-like viruses in pigs with diarrhoea in China. Transbound. Emerg. Dis. 2020, 68, 289–295. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, D.; Gao, F.; Jakovlić, I.; Zhou, H.; Zhang, J.; Li, W.X.; Wang, G.T. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour. 2019, 20, 348–355. [Google Scholar] [CrossRef]
  16. Minh, B.Q.; Nguyen, M.A.T.; Von Haeseler, A. Ultrafast approximation for phylogenetic bootstrap. Mol. Biol. Evol. 2013, 30, 1188–1195. [Google Scholar] [CrossRef] [PubMed]
  17. Nguyen, L.-T.; Schmidt, H.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2014, 32, 268–274. [Google Scholar] [CrossRef]
  18. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Pond, S.; Frost, S.D.W. Not so different after all: A comparison of methods for detecting amino acid sites under selection. Mol. Biol. Evol. 2005, 22, 1208–1222. [Google Scholar] [CrossRef] [Green Version]
  20. Murrell, B.; Moola, S.; Mabona, A.; Weighill, T.; Sheward, D.; Kosakovsky Pond, S.L.; Scheffler, K. FUBAR: A fast, unconstrained bayesian approximation for inferring selection. Mol. Biol. Evol. 2013, 30, 1196–1205. [Google Scholar] [CrossRef] [Green Version]
  21. Smith, D.; Wertheim, J.O.; Weaver, S.; Murrell, B.; Scheffler, K.; Pond, S.L.K. Less is more: An adaptive branch-site random effects model for efficient detection of episodic diversifying selection. Mol. Biol. Evol. 2015, 32, 1342–1353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Yao, B.; Zhang, L.; Liang, S.; Zhang, C. SVMTriP: A method to predict antigenic epitopes using support vector machine to integrate tri-peptide similarity and propensity. PLoS ONE 2012, 7, e45152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sueoka, N. Intrastrand parity rules of DNA base composition and usage biases of synonymous codons. J. Mol. Evol. 1995, 40, 318–325. [Google Scholar] [CrossRef] [PubMed]
  24. Sueoka, N. Translation-coupled violation of parity rule 2 in human genes is not the cause of heterogeneity of the DNA G+C content of third codon position. Gene 1999, 238, 53–58. [Google Scholar] [CrossRef]
  25. Tian, H.-F.; Hu, Q.-M.; Xiao, H.-B.; Zeng, L.-B.; Meng, Y.; Li, Z. Genetic and codon usage bias analyses of major capsid protein gene in Ranavirus. Infect. Genet. Evol. 2020, 84, 104379. [Google Scholar] [CrossRef]
  26. Wang, L.Y.; Xing, H.X.; Yuan, Y.C.; Wang, X.L.; Saeed, M.; Tao, J.C.; Feng, W.; Zhang, G.H.; Song, X.L.; Sun, X.Z. Genome-wide analysis of codon usage bias in four sequenced cotton species. PLoS ONE 2018, 13, e0194372. [Google Scholar]
  27. Sun, J.; Zhao, W.; Wang, R.; Zhang, W.; Li, G.; Lu, M.; Shao, Y.; Yang, Y.; Wang, N.; Gao, Q.; et al. Analysis of the codon usage pattern of HA and NA genes of H7N9 influenza a virus. Int. J. Mol. Sci. 2020, 21, 7129. [Google Scholar] [CrossRef]
  28. Sibila, M.; Rocco, C.; Franzo, G.; Huerta, E.; Domingo, M.; Núñez, J.I.; Segalés, J. Genotyping of porcine circovirus 2 (PCV-2) in vaccinated pigs suffering from PCV-2-systemic disease between 2009 and 2020 in Spain. Pathogens 2021, 10, 1016. [Google Scholar] [CrossRef]
  29. Ghosh, S.; Chakraborty, S. Phylogenomics analysis of SARS-CoV2 Genomes reveals distinct selection pressure on different viral strains. Biomed Res. Int. 2020, 2020, 5746461. [Google Scholar] [CrossRef]
  30. Wang, Y.; Guo, X.; Li, W.; Cui, Y.Q.; Zhang, D.; Xu, F.Z.; Jiang, S.D.; Zhou, T.H. Phylogenetic analysis and evolution of feline bocavirus in Anhui Province, eastern China. Comp. Immunol. Microb. 2021, 77, 101676. [Google Scholar] [CrossRef]
  31. Nguyen, V.; Do, H.; Huynh, T.; Park, Y.; Park, B.; Chung, H. Molecular-based detection, genetic characterization and phylogenetic analysis of porcine circovirus 4 from Korean domestic swine farms. Transbound. Emerg. Dis. 2021. [Google Scholar] [CrossRef] [PubMed]
  32. Wu, X.; Wang, S.; Xin, C.; Li, C.; Shi, J.; Peng, Z.; Liu, C.; Han, H.; Yang, Y.; Tian, Y.; et al. Mutation of the sixth amino acid of the Rep protein has no effect on porcine circovirus 2b but enhances porcine circovirus 2d replication in vitro. Arch. Virol. 2021, 166, 3189–3192. [Google Scholar] [CrossRef] [PubMed]
  33. Hu, Y.; Cai, X.; Zhan, Y.; Yuan, X.; Liu, T.; Tan, L.; Li, Y.; Zhang, L.; Yang, L.; Liu, W.; et al. Truncated rep protein of porcine circovirus 2 (PCV2) caused by a naturally occurring mutation reduced virus replication in PK15 cells. BMC Vet. Res. 2019, 15, 248. [Google Scholar] [CrossRef] [Green Version]
  34. Shi, J.; Peng, Z.; Fu, F.; Xu, S.; Xu, S.; Cong, X.; Yuan, X.; Yu, J.; Wu, J.; Sun, W.; et al. Mutant Rep protein of the porcine circovirus type 2 N-glycosylation: 23-25aa, 256-258aa mutation reduced virus replication but 286-288aa mutation enhanced virus replication in PK-15 cells. Vet. Microbiol. 2015, 177, 370–372. [Google Scholar] [CrossRef] [PubMed]
  35. Ji, Y.; Liu, T.; Du, Y.; Cui, X.; Yu, Q.; Wang, Z.; Zhang, J.; Li, Y.; Zhu, Q. A novel genotype VII Newcastle disease virus vaccine candidate generated by mutation in the L and F genes confers improved protection in chickens. Veter.-Microbiol. 2018, 216, 99–106. [Google Scholar] [CrossRef]
  36. Pang, H.; Li, L.; Liu, H.; Pan, Z. Proline to threonine mutation at position 162 of NS5B of classical swine fever virus vaccine C strain promoted genome replication and infectious virus production by facilitating initiation of RNA synthesis. Viruses 2021, 13, 1523. [Google Scholar] [CrossRef]
  37. Wang, L.; Han, F.; Duan, H.; Ji, F.; Yan, X.; Fan, Y.; Wang, K. Hepatitis B virus pre-existing drug resistant mutation is related to the genotype and disease progression. J. Infect. Dev. Ctries. 2017, 11, 727–732. [Google Scholar] [CrossRef] [Green Version]
  38. Franzo, G.; Delwart, E.; Fux, R.; Hause, B.; Su, S.; Zhou, J.Y.; Segales, J. Genotyping porcine circovirus 3 (PCV-3) nowadays: Does it make sense? Viruses 2020, 12, 265. [Google Scholar] [CrossRef] [Green Version]
  39. Franzo, G.; Segalés, J. Porcine circovirus 2 (PCV-2) genotype update and proposal of a new genotyping methodology. PLoS ONE 2018, 13, e0208585. [Google Scholar] [CrossRef] [Green Version]
  40. Saporiti, V.; Cruz, T.F.; Correa-Fiz, F.; Núñez, J.I.; Sibila, M.; Segalés, J. Similar frequency of Porcine circovirus 3 (PCV-3) detection in serum samples of pigs affected by digestive or respiratory disorders and age-matched clinically healthy pigs. Transbound. Emerg. Dis. 2019, 67, 199–205. [Google Scholar] [CrossRef]
  41. Ji, J.; Chen, Q.; Sui, C.; Yu, Z.; Xu, X.; Yao, L.; Kan, Y.; Bi, Y.; Xie, Q. Novel genotype definition and genome characteristics of duck circovirus in central and Eastern China. Transbound. Emerg. Dis. 2020, 67, 2993–3004. [Google Scholar] [CrossRef]
  42. Wu, X.; Wang, X.; Shi, T.; Luo, L.; Qiao, D.; Wang, Z.; Han, C.; Du, Q.; Tong, D.; Huang, Y. Porcine circovirus type 2 rep enhances IL-10 production in macrophages via activation of p38-MAPK pathway. Viruses 2019, 11, 1141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Khalifeh, A.; Kraberger, S.; Dziewulska, D.; Varsani, A.; Stenzel, T. A pilot study investigating the dynamics of pigeon circovirus recombination in domesticated pigeons housed in a single loft. Viruses 2021, 13, 964. [Google Scholar] [CrossRef]
  44. Piewbang, C.; Jo, W.K.; Puff, C.; Van Der Vries, E.; Kesdangsakonwut, S.; Rungsipipat, A.; Kruppa, J.; Jung, K.; Baumgärtner, W.; Techangamsuwan, S.; et al. Novel canine circovirus strains from Thailand: Evidence for genetic recombination. Sci. Rep. 2018, 8, 7524. [Google Scholar] [CrossRef]
  45. Rajkhowa, T.K.; Lalnunthanga, P.; Rao, P.L.; Subbiah, M.; Lalrohlua, B. Emergence of porcine circovirus 2g (PCV2g) and evidence for recombination between genotypes 2g, 2b and 2d among field isolates from non-vaccinated pigs in Mizoram, India. Infect. Genet. Evol. 2021, 90, 104775. [Google Scholar] [CrossRef]
  46. Stenzel, T.; Dziewulska, D.; Muhire, B.M.; Hartnady, P.; Kraberger, S.; Martin, D.P.; Varsani, A. Recombinant goose circoviruses circulating in domesticated and wild geese in Poland. Viruses 2018, 10, 107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Li, G.; He, W.; Zhu, H.; Bi, Y.; Wang, R.; Xing, G.; Zhang, C.; Zhou, J.; Yuen, K.-Y.; Gao, G.F.; et al. Origin, genetic diversity, and evolutionary dynamics of novel porcine circovirus 3. Adv. Sci. 2018, 5, 1800275. [Google Scholar] [CrossRef] [PubMed]
  48. Yang, S.; Yin, S.; Shang, Y.; Liu, B.; Yuan, L.; Khan, M.U.Z.; Liu, X.; Cai, J. Phylogenetic and genetic variation analyses of porcine circovirus type 2 isolated from China. Transbound. Emerg. Dis. 2017, 65, e383–e392. [Google Scholar] [CrossRef]
  49. Yao, X.; Fan, Q.; Yao, B.; Lu, P.; Rahman, S.U.; Chen, D.; Tao, S. Codon usage bias analysis of bluetongue virus causing livestock infection. Front. Microbiol. 2020, 11, 655. [Google Scholar] [CrossRef]
  50. Wang, L.; Li, Y.; Guo, Z.; Yi, Y.; Zhang, H.; Shangguan, H.; Huang, C.; Ge, J. Genetic changes and evolutionary analysis of canine circovirus. Arch. Virol. 2021, 166, 2235–2247. [Google Scholar] [CrossRef] [PubMed]
  51. Yu, X.; Gao, K.; Pi, M.; Li, H.; Zhong, W.; Li, B.; Ning, Z. Phylogenetic and codon usage analysis for replicase and capsid genes of porcine circovirus 3. Veter.-Res. Commun. 2021, 9, 1–9. [Google Scholar] [CrossRef] [PubMed]
  52. Xu, Y.; Jia, R.Y.; Zhang, Z.L.; Lu, Y.Y.; Wang, M.S.; Zhu, D.K.; Chen, S.; Liu, M.F.; Yin, Z.Q.; Cheng, A.C. Analysis of syn-onymous codon usage pattern in duck circovirus. Gene 2015, 557, 138–145. [Google Scholar] [CrossRef] [PubMed]
Viruses 13 02282 g001 550
Figure 1. Geographic distribution of positive samples.
Figure 1. Geographic distribution of positive samples.
Viruses 13 02282 g001
Viruses 13 02282 g002 550
Figure 2. Predicted genomic structural characteristics of PCLV-AH-23, PCLV-AH-25 and PCLV-AH-HB-2021. PCLV-AH-25 only four potential ORFs were present, whereas in strain PCLV-AH-23 and PCLV-AH-HB-2021, the locations of five potential ORFs were identified,, one of which contained the opposite transcriptional direction. A significant T7G mutation was found in the stem loop of PCLV-AH-HB-2021.
Figure 2. Predicted genomic structural characteristics of PCLV-AH-23, PCLV-AH-25 and PCLV-AH-HB-2021. PCLV-AH-25 only four potential ORFs were present, whereas in strain PCLV-AH-23 and PCLV-AH-HB-2021, the locations of five potential ORFs were identified,, one of which contained the opposite transcriptional direction. A significant T7G mutation was found in the stem loop of PCLV-AH-HB-2021.
Viruses 13 02282 g002
Viruses 13 02282 g003 550
Figure 3. Amino acid mutation analysis of the PCLV Rep. Strains obtained in this study are indicated in pink. The conserved regions of the amino acid sequences are covered in red.
Figure 3. Amino acid mutation analysis of the PCLV Rep. Strains obtained in this study are indicated in pink. The conserved regions of the amino acid sequences are covered in red.
Viruses 13 02282 g003
Viruses 13 02282 g004 550
Figure 4. Phylogenetic analysis based on complete sequences and Rep. Strains from this study are marked in red. (A) Rep: The best substitution model was TIM2 + F + G4. Branches marked in yellow are those with recombination events, and Mutation Groups 1 and 2 are strains with mutations. (B) Complete sequences: The best substitution model was TPM3 + F + G4. The tree is divided into three branches, each underlined with a different color.
Figure 4. Phylogenetic analysis based on complete sequences and Rep. Strains from this study are marked in red. (A) Rep: The best substitution model was TIM2 + F + G4. Branches marked in yellow are those with recombination events, and Mutation Groups 1 and 2 are strains with mutations. (B) Complete sequences: The best substitution model was TPM3 + F + G4. The tree is divided into three branches, each underlined with a different color.
Viruses 13 02282 g004
Viruses 13 02282 g005 550
Figure 5. Genetic distance analysis between PCLV and PCV2. The genotypes of PCV refer to the classification of Marina Sibila et al. [28], including 2a, 2b, 2c, 2d, 2e, 2f, 2g and 2h. Genetic distance was calculated using MEGA X with a p-distance model and 1000 bootstrap replicates.
Figure 5. Genetic distance analysis between PCLV and PCV2. The genotypes of PCV refer to the classification of Marina Sibila et al. [28], including 2a, 2b, 2c, 2d, 2e, 2f, 2g and 2h. Genetic distance was calculated using MEGA X with a p-distance model and 1000 bootstrap replicates.
Viruses 13 02282 g005
Viruses 13 02282 g006 550
Figure 6. Recombination analysis of PCLV. PCLV-AH-23 was the query sequence, and the major and minor parents were strain PCLV-AH-25 (accession number: MZ773068/China) and Po-Circo-like virus 21 (accession number: JF713716.1/USA). The red curve represents the similarity of strain Po-Circo-like virus 21 to strain PCLV-AH-23. The green curve represents the similarity of strain PCLV-AH-25 to strain PCLV-AH-23. Rep is the potential recombination region. Recombination events were identified using seven methods in the Recombination Detection Program v.4.39 (RDP 4.39), namely, RDP, GENECONV, Chimaera, MaxChi, BootScan, SiScan and 3Seq. The SimPlot software v.3.5.1 was used for the analysis with visual similarity plots using a window size of 200 nucleotides that was moved along in 20 nt steps.
Figure 6. Recombination analysis of PCLV. PCLV-AH-23 was the query sequence, and the major and minor parents were strain PCLV-AH-25 (accession number: MZ773068/China) and Po-Circo-like virus 21 (accession number: JF713716.1/USA). The red curve represents the similarity of strain Po-Circo-like virus 21 to strain PCLV-AH-23. The green curve represents the similarity of strain PCLV-AH-25 to strain PCLV-AH-23. Rep is the potential recombination region. Recombination events were identified using seven methods in the Recombination Detection Program v.4.39 (RDP 4.39), namely, RDP, GENECONV, Chimaera, MaxChi, BootScan, SiScan and 3Seq. The SimPlot software v.3.5.1 was used for the analysis with visual similarity plots using a window size of 200 nucleotides that was moved along in 20 nt steps.
Viruses 13 02282 g006
Viruses 13 02282 g007 550
Figure 7. Map of B-cell epitope prediction and pressure selection sites of Rep. The predicted epitopes in the Rep protein are indicated in blue. The selective pressure sites are indicated in red. Protein structure visualization was performed using PyMol.
Figure 7. Map of B-cell epitope prediction and pressure selection sites of Rep. The predicted epitopes in the Rep protein are indicated in blue. The selective pressure sites are indicated in red. Protein structure visualization was performed using PyMol.
Viruses 13 02282 g007
Viruses 13 02282 g008 550
Figure 8. Analysis of factors mediating codon usage bias. (A) Parity Rule 2 (PR2) plot was used to analyze whether the effects of mutation pressure and natural selection on the use of the Rep codon were consistent. Figure 8A shows that mutation pressure and natural selection have unequal effects on the use of codons during evolution. (B) The ENC-GC3s plot was analyzed, with the ENC values plotted against the GC3S values to determine the main factors affecting codon use bias. The expected ENC values for each GC3S were estimated. Figure 8B shows that main factor affecting codon use bias was selection pressure. (C) Neutrality analysis was used to evaluate the effects of mutation pressure and natural selection on codon usage patterns. Figure 8C shows that the fitting curve (red line) deviates from the curve y = x (black line), indicating that selection pressure plays a dominant role, which is consistent with Figure 8B.
Figure 8. Analysis of factors mediating codon usage bias. (A) Parity Rule 2 (PR2) plot was used to analyze whether the effects of mutation pressure and natural selection on the use of the Rep codon were consistent. Figure 8A shows that mutation pressure and natural selection have unequal effects on the use of codons during evolution. (B) The ENC-GC3s plot was analyzed, with the ENC values plotted against the GC3S values to determine the main factors affecting codon use bias. The expected ENC values for each GC3S were estimated. Figure 8B shows that main factor affecting codon use bias was selection pressure. (C) Neutrality analysis was used to evaluate the effects of mutation pressure and natural selection on codon usage patterns. Figure 8C shows that the fitting curve (red line) deviates from the curve y = x (black line), indicating that selection pressure plays a dominant role, which is consistent with Figure 8B.
Viruses 13 02282 g008
Table
Table 1. Positive selection site for three genotypes of PCLVs.
Table 1. Positive selection site for three genotypes of PCLVs.
MethodsPositive Selection SitesThreshold
FUBAR29, 142p < 0.1
FEL29
MEME47, 81, 142, 211, 290, 8, 309, 72, 29
SLAC/Posterior probability of 0.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, K.; Zhang, M.; Liu, Q.; Cao, Y.; Zhang, W.; Liang, Y.; Song, X.; Ji, K.; Shao, Y.; Qi, K.; et al. Epidemiology and Evolution of Emerging Porcine Circovirus-like Viruses in Pigs with Hemorrhagic Dysentery and Diarrhea Symptoms in Central China from 2018 to 2021. Viruses 2021, 13, 2282. https://doi.org/10.3390/v13112282

AMA Style

Yang K, Zhang M, Liu Q, Cao Y, Zhang W, Liang Y, Song X, Ji K, Shao Y, Qi K, et al. Epidemiology and Evolution of Emerging Porcine Circovirus-like Viruses in Pigs with Hemorrhagic Dysentery and Diarrhea Symptoms in Central China from 2018 to 2021. Viruses. 2021; 13(11):2282. https://doi.org/10.3390/v13112282

Chicago/Turabian Style

Yang, Kankan, Menghuan Zhang, Qi Liu, Yingli Cao, Wuyin Zhang, Yueqiao Liang, Xiangjun Song, Kaiyuan Ji, Ying Shao, Kezong Qi, and et al. 2021. "Epidemiology and Evolution of Emerging Porcine Circovirus-like Viruses in Pigs with Hemorrhagic Dysentery and Diarrhea Symptoms in Central China from 2018 to 2021" Viruses 13, no. 11: 2282. https://doi.org/10.3390/v13112282

APA Style

Yang, K., Zhang, M., Liu, Q., Cao, Y., Zhang, W., Liang, Y., Song, X., Ji, K., Shao, Y., Qi, K., & Tu, J. (2021). Epidemiology and Evolution of Emerging Porcine Circovirus-like Viruses in Pigs with Hemorrhagic Dysentery and Diarrhea Symptoms in Central China from 2018 to 2021. Viruses, 13(11), 2282. https://doi.org/10.3390/v13112282

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