Skip to main content
  • Research article
  • Open access
  • Published:

Systematic identification and characterization of repeat sequences in African swine fever virus genomes

Abstract

African swine fever virus (ASFV) is a large DNA virus that infects domestic pigs with high morbidity and mortality rates. Repeat sequences, which are DNA sequence elements that are repeated more than twice in the genome, play an important role in the ASFV genome. The majority of repeat sequences, however, have not been identified and characterized in a systematic manner. In this study, three types of repeat sequences, including microsatellites, minisatellites and short interspersed nuclear elements (SINEs), were identified in the ASFV genome, and their distribution, structure, function, and evolutionary history were investigated. Most repeat sequences were observed in noncoding regions and at the 5’ end of the genome. Noncoding repeat sequences tended to form enhancers, whereas coding repeat sequences had a lower ratio of alpha-helix and beta-sheet and a higher ratio of loop structure and surface amino acids than nonrepeat sequences. In addition, the repeat sequences tended to encode penetrating and antimicrobial peptides. Further analysis of the evolution of repeat sequences revealed that the pan-repeat sequences presented an open state, showing the diversity of repeat sequences. Finally, CpG islands were observed to be negatively correlated with repeat sequence occurrences, suggesting that they may affect the generation of repeat sequences. Overall, this study emphasizes the importance of repeat sequences in ASFVs, and these results can aid in understanding the virus's function and evolution.

Introduction

African swine fever (ASF) is a haemorrhagic and devastating infectious disease of pigs caused by the African swine fever virus (ASFV), with a mortality rate of up to 100% [1, 2]. ASFV has a significant impact on the pig industry in many countries [3]. Over fifty countries in Africa, Europe, and Asia have experienced ASF outbreaks [3, 4]. ASF first appeared in China in August 2018 and spread quickly throughout nearly all of the country's provinces [5]. The viral genome is double-stranded DNA, with sizes ranging from 170 to 194 kb [6]. The ASFV genome encodes more than 150 open reading frames (ORFs), which form a mature ASFV virion with a large, enveloped, and complicated architecture, making the development of an effective vaccine and drug difficult due to an incomplete understanding of the virus [6,7,8].

Repeat sequences are DNA sequences that appear repeatedly in the genome and are crucial for viral infection and variation [9,10,11]. According to the current findings, repeat sequences are widely distributed in ASFV genomes and may promote the occurrence of homologous recombination [12]. A few repeat sequences have been found to play an important role in ASFV infection in host cells [11, 13]. For example, the repeat peptide ([KPCPPP]3) acts as a cell-penetrating peptide, potentially assisting the virus in entering the cell [11]. In addition, variation in the number of repeat sequences within genes or intergenic regions contributes to the genetic diversity of ASFV [12, 14, 15]. Although the functions of a small number of repeat sequences in the virus have been identified, the structure, function, and evolution of the vast majority of repeat sequences have not been studied [13, 15]. A systematic study of ASFV repeat sequences is required to bridge this gap, which will aid in virus prevention and control.

Tandem repeats and SINEs are the most common types of repeat sequences [16, 17]. They not only have an impact on evolution, inheritance, and variation but also play an important role in gene expression, transcription, and regulation [18, 19]. For example, tandem repeats at the ends of human chromosomes play an important regulatory role in cell growth, proliferation, and apoptosis processes, such as maintaining chromosome integrity and stability, preventing chromosome degradation by nucleases, avoiding chromosome end fusion and recombination, providing telomerase substrates, and regulating cell lifecycles [20, 21]. SINEs are found in many vertebrate and invertebrate species and have lineage specificity, which plays a particularly important role in the regulation of gene expression and the creation of RNA genes [22, 23]. For example, SINEs containing CTCF binding sites can be used as boundary elements to influence chromatin structure and transcription [22, 24]. In this study, the ASFV repeat sequences were systematically identified and classified, and their structure, function, and evolutionary characteristics were further investigated. The findings demonstrated that repeat sequences are essential for ASFV expression and evolution. This work will contribute to a better understanding of ASFV function and evolution and pave the way for future research on repeat sequences in other viruses.

Materials and methods

Identification of repeat sequences in ASFV genomes

A total of 86 ASFV genome sequences were obtained from the NCBI GenBank [25] database on August 15, 2021. The gaps and nonstandard bases were removed from the ASFV genomes. REPuter (version 1) [26] and TRF (version 1) [27] were used to identify repeat sequences in ASFV genomes with the default parameters.

Annotation of ASFV genomes

Genes encoded in ASFV genomes and their relevant proteins were predicted using GeneMarkS (version 4.28) [28] with the default parameters. To infer the function of the ASFV proteins, they were queried against all ASFV proteins downloaded from the NCBI RefSeq and ASFVdb [29] databases using the NCBI BLAST program (version 2.9.1) [30].

Structure‒function analysis of repeat sequences in ASFV genomes

As shown in Table 1, the structure and function of repeat sequences were investigated using a variety of publicly available tools.

