Thank you for visiting nature.com. You are using a browser version with
limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off
compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site
without styles and JavaScript.

Abstract

Genetically distinct hantaviruses have been identified in five species of fossorial moles (order Eulipotyphla, family Talpidae) from Eurasia and North America. Here, we report the isolation and partial characterization of a highly divergent hantavirus, named Nova virus (NVAV), from lung tissue of a European mole (Talpa europaea), captured in central Poland in August 2013. Typical hantavirus-like particles, measuring 80–120 nm in diameter, were found in NVAV-infected Vero E6 cells by transmission electron microscopy. Whole-genome sequences of the isolate, designated NVAV strain Te34, were identical to that amplified from the original lung tissue, and phylogenetic analysis of the full-length L, M and S segments, using maximum-likelihood and Bayesian methods, showed that NVAV was most closely related to hantaviruses harbored by insectivorous bats, consistent with an ancient evolutionary origin. Infant Swiss Webster mice, inoculated with NVAV by the intraperitoneal route, developed weight loss and hyperactivity, beginning at 16 days, followed by hind-limb paralysis and death. High NVAV RNA copies were detected in lung, liver, kidney, spleen and brain by quantitative real-time RT-PCR. Neuropathological examination showed astrocytic and microglial activation and neuronal loss. The first mole-borne hantavirus isolate will facilitate long-overdue studies on its infectivity and pathogenic potential in humans.

Introduction

A paradigm-disruptive shift in hantavirus host diversity has emerged from the recent discovery of genetically distinct hantaviruses (family Bunyaviridae, genus Hantavirus) in multiple species of shrews and moles (order Eulipotyphla)1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19 and insectivorous bats (order Chiroptera)20,21,22,23,24 in widely separated geographic regions in Europe, Asia, Africa and North America. Moreover, phylogenetic analyses, based on partial and full-length hantavirus genomes, suggest that ancestral shrews, moles and/or bats probably predated rodents as the reservoir hosts of primordial hantaviruses25,26,27. However, despite being referred to as novel viruses, nearly all of the 36 newfound non-rodent-borne hantaviruses have not been isolated in cell culture or animals and exist only as viral sequences. As such, almost nothing is known about their biology and transmission dynamics, as well as infectivity and pathogenic potential in humans.

Hantaviruses possess a negative-sense, tripartite RNA genome, designated as large (L), medium (M) and small (S) segments, which encode an RNA-dependent RNA polymerase (RdRp), envelope glycoproteins (Gn and Gc) and a nucleocapsid (N) protein, respectively28. Some hantaviruses cause diseases of varying clinical severity, known as hemorrhagic fever with renal syndrome (HFRS) and hantavirus pulmonary syndrome (HPS)29. Until recently, rodents (order Rodentia) of the Muridae and Cricetidae families were believed to serve as the exclusive reservoir hosts. The recent flurry of descriptions of newfound non-rodent-borne hantaviruses was presaged by reports, published two to three decades ago which went largely unnoticed, of HFRS antigens in tissues of the Eurasian common shrew (Sorex araneus), Eurasian water shrew (Neomys fodiens) and European mole (Talpa europaea), captured in Russia and the former Yugoslavia30,31,32.

Among the five hantaviruses discovered to date in moles (family Talpidae) from Europe, Asia and North America, Nova virus (NVAV), originally detected in archival liver tissue of a single European mole (Talpa europaea) captured in 1999 in Zala County, Hungary18, and subsequently shown to be widespread in France33 and Poland34, is the most divergent25,26. Here, we report the isolation and partial characterization of NVAV. Since hantaviruses are notoriously difficult to isolate, the long-awaited isolation of NVAV in cell culture is a noteworthy milestone that will accelerate acquisition of new knowledge about its transmission dynamics and facilitate long-overdue studies to ascertain its infectivity and pathogenic potential in humans.

Results

Prevalence of NVAV infection

