Project description:BackgroundHuman granulocytic anaplasmosis is a tick-borne zoonotic disease caused by Anaplasma phagocytophilum. Coinfections with A. phagocytophilum and other tick-borne pathogens are reported frequently, whereas the relationship between A. phagocytophilum and flea-borne Yersnia pestis is rarely concerned.ResultsA. phagocytophilum and Yersnia pestis were discovered within a Marmota himalayana found dead in the environment, as determined by 16S ribosomal rRNA sequencing. Comparative genomic analyses of marmot-derived A. phagocytophilum isolate demonstrated its similarities and a geographic isolation from other global strains. The 16S rRNA gene and GroEL amino acid sequence identity rates between marmot-derived A. phagocytophilum (JAHLEX000000000) and reference strain HZ (CP000235.1) are 99.73% (1490/1494) and 99.82% (549/550), respectively. 16S rRNA and groESL gene screenings show that A. phagocytophilum is widely distributed in marmots; the bacterium was more common in marmots found dead (24.59%, 15/61) than in captured marmots (19.21%, 29/151). We found a higher Y. pestis isolation rate in dead marmots harboring A. phagocytophilum than in those without it (2 = 4.047, p < 0.05). Marmot-derived A. phagocytophilum was able to live in L929 cells and BALB/c mice but did not propagate well.ConclusionsIn this study, A. phagocytophilum was identified for the first time in Marmota himalayana, a predominant Yersinia pestis host. Our results provide initial evidence for M. himalayana being a reservoir for A. phagocytophilum; moreover, we found with the presence of A. phagocytophilum, marmots may be more vulnerable to plague. Humans are at risk for co-infection with both pathogens by exposure to such marmots.
Project description:The Eastern woodchuck (Marmota monax) is a classical animal model for studying hepatitis B virus (HBV) infection and hepatocellular carcinoma (HCC) in humans. Recently, we found that Marmota himalayana, an Asian animal species closely related to Marmota monax, is susceptible to woodchuck hepatitis virus (WHV) infection and can be used as a new mammalian model for HBV infection. However, the lack of genomic sequence information of both Marmota models strongly limited their application breadth and depth. To address this major obstacle of the Marmota models, we utilized Illumina RNA-Seq technology to sequence the cDNA libraries of liver and spleen samples of two Marmota monax and four Marmota himalayana. In total, over 13 billion nucleotide bases were sequenced and approximately 1.5 billion clean reads were obtained. Following assembly, 106,496 consensus sequences of Marmota monax and 78,483 consensus sequences of Marmota himalayana were detected. For functional annotation, in total 73,603 Unigenes of Marmota monax and 78,483 Unigenes of Marmota himalayana were identified using different databases (NR, NT, Swiss-Prot, KEGG, COG, GO). The Unigenes were aligned by blastx to protein databases to decide the coding DNA sequences (CDS) and in total 41,247 CDS of Marmota monax and 34,033 CDS of Marmota himalayana were predicted. The single nucleotide polymorphisms (SNPs) and the simple sequence repeats (SSRs) were also analyzed for all Unigenes obtained. Moreover, a large-scale transcriptome comparison was performed and revealed a high similarity in transcriptome sequences between the two marmota species. Our study provides an extensive amount of novel sequence information for Marmota monax and Marmota himalayana. This information may serve as a valuable genomics resource for further molecular, developmental and comparative evolutionary studies, as well as for the identification and characterization of functional genes that are involved in WHV infection and HCC development in the woodchuck model.