Table 1 Software tools for predicting the structure and function of repeat sequences

Inference of phylogenetic trees and determination of genotypes

Previously published research identified 24 genotypes based on the C-terminus of the p72 gene [37]. To determine the genotype of all analysed ASFV strains, we used MAFFT (version 7) [38] to align the genomic sequences in the 415-bp C-terminal region of the p72 gene; then, we inferred phylogenies using the maximum-likelihood algorithm in MEGAX [39] with 100 bootstrap replicates. The genotype of each ASFV was assigned based on previous studies [37, 40].

Pan-genomic analysis of the ASFV repeat sequence

Repeat sequences were clustered by treating identical repeat sequences as a cluster. The characteristic curves of the ASFV pan-repeat sequences and the core-repeat sequences were portrayed using PanGP with DG sampling algorithms [41]. The flower plot depicting the number of core and dispensable repeat sequence clusters, including unique repeat sequences in ASFV, was generated by the plotrix() function in R (version 4.0.3) [42].

Identification of CpG islands

The CpG islands in the ASFV genomes were detected using Cpgplot (version 1) [43].

Statistical analysis

All statistical analyses were conducted in R (version 4.0.3) [42] and Python (version 3.6) [44]. The Wilcoxon rank-sum test was conducted using the wilcox.test() function in R. The correlation coefficient was calculated using the cor.test() function in R.

Results

Identification, classification and distribution of repeat sequences in ASFV genomes

A total of 86 genome sequences of ASFVs were obtained from the NCBI GenBank database; these are listed in Additional file 1. The ASFVs were found to contain an abundance of repeat sequences (Figure 1). These repeat sequences were divided into three types: (1) microsatellites, which are defined as small tandem repeat sequences with less than ten bp in the repeat unit, (2) minisatellites, which are defined as small tandem repeat sequences with 10–300 bp in the repeat unit, and (3) short interspersed nuclear elements (SINEs), which are defined as a discontinuous distribution of repeat units of 100–500 bp. The analysis of the location of repeat sequences in ASFV genomes revealed that the bulk of repeat sequences were located in either the 5′ or 3′ ends of genomes (Figure 1A). Specifically, more than 90% of microsatellites and SINEs were located at the 5′ end of the ASFV genome, while minisatellites were irregularly distributed throughout the genome.

Figure 1
figure 1

The number, position, and distribution of repeat sequences in the genome. A The number of repeat sequences identified in 86 genomes (per kb). B The phylogenetic tree on the left side was built with MEGA based on the C-terminal sequences (415 bp) of the p72 gene (see Materials and methods), and the genotypes of ASFV isolates are labelled with Greek numerals. The three types of repeat sequences in ASFV genomes are indicated by different colours. The black and red arrows refer to examples of genotype-specific repeat sequences in genotypes II and I, respectively.

The vast majority of the repeat sequences were found to be genotype specific. The retrieved genomes contained nine genotypes. Most of the ASFV genomes, including those recently found in Europe and China, were genotype I or II. Genotype I included 26 ASFV genomes, with more than 30% of repeat sequences being genotype specific (see examples marked by the red arrow in Figure 1B). Genotype II encompassed 38 ASFV genomes, with more than half of the repeat sequences being genotype specific (see examples marked by the black arrow in Figure 1B).

The number of repeat sequences and ratio of repeat sequences in 86 ASFV genomes were analysed, and the ASFVs were identified to have 2.2 to 5.3 repeat sequences per 10 kb of genomic sequence, with a median of 3.7 (Figure 2A). More than 60% of these repeat sequences were minisatellites, with the number ranging from 1.9 to 4.9 per 10 kb of genomic sequence. The minisatellites covered a median of 1.8% of the ASFV genomes (Figure 2B). Approximately 30% of repeat sequences were made up of SINEs, with numbers ranging from 0.05 to 1.9 per 10 kb of genomic sequence, covering a median of 0.9% of the ASFV genomes. The remaining repeat sequences belonged to microsatellites, which had fewer than 0.5 per 10 kb of genomic sequence.

Figure 2
figure 2

The distribution of three types of repeat sequences in the ASFV genome. A The number of different types of repeat sequences per 10 kb in 86 ASFV genomes. B The ratio of genomic regions covered by three types of repeat sequences in 86 ASFV genomes. C The number of repeat units for the three repeat sequence types in 86 ASFV genomes. D The number of repeat units of different lengths. The microsatellites, minisatellites, and SINEs are coloured red, green, and blue, respectively.

As the repeat units served as the basic units of repeat sequences, the number of repeat units was analysed in ASFV genomes (Figure 2C). The number of repeat units ranged from 2 to 66, with a median of 3. In approximately 80% of the repeat sequences, the number of repeat units was less than four times; in only 2% of the repeat sequences, the number of repeat units was more than ten times. The repeat unit was only repeated two times for nearly all SINEs and three times for more than 40% of minisatellites.