Of 22 European moles (Fig. 1A), captured in Huta Dłutowska between June 21, 2013 and August 26, 2013, NVAV RNA was detected by RT-PCR in lung tissue from 11 (50%). Based on the intensity of the RT-PCR amplicon on ethidium bromide-stained gels, lung tissue homogenates were prepared from four European moles and used for virus isolation attempts.

Isolation and propagation of NVAV

Initially, NVAV RNA was detected in Vero E6 cells at 17 days following inoculation with a 1% lung tissue homogenate from one of four European moles. Viral RNA was also detected at day 14 of the first passage, but the virus isolate was lost on subsequent passages. After several other failed attempts, NVAV RNA was detected in cells and culture media at 34 days and on subsequent inoculation of fresh Vero E6 cells with culture media, indicating virus replication. At 120 days, the NVAV infectivity in culture media was 8 × 103 plaque-forming units (pfu)/mL. Like all other rodent- and shrew-borne hantaviruses isolated to date, the NVAV isolate did not produce cytopathic effect (or CPE) in Vero E6 cells. Typical hantavirus-like particles, measuring 80–120 nm in diameter, were found by thin-section transmission electron microscopy (Fig. 1B). No other virus-like particles were observed.

Sequence analysis of NVAV

The full-length S-, M- and L-genomic segments of the cell culture grown NVAV isolate, designated NVAV strain Te34, and the full-length S segment of the NVAV from the original lung tissue, were sequenced. While the overall genomic organization of NVAV was indistinguishable from other hantaviruses, the deduced-amino acid sequence similarities were significantly different from that of rodent-, shrew- and mole-borne hantaviruses, differing by more than 50% for the N protein and Gn/Gc envelope glycoproteins and 40% for the RdRp (Table 1).

The full-length, 1,825-nucleotide S-segment sequence of the NVAV isolate and the NVAV in the original lung tissue were identical, containing a single open reading frame with a predicted N protein of 428 amino acids (nucleotide positions 53 to 1,339), and 52- and 486-nucleotide 5′– and 3′–noncoding regions, respectively. Comparison of the S segment between the NVAV isolate and the prototype NVAV strain MSB95703 from Hungary showed nucleotide and amino acid sequence similarities of 85.6% and 96.5%, respectively. The secondary structure of the predicted N protein of NVAV strain Te34 was virtually identical to that reported earlier for the prototype strain (data not shown).

The entire M-genomic segment of NVAV strain Te34 was 3,590 nucleotides, with a predicted glycoprotein precursor of 1,127 amino acids. The highly conserved WAASA amino acid motif of the M segment was found at amino acid positions 641–645, and the NVAV glycoprotein had six potential N-linked glycosylation sites (five on Gn at amino acid positions 101, 133, 344, 396 and 539, and one on Gc at position 924). The very low degree of amino acid sequence similarity between the GnGc of NVAV and other hantaviruses would predict the absence of antibody cross-neutralization between NVAV and hantaviruses harbored by rodents, shrews, moles and bats (Table 1).

The full-length L-genomic segment of NVAV strain Te34 was 6,563 nucleotides, with a predicted RdRp of 2,157 amino acids, starting at nucleotide position 34 and including 56 nucleotides of the 3′–noncoding region. As with all other hantaviruses, the six conserved motifs (premotif A and motifs A, B, C, D and E) were found. Premotif A had a conserved lysine and two arginine residues. Motifs A, B and D had conserved aspartate, glycine and lysine, respectively. In motif C, there were two conserved aspartic acid residues. The XDD motif, essential for catalytic activity, and motif E, containing the E(F/Y)XS site, were also present. The higher sequence similarity of approximately 60% or more in the L segment probably signified the functional constraints of the RdRp (Table 1).

Phylogenetic analysis of NVAV

