Skip to main content

Origin and evolution of emerging Liao ning Virus (genus Seadornavirus, family Reoviridae)

Abstract

Background

Liao ning virus (LNV) is a member of the genus Seadornavirus, family Reoviridae and has been isolated from kinds of vectors in Asia and Australia. However, there are no systematic studies describe the molecular genetic evolution and migration of LNVs. With the development of bioinformatics, viral genetic data combining the information of virus isolation time and locations could be integrated to infer the virus evolution and spread in nature.

Methods

Here, a phylogenetic and phylogeographic analysis using Bayesian Markov chain Monte Carlo simulations was conducted on the LNVs isolated from a variety of vectors during 1990–2014 to identify the evolution and migration patterns of LNVs.

Results

The results demonstrated that the LNV could be divided into 3 genotypes, of which genotype 1 mainly composed of LNVs isolated from Australia during 1990 to 2014 and the original LNV strain (LNV-NE97–31) isolated from Liaoning province in northern China in 1997, genotype 2 comprised of the isolates all from Xinjiang province in western China and genotype 3 consisted the isolates from Qinghai and Shanxi province of central China. LNVs emerged about 272 years ago and gradually evolved into three lineages in the order genotype 1, genotype 2 and genotype 3. Following phylogeographic analysis, it shows genotype 1 LNVs transmitted from Australia (113°E-153°E,10°S-42°S) to Liaoning province (118°E-125°E,38°N-43°N) in Northeast Asian continent then further spread across the central part of China to western China (75°E-95°E,35°N-50°N).

Conclusion

LNVs were initially isolated from Liaoning province of China in the Northeast Asia, however, the present study revealed that LNVs were first appeared in Australia in the South Pacific region and transmitted to mainland China then rapidly spread across China and evolved three different genotypes. The above results suggested that LNV had the characteristics of long-distance transmission and there were great genetic diversity existed in the LNV population. Notably, current information of 80 strains of LNVs are limited. It is of great importance to strengthen the surveillance of LNVs to explore its real origin in nature and monitoring of the LNVs’ population variation and maintain vigilance to avoid LNV breaking through the species barrier and further clarify its relationship to human and animal infection.

Introduction

Liao ning virus (LNV) belongs to the genus Seadornavirus of the family Reoviridae [1]. The genus Seadornavirus comprises 3 species, Banna virus (BAV), Kadipiro virus (KDV) and LNV [1,2,3]. The Seadornaviruses have a genome composed of 12 segments of double-stranded RNA, which decreases with the relative molecular weight during gel electrophoresis, and is named 1–12 segments [2]. The full length of the LNV genome is about 21,000 bp, with segment lengths that range from 3747 bp (segment 1) to 759 bp (segment 12). Each segment contains a open reading frame encodes a viral protein (Viral protein, VP), of which 5 are non-structural proteins (VP5, VP6, VP7, VP11, VP12), 7 are structural proteins (VP1, VP2, VP3, VP4, VP8, VP9, VP10).VP10 protein is an outer capsid structural protein which directly interacts with the receptor of host cells and is the main region that determines the antigenicity of the virus [3].

It is reported that the Seadornaviruses could be carried by a variety of blood-sucking insects including mosquitoes, ticks, and midges and even have been isolated from pigs, cattle and patients with fever and encephalitis [4,5,6,7,8,9,10]. Studies have also shown that BAV, which is the prototype of genus Seadornavirus, is an emerging pathogen that causes human viral encephalitis [2]. Therefore, Seadornaviruses may well be a group of newly discovered viruses close related to human and animal diseases [2, 4].

LNV was initially isolated from mosquitoes (Aedes dorsalis) in Liaoning province, northeast China in 1997 [10], since then additional isolates have been isolated from mosquitoes (genus culex and Aedes) in Shanxi province and Qinghai province in China [5, 7, 8].Interestingly,it is observed that all of the LNVs have only been isolated from the northern part of China and no LNV has been reported in other provinces in China and abroad [6].Once, the LNV has been considered to be a virus specie only restricted in northern part of China. Until a recent research reported that number of LNV strains have been isolated from mosquitoes of four genera (Culex, Anopheles, Mansonia and Aedes) in Australia in 2016 [11, 12]. This result demonstrated that LNV was not only limited in mainland China but also widely distributed in Australia in the South Pacific region. All of these research suggested that LNV was a widely distributed virus that could be transmitted by kinds of blood-sucking vectors and its geographical distribution even exceeds that of the other two Seadornaviruses (BAV, KDV) [1, 4].

