The nucleotide sequence data of the N gene of Lao rabies viruses appear in the DDBJ/EMBL/GenBank nucleotide sequence databases with the accession numbers: AB981665 (strain Lao1), AB981666 (strain Lao3), AB981667 (strain Lao5), AB981668 (strain Lao7), AB981669 (strain Lao8), AB981670 (strain Lao10), AB981671 (strain Lao11), AB981672 (strain Lao12), AB981673 (strain Lao13), AB981674 (strain Lao 17), AB981675 (strain Lao21), AB981676 (strain Lao22) and AB981677 (strain Lao23). The DNA Data Bank of Japan/European Molecular Biology Laboratory/GenBank Accession numbers of each strain described following, Lao2 and Lao4: AB981663 and AB981664.

Background {#sec005}
==========

Globally, an estimated 60,000 people die of rabies annually, and more than 31,000 of these deaths occur in Asia \[[@pntd.0003645.ref001],[@pntd.0003645.ref002]\]. Among various Asian regions, the countries of the Association of South East Asian Nations (ASEAN) have been working towards substantial economic development. A call for elimination of rabies by 2020 highlights the political importance of rabies control \[[@pntd.0003645.ref003]\]. Seven out of the 10 ASEAN member states are rabies-endemic: Cambodia, Indonesia, Laos, Myanmar, the Philippines, Thailand and Vietnam \[[@pntd.0003645.ref003]\]. Among these countries, Laos has the lowest reported number of human deaths from rabies. However, the precise morbidity, mortality and molecular epidemiology of rabies in Laos is largely unknown. This is due to difficulty in collecting data and samples from remote areas of the country, and to the country's modest data collection system at the centralized diagnostic facility. Survey of prevalence of canine rabies is still in its initial stages and limited to a small central part of Laos. The first reported phylogenetic study showed that Laotian rabies viruses are grouped with viruses from Thailand, Myanmar, Cambodia and Vietnam \[[@pntd.0003645.ref004]\]. Combined analyses from geographic information system data showed that China is the likely source of the Asian rabies viruses, whereas individual migration event suggested that Cambodia may be a source of Asian rabies viruses transmission to China, Laos and Thailand \[[@pntd.0003645.ref005]\].

The Laotian government vision is to convert the country from a 'land-locked' to 'land-linked' due to its strategic position in the ASEAN region. As a result, a number of bridges and new roads have been built to connect Laos to neighboring countries. Since ancient times, Laos has relied on its neighbors for trading. For historical, travel convenience and cultural reasons, Thailand has been a primary trading partner. The Mekong River region is well known for its culture of dog-meat consumption \[[@pntd.0003645.ref006]\]. With growing economic prosperity and improved roads across this region, movement of people and allegedly, dogs for meat consumption has been increasing substantially. These factors may have influenced the dissemination and evolution of rabies viruses.

That rabies is a neglected tropical disease is exemplified by the fact that with the exception of Thailand, a full viral genome sequence has been unavailable across the rabies-endemic South East Asian countries. Rabies belongs to the genus *Lyssavirus* of the family *Rhabdoviridae*. It is a single-stranded, negative-sense RNA of approximately 12 kb that encodes five structural proteins: nucleoprotein (N), phosphoprotein (P), matric protein (M), glycoprotein (G), and RNA-dependent RNA polymerase (L). These genes are separated by intergenic regions of variable length. Phylogenetic classification of rabies viruses is based on specific genes, mostly from partial gene sequences, and thus may lead to inconsistent results. It is possible that phylogenetic classification using complete genome sequences, rather than from a partial sequence of a specific gene, would offer a more robust and comprehensive means of addressing the evolution, spread and genome-wide heterogeneity of viruses \[[@pntd.0003645.ref007]\]. Complete genome sequencing of rabies viruses from these countries will also help explain the evolution of the viruses in Laos and the Mekong River region. Because only complete genome sequence can reveal the length, nucleotide substitution patterns, and variations in the start signal, stop signal and other motifs across the genes. As a result it can provide more information to compare and can generate reliable results on genome evolution. Therefore, this study was undertaken to quantify the current animal rabies occurrence in Laos and to complete a molecular characterization of the viruses in current circulation.

Methods {#sec006}
=======

Ethics statement {#sec007}
----------------

The study was approved (No. 38/NECHR) by the National Ethics Committee for Health Research, Ministry of Health, Laos.

Collection of data {#sec008}
------------------

To determine the prevalence of rabies, we analyzed data from samples submitted from various regions of Laos to the rabies laboratory of the National Animal Health Centre, Department of Fisheries and Livestock, Ministry of Agriculture and Forestry from January 2004 through December 2011. This is the only central animal rabies diagnostic laboratory in Laos. In Laos, dog and cat rabies are notifiable, but not wildlife rabies. For animal samples, the head of the suspected animal was submitted to the laboratory, and the brain was dissected by a trained technician. These animals had signs of rabies such as aggressive behavior, a tendency to bite and excessive salivation, with or without a history of biting humans and/or animals. There are no community veterinarians in Laos; samples from suspect animals are submitted mainly by the general public.

Brain sample tests are performed using a fluorescent antibody test (FAT) on crushed smears of hippocampus for a fee of 30,000 Kip (approximately US\$ 3.70) and test results are provided on paper to the individual who submitted the sample. In emergency cases, such as when postexposure prophylaxis (PEP) is required, the person is informed of the results by phone.

Sample collection for molecular characterization of rabies virus {#sec009}
----------------------------------------------------------------

The samples submitted from Vientiane capital, Vientiane and Champasak provinces ([Fig 1](#pntd.0003645.g001){ref-type="fig"}) to the rabies laboratory of National Animal Health Centre, during February 2011 through March 2012, were used in this study for molecular characterization ([S1 Table](#pntd.0003645.s001){ref-type="supplementary-material"}).

![Map of Lao Peoples Democratic Republic.\
Map showing three provinces of Lao from where florescent antibody test positive rabies samples were collected. Samples collected from the districts of Vientiane capital (blue) are as follow: 7 from Xaythany (code VT-3), 3 from Sikhottabong (VT-5), 2 from Sisattanak (VT-7) and Saysettha (VT-8), and 1 from Hatsayfong (VT-9) and Naxaythong (VT-2). One sample was collected from each of the following districts of Champasak (light green): Sanasomeboune (CH-1), Bachieng (CH-2), Pakse (CH-3), Kong (CH-9), Phonthong (CH-5), and Pathoumphone (CH-10). One sample was collected from Meuan district of Vientiane province (green). Each red dot represents one rabies-positive sample. See [S1 Table](#pntd.0003645.s001){ref-type="supplementary-material"} for more detailed information of the samples.](pntd.0003645.g001){#pntd.0003645.g001}

RNA extraction, RT-PCR and nucleotide sequencing {#sec010}
------------------------------------------------

Total RNA was extracted from about 1 g of brain homogenate by using Trizol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instruction. Extracted RNA was stored at -30°C until further analyses.

Using random hexamer primers cDNA was synthesized by SuperScript III First-Strand Synthesis System (Invitrogen) according to the instructions of the manufacturer. The synthesized cDNA was diluted with DNase/RNase free water (Invitrogen) as a template. PCR was performed using TaKaRa ExTaq (Takara Bio Inc., Shiga, Japan). Nested PCR for nucleoprotein (N) gene was performed to detect the presence of rabies virus genome in the samples \[[@pntd.0003645.ref008]\]. The PCR for N and glycoprotein (G) genes, and G-L intergenic region was performed \[[@pntd.0003645.ref007]\]. Whole genome sequence of two strains, arbitrarily selected from Vientiane capital and Champasak samples was performed \[[@pntd.0003645.ref007]\]. For 13 other strains only full-length N gene sequence could be performed. Cycle sequencing of the amplified product was performed using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster city, CA, USA). The purified amplicons were sequenced using ABI-3130 Genetic sequencer (Applied Biosystems). All steps were done according to the manufacturer's instruction. The 5′ or 3′-terminal end of the genome was amplified using SMART RACE cDNA Amplification Kit (Clontech Laboratories, Inc., Mountain View, CA, USA) according to the instructions of the manufacturer.

Time-line evolutionary analysis {#sec011}
-------------------------------

Evolutionary analysis was done using the full-length N gene. We inferred a Maximum Clade Credibility phylogenetic tree using the Bayesian Markov Chain Monte Carlo method available in the BEAST package, v1.6.1 \[[@pntd.0003645.ref009]\]. The analysis utilized a relaxed (uncorrelated lognormal) molecular clock and GTR+Γ+I model of nucleotide substitution. The model was selected on the basis of Akaike Information Criterion using jModelTest software \[[@pntd.0003645.ref010]\]. All chains were run for 60 million generations and sampled every 3000 steps. This resulted in an effective sample size of \>521 for all estimated parameters. The posterior densities were calculated using with 10% burn-in and checked for convergence using Tracer, v.1.5.

Sequence analyses and construction of whole genome phylogenetic tree {#sec012}
--------------------------------------------------------------------

The nucleotide and amino acid sequences of genes and intergenic regions were compared among rabies viruses from Laos. The nucleotide sequences were used to construct the whole genome phylogenetic tree \[[@pntd.0003645.ref007]\]. Multiple sequence alignment was done by ClustalW ver.2 then phytogenic tree was constructed with MEGA ver. 5 using the neighbor joining (NJ) method and the branching pattern was statistically evaluated by bootstrap analysis of 1000 replicates.

Results {#sec013}
=======

Epidemiology of rabies {#sec014}
----------------------

During this eight years period, number of submitted samples peaked in 2008 and gradually decreased ([Table 1](#pntd.0003645.t001){ref-type="table"}). Whereas number of rabies positive samples varied annually, peaked in 2009 and gradually decreased. However the percentage of rabies positive samples increased significantly from 40.5% in 2004 to 60.2% in 2009 and continued at this level.

10.1371/journal.pntd.0003645.t001

###### Yearly distribution of samples submitted from different animals for the laboratory diagnosis of rabies and their proportion positive for rabies.

![](pntd.0003645.t001){#pntd.0003645.t001g}

  **Year**    **Dog**    **Cat**   **Rat**   **Monkey**   **Rabbit**   **Total**
  ----------- ---------- --------- --------- ------------ ------------ ------------------
  **2004**    66/162     0         0/1       0            0            66/163 (40.5%)
  **2005**    84/167     0         0         0            0            84/167 (50.3%)
  **2006**    87/179     0         0/1       0            0            87/180 (48.3%)
  **2007**    64/151     0/1       0         0            0            64/152 (42.1%)
  **2008**    102/198    1/1       0         0            0            103/199 (51.7%)
  **2009**    105/174    3/5       0         1/1          0/1          109/181 (60.2%)
  **2010**    72/117     0/1       0         0            0            72/118 (61.0%)
  **2011**    60/101     0         0         0            0            60/101 (59.4%)
  **Total**   640/1249   4/8       0/2       1/1          0/1          645/1261 (51.1%)

Annually, mean 157.6 {95% confidence interval (CI) 129.9--185.3} samples were submitted for rabies diagnosis, among them 80.6 (95% CI 65.3--95.9) samples were rabies positive i.e. 51.7% (95% CI 45.0--58.4%) of the submitted samples were rabies positive. Most of the submitted samples were from dogs 155.1 (95%CI 129.0--183.2), followed by cats 1.0 (95% CI 0--2.4), and monkeys 0.5 (95% CI 0--1.1). Likewise most of the rabies positive samples were also from dogs 80.0 (95%CI 65.5--94.4), followed by cats 0.5 (95%CI 0--1.4), and other animals 0.13 (95% CI 0--0.4). These represent that 99.2% of the rabies samples were from dogs followed by cats and monkeys (0.6% and 0.3%, respectively).

Province-wise ([Table 2](#pntd.0003645.t002){ref-type="table"}) data analyses showed that yearly submitted samples were mainly from Vientiane capital 115.1 (95% CI 90.2--140.0), followed by Champasak 25.4 (95% CI 14.3--36.4), Vientiane province 11.7 (95% CI 6.6--16.9), and other provinces 5.4 (95% CI 3.1--7.6). Samples from Vientiane capital represents 73.0% of the total submitted samples followed by Champasak (16.1%), Vientiane province (11.7%), and other provinces (3.4%). From other provinces sample number never exceeded more than six and usually only one per year. Samples were never received during this period from six of 18 provinces of Laos.

10.1371/journal.pntd.0003645.t002

###### Yearly distribution of province wise submitted samples for the laboratory diagnosis of rabies.

![](pntd.0003645.t002){#pntd.0003645.t002g}

          Vientiane capital   Vientiane province   Chamapasak   Salavanh   Bolikhamxay   Xekong   Xiengkhouang   Savanhakhet   Luangprabang   Attapeua   Oudomxay   Luangnamtha   Total
  ------- ------------------- -------------------- ------------ ---------- ------------- -------- -------------- ------------- -------------- ---------- ---------- ------------- -------
  2004    104                 24                   24           3          6             1        1                                                                               163
  2005    135                 14                   16                                    1                       1                                                                167
  2006    148                 11                   15                      2             1        1                            1              1                                   180
  2007    122                 13                   14                                    1        1                                                      1                        152
  2008    136                 6                    51           3          1                      1                            1                                                  199
  2009    133                 4                    39                      2             1        2                                                                               181
  2010    74                  14                   26                                    2        2                                                                               118
  2011    69                  8                    18                                    1        1                                                                 4             101
  Total   921                 94                   203          6          11            8        9              1             2              1          1          4             1261

Time-line evolutionary analysis {#sec015}
-------------------------------

The presence of rabies virus was confirmed in 21 of the 23 FAT positive samples by RT-PCR and partial nucleotide sequence of these amplicons. Lao9 and Lao19 samples from Vientiane capital were negative by RT-PCR. The full-length N gene sequence could be done in a total of 15 samples that were used in time-line evolutionary analysis. These samples were from Vientiane capital and Champasak provinces, the full-length N gene of the only sample from Vientiane province could not be amplified. The mean rate of nucleotide substitution estimated for the N gene was 2.5×10^4^ substitutions/site/year (95% HDP values 1.6--3.4×10^4^ substitutions/site/year). This rate is in agreement with previous studies \[[@pntd.0003645.ref011]\].

Phylogenetic tree ([Fig 2](#pntd.0003645.g002){ref-type="fig"}) revealed that rabies strains in Laos belong to the Southeast Asian (SEA) cluster that diverged from the Chinese strains approximately 331.5 years (95% HPD 215.1--481.9 years) ago. From the most recent common ancestor (TMRCA) of SEA cluster, approximately 84.5 years (95% HPD 58.8--120.7 years) ago, i.e., around 1928 (95% HPD range 1891 to 1953), rabies virus of Myanmar segregated. Approximately 62.6 years (95% HPD 45.1--87.4 years) ago, i.e., around 1949 (95% HPD range 1925 to 1967), a cluster only containing Laos rabies viruses segregated from TMRCA of the rest of the SEA cluster which contained viruses from Laos, Cambodia, Vietnam and Thailand. The Lao rabies viruses segregated approximately 44.2 years (95% HPD 29.0--64.1 years) ago, i.e., around 1968 (95% HPD range 1948 to 1983), into two lineages, one containing viruses circulated during 1999--2002 designated as Laos I, and the other contained viruses we identified in 2011 and 2012 designated as Laos II. This lineage had 11 viruses only 2 were from Champasak district and rest from Vientiane capital. Approximately 55.3 years (95% HPD 40.7--76.7 years) ago i. e. in 1957 (95% HPD 1935--1971) Laos III, containing one rabies virus from Vientiane capital and two from Champasak, diverged from rest of SEA cluster which still contained a rabies viruses from Laos. We designated this lineage as Laos IV which segregated from the other lineages of the SEA cluster, only containing rabies viruses from Cambodia, Vietnam and Thailand, approximately 46.4 years (95% HPD 35.6--62.7 years) ago, i.e. in 1966 (95% HPD 1949--1976). All lineages were supported by significant posterior value.

![Time-line evolutionary analysis of rabies virus.\
Bayesian maximum credibility tree representing the genealogy of rabies virus obtained by analyzing the nucleotide sequences of the full N gene sequences (1350 nt). Rabies viruses in Laos may be divided into four lineages Lao I, II, III and IV. Nodes correspond to the mean age at which they are separated from the most recent common ancestor and blue horizontal bars at nodes represent the 95% HPDs of the most recent common ancestor. Numbers at the main nodes represent posterior values. The horizontal axis at the bottom represents time scales in years starting from the year 2012. Each strain name is followed by country of origin, host, year of detection and GenBank accession number. The nucleotide sequence data of the N gene of Lao rabies viruses appear in the DDBJ/EMBL/GenBank nucleotide sequence databases with the accession numbers: AB981665 (strain Lao1), AB981666 (strain Lao3), AB981667 (strain Lao5), AB981668 (strain Lao7), AB981669 (strain Lao8), AB981670 (strain Lao10), AB981671 (strain Lao11), AB981672 (strain Lao12), AB981673 (strain Lao13), AB981674 (strain Lao17), AB981675 (strain Lao21), AB981676 (strain Lao22) and AB981677 (strain Lao23).](pntd.0003645.g002){#pntd.0003645.g002}

Molecular analysis of rabies virus strains {#sec016}
------------------------------------------

We completed the whole genome sequence of two strains, strains Lao2 and Lao4 belong to Lao II and Lao III lineage, respectively. The full genome of both strains was 11924 nt long. At whole genome level there was 97.2% identity between the strains. The N, P, M, G and L genes, and G--L intergenic region were of same length ([Table 3](#pntd.0003645.t003){ref-type="table"}). The nucleotide identity of the genes and their deduced amino acid identity were 96.2--97.6% and 97.6--99.6%, respectively ([Table 3](#pntd.0003645.t003){ref-type="table"}). To further identify the different between Lao II and Lao III lineages, strain Lao2 were compared with strain Lao4 to detect the substitutions in the deduced amino acid sequences ([Table 4](#pntd.0003645.t004){ref-type="table"}). A total of 33 substitutions were identified among them 14 were in various domains of the proteins and 19 were in portion of the protein where no site/domain/region has been identified.

10.1371/journal.pntd.0003645.t003

###### The nucleotide lengths of all genes and their deduced amino acid sequences of rabies viruses from Laos are shown.

![](pntd.0003645.t003){#pntd.0003645.t003g}

  Gene         Nucleotide Length   \% nucleotide identity   Amino Acid Length   \% Amino Acid Identity
  ------------ ------------------- ------------------------ ------------------- ------------------------
  Genome       11924               97.2                                         
  N            1353                97.6                     450                 99.6
  P            894                 97.1                     297                 98.0
  M            639                 96.2                     212                 97.6
  G            1575                97.6                     524                 98.7
  G-L region   516                 95.9                                         
  L            6387                97.4                     2128                99.4

The nucleotide and amino acid identities of different genes between Lao2 and Lao4 are also shown.

10.1371/journal.pntd.0003645.t004

###### Substitutions in genome sequence of rabies virus strains Lao2 and Lao4 from Laos.

![](pntd.0003645.t004){#pntd.0003645.t004g}

  Protein   Lao2        Lao4        Site/domain/region of protein
  --------- ----------- ----------- --------------------------------------------
  N         Ser~84~     Ala~84~     
            Val~281~    Ile~281~    
  P         Val~73~     Ala~73~     N protein binging site, Variable domain I
            Lys~76~     Gln~76~     N protein binging site, Variable domain I
            Ser~141~    Leu~141~    N protein binging site, Variable domain II
            Ala~149~    Ser~149~    N protein binging site, Variable domain II
            Ser~157~    Leu~157~    N protein binging site, Variable domain II
            Ile~257~    Leu~257~    
  M         Gly~9~      Asp~9~      
            Glu~10~     Lys~10~     
            Ser~30~     Phe~30~     
            Ala~34~     Thr~34~     
            Val~105~    Ile~105~    
  G         Ala~249~    Val~249~    Transmembrane domain
            Thr~293~    Ser~293~    Transmembrane domain
            Thr~454~    Met~454~    Cytoplasmic domain
            Thr~458~    Met~458~    Cytoplasmic domain
            Thr~464~    Ala~464~    Cytoplasmic domain
            Gly~468~    Glu~468~    
            Pro~474~    Leu~474~    
  L         Ser~25~     Ala~25~     Conserved domain I
            Lys~381~    Arg~381~    Conserved domain III
            Lys~432~    Arg~432~    Conserved domain III
            Asn~675~    Asp~675~    Conserved domain V
            Asp~684~    Asn~684~    
            Pro~1093~   Ser~1093~   
            Arg~1626~   Lys~1626~   
            Ile~1659~   Val~1659~   
            Glu~1840~   Asp~1840~   
            Val~1861~   Met~1861~   
            Ala~1883~   Thr~1883~   
            Thr~2090~   Ile~2090~   
            Gln~2105~   Arg~2105~   

Whole genome phylogeny {#sec017}
----------------------

Phylogenetic tree constructed by using whole genome showed that the bat origin rabies viruses from Americas formed a separate cluster from rest of the strains ([Fig 3](#pntd.0003645.g003){ref-type="fig"}). Rest of the strains was divided into China, Southeast Asia, arctic/arctic like, Europe-America and Sri Lanka clusters. The China cluster consists of two major sub-clusters China I and China II. All clusters and sub-clusters were supported by bootstrap value of 100. Phylogenetic tree revealed that rabies strains from Laos belong to the SEA cluster with Thai strains. Strain Lao4 from Champasak province was closer to Thai strain than strain Lao2 from Vientiane capital.

![Whole genome phylogenetic tree of rabies virus.\
Phylogenic tree constructed with the nucleotide sequence of the whole genome of strains Lao2 belongs to lineage III and Lao4 belongs to lineage II from Laos PDR and strains from different countries available in the GenBank. Mokola was used as an out-group. Based on the whole genome phylogenetic analysis rabies viruses may be divided into China, Southeast Asia, Arctic/acrtic-like. Euro-America, Sri Lanka and Bat origin clusters. Laos rabies viruses belongs to the Southeast Asian cluster with strains from Thailand. The number adjacent to the node represents the bootstrap value and values lower than 50% have not been indicated. Scale bar shows genetic distance expressed as amino acid substitutions per site. The DNA Data Bank of Japan/European Molecular Biology Laboratory/GenBank Accession numbers of each strain described following, Lao2 and Lao4: AB981663 and AB981664.](pntd.0003645.g003){#pntd.0003645.g003}

Discussion {#sec018}
==========

In Laos, a National Committee for Rabies Control was formed in 1999 and ASEAN has also called for action to eliminate rabies by 2020 from South East Asia \[[@pntd.0003645.ref012]\]. Therefore there have been continued activities such as attempt to increase rabies awareness among people, campaign for animal and human vaccines, and campaign for sterilizing roaming dogs. Whether the gradual increase of rabies samples as noted in this study indicating an increase of rabies in dogs or rabies awareness among the people is not known. If this indicates an increase of rabies then it needs urgent measures to contain the situation before it become big. In Laos annual number of humans with animal bite has increased from 8,277 in 2008 to 14,156 in 2011, (Expanded Program on Immunization unit, Mother and Child Health Centre, Hygiene and Health Promotion Department, Ministry of Health, Lao PDR). Vientiane Capital accounts for 73% of the national data for bites and rabies vaccination in humans, 99% of the bite victims are given postexposure treatment but in the form of vaccination only since rabies immunoglobulin is not available and only 30% receive a full course of vaccine \[[@pntd.0003645.ref012]\]. Dog vaccination coverage was also found to be low, only 11% of the owned dogs were vaccinated for rabies \[[@pntd.0003645.ref013]\]. Other study showed that 71--83% of dog vaccination occurs in the capital, covering 6% of the estimated dog population \[[@pntd.0003645.ref012]\]. All samples used in the present study for genetic characterizations of rabies viruses were from owned dogs but none of them were vaccinated, over that no owner restricted the dog movements at all times through fencing, leashing or confined to indoors. The low rate of vaccination might be multifactorial and one of them may be cost. Government sponsored campaign is done once a year and mainly in the big cities therefore it does not have an impact on the rural population where the burden of rabies is possibly more. Although on the campaign day the vaccine is free but other than that the dog owner has to pay for the vaccination which may cost between 20,000 and 50,000 Kip (approximately between 2.5 and 6.2 US\$). Despite the inadequate level of human and animal vaccination why the reported number of human death from rabies remains at the level of 0--1 per year deserve further investigation. Rabies epidemiology depends on human- and dog-population densities as well as to the cultural and socioeconomic environments \[[@pntd.0003645.ref014]\]. For rabies control public education on dog vaccination and responsible dog ownership is essential.

We found in this study that more than 99% of the submitted samples were from dogs and equal proportion of dog samples were likewise positive for rabies virus indicating that similar to any rabies endemic countries of Asia and Africa, in Laos also the main transmitter of rabies is dogs. However contrary to several studies which found that rabid cats pose a significant public health risk next to dogs \[[@pntd.0003645.ref015]--[@pntd.0003645.ref017]\], in Laos it seems that a relatively lower number of cats are affected by rabies. More importantly the number of samples submitted from other animals and their proportion of rabies virus positive samples were even less than cats. Why such a small number of samples are submitted from animals other than dogs is unknown, as a whole it is thought that less number of samples is received for animal rabies diagnosis in Laos. Multiple factors might be responsible for this, such as the inconveniences associated with transportation of samples from remote areas to district livestock and fisheries offices, people are unaware of rabies in other animals, and the cost of diagnosis. The provincial livestock and fisheries offices receive samples from respective district offices and send samples to the central laboratory in Vientiane capital for the diagnosis of rabies. Any disruption of this transportation process might have an adverse effect on the quality of sample. Therefore for surveillance which can reveal the real situation of rabies in Laos these hindrances should be removed.

Our data analysis revealed that Vientiane capital and Champasak province have the major share of submitted samples. Therefore we performed the genetic characterizations of rabies viruses present in the samples from these areas. It is not known why these two provinces have the major share of submitted samples, we speculate that high population density of these two provinces and ease of transport of submitted samples might be responsible for this. Among all the provinces in Laos the population density is highest in Vientiane capital followed by Champasak province\[[@pntd.0003645.ref018]\], that might be responsible for higher number of dogs in these provinces \[[@pntd.0003645.ref019]\] as a result higher number of samples. The phylogenetic tree of N genes showed that the epidemiology of Laos is complex, a particular province of Laos cannot be represented by a lineage, possibly the geographical distribution of different lineages of rabies viruses in this country is mosaic in nature. This is not unusual, co-circulation of two types of canine rabies virus has been documented in countries of Africa \[[@pntd.0003645.ref004],[@pntd.0003645.ref020]\], six lineages in China \[[@pntd.0003645.ref021],[@pntd.0003645.ref022]\] and Romania \[[@pntd.0003645.ref023]\], and up to 30 different lineages of bat rabies virus in the United States \[[@pntd.0003645.ref024]\]. In smaller countries possibly circulation of one lineage of rabies virus is usual \[[@pntd.0003645.ref017],[@pntd.0003645.ref025]\]. Considering the size of the country and movement of people and dogs to and from the rabies endemic countries, it is likely that that multiple lineages will co-circulate in Laos. Similar conditions exist in neighboring Thailand where six lineages have been identified \[[@pntd.0003645.ref026]\].

Our phylogenetic analyses were in agreement with previous studies \[[@pntd.0003645.ref007],[@pntd.0003645.ref011],[@pntd.0003645.ref022],[@pntd.0003645.ref025]\] and revealed that rabies viruses from Thailand, Malaysia, Cambodia, Laos and Vietnam had originated from a common ancestor. The lineages circulating in Laos are closely related but are not same as in other SEA countries, which are indicating that invasion of rabies viruses in Laos from these countries possibly did not occur recently. The time-line evolutionary analysis showed that among the currently circulating lineages, Laos II lineage is the oldest followed by Laos III and Laos IV. Possibly there were several clonal expansions of the rabies viruses in Laos leading to the formation of lineages and through periodic sweeps, such as Laos I lineage is no longer seen now, shaping up the circulating rabies viruses in Laos. Several amino acid substitutions identified in strains of Lao lineages also supports that genomic modifications occurred after the introduction of rabies virus strains in Laos. The selective advantages resulted by these genomic modifications needs to be addressed in future studies.

A recent study showed that in Vietnam there has been a constant influx of rabies viruses from China or vice versa \[[@pntd.0003645.ref027]\]. Both countries are neighbor of Laos and consume dog meat. Periodic dissemination of viruses from one region to other by regional movement of dogs possibly for meat consumption and subsequent evolution is possibly creating the related lineages of rabies viruses in SEA cluster. Exploring how dog meat consumption is affecting the movements of dogs might shed light on the evolution of rabies viruses in this region.

We did not find any substitutions in the G genes that are responsible for increase the virulence of rabies virus \[[@pntd.0003645.ref028],[@pntd.0003645.ref029]\], though this analysis could be done in two strains only where full genome sequence were done. The phylogenetic tree generated from whole genome sequences clearly showed the distinctness of SEA branch from other branches. Complete genome sequences of rabies viruses from other SEA countries will help to understand the complexity and dissemination of rabies virus in this region.

The limitation of the study is that samples could not be collected from all provinces of Laos yet it includes two of the provinces with highest rabies burden. This work provides insight into the evolution of rabies viruses in a geographic area where genetic analyses are sparsely reported. Small regional studies of the molecular epidemiology of rabies in Laos are of limited interest and will have to be followed after more is known of the prevalence of this disease throughout the country. Data on submitted samples clearly showed that it would take substantial amount of time to collect samples from other parts of Laos. We are collecting samples from other provinces of Laos; in future it might provide more information on the distribution, spread and evolution of rabies viruses in this region.

Conclusions {#sec019}
-----------

There is an increasing trend of rabies virus detection in animal samples in Laos, which may indicate increase of rabies since dog immunization rate is low. An integrated approach is needed such as public education, access to vaccine and rabies test at affordable cost and time. Since dogs are main animal reservoir of rabies virus in Laos, therefore it might be easier to control rabies if appropriate program is introduced. Circulating rabies viruses in Laos are closely related with the rabies viruses from neighboring countries. Possibly there is no recent invasion of rabies viruses in Laos but in the past there was periodic exchange of rabies viruses with neighboring countries. Dog movements through Laos for regional dog meat consumption did not cause any spillover recently. More regional laboratories should be established in different parts of the country for proper surveillance of rabies in Laos.

Supporting Information {#sec020}
======================

###### Information on dog samples collected for the molecular analysis of rabies virus.

Codes of corresponding districts in [Fig 1](#pntd.0003645.g001){ref-type="fig"} are also shown in the district column.

(DOCX)

###### 

Click here for additional data file.

We thank Phetlamphone Meeduangchanh and Khamkhoun Soundara for their technical assistance.

[^1]: The authors have declared that no competing interests exist.

[^2]: Conceived and designed the experiments: KA PL. Performed the experiments: KA PP PV TM PL DM MT BD BK AN. Analyzed the data: KA PP PV TM PL DM MT BD BK AN. Contributed reagents/materials/analysis tools: KA PP PV TM PL DM MT BD BK AN. Wrote the paper: KA PP PV TM PL DM MT BD BK AN.