The length of repeat units in repeat sequences was then analysed (Figure 2D). The repeat unit lengths varied from 1 to 281 bp, with a median of 24 bp. Among them, more than half of the repeat units were no longer than 30 bp, and approximately 9% of the repeat units were more than 50 bp. Interestingly, the most common repeat unit length was 10–20 bp.

A high incidence of repeat sequences in noncoding regions

We evaluated the ratios of repeat sequences in the noncoding and coding regions and discovered that approximately 50% of repeat sequences were located in noncoding regions (Figure 3A). Specifically, approximately 70% of minisatellites and 95% of microsatellites were located in noncoding regions, while less than 20% of SINEs were located in noncoding regions. The analysis of the SINEs in coding regions revealed that 70% of them were located in the multigene family (MGF) proteins, with the remaining SINEs found in other proteins and unknown proteins (Additional file 2). This indicated that SINEs were likely to cause amino acid insertions or deletions, thereby influencing the functions of the ASFV proteins.

Figure 3
figure 3

The function of repeat sequences in noncoding regions. A The distribution of the noncoding ratio for three types of repeat sequences in the ASFV genome. B Comparison of the ratio of enhancers in the repeat and nonrepeat sequences in ASFV genomes. C Comparison of the ratio of promoters in repeat and nonrepeat sequences in ASFV genomes. “**”, p value < 0.01.

Characterization of the structure and function of repeat sequences in noncoding regions

The role of repeat sequences in transcription regulation was investigated because most repeat sequences were located in noncoding regions (Figure 3A). In the repeat sequences located in noncoding regions, two kinds of functional elements were identified: enhancers and promoters. For comparison, the same number of nonrepeat sequences with the same size as repeat sequences (control group) were randomly selected from the noncoding regions in each ASFV genome. The ratio of enhancers in repeat sequences was significantly higher than that in nonrepeat sequences (0.08 vs. 0.04, p value < 0.01) (Figure 3B). Surprisingly, the median ratio of promoters in both repeat sequences and the control group was greater than 0.9 (Figure 3C), and no significant difference was observed between them.

Characterization of the sequence, structure and function of repeat sequences in coding regions

Then, the repeat sequences in coding regions were analysed. For each repeat sequence, the protein sequence encoded by the repeat sequence (defined as repeat protein sequences for clarity) and other protein sequences in the same protein (defined as nonrepeat protein sequences) were used in the analysis. The amino acid ratios in repeat protein sequences were calculated and compared to those in the nonrepeat protein sequences. Interestingly, large differences in amino acid ratios were observed in the two types of protein sequences (Figure 4A). The ratios of several hydrophobic amino acids, including phenylalanine (F), methionine (M), leucine (L), and isoleucine (I), were significantly lower in the repeat protein sequences than in the nonrepeat protein sequences, while the ratios of three hydrophilic amino acids, including cysteine (C) and threonine (T), were significantly higher in the repeat protein sequences than in the nonrepeat protein sequences (p values < 0.01).

Figure 4
figure 4

The characteristics of repeat sequences in coding regions. A The ratio of amino acids in repeat protein sequences versus nonrepeat protein sequences in ASFV strains. B The ratio of secondary structures in repeat protein sequences versus nonrepeat protein sequences. C The ratio of amino acids exposed in repeat protein sequences versus nonrepeat protein sequences. An exposed amino acid was defined as a residue with a relative solvent accessible surface area (rASA) greater than 0.25. DF The ratio of disordered, penetrating peptides and antimicrobial peptides in repeat protein sequences versus nonrepeat protein sequences. For subfigures B to D, only nonredundant protein sequences were used. “*”, p value < 0.05; “**”, p value < 0.01.

We further investigated the structural characteristics of the repeat protein sequences. The ratios of three kinds of secondary structure elements were calculated (alpha-helix, beta-sheet, and loop). The median ratio of alpha-helix in the repeat protein sequences was 0.35, which was significantly lower than that in nonrepeat protein sequences (0.35 vs. 0.43, p value < 0.01) (Figure 4B), and the median ratio of beta-sheet in the repeat protein sequences was only 0.0, which was significantly lower than that in nonrepeat protein sequences (0.0 vs. 0.06, p value < 0.01). Interestingly, the median ratio of loops in the repeat protein sequences was 0.52, which was significantly higher than that in nonrepeat protein sequences (0.52 vs. 0.45, p value < 0.05). Then, the ratio of amino acids exposed was calculated, and the repeat protein sequences were observed to have a higher ratio of amino acids exposed than the nonrepeat protein sequences (Figure 4C). In addition, the ratio of amino acids in the disordered region, which is defined as the region unlikely to form a defined protein three-dimensional structure, was calculated, and the repeat protein sequences were observed to have a lower ratio of amino acids in the disordered region than the nonrepeat protein sequences (Figure 4D).