Previous researches have reported on the molecular genetic evolution of LNV isolated in China and Australia [5, 7, 11]. However, the data used for analysis were only limited to the local LNV isolates [5, 7, 11]. China in Asia and Australia in the South Pacific are located in the northern and southern hemispheres, respectively. The natural environment in the two regions were significant different from each other. These observations raise several questions. What is the genetic relationship between the LNVs isolated from different regions? What are the characteristics of the molecular evolution of the entire LNV population? Where is the origin of LNVs? Recent revolutionary developments in virology bioinformatics provide unprecedented opportunities for analyzing the viral genetic data (nucleotide or amino acid) to model the evolutionary relationships between virus samples and help to explain epidemiological patterns and uncover processes of transmission [13,14,15]. Examples of such analyses include the estimation of origin and divergence time of Japanese encephalitis virus (JEV) [16, 17],identification of phylogenetic relationships among viral lineages in Zika virus (ZIKV) [18], inference of phylogeographic history and migration patterns of recent epidemics of Middle East respiratory syndrome coronavirus (MERS-CoV) [13] and Ebola virus (EBOV) [19]. Therefore, in the present study a comprehensive phylogeographic study of all the LNVs isolated worldwide were conducted to explore the molecular genetic evolution and migration patterns of LNVs isolated from different vectors since 1990 in China and Australia.

Materials and methods

Data set construction of the 10th segment gene of LNV

LNV is a 12-segment double-stranded RNA virus. The 10th segment gene encodes the viral cell attachment protein, which directly interacts with the receptor on host cells [3]. Thus, the phylogenetic analysis based on this segment can best demonstrated the evolutionary and dispersal trends of viral strains from different isolation sites and vectors. Thus,we downloaded all the 10th segment gene sequences of LNV from GenBank as of June 2019 (Table 1). The data set contained 80 sequences, representing the samples isolated from kinds of mosquitoes:genus Culex (n = 64), genus Aedes (n = 10), genus Ochlerotatus (n = 3),unidentified mosquitoes (n = 3). The LNV isolation sites included the Liaoning province, Shanxi province, Qinghai province and Xinjiang province in China as well as New South Wales, Northern Territory, Queensland and Western Australia in Australia. The geographical distribution extends from Australia (113°E-153°E, 10°S-42°S) to mainland China (73°E-135°E, 3°N-53°N) (Fig. 1).

Table 1 The background information of LNVs used in this study
Fig. 1
figure 1

The distribution of LNV worldwide. LN: Liaoning province, SX: Shanxi province, QH: Qinghai province, XJ: Xinjiang province, WA: Western Australia, BB: Brisbane, SY: Sydney, NSW: New South Wales, BW: Bradshaw, PR: Peel region

Time scaled phylogenetic and phylogeographic analysis of LNV

