BACKGROUND
hantaviruses, which belong to the hantavirus genus in the bunyaviridae family, are tripartite negative-sense rna viruses.
hantaviruses possess a tripartite negative-sense rna genome consisting of the following three segments: the large  segment encodes a viral rna-dependent rna polymerase; the medium-sized  segment encodes two viral glycoproteins ; and the small  segment encodes the viral nucleocapsid protein  .
the gps along with the np determine the virulence and pathogenicity of the hantavirus.
unlike other viruses of the bunyaviridae family, hantaviruses are not transmitted by arthropods; rather, they infect people though the urine, saliva and feces excreted by rodent hosts, especially muroids  <cit> .
hantaviruses only generate transient pathology in rodents, and they do not affect the life span and reproduction of their hosts.
in contrast, they can cause two severe clinical manifestations in humans: hemorrhagic fever with renal syndrome  in the old world and hantavirus cardiopulmonary syndrome  in the new world .
previous studies have indicated that at least  <dig> species and  <dig> genotypes belonging to the hantavirus genus have been isolated worldwide  <cit> .
etiological studies have shown that hfrs that has spread around the world, resulting in the production of variant hantaviruses.
in asia and europe, five types of hantaviruses can cause hfrs: hantaan virus , seoul virus , dobrava virus , saaremaa virus , and puumala virus   <cit> .
in the usa, the sin nombre virus  and the andes virus  are stable viruses that can cause hps.
most hfrs cases occur in europe and east asia   <cit> .
china is the country that is most seriously affected by hantavirus infection worldwide.
previous investigations have shown that hfrs-infected patients in east asian countries, including china, russia and korea, account for at least 90 % of hfrs patients around the world.
at least  <dig>  cases of hfrs are reported annually in china, and more than  <dig> cases are reported in korea and the eastern part of russia  <cit> .
hfrs is caused by different types of hantaviruses in different countries, and phylogenetic analysis suggests that hantaviruses and rodent hosts have coevolved  <cit> .
in china, there are two common types of hantavirus, htnv and seov, which are reportedly carried by a. agrarius and r. norvegicus, respectively, and some hantaviruses isolated from r. norvegicus seemed to be htnv  <cit> .
htnv and seov could cause serious public health problems in china, especially in the heilongjiang area of china.
heilongjiang is located in northeastern china, and it is adjacent to northern russia and the nearby jilin area to the south.
it was the first area in which the etiological agent of hfrs was isolated in china,  <cit>  and the heilongjiang area has remained a high-incidence region  <cit> .
previous studies have shown that seov and htnv are circulating in the heilongjiang area  <cit> .
the port of heilongjiang serves as bridge between countries.
a single commercial port is used for economic trade, the exchange of technology and culture, tourism, immigration and so on.
given that many anthropozoonoses are transmitted by vectors , and these vectors might be introduced into china through ports, it is essential to investigate the molecular epidemiology of hantaviruses by monitoring the rodents that are the natural hosts of hantavirus in the port cities of the heilongjiang area.
in this study,  <dig> port cities have been selected as trapping sites for our hantavirus investigation.
we captured  <dig> rodents and characterized  <dig> hantaviruses isolated from various rodent species in the port cities.
a phylogenetic analysis of the partial m segments indicated that  <dig> of  <dig> viruses belonged to seov, and  <dig> of  <dig> viruses belong to htnv.
entire m segment of htnv  which isolated from r. norvegicus have been sequenced and compared with the reference htnvs.
the results indicated that the spillover of htnv has occurred under natural condition in port cities of heilongjiang area.
methods
rodent collection
in  <dig>  rodents were trapped with rat traps or baited cages in the  <dig> port cities of luobei, jiayin, jiamusi, dongning, hulin, mudanjiang, mishan, harbin, fujin, tongjiang, suifenhe, and raohe .
the trap sites included tussocks, brushwood, residential areas, canals and fields, among other areas.
all of the trapped rodents were identified by zootaxy experts, and they were dissected as soon as they were authenticated.
lung tissues were collected from the captured rodents and transported to the heilongjiang entry-exit inspection and quarantine bureau at –80 °c and stored in a liquid nitrogen tank until further processing.
all procedures used in this study were approved by the institutional animal care and use committee of heilongjiang international travel healthcare center  .fig.
<dig> map showing the trapping sites for rodents in the heilongjiang area, china.
almost all of the port cities are located along the border with russia, except for harbin and jiamusi.
<dig> rodents were captured in luobei,  <dig> in jiayin,  <dig> in jiamusi,  <dig> in dongning,  <dig> in hulin,  <dig> in mudanjiang,  <dig> in mishan,  <dig> in harbin,  <dig> in fujin,  <dig> in tongjiang,  <dig> in suifenhe, and  <dig> in raohe
rna extraction and rt-pcr
the rnaprep pure tissue kit  was used to extract total rna from lung tissues, after the lung tissues were ground into a homogeneous mixture with a tissue lyser .
purified total rna was reverse-transcribed using the superscript first-strand synthesis system  and random primers  to obtain cdna.
the htnv and seov genotypes were amplified via nested pcr.
the initial round of pcr product amplification employed the primer pairs htm1f  and htm1r , and the second round of pcr production used the primer pairs htm2f /htm2r  and sem2f /sem2r .
the initial round of pcr cycling consisted of an initial denaturation at 95 °c for 5 min, followed by 35 cycles of denaturation at 94 °c for 30 s, annealing at 52 °c for 30 s, and elongation at 72 °c for 30 s, with a final extension at 72 °c for 10 min.
the second round of pcr cycling consisted of an initial denaturation at 95 °c for 5 min, followed by 30 cycles of denaturation at 94 °c for 30 s, annealing at 55 °c for 30 s, and elongation at 72 °c for 30 s, with a final extension at 72 °c for 10 min.
the amplified products were approximately 330 bps.
for avoiding possible false positive results caused by rt-pcr assays, a commercial hantavirus detection kit  was also used to confirm the positive results determined by rt-pcr.
the whole m gene of htnv isolate  was further amplified by rt-pcr method as following: the 5′ fragement of m gene was amplified by the primer pairs as : hm1f /hm1r , and the 3′ fragement of m gene was amplified by the primer pairs as : hm2f /hm2r .
the pcr cycling consisted of an initial denaturation at 95 °c for 5 min, followed by 35 cycles of denaturation at 94 °c for 30 s, annealing at 55 °c for 30 s, and elongation at 72 °c for 90 s, with a final extension at 72 °c for 10 min.
the sequence of entire m gene with approximately 1500 bps was obtained by assembling the above two fragements sequences by using dnastar software.
phylogenetic analysis of m gene sequences
the pcr products were purified with the ezna tm gel extraction kit  according to the manufacturer’s instructions and sequenced with the same primers used for pcr amplification.
the nucleic acid sequences of the hantaviruses and those downloaded from genbank were edited and analyzed with the dnastar program .
the identities of the isolated and downloaded hantavirus sequence were also calculated with the dnastar program.
the seqman program was used to edit the nucleic acid sequences.
the mega  <dig>  program was employed to generate a viral phylogenetic tree via the neighbor-joining  method.
the alignment program of dnastar software and gendoc software were used to complete the alignment analyses of amino acid sequence of m segment
RESULTS
descriptions of the sampled rodents
in  <dig>   <dig> rodents were captured from the  <dig> port cities in the heilongjiang area of china.
the rodents included  <dig> specimens of a. agrarius,  <dig> r. norvegicus,  <dig> clethrionomys rufocanus sundevall,  <dig> a. peninsulae,  <dig> clethrionomys rutilus pallas,  <dig> mus musculus,  <dig> microtus fortis buchner,  <dig> r. rattus,  <dig> e. sibiricus,  <dig> c. triton, 3 s. dauricus, and  <dig> sorex araneus linnaeus.
the heilongjiang area harbors a diversity of species, with a. agrarius and r. norvegicus representing the dominant species.
nested pcr targeting partial m segment sequences could be applied to screen the infected rodents in all samples.
the results showed that the infection rate of a. agrarius was  <dig>  %, and the infection rate of r. norvegicus was  <dig>  % .table  <dig> detection of hantavirus in rodent species captured at trap sites in various locations
a. agrarius
r. norvegicus
a. peninsulae
mus musculus
r. rattus
e. sibiricus
c. triton
s. dauricus
number of positive samples/number of tested samples
positive rate
no samples
genetic diversity of seov isolates
phylogenetic tree was constructed with the partial m segment sequences of the viruses .
in the phylogenetic tree analyses of the seov isolates, the calculations showed that all of isolates formed six lineages.
lineages  <dig> to  <dig> were quite closely related to one another.
all of the isolates in lineages  <dig>   <dig>  and  <dig> came from china, while the viruses in lineage  <dig> were isolated from china and neighboring countries  as well as the united states.
in lineage  <dig>  the isolates came from jiamusi, mudanjiang, harbin and raohe.
the isolates show large sequence distances from strains z <dig> and cp <dig>  which were isolated from zhejiang and beijing.
lineage  <dig> consisted of the viruses isolated from dongning and raohe, and their virus sequences were very similar.
all of the isolates were closely related to the vietnam5csg strain previously isolated from small mammals in vietnam, which is far from the heilongjiang area.
within lineage  <dig>  strain l <dig> was isolated from the jiangxi area, while strain k24-v <dig> came from the zhejiang area, and these two areas share borders.
the district from which strains r <dig> and hb <dig> were isolated, in the henan area, is located far from the zhejiang area.
lineage  <dig> consisted of gou <dig>  which was isolated from the zhejiang area, and it appeared to be distinct from the other seo viruses.
lineage  <dig> consisted of uk viruses, with the exception of one brazil strain.fig.
<dig> phylogenetic tree of seoul virus isolates based on the partial sequences of m segments.
the tree was constructed by the neighbour-joining method  of mega software.
the htnvs isolated from r. norvegicus or a. agrarius that were trapped in harbin were designated sjc <dig>  ji <dig> and sjc12; those iaolated from jiamusi, sjms <dig>  sjms <dig>  sjms <dig>  and sjms30; those isolated fromsuifenhe, srh <dig>  srh <dig>  srh <dig> and srha1; that isolated from mudanjiang, smdj24; that iaolated from tongjiang, tjf11; and those isolated from dongning, hdn <dig>  sdn <dig>  sdn <dig>  sdn <dig> and sdn <dig>  the genbank accession numbers of the viruses isolated in this study and the reference htnv stains were listed in table  <dig> in detailedly
genetic diversity of htnv
the genetic diversity of the htn isolates was lower than that of the seo isolates based on the partial m segment sequences.
the analysis of the obtained htnv phylogenetic tree revealed nine lineages from all of the isolated genes .
the partial m segment sequences of htnv that were isolated in this study belonged to only one lineage, lineage  <dig>  this lineage has been designated the far east  lineage, which was identified in patients with serious hfrs.
the viruses isolated from the port cities  in the heilongjiang area appeared to be very closely related to one another.
the bao <dig> strain showed a similar sequence to the viruses isolated from the heilongjiang area and the cjap <dig> strain isolated from the jilin area, which adjoins heilongjiang.
some of the partial m segment sequences of the isolates from port cities were also closely related to htnv strains that were isolated far from heilongjiang, such as the cghu <dig> strain originating from the guizhou area.
lineage  <dig> consisted of several korean viruses and shared a common evolutionary source with lineage  <dig>  lineage  <dig> consisted of the novel hantaan virus type nc <dig>  which was isolated from anhui.
it forms a separate branch in the phylogenetic tree.fig.
<dig> phylogenetic tree  of hantaviruses based on partial m segment sequences of htnvs from china.
the htnv that isolated from r. norvegicus or a. agrarius that were trapped in harbin was designated hji14; those isolated from tongjiang, tjf <dig>  tjf <dig> and tjf3; that isolated from hulin, hhl11; that isolated from mishan, hms1; that isolated from jiayin, hjy10; that isolated from mudanjiang, hmdj6; and those isolated from raohe were designated hrh <dig> and hrh <dig>  the genbank accession numbers of the viruses isolated in this study and the referencial seov stains were showed in table  <dig> in detailedly
phylogenetic analysis of htn, seo and prototype hantaviruses
to understand more about the htnv and seov genetic variants circulating in rodent hosts in the port cities of the heilongjiang area, we employed  <dig> partial m segment sequences of seov  and  <dig> partial m segment sequences of htnv  to construct a phylogenetic tree with the reference sequences hantaan 76- <dig> and seoul 80- <dig> .
the phylogenetic analysis showed that the mutation rates of htnv isolates were higher than the mutation rates of seov isolates, based on the reference hantavirus sequences.fig.
<dig> phylogenetic tree  of hantaviruses based on partial m segment sequences of all hantaviruses isolated in this study.
the analysis was performed with mega software.
the numbers at the nodes are bootstrap confidence levels for  <dig> replicates.
<dig> partial m segment sequences of seov isolates and  <dig> partial m segment sequences of htnv isolates were used to construct this phylogenetic tree with the reference sequences hantaan 76- <dig> and seoul 80- <dig>  the genbank accession numbers of the reference seov and htnv genes were s <dig> and y <dig> as showed in table 2
comparison of amino acid sequences of m segment among htnvs
the nucleotide sequence for m segment of tj strain which isolated from r. norvegicus is 3408 nt, and a  <dig> amino acid sequence could be translated by the nucleotide sequence.
the alignment result in fig.
<dig> showed that there was no deletion or insertion found in the deduced amino acid sequence based on m gene sequence of tj strain.
the identity of the deduced amino acid sequence between tj strain and other htnvs  were more than 95 %.
among all the htnvs, tj strain showed a highest identity with bao  <dig> strain , which was isolated from heilongjiang too.
the previous studies indicated that five n- glycosylation sites in gn  and one n- glycosylation sites in gc  were related with the function of m gene coding proteins  <cit> .
however, all of the n- glycosylation sites were conserved among these different htnvs .fig.
<dig> multiple alignment of the deduced amino acid sequences based on m gene of hantaviruses isolate tj and other reference stains.
the amino acid differences between tj isolate and other strains were shown.
the function related n- glycosylation sites were framed and labled with filled triangles
discussion
the heilongjiang area is the most seriously affected endemic area of hfrs in china  <cit> .
to gain a better understanding of the genetic characteristics of the hantaviruses present in the port cities in the heilongjiang area, we collected  <dig> rodent specimens, and we detected hantaviruses in a. agrarius and r. norvegicus.
the primary rodent hosts in the heilongjiang area were a. agrarius and r. norvegicus for all of  <dig>  among all of the evaluated rodents,  <dig>  % were positive for hantavirus infection.
in addition, we analyzed the genetic evolution of hantaviruses based on all of the positive isolates obtained during this study.
the rodent hosts in a specific area and the epidemic genotypes of hantaviruses exhibit a close relationship to epidemics of hfrs.
this study showed that two hantaviruses  are circulating in the heilongjiang frontier area of northeastern china at present.
the epidemiological investigation showed that a. agrarius and r. norvegicus are the dominant rodent species and carriers of hantavirus in the suburban district and residential areas of the port cities in the heilongjiang area .
although the number of r. norvegicus specimens captured in this study was smaller than the number of a. agrarius, the r. norvegicus infection rate was higher than that of a. agrarius.
r. norvegicus has become an advantageous rodent host for hantavirus, and the seov genotype carried by r. norvegicus has become the primary genotype of the virus.
these results showed that hfrs occurring in the heilongjiang frontier area might be caused primarily by seov.
however, the possibility that other hantaviruses also play important roles in causing hfrs in this area could not be eliminated.
therefore, further studies will be necessary to determine whether other viruses are present in the heilongjiang frontier area, and clinical samples will be needed to clarify the true etiological agents of hfrs in these regions.
molecular epidemiological analysis has indicated that the seovs obtained in this study can be divided into six lineages in phylogenetic trees based on partial m segment sequences  <cit> .
previous studies showed that lineages  <dig>   <dig> and  <dig> were present in northeastern china  <cit> .
however, we only isolated the viruses belonging to lineages  <dig> and  <dig> in the heilongjiang frontier area in this study.
notably, these two virus lineages are distributed widely in china, especially in northeastern china and the middle and lower reaches of the yangtze river  <cit> .
lineage  <dig> includes strains from jilin, zhejiang and jiangxi, while lineage  <dig> includes strains from japan, vietnam and south korea.
we surmised that hantaviruses have spread around the world along with their primary natural reservoir, harbored by r. norvegicus.
moreover, according to fig.
<dig>  we observed that lineages 1—  <dig> exhibit a very close evolutionary relationship to lineage  <dig>  these results indicated that most of the seovs isolated from around the world might have evolved from a single ancestor.
however, a few variants, such as gou <dig>  which was isolated from zhejiang and belongs to lineage  <dig>  show a high degree of variability from most of the other seovs, which may suggest the evolution of an ancient hantavirus.
the first time that htnv was reported was in korea, and it was carried by a. agrarius  <cit> .
phylogenetic analysis showed that htnv can be divided into nine lineages  <cit> .
previous work demonstrated the presence of lineage  <dig> in heilongjiang province, and lineage  <dig>  which originated from korea, is found widely throughout china  <cit> .
several of the htnvs that we detected are similar to the representative strain bao <dig>  strains bao <dig> and cjap <dig>  which belong to lineage  <dig>  and were isolated from heilongjiang and jilin, respectively, are the main pathogens responsible for hfrs in the northeast region of china.
this finding may indicate that this htnv variant is one of the pathogens associated with hfrs in northeastern china at present.
however, we have not found any isolates that are similar to strain 76- <dig>  which belongs to lineage  <dig>  <cit> .
previous studies have indicated that each dominant hantavirus genotype is associated with a specific rodent host, as htnv is carried by a. agrarius, whereas the virus isolated from r. norvegicus has been identified as seov  <cit> .
however, we isolated htnvs from r. norvegicus in this study.
this finding may be explained as a spillover infection, though we have not identified spillover of seov from r. norvegicus to a. agrarius.
furthermore, the entire m segment sequence of htnv isolate  shared high identity with bao  <dig> strain .
it indicated that both tj and bao <dig> iaolates may evolve from the same ancestry.
and we predicted that the transfer of htnv infection host might be the approach employed by hantavirus to adapt to environmental change.
in addition, all of the isolates  obtained in the present study belong to known lineages, and no new lineage of isolates has emerged in the heilongjiang area.
CONCLUSIONS
in this study, we found that htnv and seov are circulating in the heilongjiang frontier area and that seovs are the dominant common strains detected in rodent hosts.
the htnv isolates belong to lineage  <dig>  and the seov isolates belong to lineages  <dig> and  <dig>  all of the lineages are common in the northeast region of china.
although we have not identified a new lineage of hantavirus, regular epidemiological surveillance of local murine colonies is necessary and should be performed intensively.
abbreviations
andvandes virus
dobvdobrava virus
hfrshemorrhagic fever with renal syndrome
hpshantavirus cardiopulmonary syndrome
htnvhantaan virus
puuvpuumala virus
saavsaaremaa virus
seovseoul virus
snvsin nombre virus
suya cao and jian ma are co-first authors.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
ylw and jm proposed the conception of this study and supervised the project.
syc, cc and wdj performed the experiments, syc and jm analyzed the results and wrote this manuscript.
all authors read and approved the final manuscript.
this study was supported by the grant from the fundamental research funds for the central universities , the grant from the heilongjiang provincial natural science foundation , and a grant from the general administration of quality supervision, inspection and quarantine of the people’s republic of china  for scientific research subjects .

 hantaviruses, which belong to the hantavirus genus in the bunyaviridae family, are tripartite negative-sense rna viruses. hantaviruses possess a tripartite negative-sense rna genome consisting of the following three segments: the large segment encodes a viral rna-dependent rna polymerase; the medium-sized segment encodes two viral glycoproteins; and the small segment encodes the viral nucleocapsid protein. the gps along with the np determine the virulence and pathogenicity of the hantavirus. unlike other viruses of the bunyaviridae family, hantaviruses are not transmitted by arthropods; rather, they infect people though the urine, saliva and feces excreted by rodent hosts, especially muroids. hantaviruses only generate transient pathology in rodents, and they do not affect the life span and reproduction of their hosts. in contrast, they can cause two severe clinical manifestations in humans: hemorrhagic fever with renal syndrome in the old world and hantavirus cardiopulmonary syndrome in the new world. previous studies have indicated that at least species and genotypes belonging to the hantavirus genus have been isolated worldwide. etiological studies have shown that hfrs that has spread around the world, resulting in the production of variant hantaviruses. in asia and europe, five types of hantaviruses can cause hfrs: hantaan virus, seoul virus, dobrava virus, saaremaa virus, and puumala virus. in the usa, the sin nombre virus and the andes virus are stable viruses that can cause hps. most hfrs cases occur in europe and east asia. china is the country that is most seriously affected by hantavirus infection worldwide. previous investigations have shown that hfrs-infected patients in east asian countries, including china, russia and korea, account for at least 90 % of hfrs patients around the world. at least cases of hfrs are reported annually in china, and more than cases are reported in korea and the eastern part of russia. hfrs is caused by different types of hantaviruses in different countries, and phylogenetic analysis suggests that hantaviruses and rodent hosts have coevolved. in china, there are two common types of hantavirus, htnv and seov, which are reportedly carried by a. agrarius and r. norvegicus, respectively, and some hantaviruses isolated from r. norvegicus seemed to be htnv. htnv and seov could cause serious public health problems in china, especially in the heilongjiang area of china. heilongjiang is located in northeastern china, and it is adjacent to northern russia and the nearby jilin area to the south. it was the first area in which the etiological agent of hfrs was isolated in china, and the heilongjiang area has remained a high-incidence region. previous studies have shown that seov and htnv are circulating in the heilongjiang area. the port of heilongjiang serves as bridge between countries. a single commercial port is used for economic trade, the exchange of technology and culture, tourism, immigration and so on. given that many anthropozoonoses are transmitted by vectors, and these vectors might be introduced into china through ports, it is essential to investigate the molecular epidemiology of hantaviruses by monitoring the rodents that are the natural hosts of hantavirus in the port cities of the heilongjiang area. in this study, port cities have been selected as trapping sites for our hantavirus investigation. we captured rodents and characterized hantaviruses isolated from various rodent species in the port cities. a phylogenetic analysis of the partial m segments indicated that of viruses belonged to seov, and of viruses belong to htnv. entire m segment of htnv which isolated from r. norvegicus have been sequenced and compared with the reference htnvs. the results indicated that the spillover of htnv has occurred under natural condition in port cities of heilongjiang area. methods rodent collection in rodents were trapped with rat traps or baited cages in the port cities of luobei, jiayin, jiamusi, dongning, hulin, mudanjiang, mishan, harbin, fujin, tongjiang, suifenhe, and raohe. the trap sites included tussocks, brushwood, residential areas, canals and fields, among other areas. all of the trapped rodents were identified by zootaxy experts, and they were dissected as soon as they were authenticated. lung tissues were collected from the captured rodents and transported to the heilongjiang entry-exit inspection and quarantine bureau at –80 °c and stored in a liquid nitrogen tank until further processing. all procedures used in this study were approved by the institutional animal care and use committee of heilongjiang international travel healthcare center.fig. map showing the trapping sites for rodents in the heilongjiang area, china. almost all of the port cities are located along the border with russia, except for harbin and jiamusi. rodents were captured in luobei, in jiayin, in jiamusi, in dongning, in hulin, in mudanjiang, in mishan, in harbin, in fujin, in tongjiang, in suifenhe, and in raohe rna extraction and rt-pcr the rnaprep pure tissue kit was used to extract total rna from lung tissues, after the lung tissues were ground into a homogeneous mixture with a tissue lyser. purified total rna was reverse-transcribed using the superscript first-strand synthesis system and random primers to obtain cdna. the htnv and seov genotypes were amplified via nested pcr. the initial round of pcr product amplification employed the primer pairs htm1f and htm1r, and the second round of pcr production used the primer pairs htm2f /htm2r and sem2f /sem2r. the initial round of pcr cycling consisted of an initial denaturation at 95 °c for 5 min, followed by 35 cycles of denaturation at 94 °c for 30 s, annealing at 52 °c for 30 s, and elongation at 72 °c for 30 s, with a final extension at 72 °c for 10 min. the second round of pcr cycling consisted of an initial denaturation at 95 °c for 5 min, followed by 30 cycles of denaturation at 94 °c for 30 s, annealing at 55 °c for 30 s, and elongation at 72 °c for 30 s, with a final extension at 72 °c for 10 min. the amplified products were approximately 330 bps. for avoiding possible false positive results caused by rt-pcr assays, a commercial hantavirus detection kit was also used to confirm the positive results determined by rt-pcr. the whole m gene of htnv isolate was further amplified by rt-pcr method as following: the 5′ fragement of m gene was amplified by the primer pairs as: hm1f /hm1r, and the 3′ fragement of m gene was amplified by the primer pairs as: hm2f /hm2r. the pcr cycling consisted of an initial denaturation at 95 °c for 5 min, followed by 35 cycles of denaturation at 94 °c for 30 s, annealing at 55 °c for 30 s, and elongation at 72 °c for 90 s, with a final extension at 72 °c for 10 min. the sequence of entire m gene with approximately 1500 bps was obtained by assembling the above two fragements sequences by using dnastar software. phylogenetic analysis of m gene sequences the pcr products were purified with the ezna tm gel extraction kit according to the manufacturer’s instructions and sequenced with the same primers used for pcr amplification. the nucleic acid sequences of the hantaviruses and those downloaded from genbank were edited and analyzed with the dnastar program. the identities of the isolated and downloaded hantavirus sequence were also calculated with the dnastar program. the seqman program was used to edit the nucleic acid sequences. the mega program was employed to generate a viral phylogenetic tree via the neighbor-joining method. the alignment program of dnastar software and gendoc software were used to complete the alignment analyses of amino acid sequence of m segment  descriptions of the sampled rodents in rodents were captured from the port cities in the heilongjiang area of china. the rodents included specimens of a. agrarius, r. norvegicus, clethrionomys rufocanus sundevall, a. peninsulae, clethrionomys rutilus pallas, mus musculus, microtus fortis buchner, r. rattus, e. sibiricus, c. triton, 3 s. dauricus, and sorex araneus linnaeus. the heilongjiang area harbors a diversity of species, with a. agrarius and r. norvegicus representing the dominant species. nested pcr targeting partial m segment sequences could be applied to screen the infected rodents in all samples. the results showed that the infection rate of a. agrarius was %, and the infection rate of r. norvegicus was %.table detection of hantavirus in rodent species captured at trap sites in various locations a. agrarius r. norvegicus a. peninsulae mus musculus r. rattus e. sibiricus c. triton s. dauricus number of positive samples/number of tested samples positive rate no samples genetic diversity of seov isolates phylogenetic tree was constructed with the partial m segment sequences of the viruses. in the phylogenetic tree analyses of the seov isolates, the calculations showed that all of isolates formed six lineages. lineages to were quite closely related to one another. all of the isolates in lineages and came from china, while the viruses in lineage were isolated from china and neighboring countries as well as the united states. in lineage the isolates came from jiamusi, mudanjiang, harbin and raohe. the isolates show large sequence distances from strains z and cp which were isolated from zhejiang and beijing. lineage consisted of the viruses isolated from dongning and raohe, and their virus sequences were very similar. all of the isolates were closely related to the vietnam5csg strain previously isolated from small mammals in vietnam, which is far from the heilongjiang area. within lineage strain l was isolated from the jiangxi area, while strain k24-v came from the zhejiang area, and these two areas share borders. the district from which strains r and hb were isolated, in the henan area, is located far from the zhejiang area. lineage consisted of gou which was isolated from the zhejiang area, and it appeared to be distinct from the other seo viruses. lineage consisted of uk viruses, with the exception of one brazil strain.fig. phylogenetic tree of seoul virus isolates based on the partial sequences of m segments. the tree was constructed by the neighbour-joining method of mega software. the htnvs isolated from r. norvegicus or a. agrarius that were trapped in harbin were designated sjc ji and sjc12; those iaolated from jiamusi, sjms sjms sjms and sjms30; those isolated fromsuifenhe, srh srh srh and srha1; that isolated from mudanjiang, smdj24; that iaolated from tongjiang, tjf11; and those isolated from dongning, hdn sdn sdn sdn and sdn the genbank accession numbers of the viruses isolated in this study and the reference htnv stains were listed in table  in detailedly genetic diversity of htnv the genetic diversity of the htn isolates was lower than that of the seo isolates based on the partial m segment sequences. the analysis of the obtained htnv phylogenetic tree revealed nine lineages from all of the isolated genes. the partial m segment sequences of htnv that were isolated in this study belonged to only one lineage, lineage this lineage has been designated the far east lineage, which was identified in patients with serious hfrs. the viruses isolated from the port cities in the heilongjiang area appeared to be very closely related to one another. the bao strain showed a similar sequence to the viruses isolated from the heilongjiang area and the cjap strain isolated from the jilin area, which adjoins heilongjiang. some of the partial m segment sequences of the isolates from port cities were also closely related to htnv strains that were isolated far from heilongjiang, such as the cghu strain originating from the guizhou area. lineage consisted of several korean viruses and shared a common evolutionary source with lineage lineage consisted of the novel hantaan virus type nc which was isolated from anhui. it forms a separate branch in the phylogenetic tree.fig. phylogenetic tree of hantaviruses based on partial m segment sequences of htnvs from china. the htnv that isolated from r. norvegicus or a. agrarius that were trapped in harbin was designated hji14; those isolated from tongjiang, tjf tjf and tjf3; that isolated from hulin, hhl11; that isolated from mishan, hms1; that isolated from jiayin, hjy10; that isolated from mudanjiang, hmdj6; and those isolated from raohe were designated hrh and hrh the genbank accession numbers of the viruses isolated in this study and the referencial seov stains were showed in table  in detailedly phylogenetic analysis of htn, seo and prototype hantaviruses to understand more about the htnv and seov genetic variants circulating in rodent hosts in the port cities of the heilongjiang area, we employed partial m segment sequences of seov and partial m segment sequences of htnv to construct a phylogenetic tree with the reference sequences hantaan 76- and seoul 80-. the phylogenetic analysis showed that the mutation rates of htnv isolates were higher than the mutation rates of seov isolates, based on the reference hantavirus sequences.fig. phylogenetic tree of hantaviruses based on partial m segment sequences of all hantaviruses isolated in this study. the analysis was performed with mega software. the numbers at the nodes are bootstrap confidence levels for replicates. partial m segment sequences of seov isolates and partial m segment sequences of htnv isolates were used to construct this phylogenetic tree with the reference sequences hantaan 76- and seoul 80- the genbank accession numbers of the reference seov and htnv genes were s and y as showed in table 2 comparison of amino acid sequences of m segment among htnvs the nucleotide sequence for m segment of tj strain which isolated from r. norvegicus is 3408 nt, and a amino acid sequence could be translated by the nucleotide sequence. the alignment result in fig.  showed that there was no deletion or insertion found in the deduced amino acid sequence based on m gene sequence of tj strain. the identity of the deduced amino acid sequence between tj strain and other htnvs were more than 95 %. among all the htnvs, tj strain showed a highest identity with bao strain, which was isolated from heilongjiang too. the previous studies indicated that five n- glycosylation sites in gn and one n- glycosylation sites in gc were related with the function of m gene coding proteins. however, all of the n- glycosylation sites were conserved among these different htnvs.fig. multiple alignment of the deduced amino acid sequences based on m gene of hantaviruses isolate tj and other reference stains. the amino acid differences between tj isolate and other strains were shown. the function related n- glycosylation sites were framed and labled with filled triangles discussion the heilongjiang area is the most seriously affected endemic area of hfrs in china. to gain a better understanding of the genetic characteristics of the hantaviruses present in the port cities in the heilongjiang area, we collected rodent specimens, and we detected hantaviruses in a. agrarius and r. norvegicus. the primary rodent hosts in the heilongjiang area were a. agrarius and r. norvegicus for all of among all of the evaluated rodents, % were positive for hantavirus infection. in addition, we analyzed the genetic evolution of hantaviruses based on all of the positive isolates obtained during this study. the rodent hosts in a specific area and the epidemic genotypes of hantaviruses exhibit a close relationship to epidemics of hfrs. this study showed that two hantaviruses are circulating in the heilongjiang frontier area of northeastern china at present. the epidemiological investigation showed that a. agrarius and r. norvegicus are the dominant rodent species and carriers of hantavirus in the suburban district and residential areas of the port cities in the heilongjiang area. although the number of r. norvegicus specimens captured in this study was smaller than the number of a. agrarius, the r. norvegicus infection rate was higher than that of a. agrarius. r. norvegicus has become an advantageous rodent host for hantavirus, and the seov genotype carried by r. norvegicus has become the primary genotype of the virus. these results showed that hfrs occurring in the heilongjiang frontier area might be caused primarily by seov. however, the possibility that other hantaviruses also play important roles in causing hfrs in this area could not be eliminated. therefore, further studies will be necessary to determine whether other viruses are present in the heilongjiang frontier area, and clinical samples will be needed to clarify the true etiological agents of hfrs in these regions. molecular epidemiological analysis has indicated that the seovs obtained in this study can be divided into six lineages in phylogenetic trees based on partial m segment sequences. previous studies showed that lineages and were present in northeastern china. however, we only isolated the viruses belonging to lineages and in the heilongjiang frontier area in this study. notably, these two virus lineages are distributed widely in china, especially in northeastern china and the middle and lower reaches of the yangtze river. lineage includes strains from jilin, zhejiang and jiangxi, while lineage includes strains from japan, vietnam and south korea. we surmised that hantaviruses have spread around the world along with their primary natural reservoir, harbored by r. norvegicus. moreover, according to fig.  we observed that lineages 1—  exhibit a very close evolutionary relationship to lineage these results indicated that most of the seovs isolated from around the world might have evolved from a single ancestor. however, a few variants, such as gou which was isolated from zhejiang and belongs to lineage show a high degree of variability from most of the other seovs, which may suggest the evolution of an ancient hantavirus. the first time that htnv was reported was in korea, and it was carried by a. agrarius. phylogenetic analysis showed that htnv can be divided into nine lineages. previous work demonstrated the presence of lineage in heilongjiang province, and lineage which originated from korea, is found widely throughout china. several of the htnvs that we detected are similar to the representative strain bao strains bao and cjap which belong to lineage and were isolated from heilongjiang and jilin, respectively, are the main pathogens responsible for hfrs in the northeast region of china. this finding may indicate that this htnv variant is one of the pathogens associated with hfrs in northeastern china at present. however, we have not found any isolates that are similar to strain 76- which belongs to lineage. previous studies have indicated that each dominant hantavirus genotype is associated with a specific rodent host, as htnv is carried by a. agrarius, whereas the virus isolated from r. norvegicus has been identified as seov. however, we isolated htnvs from r. norvegicus in this study. this finding may be explained as a spillover infection, though we have not identified spillover of seov from r. norvegicus to a. agrarius. furthermore, the entire m segment sequence of htnv isolate shared high identity with bao strain. it indicated that both tj and bao iaolates may evolve from the same ancestry. and we predicted that the transfer of htnv infection host might be the approach employed by hantavirus to adapt to environmental change. in addition, all of the isolates obtained in the present study belong to known lineages, and no new lineage of isolates has emerged in the heilongjiang area.  in this study, we found that htnv and seov are circulating in the heilongjiang frontier area and that seovs are the dominant common strains detected in rodent hosts. the htnv isolates belong to lineage and the seov isolates belong to lineages and all of the lineages are common in the northeast region of china. although we have not identified a new lineage of hantavirus, regular epidemiological surveillance of local murine colonies is necessary and should be performed intensively. abbreviations andvandes virus dobvdobrava virus hfrshemorrhagic fever with renal syndrome hpshantavirus cardiopulmonary syndrome htnvhantaan virus puuvpuumala virus saavsaaremaa virus seovseoul virus snvsin nombre virus suya cao and jian ma are co-first authors. competing interests the authors declare that they have no competing interests. authors’ contributions ylw and jm proposed the conception of this study and supervised the project. syc, cc and wdj performed the experiments, syc and jm analyzed the results and wrote this manuscript. all authors read and approved the final manuscript. this study was supported by the grant from the fundamental research funds for the central universities, the grant from the heilongjiang provincial natural science foundation, and a grant from the general administration of quality supervision, inspection and quarantine of the people’s republic of china for scientific research subjects.