The potential functions of the repeat protein sequences were further analysed. Previous studies have shown that repeat protein sequences can act as penetrating peptides or antimicrobial peptides in herpes simplex and human respiratory syncytial viruses [11, 45]. Thus, penetrating and antimicrobial peptides were identified in repeat protein sequences. A median of 18% and 21% of repeat protein sequences were observed to have the function of penetrating peptides and antimicrobial peptides, respectively (Figures 4E and F), which were significantly higher than those in nonrepeat protein sequences. Taken together, more than 30% of the repeat protein sequences functioned as antibacterial peptides or penetrating peptides (Additional file 3), suggesting that the repeat protein sequences may play an important role in ASFV infection.

Pan-genomic analysis of repeat sequences in the ASFV genome

The pan-genomic analysis was conducted for the repeat sequences in ASFV genomes. All the repeat sequences identified in the ASFV genomes were defined as pan-repeat sequences. They were classified into 1141 clusters based on sequence similarity at the 100% level. Two kinds of clusters were obtained: core repeat clusters, which appeared in all ASFV genomes, and dispensable repeat clusters, which appeared in one or more ASFV genomes. Among the dispensable repeat clusters, those that only appeared in one genome were called unique repeat sequences. The analysis of pan-repeat sequences showed that no core repeat clusters were found, despite the presence of 1141 dispensable repeat clusters, including 748 unique repeat sequences (Figure 5A).

Figure 5
figure 5

Pan-repeat sequence analysis of ASFVs. A Flower plot showing the core repeat and dispensable repeat sequence clusters of the 86 ASFV strains. The diagram depicts the core repeat sequence cluster number (orange) and the dispensable repeat sequence cluster number (blue) for the 86 ASFV strains. The number of unique repeat sequences in each strain is shown beside the accession number of the ASFV genomes. For clarity, only the strains with unique repeat sequences are shown. B The relationship between genome number and dispensable repeat and core repeat sequence profiles. The black and blue curves refer to the least-squares fitting for the average number of dispensable repeat clusters and core repeat clusters versus the number of genomes, respectively. The mathematical functions and the R-squared values of both models are delineated on the graph.

With the increase in the number of genomes, the number of dispensable repeat clusters was fitted by the power law regression model (ypan = Apan*xBpan + Cpan), while the number of core repeat clusters was fitted by the exponential model (ycore = Acore*eBcore*x + Ccore), according to Zhao’s study [41]. The mathematical functions of both models were delineated on the graph. The Bpan in the power law regression model was greater than 0.5, suggesting that the pan-repeat sequences of ASFVs are in an open state (see the black curve in Figure 5B). This trend reflects that ASFVs have flexible repeat sequences, and the size of the pan-repeat may expand with each added genome, which contributes to new repeat sequences. As the number of analysed genomes increased, the number of core repeat sequence clusters presented a significant decrease and became zero (see the blue curve in Figure 5B). This trend reflects the diversity of repeat sequences in the ASFV genome.

Considering the genotype specificity of repeat sequences, we further analysed the pan-repeat sequences of genotypes I and II. Overall, 5 and 32 core repeat clusters and 251 and 73 dispensable repeat clusters, including 138 and 22 unique repeat sequences, respectively, were observed, which constituted the ASFV genotype I and II pan-repeats (Additional file 4). The Bpan in the power law regression model of genotypes I and II was greater than 0.5 (Additional file 4), suggesting that the pan-repeat sequences within the genotype are also in an open state. As the number of analysed genomes increased, the core-repeat cluster curve presented a converging trend. This trend reflects that repeat sequences within genotypes are conserved during evolution.

The relationship between CpG islands and repeat sequences in ASFV genomes

Previous studies have shown that CpG islands help stabilize the genome and may prevent the formation of repeat sequences. Thus, the influence of CpG islands on repeat sequences was investigated. The number of CpG islands identified in ASFV genomes ranged from 1.2 to 2.1 per 10 kb of genomic sequence (Figure 6A), with a median of 1.57. Interestingly, the number of CpG islands identified per 10 kb window of ASFV genomes was observed to have a negative correlation with the number of repeat sequences identified in the same region (Pearson correlation coefficient = −0.48) (Figure 6B). Specifically, the windows with CpG islands had a median of only one repeat sequence, while those without CpG islands had a median of five repeat sequences (Figure 6C).

Figure 6
figure 6

Distribution of the CpG islands and repeat sequences in each genome. A The number of CpG islands per 10 kb in each ASFV genome. B The number of repeat sequences in a sliding window of 10 kb (with a step size of 1 kb) versus the number of CpG islands. C Comparison of the number of CpG islands in a sliding window of 10 kb (with a step size of 1 kb) with and without CpG islands in ASFV genomes. “**”, p value < 0.01.

Discussion

