- Open Access
The isolation and molecular characterization of Mycobacterium avium subsp. paratuberculosis in Shandong province, China
Gut Pathogensvolume 8, Article number: 9 (2016)
Mycobacterium avium subspecies paratuberculosis (Map) causes Johne’s disease in domestic and wild ruminants. It has been a debate that whether Map can cause Crohn’s disease in human. To our knowledge there is no report about molecular characterization of Map in China, although several Map strains have been reported in other country. The objectives of this study was to know the recent prevalence of Johne’s disease in dairy farms in Shandong province, and have a better understanding of genotypic distribution of Map in China.
Johne’s disease was detected from 1038 individuals in 19 dairy farms by ELISA. Map in fecal and milk specimens was identified by Ziehl-Neelsen staining and confirmed using PCR-REA. In addition, frozen sections of ileum and mesenteric lymph nodes from two Map shedding cows were performed to observe the histopathological changes. Next-generation sequencing technology was performed to get whole genome sequences.
A total of 121 (11.7 %) animals were positive for Map antibody from 1038 sera tested, and 11 (57.9 %) dairy herds were positive for Map antibody. Typically histopathologic changes were observed in mesenteric lymph nodes. We have successfully isolated two Map strains, which both were Map-C. The current genome-wide analysis showed that the genome size of our isolates are respectively 4,750,273 and 4,727,050 bp with a same G + C content of 69.3 %, and the numbers of single nucleotide polymorphisms (SNPs) against Map K-10 are respectively 292 and 296.
Map is a prevalent pathogen among dairy cattle in China. This study successfully isolated two Map strains from one Chinese dairy herd with signs of diarrhoea, and identified that the two isolates were both Map-C. Furthermore, these isolates were most closely related to Map K-10.
Paratuberculosis, also known as Johne’s disease, caused by Mycobacterium avium subspecies paratuberculosis (Map), is an important and highly prevalent disease of domestic and wild ruminants manifest as a chronic granulomatous enteritis with decreased milk production and in more serious cases leading to progressive emaciation and death. The disease is transmitted via milk and colostrum to calves and by the fecal-oral route to animals of all ages. Intra-uterine transmission can also occur and Map can also be detected in the saliva of cows, indicating this as a potential further mode of transmission [1, 2]. The disease is endemic in many parts of the world with serious economic losses which in the USA, for example, are estimated at between $200 and $250 million annually .
Map belongs to the M. avium complex  and consists of three important types designated “Sheep-type” (type S) and “Cattle-type” (type C) which were originally isolated from sheep or cattle respectively  and “Bison-type” (B-type)  all of which can be differentiated by IS1311 restriction endonuclease genotyping. Despite designation of strains as B, S or C there is no host-species specificity  and strains isolated from sheep or cattle can be either C or S type . Most of the ruminant population in India are infected by Map ‘B type’ strains [9, 10].
The importance of Map as a zoonotic pathogen has not yet been clarified [11, 12]. The clinical presentation and pathology of Crohn’s disease (CD) bears more than a superficial resemblance to paratuberculosis. Despite some studies indicating no absolute relationship between Map and CD [13, 14] a proportion of individuals with CD contain Map or Map DNA in their blood or infected intestinal tissue . However, Map can be isolated from stool samples of individuals who are healthy or presenting with unrelated diseases . This latter observation indicates the general susceptibility of humans to Map infection. Map can be readily isolated from milk , including pasteurized milk  and may also be isolated from the muscles of infected cattle, consumption of which, may be an additional potential route of transmission to humans .
Complete eradication of the disease from cattle, although desirable, is extremely difficult . Vaccination can be effective to control clinical disease, reducing faecal shedding and increasing productivity but it does not completely eliminate infection [21, 22]. Various quarantine measures can be adopted at dairy farms for controlling this disease including the isolation of young calves from infected cattle. However, young stock may also excrete Map in their feces emphasizing current difficulties in disease control.
For molecular identification of Map by PCR a number of different targets have been used including f57 , HspX , genes 251 and 255 , ISMav2 , ISMpa1  and ISMAP02 , some of which (f57, HspX and ISMAP02) are unique to Map but which are present in low (1–6) copy number with associated limited potential for differentiation . For this reason IS900, with 14–18 copies is a preferable PCR target .
Additional PCR targets may be identified from Map genome sequences. Currently, the only fully annotated genome is that of Map strain K-10  although additional draft sequences from strains ATCC 19698, Map s397 and Map S5 from different countries have been published on NCBI. Despite the likely high prevalence of paratuberculosis in China, no genome sequences are available for any Chinese strains. Although Map has been isolated from Chinese farms in the past , no information is available on the strain types and molecular characterisation present. The present study thus had the following objectives: (1) investigating the prevalence rate of Map infection in Chinese farms with clinical signs of diarrhoea, (2) typing of Map strains by PCR-restriction endonuclease analysis based on polymorphisms in IS1311 and, (3) to sequence the whole genome to contribute to a better understanding of the evolution of Map.
1038 individual bovine serum samples were collected between October 2014 and August 2015 from 19 different dairy farms located in nine distinct geographic regions of Shandong province. The herds were selected at random and varied between 300–800 animals. Individual cows showing diarrhoea were selected in addition to some animals showing no clinical signs. Approximately 10 % of the animals were selected all of them being adult cows over 24 months of age. Milk, faeces and tissues were sampled from the herd showing the highest frequency of positive sera. All the faecal samples were obtained via per-rectal drags, and were kept at 4 °C up to 48 h prior to processing for culture and DNA extraction. Milk samples were collected in 50-ml tubes after four initial strippings were discarded. Ileum and mesenteric lymph nodes were collected from two euthanized cows with positive antibody for Map. All protocols and procedures were performed according to the Chinese Regulations of Laboratory Animals—The Guidelines for the Care of Laboratory Animals (Ministry of Science and Technology of People’s Republic of China) and Laboratory Animal Requirements of Environment and Housing Facilities (GB 14925–2010, National Laboratory Animal Standardization Technical Committee). The license number associated with their research protocol was 20110611–01 and the animal study proposal was approved by The Laboratory Animal Ethical Committee of China Agricultural University.
Detection of Map antibody
In order to identify the herd with highest frequency of positive sera the 1038 sera were tested for the presence of antibodies to Map using a commercial antibody test kit (IDEXX). Sera with S/P ratios ≤0.45 and ≥0.55 were considered negative and positive, respectively. Intermediate S/P values were considered “suspect” and regarded as negative for the data analysis.
Detection of bacteria
10 mg feces without any treatment was taken to spread evenly on to the surface of a clear glass slide with 20 ml water drop. For milk samples, there was a little difference. Centrifugation was needed to get more bacterium. After centrifugation at a speed of 12000 rpm for 5 min, the bottom precipitate was collected as the samples for smear, and the same procedure was used as that of faecal samples. Tissues were used as frozen sections. Slides were stained using the Ziehl–Neelsen method and examined by light microscopy for the presence of acid-fast bacilli. Map presence was confirmed by PCR amplification of IS900 as described previously .
Sample preparation and culture
A modification of the double centrifugation methods was used to cultivate Map . Briefly, two grams of faeces were added to 35 ml of sterile distilled water. Samples tubes were vortexed for 1 min and then allowed to stand undisturbed for 45 min. Five ml of the supernatant were added to 25 ml of 0.9 % hexadecylpyridinium chloride (HPC) and then vortexed for 2 min and allowed to stand undisturbed for 12 h at room temperature. Tubes were centrifuged at 1500g for 30 min, the supernatant was discarded and the pellet was resuspended in 1 ml of antibiotic mixture (50 µg/ml amphotericin B, and 100 µg/ml nalidixic acid, and 100 µg/ml vancomycin). This was mixed well and incubated for 12 h at 37 °C. Herrold’s Egg Yolk Agar with mycobactin J and ANV (BD, USA) were inoculated with 0.2 ml of the suspension and incubated at 37 °C for 6 weeks.
Isolation and identification of Map and strain typing of isolates
Colonies that appeared within 1 week were discarded and slower-growing colonies with typical colonial morphology were checked for acid fast bacteria, and then detected by PCR based on IS900 as mentioned above. In order to generate pure and sufficient biomass, a single colony of each confirmed Map isolate was subcultured into 7H9 broth supplemented with OADC and Mycobactin J and incubated at 37 °C. The culture was reconfirmed as Map by PCR and the purity checked by light microscopy with Ziehl–Neelsen method staining. The strain type was identified by IS1311 PCR and restriction endonuclease analysis (REA) using the primers and method described previously . The identification as M. avium subsp. paratuberculosis was confirmed by PCR using IS900 and REA using AlwI digestion. The IS1311 PCR-REA was used to distinguish C strains from S strains using HinfI and MseI (details shown in Table 1). PCR primers are detailed in Table 2.
Next-generation whole genome sequencing and analysis
Massive parallel sequencing using the Solexa-technology (Illumina) was performed at the Chinese Academy of Sciences. Two isolates were sequenced using the Illumina HiSeq 2500 platform using one flow cell lane with 125-cycle paired end chemistry. Short reads were assembled using SOAPdenovo (http://www.soap.genomics.org.cn), a genome assembler developed specifically for next-generation short-read sequences . The SOAP GapCloser was also used to close gaps. Reads were mapped by using BWA , and SNPs and indels were called and filtered by SAM tools . A rooted tree was computed using MEGA. Phylogenetic tree was calculated using the maximum likelihood method by whole genome alignment and performing 1000 bootstraps replicates. Several references are selected for comparation, include type C, type B, and type S (Table 3).
Statistical comparisons between milk detection and fecal detection were determined using the Chi square test with p < 0.05 denoting statistical significance. SPSS software (Chicago, IL, USA) was used.
The survey of the sera
Screening of 1038 serum samples from 19 dairy farms by the Mycobacterium paratuberculosis antibody test kit detected 121 (11.7 %) animals positive for Map antibody, and 11 (57.9 %) dairy farms were positive for Map antibody. Of 11 dairy farms with signs of diarrhoea, 7 (63.6 %) were positive and of 8 dairy farms with no signs of diarrhea, 3 (37.5 %) were positive. The highest positive rate within these herds was 78.9 % and this herd was selected for isolation and identification of Map (Table 4).
The morphology of bacteria and pathological changes in the tissues
Nineteen faecal and milk samples were sampled for Ziehl–Neelsen staining to detect individual cattle positive for acid fast bacilli. Of these 13 (68.4 %), 6 (31.6 %) faecal samples and milk samples respectively were positive for acid-fast bacilli (Fig. 1a, b). We euthanized two cattle with diarrhoea, emaciation and reduced milk production. Postmortem examination showed gross lesions in the mesenteric lymph nodes showing 5–6 fold enlargement compared with normal mesenteric lymph nodes (date not shown). However, there were no obvious gross lesions in the gastro-intestinal tract. Histopathological examination of mesenteric lymph nodes showed the characteristic diffuse granulomas with epithelioid cells, multinucleated giant cells and macrophages (Fig. 2a, b) compared with normal mesenteric lymph nodes(Fig. 2c, d). However, no granulomas were observed in the ileum (date not shown). Ziehl–Neelsen stained impression smears of the mesenteric lymph nodes showed large numbers of acid-fast bacilli (Fig. 1c). Large numbers of clustered acid-fast bacilli were observed in the smears taken from a purified colony (Fig. 1d).
The identification and strain typing of MAP
The PCR products are shown Fig. 3. Two isolates were both identified as positive using the IS900 PCR. We obtained the target fragment, whose size was 413 bp (Fig. 3a) which was digested using restriction endonuclease AlwI which resulted in two bands, one at 147 bp, and other at 266 bp (Fig. 3b). For further identification of strain type, IS1311 PCR was conducted and followed by digestion by two restriction endonucleases (HinfI/MseI). The results showed that both isolates were identified as positive with the target band of 608 bp (Fig. 3c). Restriction endonuclease analysis of the IS1311 PCR reaction products confirmed the presence of bands of 67, 218, 285 and 323 bp (Fig. 3d), which indicated that both isolates were type C.
Whole genome sequencing and phylogenetic analysis
Although we grouped the two isolates, it was not enough to deeper understand the genomic differences between our isolates and reference strains. Here Next-generation sequencing technology (Illumina, HiSeq 2500) was used to sequence the whole genome of the two isolates from cattles in China. As a result, we successfully got two draft sequences named 2015WD-1 and 2015WD-2, the accession numbers are LKUS00000000 and LKUT00000000, respectively. The results showed that the genome size of 2015WD-1 and 2015WD-2 are respectively 475,0273 and 472,7050 bp with a same G + C content of 69.3 %. The numbers of single nucleotide polymorphisms (SNPs) in our Map isolates against Map K-10 are respectively 292 and 296. Also, 84 and 96 indels are found in our isolates against Map K-10. More information between our isolates and reference strains are presented in Table 5. In the phylogenetic tree (Fig. 4) based on the nucleotide sequences, we can see our two isolates are most closely relate to type C representative strains Map K-10 and Map ATCC19698, and farthest from type S representive strains Map JQ5 and Map s397.
The first study on paratuberculosis in China was recorded in Inner Mongolia in 1953, and then it was reported consistently in China . In Shandong province, there has been one previous study of the prevalence of paratuberculosis. This indicated a prevalence of 29.34 % in 2011 and 14.93 % in 2012, detected by the interferon gamma release assay . Detection by ELISA has indicated prevalences of between 2 and 4 % in Tibet, Shanghai and Guangxi provinces and highly variable values of 0–73.4 % in Inner Mongolia [40–42].
By contrast levels of infection in Europe and other western countries appear to be equally high with individual positive rates in cattle in estimated at 20 % and the herd-level prevalence at >50 % according to the summary of Nielsen and Toft . Other studies showed seroprevalence rates of between 1.99 % in Asturias, Spain , 16.1 % in the east of Canada  and 29.8 % in Northern India . In Egypt in 26 Holstein dairy herds in 2010 and 2013 herd and individual level prevalence were recorded as 65.4 % and 13.6 % respectively. . No country claims that it is free on MAP.
Our studies showed a seroprevalence of 11.7 % in individual level in Shandong province of China, which is a little higher than most of provinces of China and a little lower than Europe, Indian and Canada. And the herd level prevalence was 57.9 % in 19 dairy herds we have investigated, indicated that Map has become a common pathogen in dairy farms in China. Among the 19 herds, there were more or less individuals behaved diarrhea, but the Map antibody positive rate were zero. Maybe these symptoms of diarrhea were caused by any other organisms than Map.
Two strains of Map were isolated from a farm showing heavy losses in their dairy cattle from diarrhoea resulting in weakness and reductions of milk production.
We successfully isolated two MAP from one serious loss of dairy farm in which many cattle behaved diarrhoea, weakness, reduction of milk production. Environmental bacteria such as streptococci were observed to be present in the samples (data not shown) but decontamination using an adaptation off two HPC methods (; Report, 2011) was clearly successful.
Several methods were used to detect the pathogen; fast-acid staining showed that the isolation rate of milk (9/19) was less than in faeces (14/19), but the difference is not significant (p = 0.184). The histological picture and bacillary morphology were consistent with paratuberculosis.
IS1311 PCR-REA provides similar information as IS900 RFLP analysis but is more useful than RFLP analysis where DNA is degraded or present in low concentrations . To further confirm identification we used the IS900 PCR-REA to confirm the identity of isolates as Map and IS1311 PCR-REA for strain typing use of which indicated that the isolates we obtained were Map and were C strains.
Most previous studies have indicated that C strains were the dominant strains in cattle in countries as different as Australia, New Zealand, South America, North America and Europe [5, 48–50]. The exception appears to be that B strains are the most prevalent (82 %) genotype of Map in all domestic ruminants in India with C strains present as a minority strain . It was thus important to know which types existed in China. Although the present study involved two strains only these were both C strains, which is at least consistent with most other previous studies.
In terms of whole genome sequencing, our two isolates have less base pairs than Map K-10, but the same G + C content with Map K-10. Our two isolates have least number of SNPs (N = 292 for 2015WD-1 and 296 for 2015WD-2) and indels (N = 84 2015WD-1 and 96 for 2015WD-2) against Map K-10 when compared with Map s397 and Map s5, indicating that Chinese isolates are most likely originated from Map K-10 with few changes. From phylogenetic tree, we can also find that 2015WD-1and 2015WD-2 are most closely related to Map K-10, which further confirmed the identified strain type of our isolates.
Although there is no definite conclusion on the relationship between Map and CD, Map is detected sometimes in normal individuals, which less often than in Crohn’s patients . The fact that Map can be found in normal individuals does not exclude it as the cause of the disease, because the normal individuals maybe subclinically infected. Perhaps CD is caused by complex factors. Map is one of them and other element such as environmental stress, immunosuppressive conditions and so on are involved.
Map has become a common pathogen in dairy farms in China. Present study successfully isolated two Map strains from one Chinese dairy farm with signs of diarrhoea, and identified that the two isolates were both Map-C. Farther more, we found that our isolates were most closely relate to Map K-10.
Stevenson K. Genetic diversity of Mycobacterium avium subspecies paratuberculosis and the influence of strain type on infection and pathogenesis: a review. Vet Res. 2015;46:64.
Sorge US, Kurnick S, Sreevatsan S. Detection of Mycobacterium avium subspecies paratuberculosis in the saliva of dairy cows: a pilot study. Vet Microbiol. 2013;164:383–6.
Ott SL, Wells SJ, Wagner BA. Herd-level economic losses associated with Johne’s disease on US dairy operations. Prev Vet Med. 1999;40:179–92.
Thorel MF, Krichevsky M, Levy-Frebault VV. Numerical taxonomy of mycobactin-dependent mycobacteria, emended description of Mycobacterium avium, and description of Mycobacterium avium subsp. avium subsp. nov., Mycobacterium avium subsp. paratuberculosis subsp. nov., and Mycobacterium avium subsp. silvaticum subsp. nov. Int J Syst Bacteriol. 1990;40:254–60.
Collins DM, Gabric DM, de Lisle GW. Identification of two groups of Mycobacterium paratuberculosis strains by restriction endonuclease analysis and DNA hybridization. J Clin Microbiol. 1990;28:1591–6.
Sohal JS, Singh SV, Singh PK, Singh AV. On the evolution of ‘Indian Bison type’ strains of Mycobacterium avium subspecies paratuberculosis. Microbiol Res. 2010;165:163–71.
Ronai Z, Csivincsik A, Gyuranecz M, Kreizinger Z, Dan A, et al. Molecular analysis and MIRU-VNTR typing of Mycobacterium avium subsp. paratuberculosis strains from various sources. J Appl Microbiol. 2015;118:275–83.
Stevenson K, Hughes VM, de Juan L, Inglis NF, Wright F, et al. Molecular characterization of pigmented and nonpigmented isolates of Mycobacterium avium subsp. paratuberculosis. J Clin Microbiol. 2002;40:1798–804.
Sevilla I, Singh SV, Garrido JM, Aduriz G, Rodriguez S, et al. Molecular typing of Mycobacterium avium subspecies paratuberculosis strains from different hosts and regions. Rev Sci Tech. 2005;24:1061–6.
Yadav D, Singh SV, Singh AV, Sevilla I, Juste RA, et al. Pathogenic ‘Bison-type’ Mycobacterium avium subspecies paratuberculosis genotype characterized from riverine buffalo (Bubalus bubalis) in North India. Comp Immunol Microbiol Infect Dis. 2008;31:373–87.
Waddell LA, Rajic A, Stark KD, Mc ES. The zoonotic potential of Mycobacterium avium ssp. paratuberculosis: a systematic review and meta-analyses of the evidence. Epidemiol Infect. 2015;143:3135–57.
Rowbotham DS, Mapstone NP, Trejdosiewicz LK, Howdle PD, Quirke P. Mycobacterium paratuberculosis DNA not detected in Crohn’s disease tissue by fluorescent polymerase chain reaction. Gut. 1995;37:660–7.
Barrett JC, Hansoul S, Nicolae DL, Cho JH, Duerr RH, et al. Genome-wide association defines more than 30 distinct susceptibility loci for Crohn’s disease. Nat Genet. 2008;40:955–62.
Over K, Crandall PG, O’Bryan CA, Ricke SC. Current perspectives on Mycobacterium avium subsp. paratuberculosis, Johne’s disease, and Crohn’s disease: a review. Crit Rev Microbiol. 2011;37:141–56.
Kirkwood CD, Wagner J, Boniface K, Vaughan J, Michalski WP, et al. Mycobacterium avium subspecies paratuberculosis in children with early-onset Crohn’s disease. Inflamm Bowel Dis. 2009;15:1643–55.
Tuci A, Tonon F, Castellani L, Sartini A, Roda G, et al. Fecal detection of Mycobacterium avium paratuberculosis using the IS900 DNA sequence in Crohn’s disease and ulcerative colitis patients and healthy subjects. Dig Dis Sci. 2011;56:2957–62.
Stabel JR, Bradner L, Robbe-Austerman S, Beitz DC. Clinical disease and stage of lactation influence shedding of Mycobacterium avium subspecies paratuberculosis into milk and colostrum of naturally infected dairy cows. J Dairy Sci. 2014;97:6296–304.
Ellingson JL, Anderson JL, Koziczkowski JJ, Radcliff RP, Sloan SJ, et al. Detection of viable Mycobacterium avium subsp. paratuberculosis in retail pasteurized whole milk by two culture methods and PCR. J Food Prot. 2005;68:966–72.
Alonso-Hearn M, Molina E, Geijo M, Vazquez P, Sevilla I, et al. Isolation of Mycobacterium avium subsp. paratuberculosis from muscle tissue of naturally infected cattle. Foodborne Pathog Dis. 2009;6:513–8.
Wolf R, Orsel K, De Buck J, Barkema HW. Calves shedding Mycobacterium avium subspecies paratuberculosis are common on infected dairy farms. Vet Res. 2015;46:71.
Singh SV, Singh PK, Kumar N, Gupta S, Chaubey KK, et al. Evaluation of goat based ‘indigenous vaccine’ against bovine Johne’s disease in endemically infected native cattle herds. Indian J Exp Biol. 2015;53:16–24.
Tewari D, Hovingh E, Linscott R, Martel E, Lawrence J, et al. Mycobacterium avium subsp. paratuberculosis antibody response, fecal shedding, and antibody cross-reactivity to Mycobacterium bovis in M. avium subsp. paratuberculosis-infected cattle herds vaccinated against Johne’s disease. Clin Vaccine Immunol. 2014;21:698–703.
Vansnick E, De Rijk P, Vercammen F, Geysen D, Rigouts L, et al. Newly developed primers for the detection of Mycobacterium avium subspecies paratuberculosis. Vet Microbiol. 2004;100:197–204.
Miller JM, Jenny AL, Ellingson JL. Polymerase chain reaction identification of Mycobacterium avium in formalin-fixed, paraffin-embedded animal tissues. J Vet Diagn Invest. 1999;11:436–40.
Bannantine JP, Baechler E, Zhang Q, Li L, Kapur V. Genome scale comparison of Mycobacterium avium subsp. paratuberculosis with Mycobacterium avium subsp. avium reveals potential diagnostic sequences. J Clin Microbiol. 2002;40:1303–10.
Shin SJ, Chang YF, Huang C, Zhu J, Huang L, et al. Development of a polymerase chain reaction test to confirm Mycobacterium avium subsp. paratuberculosis in culture. J Vet Diagn Invest. 2004;16:116–20.
Olsen I, Johansen TB, Billman-Jacobe H, Nilsen SF, Djonne B. A novel IS element, ISMpa1, in Mycobacterium avium subsp. paratuberculosis. Vet Microbiol. 2004;98:297–306.
Stabel JR, Bannantine JP. Development of a nested PCR method targeting a unique multicopy element, ISMap02, for detection of Mycobacterium avium subsp. paratuberculosis in fecal samples. J Clin Microbiol. 2005;43:4744–50.
Castellanos E, de Juan L, Dominguez L, Aranaz A. Progress in molecular typing of Mycobacterium avium subspecies paratuberculosis. Res Vet Sci. 2012;92:169–79.
Bull TJ, Hermon-Taylor J, Pavlik I, El-Zaatari F, Tizard M. Characterization of IS900 loci in Mycobacterium avium subsp. paratuberculosis and development of multiplex PCR typing. Microbiology. 2000;146((Pt 9)):2185–97.
Li L, Bannantine JP, Zhang Q, Amonsin A, May BJ, et al. The complete genome sequence of Mycobacterium avium subspecies paratuberculosis. Proc Natl Acad Sci USA. 2005;102:12344–9.
Zhang J. Detecting pathogeny of bovine paratuberculosis [master degree]: Yunnan Agricultural University. 2008.
Marsh I, Whittington R, Cousins D. PCR-restriction endonuclease analysis for identification and strain typing of Mycobacterium avium subsp. paratuberculosis and Mycobacterium avium subsp. avium based on polymorphisms in IS1311. Mol Cell Probes. 1999;13:115–26.
Douarre PE, Cashman W, Buckley J, Coffey A, O’Mahony JM. Isolation and detection of Mycobacterium avium subsp. paratuberculosis (MAP) from cattle in Ireland using both traditional culture and molecular based methods. Gut Pathog. 2010;2:11.
Luo R, Liu B, Xie Y, Li Z, Huang W, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience.1:18.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, et al. The sequence alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Cui Z. Paratuberculosis. Diseases of dairy cattle. 1 ed. China Agriculture Press. 2007. p. 220–3.
Meng Han LZ, Wang Jilong, Ding Jiabo. Epidemiological investigation of large-scale dairy farm on tuberculosis and paratuberculosis in Shandong province. Chin Agric Sci Bullet. 2014;30:10–3.
Xiaoqiang Wang ZH, Liu Xinyu, Jiang Wenteng. Seroprevalence of paratuberculosis in Yaks in some regions of Tibet, China. China Dairy Cow. 2015;14:63–5.
Hui Zhai GF, Xiao gu, Dongrao Li. Serological investigation and prevention measures of cow paratuberculosis. China Livest Poult Ind. 2015:101–2.
Bao Tian YB, Jiang Yonglu, Sun Qingyu. Serological survey of paratuberculosis in cattle of some dairy farms in inner Mongolia. Anim Husb Feed Sci. 2014;35:21–2.
Nielsen SS, Toft N. A review of prevalences of paratuberculosis in farmed animals in Europe. Prev Vet Med. 2009;88:1–14.
Gasteiner J, Wenzl H, Fuchs K, Jark U, Baumgartner W. Serological cross-sectional study of paratuberculosis in cattle in Austria. Zentralbl Veterinarmed B. 1999;46:457–66.
McKenna SL, Keefe GP, Barkema HW, McClure J, Vanleeuwen JA, et al. Cow-level prevalence of paratuberculosis in culled dairy cows in Atlantic Canada and Maine. J Dairy Sci. 2004;87:3770–7.
Singh SV, Singh AV, Singh R, Sharma S, Shukla N, et al. Sero-prevalence of bovine Johne’s disease in buffaloes and cattle population of North India using indigenous ELISA kit based on native Mycobacterium avium subspecies paratuberculosis ‘Bison type’ genotype of goat origin. Comp Immunol Microbiol Infect Dis. 2008;31:419–33.
Amin AS, Hsu CY, Darwish SF, Ghosh P, AbdEl-Fatah EM, et al. Ecology and genomic features of infection with Mycobacterium avium subspecies paratuberculosis in Egypt. Microbiology. 2015;161:807–18.
Whittington RJ, Hope AF, Marshall DJ, Taragel CA, Marsh I. Molecular epidemiology of Mycobacterium avium subsp. paratuberculosis: IS900 restriction fragment length polymorphism and IS1311 polymorphism analyses of isolates from animals and a human in Australia. J Clin Microbiol. 2000;38:3240–8.
Whipple D, Kapke P, Vary C. Identification of restriction fragment length polymorphisms in DNA from Mycobacterium paratuberculosis. J Clin Microbiol. 1990;28:2561–4.
Cousins DV, Williams SN, Hope A, Eamens GJ. DNA fingerprinting of Australian isolates of Mycobacterium avium subsp paratuberculosis using IS900 RFLP. Aust Vet J. 2000;78:184–90.
Kaur P, Filia G, Singh SV, Patil PK, Ravi Kumar GV, et al. Molecular epidemiology of Mycobacterium avium subspecies paratuberculosis: IS900 PCR identification and IS1311 polymorphism analysis from ruminants in the Punjab region of India. Comp Immunol Microbiol Infect Dis. 2011;34:163–9.
RY carried out the isolation and molecular genetic studies and drafted the manuscript. CL and YC participated in the survey of the sera. PB participated in the modification of this paper. FL participated in the phylogenetic analysis. LY performed the statistical analysis. DZ and XZ conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
This work was supported by the MoSTRCUK international cooperation project (Project No. 2013DFG32500), National Natural Science Foundation of China (Project No. 31572487), Funding of State Key Lab of Agrobiotechnology (Project No. 2012SKLAB06-14). 2015 CAU Foreign Experts Major Projects (Project No: 2012z018). Hign-end Foreign Experts Recruitment Program (Project No: GDW20151100036).
The authors declare that they have no competing interests.