As determined by maximum-likelihood and Bayesian methods, phylogenetic analysis, based on full-length sequences of each genomic segment, demonstrated tree topologies, well supported by bootstrap analysis and posterior node probabilities, showing the NVAV isolate forming a distinctly divergent clade, which comprised other NVAV strains from Hungary, France and Poland (Fig. 2). Similar topologies were found for deduced amino acid sequences of the encoded proteins (data not shown). Collectively, these data strongly supported an ancient hantavirus-host relationship without evidence of host switching.

Figure 2: Phylogenetic trees generated by the Bayesian method, using the GTR+I+Γ model of evolution as estimated from the data, based on the alignment of the coding regions of the full- length S, M and L segments of NVAV strain Te34.

Experimental NVAV infection

Of 20 infant Swiss Webster mice, inoculated with 200 pfu of NVAV strain Te34 by the intraperitorneal route, all developed weight loss and hyperactivity, beginning at 16 days. In 11 mice, disease progressed to include hind-limb paralysis, during the ensuing four-day period (Fig. 3A). Moribund mice, with clinical scores of 5, were subsequently euthanized, and their sera had antibodies against NVAV, as determined by the indirect immunofluorescent antibody (IFA) test (data not shown). Using conventional RT-PCR and quantitative real-time RT-PCR, NVAV RNA was detected in lung, liver, kidney, spleen and brain from experimentally infected mice. NVAV RNA copies in the brain were 10- to 100-fold higher than in other tissues (Fig. 3B).

Moreover, brain tissue homogenates from NVAV-infected moribund mice, when inoculated intraperitoneally into new litters of infant mice, produced lethal meningoencephalitis, which was indistinguishable to that observed in infant mice experimentally infected with the NVAV cell-culture isolate. In addition, NVAV sequences amplified from tissues of infected mice were identical to that of the NVAV cell-culture isolate used for experimental infection, as well as to the NVAV originally detected in the European mole lung tissue.

As assessed by light-microscopic examination of hematoxylin and eosin-stained tissues and by immunohistochemical staining, pathological abnormalities in NVAV-infected infant mice were largely confined to the brain. Perivascular and subpial lymphocytic infiltrates were found, particularly in the cerebellum. Minimal focal loss of neurons was evident in the cerebral cortex, using a neuronal nuclei-specific antibody (NeuN) (Fig. 4A). Microglial activation and astrocytosis, as evidenced by staining for a microglia-specific calcium-binding protein (Iba1) and glial fibrillary acidic protein (GFAP), respectively, were also prominent (Fig. 4B,C). In the lung, perivascular lymphocytic infiltrates were conspicuous with mild widening of the alveolar septa and interstitial infiltration, but without evidence of hyaline membranes or intra-alveolar edema (Fig. 4D).

Figure 4: Histopathology in infant mice, 17 days after inoculation by the intraperitoneal route with 200 pfu of NVAV strain Te34.

Materials and Methods

Trapping and specimen processing

Moles were trapped using PVC pipes, equipped with flat aluminum latches at both ends, in Huta Dłutowska (N51°35′49.51, E19°22′46.80), near Łódź in central Poland, during June to August 2013. Live-caught moles were euthanized by cervical dislocation, and stored at 4 °C or −20 °C for several hours or days before harvesting of lung tissues, which were stored at −80 °C until used for RT-PCR and virus isolation.

Ethics statement

Field procedures and protocols, including trapping, euthanasia and tissue processing, were approved by the Łódź Ethical Committee on Animal Testing (14/LB/511/2010 and 29/LB/548/2011) and the Directorate General for Environmental Protection (DOP-OZGiZ.4200/N2732/10/JRO, DOP- OZGiZ.6401.05.25.2011kp.3 and DOP-OZGiZ.6401.05.28.2011kp.1), as well as by the Łódź Regional Directorate for Environmental Protection (WPN-I6631.2010.MS and WPN-L6400.59.2011.MS), in accordance to guidelines of the American Society of Mammalogists35. For experimental NVAV infection in Swiss Webster mice, approval was received from the University of Hawaii Institutional Animal Care and Use Committee (protocol 06-055), which adhered to guidelines and regulations of the U.S. Public Health Service Policy on Humane Care and Use of Laboratory Animals and the American Veterinary Medical Association. All experiments involving mice were conducted in consultation with veterinary and animal care staff of the University of Hawaii animal biosafety level-3 laboratory.