The repeat sequences in the ASFV genome were studied for their classification, distribution, structure, function, and evolution. In this study, four features of ASFV repeat sequences were revealed for the first time. First, the repeat sequences were enriched at the 5’ end of the ASFV genome. This is consistent with Campbell et al.’s study that discovered a large number of simple sequence repeats (SSRs) within the terminal region of the Fowlpox virus genome [46]. The repeat sequences tended to be located in noncoding regions because insertions and deletions of repeat sequences in the coding region may disrupt the ORF, which may cause severe damage or even death to the virus [47, 48]. Therefore, only a few repeat sequences were observed in the coding regions of viral genomes such as Dengue virus (DEN) and Yellow fever [49, 50]. However, repeat sequences in noncoding regions also play an important role in the viral life cycle [51,52,53]. Some repeat sequences can affect the replication and transcription of viruses [51, 53]. For example, the repeat sequences of the herpes simplex virus contain a promoter of a gene that can regulate virus gene transcription [51]. Thus, our results provide evidence to conclude that repeat sequences are structural elements of the ASFV genome that underpin its distinct functionality.

Second, repeat sequences were predicted to have a higher ratio of loop structures and a lower ratio of disorder regions compared to nonrepeat protein sequences. Similar to the results of previous studies, repeat protein sequences had a large ratio of loop structures and could provide more protein binding sites [54, 55]. The repeat sequences may be less likely to form disordered regions due to their flexible nature. Because the repeat sequences contain a low ratio of polar amino acids (such as tyrosine (Y) and arginine (R)), it is difficult to form polyelectrolytes, which are an important part of the disordered region [56]. Based on these results, we conclude that repeat sequences are important in the conformation of ASFV-specific structural peptides.

Third, our results demonstrated that the ASFV repeat sequences tended to encode penetrating peptides and antimicrobial peptides. A small number of repeat sequences have been identified to encode penetrating peptides and antimicrobial peptides [45, 57]. The repeat sequences encoding penetrating peptides can play an important role in ASFV entering host cells by receptor-mediated endocytosis [11]. Moreover, the repeat sequences encoding antimicrobial peptides can help in the development of new antimicrobial peptide drugs [58, 59], as antimicrobial peptides originating from viral proteins are key players in the development of innovative drug delivery systems [45, 57]. These findings imply that the ASFV repeat sequences have a specific role in functional peptides.

Fourth, our results showed that the pan-repeat sequences of ASFV presented an open state and that core repeats were non-existent in all of the analysed strains. Similar to bacteria, open pan-repeat sequences indicate that the genetic material of ASFV may be constantly mutated [60]. The core repeat sequences in ASFV were not observed, which suggests that ASFV repeat sequences are diverse. Furthermore, repeat sequences have been reported to be related to homologous recombination [12]. As a result, the diversity of repeat sequences can further promote the diversity of ASFV. Overall, our findings underscore the involvement of repeat sequences in the genomic evolution and variation of ASFV, highlighting the importance of further investigation into repeat sequences in other viruses.

This study has two limitations that need to be acknowledged. First, the location, size and number of repeat sequences observed in ASFV genomes may be affected by the repeat sequence identification algorithm. For comprehensive detection of repeat sequences, two complementary methods (TRF and REPuter) were used for identifying repeat sequences in ASFV genomes. Second, the predicted structure and functional features of repeat sequences have not been examined experimentally. More efforts are needed to validate these predictions and to clarify the roles of repeat sequences of the ASFV genome.

In summary, this study is the first to illustrate the distribution, structural, functional, and evolutionary properties of repeat sequences in ASFVs, thus advancing our understanding of the virus and facilitating the development of effective anti-ASFV medications. The methodological framework established in this study can be utilized to guide future research on repeat sequences of other viruses.