it can infect humans through contaminated rodent excreta and causes two types of fatal human diseases: hemorrhagic fever with renal syndrome and hantavirus pulmonary syndrome. to obtain additional insights into the genetic characteristics of hantaviruses in the port cities of the heilongjiang area in china, a molecular epidemiological investigation of hantaviruses isolated from rodents was performed in a total of rodents were caught in port cities in heilongjiang. among these rodents, the most common species was a. agrarius, and the second-most common was r. norvegicus. these results suggested that htnv and seov are circulating in a. agrarius and r. norvegicus in the port cities in the area of heilongjiang, but seov may be the dominant common hantavirus. a viral gene pcr assay revealed the presence of two specific genotypes of hantavirus, referred to as hantaan virus and seoul virus, and the positive seov infection rate was higher than that for htnv. hantavirusseovhtnvepidemiologyheilongjiangthe fundamental research funds for the central universitiesdl13ca04wang yulong heilongjiang provincial natural science foundationc201116wang yulong general administration of quality supervision, inspection and quarantine of the people’s republic of china for scientific research subjects2014ik407ju wendong issue-copyright-statement© the author 2016 china exhibits the highest hfrs occurrence rate in the world, and the heilongjiang area is one of the most severely infected regions. a genetic analysis based on partial m segment sequences indicated that all of the isolates belonging to seov could be assigned to two genetic lineages, whereas the isolate belonging to htnv could be assigned to only one genetic lineage.
 hantavirus is a tripartite negative-sense rna virus. it can infect humans through contaminated rodent excreta and causes two types of fatal human diseases: hemorrhagic fever with renal syndrome and hantavirus pulmonary syndrome. china exhibits the highest hfrs occurrence rate in the world, and the heilongjiang area is one of the most severely infected regions.  to obtain additional insights into the genetic characteristics of hantaviruses in the port cities of the heilongjiang area in china, a molecular epidemiological investigation of hantaviruses isolated from rodents was performed in a total of rodents were caught in port cities in heilongjiang. among these rodents, the most common species was a. agrarius, and the second-most common was r. norvegicus. a viral gene pcr assay revealed the presence of two specific genotypes of hantavirus, referred to as hantaan virus and seoul virus, and the positive seov infection rate was higher than that for htnv. a genetic analysis based on partial m segment sequences indicated that all of the isolates belonging to seov could be assigned to two genetic lineages, whereas the isolate belonging to htnv could be assigned to only one genetic lineage.  these results suggested that htnv and seov are circulating in a. agrarius and r. norvegicus in the port cities in the area of heilongjiang, but seov may be the dominant common hantavirus. keywords hantavirusseovhtnvepidemiologyheilongjiangthe fundamental research funds for the central universitiesdl13ca04wang yulong heilongjiang provincial natural science foundationc201116wang yulong general administration of quality supervision, inspection and quarantine of the people’s republic of china for scientific research subjects2014ik407ju wendong issue-copyright-statement© the author 2016 