RNA extraction, cDNA synthesis and RT-PCR amplification

As an initial screening for NVAV infection, total RNA was extracted from 20–50 mg of mole lung tissues, using the PureLink Micro-to-Midi total RNA purification kit (Invitrogen, San Diego, CA), and cDNA was synthesized using the SuperScript III First-Strand Synthesis Systems (Invitrogen), then analyzed for NVAV RNA by RT-PCR, using NVAV-specific oligonucleotide primers33,34.

Virus isolation and propagation

Attempts to isolate hantavirus from lung tissues of moles with genetic evidence of NVAV infection were made in Vero E6 cells (CRL 1586, American Type Culture Collection, Manassas, VA), using previously described methods6. Briefly, subconfluent monolayers of Vero E6 cells, grown in 25-cm2 flasks, were inoculated with 1% (w/v) suspensions of lung tissue homogenates, prepared in Dulbecco’s minimum essential medium (DMEM) containing 30% fetal bovine serum, from four NVAV RNA-positive European moles. Cells were sub-cultured at two- to five-week intervals, at which time aliquots of cells and supernatants were examined for NVAV RNA by RT-PCR. Supernatant from NVAV-infected cultures were then used to inoculate fresh Vero E6 cells to demonstrate virus replication and to determine infectivity titers by plaque assay6.

Whole genome sequencing

Oligonucleotide primers were designed using the MegAlign Clustal W program (DNASTAR Inc., Madison, WI) to obtain the whole genome of NVAV strain Te34 from cDNA prepared from total RNA, extracted from NVAV-infected Vero E6 cells and from lung tissue of the original wild-caught European mole (Table 2). The 5′– and 3′–termini of each segment were amplified using the 3′–Full RACE Core Set (Takara Bio Inc., Otsu, Japan). Nested or hemi-nested PCR was performed in 20 μL reaction mixtures, containing 250 μM dNTP, 2.5 mM MgCl2, 1 U of Takara LA Taq polymerase (Takara) and 0.25 μM of each primer (Table 2). Initial denaturation at 94 °C for 2 min was followed by two cycles each of denaturation at 94 °C for 30 sec, two-degree step-down annealing from 46 °C to 38 °C for 40 sec, and elongation at 72 °C for 1 min, then 30 cycles of denaturation at 94 °C for 30 sec, annealing at 42 °C for 40 sec, and elongation at 72 °C for 1 min, in a GeneAmp PCR 9700 thermal cycler (Perkin-Elmer, Waltham, MA)34. PCR products were separated, using MobiSpin S-400 spin columns (MoBiTec, Goettingen, Germany), and amplicons were sequenced directly using an ABI Prism 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA)34.

Phylogenetic analysis

To assess the evolutionary relationships, phylogenetic analysis was performed on the coding regions of the full-length S-, M- and L-genomic segment sequences of NVAV and other representative rodent-, shrew-, mole- and bat-borne hantaviruses, using maximum likelihood and Bayesian methods, implemented in PAUP* (Phylogenetic Analysis Using Parsimony, 4.0b10)36, RAxML Blackbox webserver37 and MrBayes 3.138, under the best-fit GTR+I+Γ model of evolution selected by hierarchical likelihood-ratio test in MrModeltest v2.339 and jModelTest version 0.140. Two replicate Bayesian Metropolis–Hastings Markov Chain Monte Carlo runs, each comprising six chains of 10 million generations sampled every 100 generations with a burn-in of 25,000 (25%), resulted in 150,000 trees overall. Sequences were aligned using Clustal W41, and each genomic segment was treated separately in phylogenetic analyses. The posterior node probabilities were based on 2 million generations and estimated sample sizes over 100 (implemented in MrBayes).