References

  1. Parker J, Plowright W, Pierce M (1969) The epizootiology of African swine fever in Africa. Vet Rec 85:668–674

    CAS  PubMed  Google Scholar 

  2. Salguero FJ (2020) Comparative pathology and pathogenesis of African swine fever infection in swine. Front Vet Sci 7:282

    Article  PubMed  PubMed Central  Google Scholar 

  3. Gaudreault NN, Madden DW, Wilson WC, Trujillo JD, Richt JA (2020) African swine fever virus: an emerging DNA arbovirus. Front Vet Sci 7:215

    Article  PubMed  PubMed Central  Google Scholar 

  4. Sun E, Zhang Z, Wang Z, He X, Zhang X, Wang L, Wang W, Huang L, Xi F, Huangfu H, Tsegay G, Huo H, Sun J, Tian Z, Xia W, Yu X, Li F, Liu R, Guan Y, Zhao D, Bu Z (2021) Emergence and prevalence of naturally occurring lower virulent African swine fever viruses in domestic pigs in China in 2020. Sci China Life Sci 64:752–765

    Article  CAS  PubMed  Google Scholar 

  5. Ge S, Li J, Fan X, Liu F, Li L, Wang Q, Ren W, Bao J, Liu C, Wang H, Liu Y, Zhang Y, Xu T, Wu X, Wang Z (2018) Molecular characterization of African swine fever virus, China, 2018. Emerg Infect Dis 24:2131–2133

    Article  PubMed  PubMed Central  Google Scholar 

  6. Dixon LK, Chapman DA, Netherton CL, Upton C (2013) African swine fever virus replication and genomics. Virus Res 173:3–14

    Article  CAS  PubMed  Google Scholar 

  7. Cackett G, Matelska D, Sýkora M, Portugal R, Malecki M, Bähler J, Dixon L, Werner F (2020) The African swine fever virus transcriptome. J Virol 94:e00119-00120

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Alejo A, Matamoros T, Guerra M, Andrés G (2018) A proteomic atlas of the African swine fever virus particle. J Virol 92:e01293

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Versalovic J, Koeuth T, Lupski R (1991) Distribution of repetitive DNA sequences in eubacteria and application to finerpriting of bacterial enomes. Nucleic Acids Res 19:6823–6831

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Yao Q, Compans RW (1996) Peptides corresponding to the heptad repeat sequence of human parainfluenza virus fusion protein are potent inhibitors of virus infection. Virology 223:103–112

    Article  CAS  PubMed  Google Scholar 

  11. Yang S, Zhang X, Cao Y, Li S, Shao J, Sun S, Guo H, Yin S (2021) Identification of a new cell-penetrating peptide derived from the african swine fever virus CD2v protein. Drug Deliv 28:957–962

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zhu Z, Xiao CT, Fan Y, Cai Z, Lu C, Zhang G, Jiang T, Tan Y, Peng Y (2019) Homologous recombination shapes the genetic diversity of African swine fever viruses. Vet Microbiol 236:108380

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Farlow J, Donduashvili M, Kokhreidze M, Kotorashvili A, Vepkhvadze NG, Kotaria N, Gulbani A (2018) Intra-epidemic genome variation in highly pathogenic African swine fever virus (ASFV) from the country of Georgia. Virol J 15:190

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Goller KV, Malogolovkin AS, Katorkin S, Kolbasov D, Titov I, Höper D, Beer M, Keil GM, Portugal R, Blome S (2015) Tandem repeat insertion in African swine fever virus, Russia, 2012. Emerg Infect Dis 21:731

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kim SH, Lee SI, Jeong HG, Yoo J, Jeong H, Choi Y, Son K, Jheong WH (2021) Rapid emergence of African swine fever virus variants with different numbers of a tandem repeat sequence in South Korea. Transbound Emerg Dis 68:1726–1730

    Article  CAS  PubMed  Google Scholar 

  16. Shang WH, Hori T, Toyoda A, Kato J, Popendorf K, Sakakibara Y, Fujiyama A, Fukagawa T (2010) Chickens possess centromeres with both extended tandem repeats and short non-tandem-repetitive sequences. Genome Res 20:1219–1228

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Shammas MA (2011) Repetitive sequences, genomic instability, and Barrett’s esophageal adenocarcinoma. Mobile Genet Elem 1:208–212

    Article  Google Scholar 

  18. Boeva V (2016) Analysis of genomic sequence motifs for deciphering transcription factor binding and transcriptional regulation in eukaryotic cells. Front Genet 7:24

    Article  PubMed  PubMed Central  Google Scholar 

  19. Lynch VJ, Leclerc RD, May G, Wagner GP (2011) Transposon-mediated rewiring of gene regulatory networks contributed to the evolution of pregnancy in mammals. Nat Genet 43:1154–1159

    Article  CAS  PubMed  Google Scholar 

  20. Gymrek M, Willems T, Guilmatre A, Zeng H, Markus B, Georgiev S, Daly MJ, Price AL, Pritchard JK, Sharp AJ, Erlich Y (2016) Abundant contribution of short tandem repeats to gene expression variation in humans. Nat Genet 48:22–29

    Article  CAS  PubMed  Google Scholar 

  21. Gymrek M, Willems T, Reich D, Erlich Y (2017) Interpreting short tandem repeat variations in humans using mutational constraint. Nat Genet 49:1495–1501

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Lunyak VV, Prefontaine GG, Núñez E, Cramer T, Ju BG, Ohgi KA, Hutt K, Roy R, García-Díaz A, Zhu X, Yung Y, Montoliu L, Glass CK, Rosenfeld MG (2007) Developmentally regulated activation of a SINE B2 repeat as a domain boundary in organogenesis. Science 317:248–251

    Article  CAS  PubMed  Google Scholar 

  23. Lu JY, Shao W, Chang L, Yin Y, Li T, Zhang H, Hong Y, Percharde M, Guo L, Wu Z, Liu L, Liu W, Yan P, Ramalho-Santos M, Sun Y, Shen X (2020) Genomic repeats categorize genes with distinct functions for orchestrated regulation. Cell Rep 30:3296-3311.e5

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Schmidt D, Schwalie PC, Wilson MD, Ballester B, Gonçalves Â, Kutter C, Brown GD, Marshall A, Flicek P, Odom DT (2012) Waves of retrotransposon expansion remodel genome organization and CTCF binding in multiple mammalian lineages. Cell 148:335–348

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Sayers EW, Cavanaugh M, Clark K, Ostell J, Pruitt KD, Karsch-Mizrachi I (2019) GenBank. Nucleic Acids Res 47:D94–D99

    Article  CAS  PubMed  Google Scholar 

  26. Kurtz S, Choudhuri JV, Ohlebusch E, Schleiermacher C, Stoye J, Giegerich R (2001) REPuter: the manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res 29:4633–4642

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Benson G (1999) Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res 27:573–580

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Besemer J, Lomsadze A, Borodovsky M (2001) GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res 29:2607–2618

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhu Z, Meng G (2020) ASFVdb: an integrative resource for genomic and proteomic analyses of African swine fever virus. https://asfvdb.popgenetics.net/. Accessed 20 Mar 2021

  30. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:3389–3402

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Liu B, Li K, Huang DS, Chou KC (2018) iEnhancer-EL: identifying enhancers and their strength with ensemble learning approach. Bioinformatics 34:3835–3842

    Article  CAS  PubMed  Google Scholar 

  32. Liu B, Yang F, Huang DS, Chou KC (2018) iPromoter-2L: a two-layer predictor for identifying promoters and their types by multi-window-based PseKNC. Bioinformatics 34:33–40

    Article  CAS  PubMed  Google Scholar 

  33. Walsh I, Martin AJ, Di Domenico T, Tosatto SC (2012) ESpritz: accurate and fast prediction of protein disorder. Bioinformatics 28:503–509

    Article  CAS  PubMed  Google Scholar 

  34. Klausen MS, Jespersen MC, Nielsen H, Jensen KK, Jurtz VI, Soenderby CK, Sommer MOA, Winther O, Nielsen M, Petersen B, Marcatili P (2019) NetSurfP-2.0: improved prediction of protein structural features by integrated deep learning. Proteins 87:520–527

    Article  CAS  PubMed  Google Scholar 

  35. Waghu FH, Barai RS, Gurung P, Idicula-Thomas S (2016) CAMPR3: a database on sequences, structures and signatures of antimicrobial peptides. Nucleic Acids Res 44:D1094–D1097

    Article  CAS  PubMed  Google Scholar 

  36. Manavalan B, Subramaniyam S, Shin TH, Kim MO, Lee G (2018) Machine-learning-based prediction of cell-penetrating peptides and their uptake efficiency with improved accuracy. J Proteome Res 17:2715–2726

    Article  PubMed  Google Scholar 

  37. Bastos AD, Penrith ML, Cruciere C, Edrich J, Hutchings G, Roger F, Couacy-Hymann E, Thomson R (2003) Genotyping field strains of African swine fever virus by partial p72 gene characterisation. Arch Virol 148:693–706

    Article  CAS  PubMed  Google Scholar 

  38. Katoh K, Standley DM (2013) MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30:772–780

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Kumar S, Stecher G, Li M, Knyaz C, Tamura K (2018) MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol 35:1547–1549

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Muangkram Y, Sukmak M, Wajjwalku W (2015) Phylogeographic analysis of African swine fever virus based on the p72 gene sequence. Genet Mol Res 14:4566–4574

    Article  CAS  PubMed  Google Scholar 

  41. Zhao Y, Jia X, Yang J, Ling Y, Zhang Z, Yu J, Wu J, Xiao J (2014) PanGP: a tool for quickly analyzing bacterial pan-genome profile. Bioinformatics 30:1297–1299

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. R Core Team. R (2019) A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/. Accessed 12 Jul 2020

  43. Rice P, Longden I, Bleasby A (2000) EMBOSS: the European molecular biology open software suite. Trends Genet 16:276–277

    Article  CAS  PubMed  Google Scholar 

  44. Sanner, M. F (1999) Python: a programming language for software integration and development. J Mol Graph Model 17:57–61. https://www.python.org/. Accessed 10 June 2020

  45. Freire JM, Almeida Dias S, Flores L, Veiga AS, Castanho MA (2015) Mining viral proteins for antimicrobial and cell-penetrating drug delivery peptides. Bioinformatics 31:2252–2256

    Article  CAS  PubMed  Google Scholar 

  46. Campbell JI, Binns MM, Tomley FM, Boursnell ME (1989) Tandem repeated sequences within the terminal region of the fowlpox virus genome. J Gen Virol 70:145–154

    Article  CAS  PubMed  Google Scholar 

  47. Singh M, Kishore A, Maity D, Sunanda P, Krishnarjuna B, Vappala S, Raghothama S, Kenyon LC, Pal D, Sarma JD (2019) A proline insertion-deletion in the spike glycoprotein fusion peptide of mouse hepatitis virus strongly alters neuropathology. J Biol Chem 294:8064–8087

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Wang Z, Pan Q, Gendron P, Zhu W, Guo F, Cen S, Wainberg MA, Liang C (2016) CRISPR/Cas9-derived mutations both inhibit HIV-1 replication and accelerate viral escape. Cell Reps 15:481–489

    Article  CAS  Google Scholar 

  49. Shurtleff AC, Beasley DW, Chen JJ, Ni H, Suderman MT, Wang H, Xu R, Wang E, Weaver SC, Watts DM, Russell KL, Barrett AD (2001) Genetic variation in the 3′ non-coding region of dengue viruses. Virology 281:75–87

    Article  CAS  PubMed  Google Scholar 

  50. Wang E, Weaver SC, Shope RE, Tesh RB, Watts DM, Barrett AD (1996) Genetic variation in yellow fever virus: duplication in the 3′ noncoding region of strains from Africa. Virology 225:274–281

    Article  CAS  PubMed  Google Scholar 

  51. Chou J, Roizman B (1986) The terminal a sequence of the herpes simplex virus genome contains the promoter of a gene located in the repeat sequences of the L component. J Virol 57:629–637

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Madhry D, Pandey KK, Kaur J, Rawat Y, Sapra L, Kumar R, Srivastava RK, Bhattacharyya S, Verma B (2021) Role of non-coding RNAs in dengue virus-host interaction. Front Biosci (Schol Ed) 13:44–55

    Article  CAS  PubMed  Google Scholar 

  53. Vychodil T, Conradie AM, Trimpert J, Aswad A, Bertzbach LD, Kaufer BB (2021) Marek’s disease virus requires both copies of the inverted repeat regions for efficient in vivo replication and pathogenesis. J Virol 95:e01256-e1320

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Weng Y, Weiss CD (1998) Mutational analysis of residues in the coiled-coil domain of human immunodeficiency virus type 1 transmembrane protein gp41. J Virol 72:9676–9682

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Qi D, DeYoung BJ, Innes RW (2012) Structure-function analysis of the coiled-coil and leucine-rich repeat domains of the RPS5 disease resistance protein. Plant Physiol 158:1819–1832

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Das RK, Ruff KM, Pappu RV (2015) Relating sequence encoded information to form and function of intrinsically disordered proteins. Curr Opin Struct Biol 32:102–112

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Fox JL (2013) Antimicrobial peptides stage a comeback: Better understanding of the mechanisms of action, modification and synthesis of antimicrobial peptides is reigniting commercial development. Nat Biotechnol 31:379–383

    Article  CAS  PubMed  Google Scholar 

  58. Huang Y, Jiang Y, Wang H, Wang J, Shin MC, Byun Y, He H, Liang Y, Yang VC (2013) Curb challenges of the “Trojan Horse” approach: smart strategies in achieving effective yet safe cell-penetrating peptide-based drug delivery. Adv Drug Deliv Rev 65:1299–1315

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Moretta A, Scieuzo C, Petrone AM, Salvia R, Manniello MD, Franco A, Lucchetti D, Vassallo A, Vogel H, Sgambato A, Falabella B (2021) Antimicrobial peptides: a new hope in biomedical and pharmaceutical fields. Front Cell Infect Microbiol 11:668632

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Wang Z, Jia L, Li J, Liu H, Liu D (2020) Pan-genomic analysis of African swine fever virus. Virol Sin 35:662–665

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Natural Science Foundation of China (32170651) and Hunan Provincial Natural Science Foundation of China (2020JJ3006).