The 10th segment sequence database of LNV was analyzed using Bayesian Markov chain Monte Carlo (MCMC) method [20]. The GTR + I + G substitution model was selected to be the optimal model by MrModelTest [21]. Bayesian time scaled phylogenetic analysis and nucleotide substitution rate and most recent common ancestor (tMRCA) of LNVs were coestimated using the BEAST software package [22]. The relaxed clock model with different demographic models was tested [23], and the best models were selected by means of a Bayes factor (BF) test using marginal likelihoods values (2lnBF > 2) and 95% highest posterior density (HPD) intervals. The chain length was 1,000,000,000 generations. Convergence of parameters was checked using TRACER (http://beast.community/tracer) and was indicated as effective sample size (ESS>200), and the maximum clade credibility (MCC) tree was constructed using TreeAnnotator (http://beast.community/treeannotator) with 10% burn-in.

In order to clarify the geographical dispersal history of LNVs, the bayesian stochastic search variable selection (BSSVS) was used to provide evidence for statistically supported diffusion between state variables under BEAST software package [24]. This method estimates the most probable state at each node in the MCC trees, allowing us to reconstruct ancestral positions for ancestral viral lineages along the tree. For phylogeographic reconstructions, each region was coded as a discrete trait. BSSVS output and surfaces representing uncertainty for continuous diffusion processes were formatted as KML using the SPREAD software [25]. Determination of each locality was coordinated and performed using Google Earth. ArcGis was finally used to display the dispersal pattern of LNV based on the phylogeographic analysis.

The sequence analysis of the 10th segment of LNV

Sequence similarity and nucleotide base composition analysis were performed using the CLUSTALX software [26], BioEdit software [27] MEGA-X [28] and MegAlign (DNASTAR, Madison, WI, USA). Hemi (version1.0) [29] software visualized the similarity comparison results of LNV and Rstudio made the nucleotide base composition analysis map of LNV.

Three-dimensional structural analysis of VP10 of LNVs

The crystal structure of the cell attachment protein VP9 of BAV (PDB 1w9z) was selected as the best template for the homology modeling. The YASARA software [30], PyMol software [31] and VMD (version 1.9.2) [32] were used to analyze the structures and surface charge density of VP10 among LNV isolates.

Result

The evolution of LNV isolated worldwide based on the time-scaled phylogenetic analysis

According to the 95% HPD intervals and bayes factor, the Bayesian skyline model with a relaxed molecular clock was selected to be the best fit model. The MCC tree was established using the Bayesian Markov chain Monte Carlo (MCMC) approach (Fig. 2). The posterior probability values of each branch node were all greater than 0.7, showing the robustness of the result. The mean nucleotide substitution rate for the entire LNV population was estimated to be 1.6 × 10− 3 s/s/y (95%HPD:2.8 × 10− 4,3.3 × 10− 3) (Table 2). At this evolution rate, the tMRCA of the LNV was calculated to be 272.6 years ago (95%HPD:58.9,790.3). The LNV evolved chronologically into three major evolutionary populations, genotype 1, genotype 2 and genotype 3.Genotype1, which included the initial isolate (LNV-NE97–31) from Liaoning province of China and seven strains from Australian, emerged 73.0 years ago (95%HPD:28.9,193.0) and demonstrated to be the oldest lineage. Genotype 2 emerged about 46.5 years ago (95%HPD:18.8,113.5) and was composed of LNV strains isolated from Xinjiang province in western China. Genotype 3, which appeared approximately 25.3 years ago (95%HPD:17.4,46.6), included isolates from Qinghai province, Shanxi province and Liaoning province (another initial isolate LNV-NE97–12), and was the youngest LNV lineage.

Fig. 2
figure 2

Evolution analysis of LNV isolated worldwide from 1990 to 2014. Maximum clade credibility tree for the 10th segment of LNVs isolated from 1990 to 2014 worldwide. The tree identified 3 distinct lineages, genotype 1, genotype 2 and genotype 3.Genotype 1 contained the initial LNV isolate (LNSV-NE97–31) and seven Australian isolates (Blue). Genotype 2 consisted of isolates all from Xingjiang province in western China (Green). Genotype 3 comprised LNV isolates from Qinghai province, Shanxi province and Liaoning province (another initial LNV isolate LNSV-NE97–12) (Red). Estimated tMRCAs of these lineages (with their 95% HPD values in parentheses) are presented. The posterior probability values of main branch nodes were indicated. The letters G1, G2 and G3 represent genotypes 1,2 and 3, respectively

Table 2 Sequence similarity and Bayesian Markov chain Monte Carlo analysis of LNVs

Population dynamics of LNV

The skyline plot of the LNV population dynamics is shown in Fig. 3. The population dynamic of LNV kept stable from 1990 to 1997.Since 1997, the population diversity began to decrease and reached the lowest point at the year of 2003.From 2003 the LNV population diversity increased, peaking around 2005.After a minor decrease during around 2005 to 2006, the virus population was stable from 2006 to 2014.

Fig. 3
figure 3

Bayesian skyline plot for LNV. Highlighted areas correspond to 95% HPD intervals

The sequence similarity analysis of LNV

The Sequence analyses of 80 strains of LNV isolated from different vectors and locations during 1990–2014 revealed that the level of nucleotide similarity varied from 70.4 to 100% with an average value of μ = 92.6% while the amino acid similarity ranged from 75.0 to 100% with an average value μ = 94.6% (Fig. 4). Genotype 1 exhibited nucleotide similarity between 87.5–100% (μ = 94.3%) and amino acid similarity range from 91.2 to 100% (μ = 96.4%).Genotype 2 exhibited nucleotide similarity from 92.5 to 100% (μ = 97.8%) and amino acid similarity from 92.5 to 100% (μ = 98.9%). Genotype3 exhibited nucleotide similarity from 97.5 to 100% (μ = 99.2%) and amino acid similarity from 98.1 to 100% (μ = 99.4%). The results demonstrated that the genotype 1 population possessed the highest degree of nucleotide and amino acid variation among all the three LNV genotypes.

Fig. 4
figure 4

The Nucleotide and amino acid similarity of the 10th segment of LNV. a: Nucleotide similarity of LNV. b: Amino acid similarity of LNV. The horizontal and vertical axes represent the names of different genotypes of LNV isolates. The names of the isolates colored blue, green and red represent the genotypes 1,2 and 3, respectively

Nucleotide composition of the 10th segment genome of LNV

The values of nucleotide contents in the entire LNV population and each genotype were analyzed (Fig. 5). The results showed that VP10 genome of LNV was rich in A and U (A + U > C + G) and the percentage of AU was reaching 59.5%. The A%, U%, G%, C% are 26.0% ±0.6 (mean ± SD), 32.3% ± 1.3, 18.7% ± 0.4, 23.1% ± 0.4, respectively. There was a common trend in the usage of A,U,G,C between the three genotypes that was U > A > 0.25 > C > G (Fig. 5). As for each genotype, the average contents of A,U,G,C for genotype 1 were 27.6% ± 0.5, 29.0% ± 0.5, 19.8% ± 0.2, 23.5% ± 0.3,for genotypes 2 were 25.8% ± 0.3, 32.8% ± 0.5,18.5% ± 0.2, 23.0% ± 0.4 and these values for genotypes 3 were 26.3% ± 0.2, 30.7% ± 0.3, 19.4% ± 0.1, 23.7% ± 0.2. The value of the A content in genotype 1 was the highest while that in genotype 2 was the lowest.

Fig. 5
figure 5

The A, U, G and C contents of different genotypes of the 10th segment of LNV. The vertical axes indicates the percentage of nucleotide contents. The cross bar represents the different isolates in the same genotype. The red bar identifies adenine (A), green identifies cytosine (C), cyan identifies guanine (G) and violet identifies uracil (U)

The dispersal route of LNV base on phylogeographic analysis

The estimated history of LNV dispersal route over time is shown in detail in Fig. 6. According to our results of the phylogeographic analysis, LNV was originated in Australia in the south pacific region and then initially introduced to Liaoning province in China in the Northeast Asia around 1980s.Subsequently, the virus was further westwarded spread to Shanxi province, Qinghai province and Xingjiang province of China during 1990s.

Fig. 6
figure 6

Dispersal route of LNVs based on phylogeographic analysis. Gray shadow presents the areas with LNV isolated. Line and arrow represents the transmission route and direction of LNV. Red line: Spread from Australia to Liaoning province of China in 1980s; Green line: Spread from Liaoning province to Shanxi province, Qinghai province and Xinjiang province in 1990s. LN: Liaoning province, SX: Shanxi province, QH: Qinghai province, XJ: Xinjiang province, AU: Australia

The amino acid comparison and the three dimensional structural analysis of the cell attachment protein (VP10) of LNV

To characterize the mutations in the cell attachment protein (VP10) between the three genotypes of LNV, we analyzed the amino acids encoded by the VP10 genes derived from the previously mentioned LNV isolates. The results revealed that the ORF of VP10 was 753 nucleotides in length and encoded 250 amino acids. The amino acid sequences of genotype 1 were significantly different from those of the rest two genotypes while the genotype 2 and 3 shared a great similarity in amino acid sequences. The detailed results were that a total of 20 common amino acid differences identified between genotype 1 to genotype 2 and 3 (Table 3). The amino acid sequences were highly similar between genotype2 and 3 and there were only 3 amino acid differences were identified (Table 3). When compared the amino acid sequences of genotype 1 to those of genotype 2 and 3, twenty common amino acid differences were identified as follows: S44A, A60G, D64G, S65N, H68V, N82A, P99A, V122T, M131S, A132G, S160N, I167S, Q171S, A181L, E207D, S209V, T217N, S219N, N223H, N238S.The amino acid sequences of genotype 2 and 3 were highly similar, there were only three amino acid mutations: S58A, A130N, L210H.

Table 3 Comparison of VP10 amino acid sequences between different LNV genotypes

In order to further explore whether these mutations affected the three shape and the surface charge density of the cell attachment protein of LNV, the three dimensional structures and electrostatic potential analysis have been conducted between different genotypes of LNV. Seven of the twenty common amino acid differences were located in the β-sheet (residues 82, 122, 131, 132, 181, 223, and 238), while the remaining thirteen were located in the loop region (residues 44, 60, 64, 65, 68, 99, 160, 167, 171, 207, 209, 217 and 219). These differences have some effect on the structure and electrostatic properties of VP10 between the different genotypes. In particular, the residue at position 64, 65, 131, 132, 167, 171, 207, 209 and 219.There were great differences existed the R group of the amino acid and the electrical polarity of the mentioned 9 mutations, resulting in significant differences in structure and charge at these sites between the LNVs (Fig. 7a’, b’, c’).

Fig. 7
figure 7

The theoretical three-dimensional structure and surface charge of different genotypes of LNV. a, b, c: The theoretical three-dimensional structure of LNV of genotype1,2,3 and the differences sites of amino acid within 3 genotypes. a’, b’, c’: The theoretical three-dimensional structure as well as the surface charge distribution of LNV of genotype 1,2,3. def and d’e’f’ are the angles of view of abc and a’b’c’ rotated the Y axis by 180 degrees, respectively

Discussion

During the arboviruses survey in China in 1997, two virus strains (LNSV-NE9731 and LNV-NE9712) were obtained from Aedes dorsalis mosquitoes in Liaoning province in North-East of China which were found to cause cytopathic effect in C6/36 cells. The virus was designated Liao ning virus (LNV) [3, 10]. Since then, additional strains of LNV were isolated from Shanxi province, Qinhai province and Xinjiang province [7, 8]. A more than 30 years national arbovirus surveillance in mainland China revealed the LNV strains had only been isolated in a long and narrow region covering 73–125 ° E and 31–48 ° N in the northwest to northeast part of the country. No strains of LNV have been reported isolated in other regions of China or abroad [6]. However, a recent research reported a total of 35 strains of LNV were isolated from mosquitoes belonging to four genera (Culex, Anopheles, Mansonia and Aedes) collected from 1988 to 2014 in Australia in the southern hemisphere. The isolation sites included New South Wales, Northern Territory, Queensland and Western Australia, almost across the entire Australian continent. Thus, it is obvious that the geographical distribution and the vectors of LNV are wide and variable in Australia. What is more, the initial LNV isolate in Australian was obtained from mosquitoes collected in 1988, predating the first Chinese LNV isolate which was obtained in 1997 [10, 11].

Several phylogenetic analysis have been conducted on the LNV’s nucleotide sequences, previously. The LNVs isolated in China could be divided into 3 evolutionary branches. The LNSV-NE9731 strain isolated in Liaoning province in 1997 formed an independent evolutionary branch while another Liaoning isolates named LNSV-NE9712 clustered together with the isolates from Qinghai province and Shanxi province that formed a branch. All of the Xinjiang isolates grouped together formed the largest evolutionary branch [5, 33]. All of the Australian LNV isolates were divided into two disparate lineage, one composed the isolates from eastern and northern Australia and the other included the isolates from western and southern Australia [11, 12]. In this study, a comprehensive molecular phylogenetic analysis of all the LNVs isolated from China and Austrila demonstrated that the LNVs could be divide into three genotypes (Fig. 2). Genotype 1 included LNSV-NE9731 (isolated in Liaoning provinve of China in 1997) and 7 strains from Austrila (the mosquito collection time was 1990,2005,2007,2013 and 2014, respectively). A total of 68 strains of LNV isolated from Xinjiang province of western China clustered together that formed an independent branch. And the strain LNV-NE9712 isolated in 1997 together with the Qinghai and Shanxi isolates formed the genotype 3 evolutionary branch. This result suggested that the Australia LNV isolates were not only restricted and circulated in Australia and even this virus population transmitted a long way to northern China in Asia then evolved into new LNV populations with local circulation characteristics.

The time-scaled evolutionary analysis of LNV revealed the branching of the lineages occurred in the following order:genotype 1 emerged 73 years ago (95%HPD:28.9, 193.0), the genotype 2 at 46.5 years ago (95%HPD:18.8, 113.5) and genotype 3 at 25.3 years ago (95%HPD:17.4, 46.6.) Thus, genotype 1 is the oldest lineage and the genotype 3 is the youngest one (Table 2). The current results lead to an estimate that the most recent common ancestor of LNV appeared about 277 years ago which similar to that of another member of seadornaviruses BAV, whose tMRCA is 105 years ago [34], indicating that LNV is an emerging virus population.

Although the LNVs were initially isolated from several inland provinces in China, the phylogeopgraphic analysis of our study showed that the LNV likely originated from Australia in the South Pacific region and the genotype 1 viruses spread northward from Australia to Liaoning province in northeast China. The reasons that LNV originated in Australia and spread to the Asian continent may be related to the following factors: First, Australia is located in the South Pacific region and contains a variety of geographic climate types, including the subtropical humid climate in the east, the savanna climate in the northwest, the tropical rain forest climate in the northeast, and Mediterranean climate in the southwest [35]. The climate types in Australia facilitate the local host vectors diverse and can breed LNV population with a strong adaptability. Second, the nucleotide similarity of genotype 1 was the most divergent between all of the three genotypes, indicating greater population variability exists in Australian LNVs.

Additionally, there was a very interesting finding that the AU content was relative high in VP10 of LNV, whereas the AU content was also high within the LNV’s whole viral genome (data not shown). These results were consistent with the prior studies wherein A and U frequencies were higher than C and G frequencies for avian rotaviruses and some flaviviruses including dengue virus (DENV), West Nile virus (WNV), yellow fever virus (YFV) and JEV [36,37,38,39]. Besides, the nucleotide composition of the BAV, which is the prototype species of genus Seadornavirus, was also observed a higher AU content [40]. It has been reported that the AU-rich genome structure facilitate the Human Immunodeficiency virus (HIV) to avoid recognition by the innate immune system of host cell [41]. However, the biological causes and the consequence for increased A and U within the LNV and other mentioned arboviruses genome are still unknown, so enhanced experimental studies are required to explore the AU rich molecular function in LNV.

This may be the reasons why LNVs have been continuously isolated from 12 species of mosquitoes belonging to four genera in different geographical and climatic regions in Australia, during the last 30 years [11, 12]. However, in China the LNV has only been isolated from a long and narrow region restricted from 31°N to 48°N, where belongs to the north temperate zone with cold climate, less rain and low species diversity [5, 7, 10]. What is more, after transmission from Australia to China, the genotype 1 LNV population gradually adapted to the local natural environment and evolved new strains with regional genetic characteristics. The newly evolved LNVs (genotypes 2 and 3) contain 20 amino acid mutations compared with the initial LNVs (genotype 1) in the cell attachment protein (VP10), which altered the structure and electrostatic presentation influencing the binding properties to host vectors. This might be one of the reasons why the newly evolved LNVs (genotyp 2 and 3) were restricted to a relative narrow range of vectors and habitat. Compared with genotype 1 and genotype 2, the genotype 3 which located in the central part of China contains 2 unique amino acids, which were identical with genotype 1 but different from genotype 2,indicating that this lineage is at evolutionary transitional position which preserved the genetic information of the original LNV population and also evolved novel genetic information sites during spread to new locations. When compared with original LNV population, the genotype 2 LNVs in Xinjiang province contained the the maximum numbers of amino acid mutations thus it was the the most divergent lineage and formed a independed evolutionary branch. The genetic informative sites of the entire LNVs population confirmed its transmission path from north to south, and then from east to west.

Conclusions

In this study, the evolution analysis of LNVs revealed that the virus belongs to the emerging virus group. In particular, it was suggested that the genotype 1 LNVs were a group of segmented double-strand RNA virus with extremely adaptability, which have a large group of vectors, a high rate of genetic variation and apparent active transmissibility that can adapt to different geographical environments over a quite long distance. Recently, the genome information of LNV was reported to be detected in Aedes aegypti of African origin, reminding us that LNVs was not only limited to China in Asia and Australia in the South Pacific region but may well be extended to Africa and even posing a high risk of spreading to new areas such as Central Asia and Europe. It is well known that genetic mutation such as recombination or reassortment can easily occur in segmented RNA virus. For example, BAV which also belongs to the Seadornavirus, was originally discovered in southern China. However, several novel variants of BAV have been found, such as the BALV isolate from Hungary [33] and the Mangshi virus from southern China [42]. Considering LNV was the only specie in Seadornavirus that can replicate in mammalian cell lines and cause fatal haemorrhagic symptoms in mice [3], it might be a pathogen that has great potential to cause disease in human and (or) animals. Therefore, to strengthen the research on genetic variation of LNV and to clarify the relationship between LNV and zoonosis is not only a research topic for virologists, but also a scientific issues for public health communities.

Availability of data and materials

All data and materials are included in this published article.

Abbreviations

LNV:

Liao ning virus

LNVs:

Liao ning viruses

BAV:

Banna virus

KDV:

Kadipiro virus

RNA:

RibonucleicAcid

bp:

basepair

VP:

Viral protein

VP10:

The tenth segnment of viral protein

MCMC:

The Bayesian Markov chain Monte Carlo

tMRCA:

The time of most recent ancestor

BEAST:

Bayesian Evolutionary Analysis Sampling Trees

BF:

Bayes factor

HPD:

Highest posterior density

MCC:

The maximum clade credibility

ESS:

The effective sample size

BSSVS:

The bayesian stochastic search variable selection

KML:

Keyhole Markup Language

SPREAD:

Spatial phylogenetic reconstruction of evolutionary dynamics

SD:

Standard Deviation

A,U,G,C:

Adenine,Uracil,Guanine,Cytosine

ORF:

Open reading frame

BALV:

Banna-like virus

References

  1. Attoui H, Billoir F, Biagini P, de Micco P, de Lamballerie X. Complete sequence determination and genetic analysis of Banna virus and Kadipiro virus: proposal for assignment to a new genus (Seadornavirus) within the family Reoviridae. J Gen Virol. 2000;81:1507–15.

    Article  CAS  Google Scholar 

  2. Attoui H, Mohd Jaafar F, de Micco P, de Lamballerie X. Coltiviruses and seadornaviruses in North America, Europe, and Asia. Emerg Infect Dis. 2005;11:1673–9.

    Article  CAS  Google Scholar 

  3. Attoui H, Mohd Jaafar F, Belhouchet M, Tao S, Chen B, Liang G, et al. Liao ning virus, a new Chinese seadornavirus that replicates in transformed and embryonic mammalian cells. J Gen Virol. 2006;87:199–208.

    Article  CAS  Google Scholar 

  4. Liu H, Li MH, Zhai YG, Meng WS, Sun XH, Cao YX, et al. Banna virus, China, 1987-2007. Emerg Infect Dis. 2010;16:514–7.

    Article  Google Scholar 

  5. Lu Z, Liu H, Fu S, Lu X, Dong Q, Zhang S, et al. Liao ning virus in China. Virol J. 2011;8:282.

    Article  Google Scholar 

  6. Liang G, Li X, Gao X, Fu S, Wang H, Li M, et al. Arboviruses and their related infections in China: a comprehensive field and laboratory investigation over the last 3 decades. Rev Med Virol. 2018;28:e1959.

    Article  Google Scholar 

  7. Lv X, Mohd Jaafar F, Sun X, Belhouchet M, Fu S, Zhang S, et al. Isolates of Liao ning virus from wild-caught mosquitoes in the Xinjiang province of China in 2005. PLoS One. 2012;7:e37732.

    Article  CAS  Google Scholar 

  8. Li WJ, Li MH, Wang JL, Lv Z, Cao YX, Wang ZL, et al. Isolation and identification of Liao ning virus from Culex pipiens mosquitoes in Qinghai Province. Chinese J Vector Biology and Control. 2014;25:12–4.

    CAS  Google Scholar 

  9. Sun XH, Meng WS, Fu SH, Feng Y, Zhai YG, Wang JL, et al. The first report of Kadipiro virus isolation in China. Bing Du Xue Bao. 2009;25:173–7.

    CAS  PubMed  Google Scholar 

  10. Tao S, Cai Z. Yang D: [new subtype of coltivirus isolated from mosquitoes in the northeast part of China]. Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi. 1999;13:228–30.

    CAS  PubMed  Google Scholar 

  11. Prow NA, Mah MG, Deerain JM, Warrilow D, Colmant AMG, O'Brien CA, et al. New genotypes of Liao ning virus (LNV) in Australia exhibit an insect-specific phenotype. J Gen Virol. 2018;99:596–609.

    Article  CAS  Google Scholar 

  12. Coffey LL, Page BL, Greninger AL, Herring BL, Russell RC, Doggett SL, et al. Enhanced arbovirus surveillance with deep sequencing: identification of novel rhabdoviruses and bunyaviruses in Australian mosquitoes. Virology. 2014;448:146–58.

    Article  CAS  Google Scholar 

  13. Ibrahim B, McMahon DP, Hufsky F, Beer M, Deng L, Mercier PL, et al. A new era of virus bioinformatics. Virus Res. 2018;251:86–90.

    Article  CAS  Google Scholar 

  14. Geoghegan JL, Holmes EC. Evolutionary virology at 40. Genetics. 2018;210:1151–62.

    Article  CAS  Google Scholar 

  15. Wolf YI, Kazlauskas D, Iranzo J, Lucia-Sanz A, Kuhn JH, Krupovic M, et al. Origins and evolution of the global RNA Virome. mBio. 2018;9:1.

    Article  CAS  Google Scholar 

  16. Pan XL, Liu H, Wang HY, Fu SH, Liu HZ, Zhang HL, et al. Emergence of genotype I of Japanese encephalitis virus as the dominant genotype in Asia. J Virol. 2011;85:9847–53.

    Article  CAS  Google Scholar 

  17. Gao X, Liu H, Li M, Fu S, Liang G. Insights into the evolutionary history of Japanese encephalitis virus (JEV) based on whole-genome sequences comprising the five genotypes. Virol J. 2015;12:43.

    Article  CAS  Google Scholar 

  18. Liu H, Shen L, Zhang XL, Li XL, Liang GD, Ji HF. From discovery to outbreak: the genetic evolution of the emerging Zika virus. Emerg Microbes Infect. 2016;5:e111.

  19. Dudas G, Carvalho LM, Bedford T, Tatem AJ, Baele G, Faria NR, Park DJ, Ladner JT, Arias A, Asogun D, et al. Virus genomes reveal factors that spread and sustained the Ebola epidemic. Nature. 2017;544:309–15.

  20. van Ravenzwaaij D, Cassey P, Brown SD. A simple introduction to Markov chain Monte-Carlo sampling. Psychon Bull Rev. 2018;25:143–54.

    Article  Google Scholar 

  21. Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.

    Article  CAS  Google Scholar 

  22. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.

    Article  CAS  Google Scholar 

  23. Ho SY, Duchene S. Molecular-clock methods for estimating evolutionary rates and timescales. Mol Ecol. 2014;23:5947–65.

    Article  Google Scholar 

  24. Lemey P, Rambaut A, Drummond AJ, Suchard MA. Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009;5:e1000520.

    Article  Google Scholar 

  25. Bielejec F, Rambaut A, Suchard MA, Lemey P. SPREAD: spatial phylogenetic reconstruction of evolutionary dynamics. Bioinformatics. 2011;27:2910–2.

    Article  CAS  Google Scholar 

  26. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.

    Article  CAS  Google Scholar 

  27. Tippmann HF. Analysis for free: comparing programs for sequence analysis. Brief Bioinform. 2004;5:82–7.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  29. Deng W, Wang Y, Liu Z, Cheng H, Xue Y. HemI: a toolkit for illustrating heatmaps. PLoS One. 2014;9:e111988.

    Article  Google Scholar 

  30. Land H, Humble MS. YASARA: a tool to obtain structural guidance in biocatalytic investigations. Methods Mol Biol. 2018;1685:43–67.

    Article  CAS  Google Scholar 

  31. Grell L, Parkin C, Slatest L, Craig PA. EZ-Viz, a tool for simplifying molecular viewing in PyMOL. Biochem Mol Biol Educ. 2006;34:402–7.

    Article  CAS  Google Scholar 

  32. Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14:33–8 27-38.

    Article  CAS  Google Scholar 

  33. Reuter G, Boros A, Delwart E, Pankovics P. Novel seadornavirus (family Reoviridae) related to Banna virus in Europe. Arch Virol. 2013;158:2163–7.

    Article  CAS  Google Scholar 

  34. Liu H, Gao XY, Fu SH, Li MH, Zhai YG, Meng WS, et al. Molecular evolution of emerging Banna virus. Infect Genet Evol. 2016;45:250–5.

    Article  CAS  Google Scholar 

  35. Haines HA, Olley JM, Kemp J, English NB. Progress in Australian dendroclimatology: identifying growth limiting factors in four climate zones. Sci Total Environ. 2016;572:412–21.

    Article  CAS  Google Scholar 

  36. Kattoor JJ, Malik YS, Sasidharan A, Rajan VM, Dhama K, Ghosh S, et al. Analysis of codon usage pattern evolution in avian rotaviruses and their preferred host. Infect Genet Evol. 2015;34:17–25.

    Article  CAS  Google Scholar 

  37. Lara-Ramirez EE, Salazar MI, Lopez-Lopez Mde J, Salas-Benito JS, Sanchez-Varela A, Guo X. Large-scale genomic analysis of codon usage in dengue virus and evaluation of its phylogenetic dependence. Biomed Res Int. 2014;2014:851425.

    Article  Google Scholar 

  38. Moratorio G, Iriarte A, Moreno P, Musto H, Cristina J. A detailed comparative analysis on the overall codon usage patterns in West Nile virus. Infect Genet Evol. 2013;14:396–400.

    Article  CAS  Google Scholar 

  39. van Hemert F, Berkhout B. Nucleotide composition of the Zika virus RNA genome and its codon usage. Virol J. 2016;13:95.

    Article  Google Scholar 

  40. Long S, Yao H, Wu Q, Li G. Analysis of compositional bias and codon usage pattern of the coding sequence in Banna virus genome. Virus Res. 2018;258:68–72.

    Article  CAS  Google Scholar 

  41. Vabret N, Bailly-Bechet M, Najburg V, Muller-Trutwin M, Verrier B, Tangy F. The biased nucleotide composition of HIV-1 triggers type I interferon response and correlates with subtype D increased pathogenicity. PLoS One. 2012;7:e33502.

    Article  CAS  Google Scholar 

  42. Wang J, Li H, He Y, Zhou Y, Meng J, Zhu W, et al. Isolation and genetic characterization of Mangshi virus: a newly discovered Seadornavirus of the Reoviridae family found in Yunnan Province, China. PLoS One. 2015;10:e0143601.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by grants from the National Natural Science Foundation of China (81301479), a Development Grant of the State Key Laboratory of Infectious Disease Prevention and Control (2014SKLID103,2018SKLID309), the Natural Science Foundation of Shandong Province, China (ZR2013HQ008, ZR2019MH020) and the Health Commission of Hubei Province [WJ2019Q060] the young scientist supportive program of Shandong University of Technology. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

GDL and HL conceived and designed the study and wrote the paper. GDL, HL and JZ revised the manuscript. JZ Contributed materials/analysis tools, performed the study, analyzed the data. JZ, JHW, JhW, JMZ, JYW and XZ performed the study. ZFD, HX, HFJ, CYZ, QZ and HL analysis and revised data. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Hong Liu or Guodong Liang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that there are no conflicts of interest.

Additional information

Publisher’s Note

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

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

Zhang, J., Liu, H., Wang, J. et al. Origin and evolution of emerging Liao ning Virus (genus Seadornavirus, family Reoviridae). Virol J 17, 105 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12985-020-01382-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12985-020-01382-2

Keywords