Experimental NVAV infection

To ascertain if experimental NVAV infection resembled that observed with hantaviruses isolated from murid rodents and crocidurine shrews, two-day old Swiss Webster mice were inoculated with 200 pfu of NVAV strain Te34, by the intraperitoneal route. Mice were monitored for weight loss and neurological signs, and the severity of their clinical condition was scored, using a standardized scale42,43. Terminal sera were tested for antibodies against NVAV by the IFA technique6, using NVAV-infected Vero E6 cells spotted onto 10-well slides (Tekdon, Inc., Myakka City, FL) and Alexa Fluor® 488 conjugated goat anti-mouse IgG (H+L) as the secondary antibody (ThermoFisher, Cat. #A-11001, Waltham, MA). Lung, liver, kidney, spleen and brain tissues were collected from moribund mice for histopathology and immunohistochemistry, as well as for quantitation of NVAV RNA by real-time PCR. In addition, brain tissue homogenates from NVAV-infected moribund mice were inoculated intraperitoneally into new litters of infant mice, and NVAV sequences in tissues were compared to that of the NVAV cell-culture isolate and original lung tissue used in the isolation attempts.

Quantitative Taqman real-time PCR

The primers and probe targeting the NVAV S segment were designed using the Primer Express® software version 3.0 (Applied Biosystems, Foster City, CA). The S-segment forward and reverse primer sequences for real-time PCR were 5′–GCGGTGTTAAGGTCCCAAAA–3′ and 5′–CTTCTGCCTTGATAGCTGATTGAG–3′, respectively, and the probe was 5′–CTTTATGTCTCACTCCCAAC–3′. The probe was labeled with the reporter dye FAM at the 5′–end and quencher dye MGB/non fluorescent at the 3′–end, respectively. Each 20-μL reaction contained 1 μL cDNA, 10 μL 2X TaqMan Gene Expression master mix (Applied Biosystems), 0.5 μL forward and reverse primers (36 μM), 0.5 μL fluorescent probe (10 μM), and 7.5 μL double deionized water. The reaction was performed at 50 °C for 2 min and 90 °C for 10 min, followed by 40 cycles at 95 °C for 15 sec and 60 °C for 1 min, in an ABI 7500 Real-Time PCR system (Applied Biosystems). The standard curve was constructed using RNA from NVAV-infected Vero E6 cells. RNA concentration was measured using a NanoDrop spectrophotometer. After RT-PCR using random primers (10 μM), 10-fold serial dilutions of the cDNA were used in duplicate to generate a standard curve.

Discussion

Prior to this report, only two non-rodent-borne hantaviruses had been isolated. One is the prototype shrew-borne hantavirus, known as Thottapalayam virus (TPMV), isolated initially in infant mice injected with spleen tissue homogenate from an Asian house shrew (Suncus murinus), captured during 1964 in southern India44. The other is Imjin virus (MJNV), isolated from the Ussuri white-toothed shrew (Crocidura lasiura), captured near the demilitarized zone in the Republic of Korea6. This report represents the isolation and partial characterization of the first mole-borne hantavirus, known as NVAV, which is widespread across the geographic distribution of its talpid host, the European mole33,34.

Because most of the newfound hantaviruses detected in shrews, moles and insectivorous bats have relied on archival tissues, poor tissue preservation, suboptimal RNA quality and limited sample quantity have all conspired to thwart intensive efforts at virus isolation and at acquiring full-length genome sequences. In addition, the high sequence diversity of non-rodent-borne hantaviruses has hampered reliable primer design. Thus, only a few of the hantaviruses harbored by shrews, moles and bats have been fully sequenced, preventing definitive phylogenetic analysis.

The availability of NVAV in cell culture has facilitated the acquisition of its whole genome, resulting in robust phylogenetic trees based on all three segments, with each topology demonstrating that NVAV occupies a highly divergent lineage. NVAV was closely related to bat-borne hantaviruses, suggesting a complex evolutionary history, presumably involving host-switching events, as shown previously for Asama virus in the Japanese shrew mole (Urotrichus talpoides)16, Oxbow virus in the American shrew mole (Neurotrichus gibbsii)17 and Rockport virus in the eastern mole (Scalopus aquaticus)19, as well as periodic episodes of host/pathogen co-divergence through deep evolutionary time19,25,26.

The Gn and Gc glycoproteins, encoded by the M segment, are antigenic and serve as determinants of the protective humoral immune response45,46. One of the traditional criteria to assign species status to hantaviruses is the absence of antibody cross neutralization, as determined by plaque-reduction neutralization tests (PRNT)47. Previously, no antibody cross neutralization could be demonstrated between TPMV or MJNV and rodent-borne hantaviruses6,48. In the present report, PRNT were not performed because the less than 50% amino acid sequence similarity between the glycoproteins of NVAV and representative hantaviruses for which isolates are available predicted no antibody cross neutralization.

The N-linked glycosylation sites on the Gn (N134, N235, N347, N399 and N609) and Gc (N928) of prototype Hantaan virus (HTNV) strain 76-118 are relatively well conserved among murid rodent-borne hantaviruses and are critical for the conformation of these proteins for transport, receptor binding and antigenicity49,50,51,52. Overall, the first, third and fourth N-linked glycosylation sites on the Gn are also generally preserved among cricetid rodent-borne hantaviruses, as well as the few shrew- and mole-borne hantaviruses for which full-length M-segment sequences are available6,10,12,16,17,19,52. By contrast, in both NVAV strain Te34, reported here, and NVAV strain BE/Namur/TE/2013/1 from Belgium (GenBank KT004446)53, the glycosylation site of N101 on the Gn is unique among all hantaviruses reported to date, including Laibin virus24, the only bat-borne hantavirus which has a whole M-segment sequence. To what extent this unique N-linked glycosylation site impacts the function of the Gn is unknown.

As mentioned above, the L-protein motifs of NVAV closely resembled that of other rodent-borne hantaviruses54, as well as the six shrew- and one mole-borne hantaviruses for which full-length L segment sequences are available6,7,10,12,14,19. Recently, antibodies directed against the middle portion of the Tula virus L protein were used to study the membrane association of the L protein in infected cells55. A similar approach would allow the detection and monitoring of L protein expression in NVAV-infected cells and tissues. These and other studies are now possible with the isolation of NVAV.

The lack of a small animal model has hampered studies aimed at elucidating the pathogenesis of HFRS. Infant mice, as well as certain strains of adult mice (such as C57BL/6, BALB/c, AKR/J, and SJL/J), develop fatal meningoencephalitis following inoculation with HTNV and other HFRS-causing murid rodent-borne hantaviruses56,57,58,59,60,61,62, as well as TPMV and MJNV from crocidurine shrews63 (Gu et al. unpublished observations). Our data indicate that experimental NVAV infection also resulted in neuroinvasive disease. By contrast, infant mice are not susceptible to infection with cricetid rodent-borne hantaviruses, such as Sin Nombre virus (SNV). The absence of convincing human correlates of acute hantaviral neurological disease and of hantavirus antigen and RNA in brains of HFRS and HPS patients argue against the utility and significance of these infant mouse models. That is, the infectivity of hantaviruses in infant and adult mice does not predict the pathogenic potential of hantaviruses in humans.

On the other hand, Syrian hamsters (Mesocricetus auratus) may serve as a useful experimental host. Because of their susceptibility to experimental hantavirus infection, Syrian hamsters have been used extensively in the testing of candidate vaccines for HFRS64,65. Moreover, a lethal HPS model has been developed in adult Syrian hamsters, which recapitulates the clinical and pathological features of HPS in humans66,67. Recently, we demonstrated that infant and juvenile Syrian hamsters develop a fatal disease, characterized by respiratory distress, weight loss, hind limb paralysis and seizures, following intraperitoneal inoculation with MJNV68. High viral loads were found in tissues, and histopathology included hepatic necrosis, widespread inflammation in lungs without pulmonary edema, and extensive MJNV antigens in microvascular endothelial cells. Adult hamsters, on the other hand, were resistant to MJNV infection. Whether or not Syrian hamsters, which are transiently immunosuppressed with dexamethasone and cyclophosphamide, will be rendered susceptible to MJNV infection, as has been shown for SNV, the HPS-causing prototype hantavirus69, awaits future study.

Similarly, if NVAV is pathogenic in immunocompetent or immunosuppressed Syrian hamsters, new insights will be gained about hantavirus pathogenesis, and the development of medical countermeasures for hantavirus disease prevention will be accelerated. Alternatively, the development of a chronic NVAV infection model in Syrian hamsters with persistent virus shedding in saliva, urine and/or feces, that presumably resemble natural NVAV infection in European moles, would serve as a convenient model system in which to study NVAV transmission dynamics and genomic reassortment.

Finally, the long-awaited isolation of a highly divergent mole-borne hantavirus will accelerate the acquisition of new knowledge about its infectivity and pathogenic potential in humans. That is, because NVAV infection is highly prevalent in European moles and because European moles often reside near human habitation, clinicians and healthcare personnel should be vigilant for possible NVAV infection in individuals, who develop febrile illnesses or unusual clinical syndromes following known or suspected exposures. The availability of the NVAV isolate will facilitate the development of robust diagnostic tests in search for human infection and disease.

Acknowledgements

This research was supported in part by U.S. Public Health Service grant R01AI075057 from the National Institute of Allergy and Infectious Diseases and P20GM103516 and P30GM114737 from the National Institute of General Medical Sciences, National Institutes of Health; N N303 604538 from the State Committee for Scientific Research in Poland; and by the Healthy Ageing Research Centre project (REGPOT-2012-2013-1, 7FP). The expert assistance provided by Tina Carvalho, of the Biological Electron Microscopy Facility, is gratefully acknowledged, as are the services of the Genomics Core Facility, funded partially by the Centers of Biomedical Research Excellence program (P30GM103341). Grant DBI-1040548 from the National Science Foundation provided funding to acquire the Hitachi HT7700 TEM system and RMC Powertome ultramicrotome. The funding agencies had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Affiliations

Departments of Pediatrics and Tropical Medicine, Medical Microbiology and Pharmacology, John A. Burns School of Medicine, University of Hawaii at Manoa, Honolulu, Hawaii, USA

Se Hun Gu

, Mukesh Kumar

& Richard Yanagihara

Department of Molecular Pathology and Neuropathology, Faculty of Medicine, Medical University of Łódź, 92-216 Łódź, Poland

Beata Sikorska

& Paweł P. Liberski

Department of Biodiversity Studies, Didactics and Bioeducation, Faculty of Biology and Environmental Protection, University of Łódź, 90-237 Łódź, Poland

Janusz Hejduk

& Janusz Markowski

Department of Experimental Zoology and Evolutionary Biology, Faculty of Biology and Environmental Protection, University of Łódź, 90-237 Łódź, Poland

Search for Paweł P. Liberski in:

Search for Richard Yanagihara in:

Contributions

S.H.G., M.K. and R.Y. conceived and designed the experiments; J.H., J.M. and M.M. conducted the trapping and field collections; S.H.G., M.K., B.S. and P.P.L. performed the experiments; S.H.G., M.K. and R.Y. contributed mouse model expertise; S.H.G., M.K., B.S., P.P.L. and R.Y. analyzed the data and prepared the figures; and all authors prepared and reviewed the manuscript.

Competing interests

Corresponding author

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/