Author information

Authors and Affiliations

Authors

Contributions

ZZ, SG and YP carried out the experiments and drafted the manuscript. YP, ZZ, XW and LM critically revised the manuscript and the experimental design. ZZ, SG, RJ, ZC, YW, CL, PF, XW and YP helped with the experiments. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Longfei Mao, Xiaodong Wu or Yousong Peng.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1. The viruses used in this study.

Additional file 2. The functional classification of ASFV proteins with short interspersed repeated sequence (SINE).

Additional file 3. The intersection of penetrating peptides and antimicrobial peptides in repeat protein sequences.

Additional file 4. Pan-repeat sequence analysis of ASFV genotypes I and II.

(A and C) Flower plot showing the core repeat and dispensable repeat sequence clusters of ASFV genotypes I and II, respectively. The diagram depicts the core repeat sequence cluster number (orange) and the dispensable repeat sequence cluster number (blue) for ASFV genotypes I and II, respectively. The number of unique repeat sequences in each strain is shown beside the accession number of the ASFV genome. (B and D) The relationship between genome number and dispensable repeat and core repeat sequence profiles for ASFV genotypes I and II, respectively. The black and blue curves refer to the least-squares fitting for the average number of dispensable repeat clusters and core repeat clusters versus the number of genomes, respectively. The mathematical functions and the R-squared values of both models are delineated on the graph.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, Z., Ge, S., Cai, Z. et al. Systematic identification and characterization of repeat sequences in African swine fever virus genomes. Vet Res 53, 101 (2022). https://doi.org/10.1186/s13567-022-01119-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13567-022-01119-9

Keywords