Advertisement
Research| Volume 104, ISSUE 1, P211-220, January 2021

Genotypic analysis of nontuberculous mycobacteria isolated from raw milk and human cases in Wisconsin

Open ArchivePublished:November 05, 2020DOI:https://doi.org/10.3168/jds.2020-18214

      ABSTRACT

      Nontuberculous mycobacteria (NTM) compose a group of mycobacteria that do not belong to the Mycobacterium tuberculosis complex group. They are frequently isolated from environmental samples such as water, soil, and, to a lesser extent, food samples. Isolates of NTM represent a major health threat to humans worldwide, especially those who have asthma or are immunocompromised. Human disease is acquired from environmental exposures and through consumption of NTM-contaminated food. The most common clinical manifestation of NTM disease in human is lung disease, but lymphatic, skin and soft tissue, and disseminated disease are also important. The main objective of the current study was to profile the farm-level contamination of cow milk with NTM by examining milk filters and bulk tank milk samples. Five different NTM species were isolated in one dairy herd in Wisconsin, with confirmed 16S rRNA genotypes including Mycobacterium fortuitum, Mycobacterium avium ssp. hominissuis, Mycobacterium abscessus, Mycobacterium simiae, and Mycobacterium avium ssp. paratuberculosis (Mycobacterium paratuberculosis). In tank milk samples, M. fortuitum was the predominant species in 48% of the samples, whereas M. chelonae/abscessus and M. fortuitum were the only 2 species obtained from 77 and 23% of the examined filters, respectively. Surprisingly, M. avium ssp. hominissuis, M. paratuberculosis, and M. simiae were isolated from 16.7, 10.4, and 4% of the examined milk samples, respectively, but not from milk filters. Interestingly, NTM isolates from human clinical cases in Wisconsin clustered very closely with those from milk samples. These findings suggest that the problem of NTM contamination is underestimated in dairy herds and could contribute to human infections with NTM. Overall, the study validates the use of bulk tank samples rather than milk filters to assess contamination of milk with NTM. Nontuberculous mycobacteria represent one type of pathogens that extensively contaminate raw milk at the farm level. The significance of our research is in evaluating the existence of NTM at the farm level and identifying a simple approach to examine the potential milk contamination with NTM members using tank milk or milk filters from dairy operations. In addition, we attempted to examine the potential link between NTM isolates found in the farm to those circulating in humans in Wisconsin.

      Key words

      INTRODUCTION

      Nontuberculous mycobacteria (NTM), also known as environmental mycobacteria due to their abundant existence in the environment, are opportunistic pathogens of public health concern, especially for immunocompromised individuals (
      • Claeys T.A.
      • Robinson R.T.
      The many lives of nontuberculous mycobacteria.
      ). Members of NTM are the causative agents of many human infections such as lymphadenitis in children and pulmonary, skin, and soft tissue infections (
      • Nishiuchi Y.
      • Iwamoto T.
      • Maruyama F.
      Infection sources of a common non-tuberculous mycobacterial pathogen, Mycobacterium avium complex.
      ). Although human-to-human transmission is rare, it is essential to identify potential sources and routes of exposure to NTM. Although the risk of NTM transmission through water and several foods is documented (
      • Argueta C.
      • Yoder S.
      • Holtzman A.E.
      • Aronson T.W.
      • Glover N.
      • Berlin O.G.
      • Stelma Jr., G.N.
      • Froman S.
      • Tomasek P.
      Isolation and identification of nontuberculous mycobacteria from foods as possible exposure sources.
      ), the prevalence and risk of NTM in animal products, such as milk, is understudied. Globally, NTM isolates have been detected in milk products from both high- and low-income countries (
      • Yoo J.W.
      • Jo K.W.
      • Kim M.N.
      • Lee S.D.
      • Kim W.S.
      • Kim D.S.
      • Shim T.S.
      Increasing trend of isolation of non-tuberculous mycobacteria in a tertiary university hospital in South Korea.
      ;
      • Kendall B.A.
      • Winthrop K.L.
      Update on the epidemiology of pulmonary nontuberculous mycobacterial infections.
      ;
      • Panagiotou M.
      • Papaioannou A.I.
      • Kostikas K.
      • Paraskeua M.
      • Velentza E.
      • Kanellopoulou M.
      • Filaditaki V.
      • Karagiannidis N.
      The epidemiology of pulmonary nontuberculous mycobacteria: Data from a general hospital in Athens, Greece, 2007–2013.
      ). Nontuberculous mycobacteria are considered emerging milk-borne pathogens worldwide. In Brazil, 21.7% of mozzarella cheese manufactured from raw buffalo milk was positive for NTM, demonstrating a risk for consumers (
      • Jordão Junior, C.M.
      • Lopes F.C.M.
      • David S.
      • Farache Filho A.
      • Leite C.Q.F.
      Detection of nontuberculous mycobacteria from water buffalo raw milk in Brazil.
      ). Another study, in Spain, reported the presence of viable Mycobacterium avium ssp. hominissuis in a powdered infant formula and mycobacterial DNA (M. avium, Mycobacterium tuberculosis complex, and other NTM mycobacteria) was detected in 15% of the examined dairy products (
      • Sevilla I.A.
      • Molina E.
      • Tello M.
      • Elguezabal N.
      • Juste R.A.
      • Garrido J.M.
      Detection of mycobacteria by culture and DNA-based methods in animal-derived food products purchased at Spanish supermarkets.
      ). Finally, Mycobacterium gordonae, Mycobacterium fortuitum, Mycobacterium senegalense, and M. avium have been isolated from fresh milk sold in Nigeria, suggesting a high prevalence of NTM in dairy food (
      • Agada C.A.
      • Adesokan H.K.
      • Igwe D.
      • Cadmus S.I.
      Mycobacterium africanum and nontuberculous mycobacteria from fresh milk of pastoral cattle and soft cheese in Oyo State—implications for public health.
      ). In Taiwan, the incidence of NTM infections increased from 2.7 to 10.2 cases per 100,000 between the years 2000 and 2008 (
      • Lai C.C.
      • Tan C.K.
      • Chou C.H.
      • Hsu H.L.
      • Liao C.H.
      • Huang Y.T.
      • Yang P.C.
      • Luh K.T.
      • Hsueh P.R.
      Increasing incidence of nontuberculous mycobacteria, Taiwan, 2000–2008.
      ). The incidence of NTM infection increased from 0.9 per 100,000 population in 1995 to 2.9 per 100,000 in 2006 in England, Wales, and Northern Ireland (
      • Moore J.E.
      • Kruijshaar M.E.
      • Ormerod L.P.
      • Drobniewski F.
      • Abubakar I.
      Increasing reports of non-tuberculous mycobacteria in England, Wales and Northern Ireland, 1995–2006.
      ). The increase in NTM prevalence highlights the need to examine potential sources of infection and develop protocols to eliminate them. With the increasing popularity of raw milk consumption throughout the United States (allowed in 30 of 50 states;
      • Costard S.
      • Espejo L.
      • Groenendaal H.
      • Zagmutt F.J.
      Outbreak-related disease burden associated with consumption of unpasteurized cow's milk and cheese, United States, 2009–2014.
      ), surveillance of NTM in milk and dairy herds is becoming of a paramount importance. This study was designed to help close this knowledge gap.
      Previously, NTM have been detected in meat and cow milk (raw and pasteurized;
      • Sgarioni S.A.
      • Hirata R.D.
      • Hirata M.H.
      • Leite C.Q.
      • de Prince K.A.
      • de Andrade Leite S.R.
      • Filho D.V.
      • Siqueira V.L.
      • Caleffi-Ferracioli K.R.
      • Cardoso R.F.
      Occurrence of Mycobacterium bovis and non-tuberculous mycobacteria (NTM) in raw and pasteurized milk in the northwestern region of Parana, Brazil.
      ), major food sources for humans. Water, soil, and food have been suggested as important NTM sources for infection transmission to humans (
      • Shitaye J.E.
      • Horvathova A.
      • Bartosova L.
      • Moravkova M.
      • Kaevska M.
      • Donnelly N.
      • Pavlik I.
      Distribution of nontuberculosis mycobacteria in environmental samples from a slaughterhouse and in raw and processed meats.
      ;
      • Sgarioni S.A.
      • Hirata R.D.
      • Hirata M.H.
      • Leite C.Q.
      • de Prince K.A.
      • de Andrade Leite S.R.
      • Filho D.V.
      • Siqueira V.L.
      • Caleffi-Ferracioli K.R.
      • Cardoso R.F.
      Occurrence of Mycobacterium bovis and non-tuberculous mycobacteria (NTM) in raw and pasteurized milk in the northwestern region of Parana, Brazil.
      ). The relationship between environmental sources of NTM and disease in humans has been studied for decades (
      • Wolinsky E.
      Nontuberculous mycobacteria and associated diseases.
      ), but mostly as association studies because of a lack of data on the genetic basis of the relationship. However, with the advancement of molecular techniques, several reports have suggested a stronger link between NTM in food samples and those circulating in patients (
      • Yoder S.
      • Argueta C.
      • Holtzman A.
      • Aronson T.
      • Berlin O.G.
      • Tomasek P.
      • Glover N.
      • Froman S.
      • Stelma Jr., G.
      PCR comparison of Mycobacterium avium isolates obtained from patients and foods.
      ). In California, 25 isolates of NTM, including M. avium, M. gordonae, and Mycobacterium simiae, were recovered from 121 supermarket food items, mainly of plant origin (
      • Argueta C.
      • Yoder S.
      • Holtzman A.E.
      • Aronson T.W.
      • Glover N.
      • Berlin O.G.
      • Stelma Jr., G.N.
      • Froman S.
      • Tomasek P.
      Isolation and identification of nontuberculous mycobacteria from foods as possible exposure sources.
      ). Unfortunately, the published data regarding the prevalence of NTM in food of animal origin are scarce in the United States. In Brazil, NTM was found in 9 to 25% of raw milk samples and related dairy products (
      • Franco M.M.
      • Paes A.C.
      • Ribeiro M.G.
      • de Figueiredo Pantoja J.C.
      • Santos A.C.
      • Miyata M.
      • Leite C.Q.
      • Motta R.G.
      • Listoni F.J.
      Occurrence of mycobacteria in bovine milk samples from both individual and collective bulk tanks at farms and informal markets in the southeast region of Sao Paulo, Brazil.
      ). Similar levels were obtained from analyses of raw and processed meat samples from the Czech Republic (
      • Shitaye J.E.
      • Horvathova A.
      • Bartosova L.
      • Moravkova M.
      • Kaevska M.
      • Donnelly N.
      • Pavlik I.
      Distribution of nontuberculosis mycobacteria in environmental samples from a slaughterhouse and in raw and processed meats.
      ). High levels of NTM contamination were also reported in raw milk samples analyzed in several countries from Europe and Africa (
      • Kazwala R.R.
      • Daborn C.J.
      • Kusiluka L.J.
      • Jiwa S.F.
      • Sharp J.M.
      • Kambarage D.M.
      Isolation of Mycobacterium species from raw milk of pastoral cattle of the Southern Highlands of Tanzania.
      ;
      • Konuk M.
      • Korcan E.
      • Dulgerbaki S.
      • Altindis M.
      Isolation and identification of Mycobacteria from raw milk samples in Afyonkarahisar district of Turkey.
      ). Both PCR and genetic fingerprinting-based approaches (e.g., variable number of tandem repeats, DNA sequencing) were suggested to accelerate the detection of the prevalence and persistence of NTM in food samples (
      • Bolaños C.A.D.
      • Paula C.L.
      • Guerra S.T.
      • Franco M.M.J.
      • Ribeiro M.G.
      Diagnosis of mycobacteria in bovine milk: An overview.
      ). Here, we hypothesized that the diversity of NTM contamination is underestimated in dairy products, which should be of concern to public health authorities. Our study was conducted to develop a simple approach to examine the potential milk contamination with NTM on the farm level using tank milk and milk filters from dairy operations. In addition, we estimated the level of diversity among NTM present in raw milk and examined any potential link to the NTM isolated from patients in Wisconsin.

      MATERIALS AND METHODS

      Milk Sample Collection

      A medium-size local dairy herd with 80 lactating cows was used for this investigation from January to May 2016. Milk filters and raw milk samples were collected from the bulk tank with a total of 85 samples for each (milk filters and tank milk samples = 170 samples). Milk filters are usually placed just before the bulk tank to remove gross foreign material from milk. The 85 raw milk samples (50 mL each) were collected using aseptic technique directly from the bulk tank that contains the milk that has passed through the milk filter. The milk filters and tank milk samples were maintained at 4°C during transport and prepared as described before (
      • Slana I.
      • Kralik P.
      • Kralova A.
      • Babak V.
      • Pavlik I.
      Short communication: Examination of milk filters by real-time PCR as a herd-level indicator of the presence of Mycobacterium avium subspecies paratuberculosis in dairy herds.
      ). Briefly, each filter was placed separately in a stomacher bag with 100 mL of PBS (pH 7.0) containing 0.05% Tween 20 (PBS-T) and homogenized in a stomacher for 2 min. Sample homogenates were centrifuged (15 min at 2,500 × g) and pellets were resuspended in 10 mL of PBS-T for culturing.

      Isolation of Mycobacteria from Filters and Milk Samples

      Ten milliliters from milk samples and 1 mL from the processed filter homogenates were centrifuged (15 min at 2,500 × g) before the decontamination step. The obtained pellets and cream layers were resuspended in 10 mL of 0.75% 1-hexadecylpyridinium chloride (Sigma-Aldrich, St. Louis, MO) and incubated for 5 h at room temperature for decontamination. Samples were centrifuged again (15 min at 2,500 × g) and pellets were resuspended in 1 mL of PBS for inoculation onto Lowenstein Jensen (LJ) medium (LJ slants) and Middlebrook 7H10 agar (BD, Franklin Lakes, NJ) as described before (
      • Dundee L.
      • Grant I.R.
      • Ball H.J.
      • Rowe M.T.
      Comparative evaluation of four decontamination protocols for the isolation of Mycobacterium avium ssp. paratuberculosis from milk.
      ). Cultures were incubated in the presence of 5 to 10% CO2 at 37°C for 90 d and inspected weekly for bacterial growth (
      • Ghosh P.
      • Hsu C.
      • Alyamani E.J.
      • Shehata M.M.
      • Al-Dubaib M.A.
      • Al-Naeem A.
      • Hashad M.
      • Mahmoud O.M.
      • Alharbi K.B.
      • Al-Busadah K.
      • Al-Swailem A.M.
      • Talaat A.M.
      Genome-wide analysis of the emerging infection with Mycobacterium avium subspecies paratuberculosis in the Arabian camels (Camelus dromedarius).
      ). One set of Middlebrook 7H10 plates was supplemented with mycobactin J (MJ, 2 mg/mL), an iron-chelating cell wall component for the isolation of M. paratuberculosis, and 2 µg of vancomycin, 30 µg of amphotericin B, and 20 µg of nalidixic acid (VAN) per milliliter. All antibiotics were from Sigma Chemical Co. (St. Louis, MO) Colonies growing on Mycobacterium-specific media were characterized by Ziehl-Neelsen staining and 16S rRNA PCR. For Ziehl-Neelsen staining, wet smears from each suspected colony were prepared on a clean and dry glass slide and heat fixed over flame until dried. Slides were placed on a 90°C adjusted heat block, and carbol-fuchsin solution was added. After 5 min, slides were washed with distilled water, followed by an additional wash with 1% acid alcohol for discoloration. The slides were counterstained with methylene blue solution for 30 s and washed with distilled water before air drying and microscopical examination.

      NTM Isolates from Human Cases

      With the help of Wisconsin State Laboratory of Hygiene (Madison), we obtained 13 de-identified isolates of NTM that were collected from human patients with diagnosed NTM infections cultured from various sources including sputum, urine, blood, or neck tissues. All isolates were cultured as detailed above and their DNA was extracted for PCR amplification as detailed below.

      Extraction of Genomic DNA

      To obtain high-quality genomic DNA (gDNA), we followed a method previously developed in our group (
      • Talaat A.M.
      • Reimschuessel R.
      • Trucksis M.
      Identification of mycobacteria infecting fish to the species level using polymerase chain reaction and restriction enzyme analysis.
      ). Briefly, mycobacterial cultures were pelleted by centrifugation, and the pellets were washed and resuspended in equal volume of Tris-EDTA buffer. Bacterial suspensions were placed at 80°C for 20 min to kill all living organisms. Tubes were allowed to cool at room temperature, and 10 µL of 100 mg/mL lysozyme was then added to each tube, followed by incubation at 37°C for 3 h with occasional mixing. A solution of 10% SDS and proteinase K (20 mg/mL) was added to each tube at a ratio of 88:12. Tubes were incubated at 65°C for 2 h followed by addition of 100 µL of 5 M NaCl at 65°C for additional 10 min. A 10% solution of cetyl trimethyl ammonium bromide (80 µL; Sigma-Aldrich) was added, followed by mixing and incubation at 65°C for another 10 min. Then, DNA was extracted with an equal volume of phenol:chloroform:isoamylalcohol (25:24:1 vol:vol:vol) and similarly with chloroform:isoamylalcohol (24:1 vol:vol) followed by precipitation with 0.6 volumes of ice-cold isopropanol and centrifugation (14,000 × g) at 4°C for 15 min. Pellets were washed with cold 75% ethanol and dried in a SpeedVac (Thermo Scientific, Wilmington, DE) for 5 min. The DNA was finally resuspended in 50 µL of sterile distilled water. Quality of the gDNA was verified by both NanoDrop (Thermo Scientific) machine and electrophoresis. To detect the presence of mycobacterial DNA in the collected samples, gDNA was extracted directly following our previously described procedure (
      • Talaat A.M.
      • Reimschuessel R.
      • Trucksis M.
      Identification of mycobacteria infecting fish to the species level using polymerase chain reaction and restriction enzyme analysis.
      ).

      Polymerase Chain Reaction

      To confirm the identity of the mycobacterial colonies and gDNA extracted directly from the examined samples, we subjected them to PCR amplification targeting the 16S rRNA gene. A partial sequence (938 bp) of the 16S rRNA gene was amplified from all isolates and samples. Primers used for this purpose were designed by our group (
      • Talaat A.M.
      • Reimschuessel R.
      • Trucksis M.
      Identification of mycobacteria infecting fish to the species level using polymerase chain reaction and restriction enzyme analysis.
      ); namely, AMT 36 (GCGAACGGGTGAGTAACACG) and AMT 37 (TGCACACAGGCCACAAGC). Each PCR amplification was 25 µL, which contained 1 M betaine, 50 mM potassium glutamate, 10 mM Tris-HCl pH 8.8, 0.1% Triton X-100, 2 mM magnesium chloride, 0.2 mM dNTPs, 0.5 mM each primer, 0.5 U of Taq DNA polymerase (Promega, Madison, WI), and 25 ng of genomic DNA. The amplification cycle consisted of an initial denaturation step of 94°C for 5 min; 35 cycles of denaturation at 94°C for 30 s, annealing at 55°C for 30 s, and extension at 72°C for 1 min; and a final extension at 72°C for 7 min. Amplicons were evaluated by electrophoresis in 2% agarose gels prestained with ethidium bromide (0.5 mg/mL). Single band products of amplicons were purified from agarose gels and extracted using Wizard SV Gel and PCR cleanup system (Promega) as detailed before (
      • Talaat A.M.
      • Reimschuessel R.
      • Trucksis M.
      Identification of mycobacteria infecting fish to the species level using polymerase chain reaction and restriction enzyme analysis.
      ).

      Sanger DNA Sequencing and Phylogenetic Analysis

      Following PCR amplification, purified PCR fragments were sequenced with BigDye Terminator v3.1 (Applied Biosystems, Foster City, CA) with AMT 36 and AMT 37 primers, according to the manufacturer's instruction. All sequences were analyzed with BLASTn algorithm on the NCBI web portal (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
      A rooted tree was computed using MEGA7 software (
      • Kumar S.
      • Stecher G.
      • Tamura K.
      MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets.
      ) and neighbor-joining method to align all sequences of 16S rRNA and build a phylogenetic tree with the presence of reference sequences from several mycobacterial species deposited in GenBank (
      • Tamura K.
      • Stecher G.
      • Peterson D.
      • Filipski A.
      • Kumar S.
      MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0.
      ).

      Whole-Genome Sequence Analysis and Phylogenetic Analysis

      Selected isolates from farm tank milk samples and de-identified human isolates from the Wisconsin State Laboratory of Hygiene were submitted to the Wisconsin Biotechnology Center for whole-genome sequencing (WGS) as detailed before by our group (
      • Hsu C.Y.
      • Wu C.W.
      • Talaat A.M.
      Genome-wide sequence variation among Mycobacterium avium subspecies paratuberculosis isolates: A better understanding of Johne's disease transmission dynamics.
      ;
      • Abdelaal H.F.M.
      • Spalink D.
      • Amer A.
      • Steinberg H.
      • Hashish E.A.
      • Nasr E.A.
      • Talaat A.M.
      Genomic polymorphism associated with the emergence of virulent isolates of Mycobacterium bovis in the Nile Delta.
      ). Genomic libraries for a total of 19 isolates were prepared and run on Illumina-MiSeq 2000 platform (Illumina Inc., San Diego, CA). Reads with an average length of 250 bp were assembled using SPAdes 3.14.1 (
      • Bankevich A.
      • Nurk S.
      • Antipov D.
      • Gurevich A.A.
      • Dvorkin M.
      • Kulikov A.S.
      • Lesin V.M.
      • Nikolenko S.I.
      • Pham S.
      • Prjibelski A.D.
      • Pyshkin A.V.
      • Sirotkin A.V.
      • Vyahhi N.
      • Tesler G.
      • Alekseyev M.A.
      • Pevzner P.A.
      SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing.
      ) for de novo genome assembly. The assembled contigs were searched against the NCBI database using BLASTn (https://blast.ncbi.nlm.nih.gov/Blast.cgi) to ensure the identity of each isolate. CLC Bio Workbench version 8.5.1 (https://digitalinsights.qiagen.com/) was used, and we used standard genome corresponding to identification from BLASTn; namely, M. avium 104 (accession number: CP000479.1), M. abscessus ATCC19977 (accession number: NC_010397.1), or M. fortuitum CT6 (accession number: CP011269.1), to identify SNP; insertions/deletions (InDels) event with minimum variation frequency was set to 50% to count SNP. Consensus sequences for mapped reads were exported from CLC Bio workbench version 8.5.1 as FASTA format and aligned using the Harvest software package (
      • Treangen T.J.
      • Ondov B.D.
      • Koren S.
      • Phillippy A.M.
      The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes.
      ) and the “parsnp” function to output a newick tree file that was visualized using MEGA X software.

      Data Availability

      The sequences of the 16S rRNA gene have been deposited in the GenBank database under the accession numbers reported in Table 1. All whole-genome sequences for mycobacterial isolates used in this study were deposited at the NCBI GenBank (BioProject no. PRJNA656902) with BioSample numbers from SAMN15806900 to SAMN15806920. See Table 2 for isolate names, BioSample numbers, and isolate sources.
      Table 1Different nontuberculous mycobacteria (NTM) isolates (M. = Mycobacterium) detected from tank milk and milk filters
      Numbers followed by parenthetical values are the total number of isolates for tank milk samples or milk filters and the percentage of each NTM isolate. The collected tank milk and milk filter samples were cultured on 3 selective media including Lowenstein-Jensen (LJ) medium, 7H10 Middlebrook medium, and 7H10 Middlebrook medium supplemented with Mycobactin J (MJ).
      SourceNTM isolate
      M. fortuitumM. aviumssp.hominissuisM. abscessusM. simiaeM. aviumssp.paratuberculosis
      Tank milk2381025
       Accession numberMN880273–MN880298MN880384–MN880391MN880315–MN880324MN880382–MN880383MN880377–MN880381
       LJ medium8 (34.8)4 (50)1 (10)
       Middlebrook medium10 (43.5)4 (50)4 (40)1 (50)
       MJ Middlebrook medium5 (21.7)5 (50)1 (50)5 (100)
      Milk filters310
       Accession numberMN880309–MN880311MN880300–MN880308 and MN880272(F39)
       LJ medium
       Middlebrook medium
       MJ Middlebrook medium3 (100)10 (100)
      1 Numbers followed by parenthetical values are the total number of isolates for tank milk samples or milk filters and the percentage of each NTM isolate. The collected tank milk and milk filter samples were cultured on 3 selective media including Lowenstein-Jensen (LJ) medium, 7H10 Middlebrook medium, and 7H10 Middlebrook medium supplemented with Mycobactin J (MJ).
      Table 2Summary data of whole-genome sequence analysis of Mycobacterium (M.) species isolates from milk tank samples and human specimens in Wisconsin
      Syn = synonymous; non-syn = nonsynonymous; MNV = multiple nucleotide variant.
      Sample IDBioSample no.ReferenceMapping %Average coverageNo. of SNPNo. of syn SNPNo. of non-syn SNPNo. of intragenic SNPNo. of MNVNo. of replacementsNo. of insertionsNo. of deletions
      W9SAMN15806900M. avium10482.59309.8136,65023,7399,1303,7811,54093473475
      W11SAMN15806902M. avium10482.10454.2732,55021,8957,5803,0751,769107482446
      W12SAMN15806903M. avium10472.86295.3423,23616,2975,4551,4841,96994398391
      W14SAMN15806906M. avium10486.50412.0837,41124,7998,9613,6511,811106538532
      W15SAMN15806905M. avium10477.23243.0232,62521,8647,7553,0061,836106491446
      T9SAMN15806917M. avium10467.1748.1537,93824,7449,1154,0791,898113566598
      T52bSAMN15806914M. avium10484.85128.338,91025,5479,4843,8792,030132578609
      T62bSAMN15806915M. avium10487.69153.5338,65625,3249,4263,9061,981133585624
      T84SAMN15806916M. avium10488.8179.0837,15024,1949,1333,8231,655120564587
      T2SAMN15806913M. avium10488.197.9947,90926,2299,1003,5692,204136631595
      W7SAMN15806907M. abscessusATCC 1997791.28378.06103,14876,85818,3047,9863,489166675781
      W8SAMN15806908M. abscessusATCC 1997898.74425.4652332103076
      T7SAMN15806920M. abscessusATCC 1997793.53128.18103,44877,10918,3887,9513,648213743983
      T29CSAMN15806918M. abscessusATCC 1997793.3975.14103,26976,95918,3637,9473,560196726962
      T84bSAMN15806919M. abscessusATCC 1997792.87162.65103,61677,22318,4427,9513,6952097621007
      W4SAMN15806909M. fortuitumCT679.30296.11152,070114,50527,13710,4288,3864031,2591,295
      W5SAMN15806910M. fortuitumCT683.60336.867,50613,59949,2144,6933,625170635594
      W6SAMN15806911M. fortuitumCT668.40238.671,97955,01813,8815,0804,271192730685
      T60aSAMN15806912M. fortuitumCT684.3051.968,12641,47113,5115,1473,973219815766
      1 Syn = synonymous; non-syn = nonsynonymous; MNV = multiple nucleotide variant.

      RESULTS

      Milk Contamination with NTM

      We hypothesized that examination of milk filters and tank milk would provide a means of detecting mycobacteria contamination on the herd level. Although milk filters trap debris, typically in the range of 100 to 150 μm, they do not prevent pathogenic bacteria (~1–10 μm) from passing through. Culture-positive milk filters likely indicate that tank milk has been contaminated, but the reverse might not always be correct. Detecting mycobacterial DNA directly from milk filters and the examined milk samples resulted in a total of 48 positive milk filter samples out of 85 (56.5%) with almost the same percentage for the tank milk (44 samples out of 85, 51.8%).
      Using culturing on mycobacteria-specific media, tank milk and milk filter samples were positive in 17 (20%) and 10 (11.7%) samples, including 48 and 13 mycobacterial isolates (Table 1), respectively, with a total of 61 independent isolates cultured from the examined dairy herd.
      A total of 5 mycobacterial species were isolated from the tank milk samples and only 2 were obtained from the milk filters (Table 1), an indication of the validity of sampling tank milk to assess herd health and milk contamination with NTM. All isolates grown on the selected media were positive for acid fast staining.

      Genotypes of NTM Present in Raw Milk

      To confirm isolate identity, we amplified the 16S rRNA gene from all isolates. As expected, specific DNA bands of approximately 938 bp were observed on agarose gel for all positive samples, indicating that all isolates were members of the genus Mycobacterium. To identify the species of mycobacterial isolates, amplicons from all isolates were sequenced using Sanger sequencing (
      • Talaat A.M.
      • Reimschuessel R.
      • Trucksis M.
      Identification of mycobacteria infecting fish to the species level using polymerase chain reaction and restriction enzyme analysis.
      ). DNA sequences from all isolates were analyzed using the BLASTn (https://blast.ncbi.nlm.nih.gov/Blast.cgi) algorithm against all bacterial genomes present in the GenBank database (
      • NCBI Resource Coordinators
      Database resources of the National Center for Biotechnology Information.
      ). From the tank milk samples cultured on LJ slants, 13 independent isolates were grown. Among these isolates, 8 belonged to the M. fortuitum group, 4 belonged to M. avium ssp. hominissuis, and 1 was identified as M. abscessus (Table 1). The other set of tank milk samples was cultured on 7H10 medium, and 19 isolates were grown. Among them, M. fortuitum was the predominant species (n = 10). Four isolates belonged to the M. avium ssp. hominissuis group, while 1 isolate was identified as M. abscessus and another was identified as M. simiae. In addition, when 7H10 medium was supplemented with MJ, 16 isolates were grown. Among them, 5 isolates belonged to each group for M. avium ssp. paratuberculosis, M. abscessus, and M. fortuitum, and 1 isolate was identified as M. simiae. For the milk filter samples, both LJ slants and 7H10 Middlebrook medium without antibiotics were heavily contaminated, a potential cause of isolation failure; however, all isolates were obtained from Middlebrook 7H10 medium supplemented with MJ and the antibiotics mixture (VAN). A total of 13 isolates were recovered from milk filters. Among them, 10 isolates belonged to the M. abscessus group, while the remaining 3 isolates belonged to the M. fortuitum group (Table 1).

      Phylogenetic Analysis of the NTM

      To better examine the relationship between NTM isolates, we examined the evolutionary distances among isolates based on 16S rRNA sequences. As expected, clear phylogenetic clusters formed between members of the M. abscessus/chelonae complex and the M. avium ssp. hominissuis/paratuberculosis complex. No significant genetic distance was found between members of the M. abscessus/chelonae complex whether they were isolated from tank milk or milk filters, suggesting a common source of contamination. All isolates of M. fortuitum occupied a different cluster between the other 2 clusters. As expected, 2 isolates of M. simiae clustered independently, another confirmation of their identity and potential different source of transmission. However, the low bootstrapping percentage (<50%) that distinguished the M. simiae from the M. avium/paratuberculosis complex could indicate their evolutionary relatedness (Figure 1).
      Figure thumbnail gr1
      Figure 1Evolutionary tree of all nontuberculous mycobacteria isolates from tank milk (identified as T samples) and milk filters (identified as F samples) isolated from a Wisconsin dairy herd. The tree is based on 16S rRNA sequences using the neighbor-joining method in MEGA7 software (https://www.megasoftware.net/). The tree is drawn to scale with a bar representing the phylogenetic distance. Numbers next to the branches represent percentages of 1,000 bootstrap test.
      In another analysis, we tested the phylogenetic relatedness between selected farm NTM isolates (n = 9) and those circulating in human patients (n = 10) from Wisconsin, using de-identified samples from the collection of the Wisconsin State Laboratory of Hygiene. In Table 2, we added a summary of the WGS of 19 NTM isolates from both milk and human samples. Single nucleotide polymorphisms for each isolate were compared with genomes representing the standard strains of M. avium ssp. hominissuis, M. fortuitum, and M. abscessus (Figure 2b). These isolates were the most commonly isolated mycobacteria from tank milk samples. Interestingly, the whole-genome sequences from all human isolates were clustered with isolates from milk samples, an indication of highly similar genotypes of NTM present in both milk samples and those circulating in humans infected with NTM (Figure 2). Moreover, the overall SNP 16S rRNA genotyping did not accurately identify 3 isolates of M. fortuitum (T2, T9, and T84), which WGS identified as M. abscessus for T2 and M. avium ssp. hominissuis for T9 and T84, respectively. On the other hand, T52b and T62b isolates that were identified by 16S rRNA as M. avium ssp. hominissuis were identified by WGS as M. abscessus and M. fortuitum, respectively. This finding further confirms the superiority of WGS for strain identification compared with other standard genotyping approaches.
      Figure thumbnail gr2
      Figure 2Evolutionary analysis of mycobacterial isolates from tank milk samples and human specimens. (a) Evolutionary tree based on whole-genome sequence analysis of 22 nontuberculous mycobacteria (NTM) isolates representing all groups detected in a Wisconsin dairy herd (designated with T) and human isolates from Wisconsin patients (designated with W). The unrooted tree is based on whole-genome sequences using Maximum Unique Match algorithm in the Harvest package and visualized by MEGA X software (https://www.megasoftware.net/). The nodes represent unique matches between bacterial genomes. The value on each branch is the posterior probability, showing the percentage support for the following node. (b) A heat map of NTM isolates displaying percentage of the total SNPs with synonymous (Syn) and nonsynonymous (Non-Syn) polymorphism when each isolate is compared with its own standard strain. The heat map was generated using GraphPad Prism version 7 (GraphPad Software Inc., La Jolla, CA).

      DISCUSSION

      Nontuberculous mycobacteria are a group of opportunistic mycobacterial species that do not belong to the M. tuberculosis complex (
      • Brode S.K.
      • Marchand-Austin A.
      • Jamieson F.B.
      • Marras T.K.
      Pulmonary versus nonpulmonary nontuberculous mycobacteria, Ontario, Canada.
      ). These mycobacterial species have been isolated from many environmental sources, such as water, soil, and food (
      • Nishiuchi Y.
      • Iwamoto T.
      • Maruyama F.
      Infection sources of a common non-tuberculous mycobacterial pathogen, Mycobacterium avium complex.
      ), implicating these sources as vehicles for NTM transmission to humans (
      • Sgarioni S.A.
      • Hirata R.D.
      • Hirata M.H.
      • Leite C.Q.
      • de Prince K.A.
      • de Andrade Leite S.R.
      • Filho D.V.
      • Siqueira V.L.
      • Caleffi-Ferracioli K.R.
      • Cardoso R.F.
      Occurrence of Mycobacterium bovis and non-tuberculous mycobacteria (NTM) in raw and pasteurized milk in the northwestern region of Parana, Brazil.
      ). Previously, analysis of raw milk and dairy products indicated the presence of high levels of NTM in some developing countries (
      • Kazwala R.R.
      • Daborn C.J.
      • Kusiluka L.J.
      • Jiwa S.F.
      • Sharp J.M.
      • Kambarage D.M.
      Isolation of Mycobacterium species from raw milk of pastoral cattle of the Southern Highlands of Tanzania.
      ;
      • Konuk M.
      • Korcan E.
      • Dulgerbaki S.
      • Altindis M.
      Isolation and identification of Mycobacteria from raw milk samples in Afyonkarahisar district of Turkey.
      ). Unfortunately, in tuberculosis-endemic countries (e.g., India, China), the problem of NTM is usually masked by problems associated with tuberculosis (
      • Gopinath K.
      • Singh S.
      Non-tuberculous mycobacteria in TB-endemic countries: Are we neglecting the danger?.
      ), suggesting an underestimation of NTM prevalence. Immunocompromised individuals are at high risk of contracting NTM infections, especially with the opportunistic members of the M. avium complex (
      • Shojaei H.
      • Heidarieh P.
      • Hashemi A.
      • Feizabadi M.M.
      • Daei Naser A.
      Species identification of neglected nontuberculous mycobacteria in a developing country.
      ). The current study was conducted to examine the extent of NTM presence at the farm level, using a dairy herd in Wisconsin, where mycobacterial prevalence is low (
      • Wisconsin Department of Health Services
      Tuberculosis statistics.
      ) compared with other locations (
      • Gopinath K.
      • Singh S.
      Non-tuberculous mycobacteria in TB-endemic countries: Are we neglecting the danger?.
      ). The milk filters are usually located just before the bulk tank and have pores of 65 to 75 μm to filter debris and fecal material from the milk (
      • McKee R.
      • Grant I.R.
      • Rowe M.T.
      • Buckley H.G.
      • Buckley J.F.
      • Fanning S.
      Examination of in-line filters to detect Mycobacterium avium ssp. paratuberculosis infection at farm level.
      ). Microbial analysis of milk filters represents a first-choice sampling strategy to assess the overall microbial quality of milk, hence we compared the utility of milk filters with that of tank milk samples to assess mycobacterial presence in unpasteurized milk. Unexpectedly, our analysis of both sampling strategies indicated the superiority of bulk tank milk sampling over milk filters in both the number of isolates and diversity of isolated species. This result is likely due to high level of contamination that could not be trapped by milk filters.
      One of the major findings obtained in this study is that M. fortuitum was the predominate NTM isolate from tank milk as shown in Table 1, while M. chelonae/abscessus was the prevailing isolate detected from milk filters. In an earlier study of NTM isolates, the M. avium complex, specifically M. avium ssp. paratuberculosis (M. paratuberculosis), represented the highest prevalence, followed by M. fortuitum (
      • Gopinath K.
      • Singh S.
      Non-tuberculous mycobacteria in TB-endemic countries: Are we neglecting the danger?.
      ). In this report, M. paratuberculosis was isolated from 5 tank milk samples but was not detected at all from any of the milk filters, most likely because of its low level in the examined herd. Similar studies reported the presence of M. paratuberculosis in raw milk samples worldwide (
      • Grant I.R.
      • Ball H.J.
      • Rowe M.T.
      Incidence of Mycobacterium paratuberculosis in bulk raw and commercially pasteurized cows' milk from approved dairy processing establishments in the United Kingdom.
      ;
      • Ayele W.Y.
      • Svastova P.
      • Roubal P.
      • Bartos M.
      • Pavlik I.
      Mycobacterium avium subspecies paratuberculosis cultured from locally and commercially pasteurized cow's milk in the Czech Republic.
      ;
      • Slana I.
      • Kralik P.
      • Kralova A.
      • Pavlik I.
      On-farm spread of Mycobacterium avium ssp. paratuberculosis in raw milk studied by IS900 and F57 competitive real time quantitative PCR and culture examination.
      ;
      • Botsaris G.
      • Slana I.
      • Liapi M.
      • Dodd C.
      • Economides C.
      • Rees C.
      • Pavlik I.
      Rapid detection methods for viable Mycobacterium avium subspecies paratuberculosis in milk and cheese.
      ;
      • Shankar H.
      • Singh S.V.
      • Singh P.K.
      • Singh A.V.
      • Sohal J.S.
      • Greenstein R.J.
      Presence, characterization, and genotype profiles of Mycobacterium avium subspecies paratuberculosis from unpasteurized individual and pooled milk, commercial pasteurized milk, and milk products in India by culture, PCR, and PCR-REA methods.
      ).
      Related studies indicated the presence of NTM in both raw and pasteurized milk in Brazil including M. fortuitum, Mycobacterium marinum, Mycobacterium kansasii, and M. gordonae (
      • Leite C.Q.
      • Anno I.S.
      • Leite S.R.
      • Roxo E.
      • Morlock G.P.
      • Cooksey R.C.
      Isolation and identification of mycobacteria from livestock specimens and milk obtained in Brazil.
      ). In other studies, M. chelonae and Mycobacterium scrofulaceum were also isolated from raw and pasteurized milk (
      • Sgarioni S.A.
      • Hirata R.D.
      • Hirata M.H.
      • Leite C.Q.
      • de Prince K.A.
      • de Andrade Leite S.R.
      • Filho D.V.
      • Siqueira V.L.
      • Caleffi-Ferracioli K.R.
      • Cardoso R.F.
      Occurrence of Mycobacterium bovis and non-tuberculous mycobacteria (NTM) in raw and pasteurized milk in the northwestern region of Parana, Brazil.
      ). In Turkey, other members of the NTM group were detected in raw milk including, Mycobacterium terrae, M. kansasii, Mycobacterium haemophilum, and Mycobacterium agri (
      • Konuk M.
      • Korcan E.
      • Dulgerbaki S.
      • Altindis M.
      Isolation and identification of Mycobacteria from raw milk samples in Afyonkarahisar district of Turkey.
      ). In our study, 5 NTM species involved in a variety of human illnesses were detected including M. fortuitum, M. chelonae/abscessus, M. paratuberculosis, M. avium ssp. hominissuis, and M. simiae.
      To better identify each isolate, we sequenced a large proportion of the 16S rRNA gene, which confirmed the identity of all mycobacterial isolates (Figure 1) and provided a tool for further analysis of such isolates in the context of other circulating mycobacterial species. Interestingly, phylogenetic analysis of milk isolate sequences and their counterparts from human cases of infection with NTM species circulating in Wisconsin patients indicated a close relationship between both groups of NTM (Figure 2) when whole genomes were analyzed. Such analysis suggested a potential common source of NTM isolates, most likely potable tap water as indicated before (
      • Honda J.R.
      • Virdi R.
      • Chan E.D.
      Global environmental nontuberculous mycobacteria and their contemporaneous man-made and natural niches.
      ) for M. avium complex but not for M. paratuberculosis, M. abscessus, and M. fortuitum complexes. It is possible that animal infections with such organisms are transmissible to human through the consumption of raw milk or dairy products made with unpasteurized milk. Results presented here should prompt examination of more milk and dairy products for mycobacterial contamination, especially with the increase of consumption of raw milk throughout the United States (
      • Costard S.
      • Espejo L.
      • Groenendaal H.
      • Zagmutt F.J.
      Outbreak-related disease burden associated with consumption of unpasteurized cow's milk and cheese, United States, 2009–2014.
      ) and in low-income countries as well. Moreover, contradictory results from 16S rRNA genotyping and WGS of members of M. fortuitum and M. avium groups suggested the need to use more targets for genotyping (e.g., rpoB and hsp65) or the need to use WGS.
      Overall, our study illustrates that examination of milk samples could help in the evaluation of mycobacterial diversity in a dairy herd, with bulk tank milk sampling being the superior strategy. The study validates the use of bulk tank samples over milk filters to assay milk contamination with NTM. Moreover, consumption of raw milk or unpasteurized dairy products constitutes a potential risk factor for contracting NTM infections in humans. The implementation of control measures that reduce or prevent contamination of raw milk with NTM during and after milking is of a paramount importance to avoid human illness associated with NTM infection.

      ACKNOWLEDGMENTS

      We also thank David Warshauer from the Wisconsin State Laboratory of Hygiene (Madison) for providing the NTM isolates from human cases. This work was partially supported by a USDA (Washington, DC) grant (NIFA-2018-67015-28243), United States Agency for International Development and National Academy of Sciences subaward no. 2000010565, and animal formula fund (no. WIS02015) awarded to AT. Any findings, conclusions, or recommendations expressed in this article are those of the authors alone, and do not necessarily reflect the views of USAID or NAS. In addition, Egyptian Culture and Educational Bureau Awards were made to ZA and MH. The authors have not stated any conflicts of interest.

      REFERENCES

        • Abdelaal H.F.M.
        • Spalink D.
        • Amer A.
        • Steinberg H.
        • Hashish E.A.
        • Nasr E.A.
        • Talaat A.M.
        Genomic polymorphism associated with the emergence of virulent isolates of Mycobacterium bovis in the Nile Delta.
        Sci. Rep. 2019; 9 (31406159)11657
        • Agada C.A.
        • Adesokan H.K.
        • Igwe D.
        • Cadmus S.I.
        Mycobacterium africanum and nontuberculous mycobacteria from fresh milk of pastoral cattle and soft cheese in Oyo State—implications for public health.
        Afr. J. Med. Med. Sci. 2014; 43 (26949776): 13-20
        • Argueta C.
        • Yoder S.
        • Holtzman A.E.
        • Aronson T.W.
        • Glover N.
        • Berlin O.G.
        • Stelma Jr., G.N.
        • Froman S.
        • Tomasek P.
        Isolation and identification of nontuberculous mycobacteria from foods as possible exposure sources.
        J. Food Prot. 2000; 63 (10914663): 930-933
        • Ayele W.Y.
        • Svastova P.
        • Roubal P.
        • Bartos M.
        • Pavlik I.
        Mycobacterium avium subspecies paratuberculosis cultured from locally and commercially pasteurized cow's milk in the Czech Republic.
        Appl. Environ. Microbiol. 2005; 71 (15746320): 1210-1214
        • Bankevich A.
        • Nurk S.
        • Antipov D.
        • Gurevich A.A.
        • Dvorkin M.
        • Kulikov A.S.
        • Lesin V.M.
        • Nikolenko S.I.
        • Pham S.
        • Prjibelski A.D.
        • Pyshkin A.V.
        • Sirotkin A.V.
        • Vyahhi N.
        • Tesler G.
        • Alekseyev M.A.
        • Pevzner P.A.
        SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing.
        J. Comput. Biol. 2012; 19: 455-477
        • Bolaños C.A.D.
        • Paula C.L.
        • Guerra S.T.
        • Franco M.M.J.
        • Ribeiro M.G.
        Diagnosis of mycobacteria in bovine milk: An overview.
        Rev. Inst. Med. Trop. São Paulo. 2017; 59 (28591268): e40
        • Botsaris G.
        • Slana I.
        • Liapi M.
        • Dodd C.
        • Economides C.
        • Rees C.
        • Pavlik I.
        Rapid detection methods for viable Mycobacterium avium subspecies paratuberculosis in milk and cheese.
        Int. J. Food Microbiol. 2010; 141 (20381185): S87-S90
        • Brode S.K.
        • Marchand-Austin A.
        • Jamieson F.B.
        • Marras T.K.
        Pulmonary versus nonpulmonary nontuberculous mycobacteria, Ontario, Canada.
        Emerg. Infect. Dis. 2017; 23 (29048292): 1898-1901
        • Claeys T.A.
        • Robinson R.T.
        The many lives of nontuberculous mycobacteria.
        J. Bacteriol. 2018; 200 (29483164)e00739-17
        • Costard S.
        • Espejo L.
        • Groenendaal H.
        • Zagmutt F.J.
        Outbreak-related disease burden associated with consumption of unpasteurized cow's milk and cheese, United States, 2009–2014.
        Emerg. Infect. Dis. 2017; 23 (28518026): 957-964
        • Dundee L.
        • Grant I.R.
        • Ball H.J.
        • Rowe M.T.
        Comparative evaluation of four decontamination protocols for the isolation of Mycobacterium avium ssp. paratuberculosis from milk.
        Lett. Appl. Microbiol. 2001; 33 (11555198): 173-177
        • Franco M.M.
        • Paes A.C.
        • Ribeiro M.G.
        • de Figueiredo Pantoja J.C.
        • Santos A.C.
        • Miyata M.
        • Leite C.Q.
        • Motta R.G.
        • Listoni F.J.
        Occurrence of mycobacteria in bovine milk samples from both individual and collective bulk tanks at farms and informal markets in the southeast region of Sao Paulo, Brazil.
        BMC Vet. Res. 2013; 9 (23618368): 85
        • Ghosh P.
        • Hsu C.
        • Alyamani E.J.
        • Shehata M.M.
        • Al-Dubaib M.A.
        • Al-Naeem A.
        • Hashad M.
        • Mahmoud O.M.
        • Alharbi K.B.
        • Al-Busadah K.
        • Al-Swailem A.M.
        • Talaat A.M.
        Genome-wide analysis of the emerging infection with Mycobacterium avium subspecies paratuberculosis in the Arabian camels (Camelus dromedarius).
        PLoS One. 2012; 7 (22393374)e31947
        • Gopinath K.
        • Singh S.
        Non-tuberculous mycobacteria in TB-endemic countries: Are we neglecting the danger?.
        PLoS Negl. Trop. Dis. 2010; 4 (20436962): e615
        • Grant I.R.
        • Ball H.J.
        • Rowe M.T.
        Incidence of Mycobacterium paratuberculosis in bulk raw and commercially pasteurized cows' milk from approved dairy processing establishments in the United Kingdom.
        Appl. Environ. Microbiol. 2002; 68 (11976118): 2428-2435
        • Honda J.R.
        • Virdi R.
        • Chan E.D.
        Global environmental nontuberculous mycobacteria and their contemporaneous man-made and natural niches.
        Front. Microbiol. 2018; 9 (30214436)2029
        • Hsu C.Y.
        • Wu C.W.
        • Talaat A.M.
        Genome-wide sequence variation among Mycobacterium avium subspecies paratuberculosis isolates: A better understanding of Johne's disease transmission dynamics.
        Front. Microbiol. 2011; 2 (22163226): 236
        • Jordão Junior, C.M.
        • Lopes F.C.M.
        • David S.
        • Farache Filho A.
        • Leite C.Q.F.
        Detection of nontuberculous mycobacteria from water buffalo raw milk in Brazil.
        Food Microbiol. 2009; 26 (19527843): 658-661
        • Kazwala R.R.
        • Daborn C.J.
        • Kusiluka L.J.
        • Jiwa S.F.
        • Sharp J.M.
        • Kambarage D.M.
        Isolation of Mycobacterium species from raw milk of pastoral cattle of the Southern Highlands of Tanzania.
        Trop. Anim. Health Prod. 1998; 30 (9760715): 233-239
        • Kendall B.A.
        • Winthrop K.L.
        Update on the epidemiology of pulmonary nontuberculous mycobacterial infections.
        Semin. Respir. Crit. Care Med. 2013; 34 (23460008): 87-94
        • Konuk M.
        • Korcan E.
        • Dulgerbaki S.
        • Altindis M.
        Isolation and identification of Mycobacteria from raw milk samples in Afyonkarahisar district of Turkey.
        Int. J. Food Microbiol. 2007; 115 (17303274): 343-347
        • Kumar S.
        • Stecher G.
        • Tamura K.
        MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets.
        Mol. Biol. Evol. 2016; 33: 1870-1874
        • Lai C.C.
        • Tan C.K.
        • Chou C.H.
        • Hsu H.L.
        • Liao C.H.
        • Huang Y.T.
        • Yang P.C.
        • Luh K.T.
        • Hsueh P.R.
        Increasing incidence of nontuberculous mycobacteria, Taiwan, 2000–2008.
        Emerg. Infect. Dis. 2010; 16 (20113563): 294-296
        • Leite C.Q.
        • Anno I.S.
        • Leite S.R.
        • Roxo E.
        • Morlock G.P.
        • Cooksey R.C.
        Isolation and identification of mycobacteria from livestock specimens and milk obtained in Brazil.
        Mem. Inst. Oswaldo Cruz. 2003; 98 (12886409): 319-323
        • McKee R.
        • Grant I.R.
        • Rowe M.T.
        • Buckley H.G.
        • Buckley J.F.
        • Fanning S.
        Examination of in-line filters to detect Mycobacterium avium ssp. paratuberculosis infection at farm level.
        in: Juste R.A. Geijo M.V. Garrido J.M. Proc. 7th Int. Colloq. Paratuberculosis. International Association for Paratuberculosis Inc., Madison, WI2002: 258-260
        • Moore J.E.
        • Kruijshaar M.E.
        • Ormerod L.P.
        • Drobniewski F.
        • Abubakar I.
        Increasing reports of non-tuberculous mycobacteria in England, Wales and Northern Ireland, 1995–2006.
        BMC Public Health. 2010; 10 (20950421): 612
        • NCBI Resource Coordinators
        Database resources of the National Center for Biotechnology Information.
        Nucleic Acids Res. 2016; 44: D7-D19
        • Nishiuchi Y.
        • Iwamoto T.
        • Maruyama F.
        Infection sources of a common non-tuberculous mycobacterial pathogen, Mycobacterium avium complex.
        Front. Med. (Lausanne). 2017; 4 (28326308): 27
        • Panagiotou M.
        • Papaioannou A.I.
        • Kostikas K.
        • Paraskeua M.
        • Velentza E.
        • Kanellopoulou M.
        • Filaditaki V.
        • Karagiannidis N.
        The epidemiology of pulmonary nontuberculous mycobacteria: Data from a general hospital in Athens, Greece, 2007–2013.
        Pulm. Med. 2014; 2014 (25132991)894976
        • Sevilla I.A.
        • Molina E.
        • Tello M.
        • Elguezabal N.
        • Juste R.A.
        • Garrido J.M.
        Detection of mycobacteria by culture and DNA-based methods in animal-derived food products purchased at Spanish supermarkets.
        Front. Microbiol. 2017; 8 (28649235)1030
        • Sgarioni S.A.
        • Hirata R.D.
        • Hirata M.H.
        • Leite C.Q.
        • de Prince K.A.
        • de Andrade Leite S.R.
        • Filho D.V.
        • Siqueira V.L.
        • Caleffi-Ferracioli K.R.
        • Cardoso R.F.
        Occurrence of Mycobacterium bovis and non-tuberculous mycobacteria (NTM) in raw and pasteurized milk in the northwestern region of Parana, Brazil.
        Braz. J. Microbiol. 2014; 45 (25242962): 707-711
        • Shankar H.
        • Singh S.V.
        • Singh P.K.
        • Singh A.V.
        • Sohal J.S.
        • Greenstein R.J.
        Presence, characterization, and genotype profiles of Mycobacterium avium subspecies paratuberculosis from unpasteurized individual and pooled milk, commercial pasteurized milk, and milk products in India by culture, PCR, and PCR-REA methods.
        Int. J. Infect. Dis. 2010; 14 (19576834): e121-e126
        • Shitaye J.E.
        • Horvathova A.
        • Bartosova L.
        • Moravkova M.
        • Kaevska M.
        • Donnelly N.
        • Pavlik I.
        Distribution of nontuberculosis mycobacteria in environmental samples from a slaughterhouse and in raw and processed meats.
        Czech J. Food Sci. 2009; 27: 194-202
        • Shojaei H.
        • Heidarieh P.
        • Hashemi A.
        • Feizabadi M.M.
        • Daei Naser A.
        Species identification of neglected nontuberculous mycobacteria in a developing country.
        Jpn. J. Infect. Dis. 2011; 64 (21788699): 265-271
        • Slana I.
        • Kralik P.
        • Kralova A.
        • Babak V.
        • Pavlik I.
        Short communication: Examination of milk filters by real-time PCR as a herd-level indicator of the presence of Mycobacterium avium subspecies paratuberculosis in dairy herds.
        J. Dairy Sci. 2012; 95 (22365199): 1162-1165
        • Slana I.
        • Kralik P.
        • Kralova A.
        • Pavlik I.
        On-farm spread of Mycobacterium avium ssp. paratuberculosis in raw milk studied by IS900 and F57 competitive real time quantitative PCR and culture examination.
        Int. J. Food Microbiol. 2008; 128 (18824269): 250-257
        • Talaat A.M.
        • Reimschuessel R.
        • Trucksis M.
        Identification of mycobacteria infecting fish to the species level using polymerase chain reaction and restriction enzyme analysis.
        Vet. Microbiol. 1997; 58 (9453133): 229-237
        • Tamura K.
        • Stecher G.
        • Peterson D.
        • Filipski A.
        • Kumar S.
        MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0.
        Mol. Biol. Evol. 2013; 30 (24132122): 2725-2729
        • Treangen T.J.
        • Ondov B.D.
        • Koren S.
        • Phillippy A.M.
        The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes.
        Genome Biol. 2014; 15 (25410596): 524
        • Wisconsin Department of Health Services
        Tuberculosis statistics.
        https://www.dhs.wisconsin.gov/tb/data.htm
        Date: 2019
        Date accessed: April , 2019
        • Wolinsky E.
        Nontuberculous mycobacteria and associated diseases.
        Am. Rev. Respir. Dis. 1979; 119 (369415): 107-159
        • Yoder S.
        • Argueta C.
        • Holtzman A.
        • Aronson T.
        • Berlin O.G.
        • Tomasek P.
        • Glover N.
        • Froman S.
        • Stelma Jr., G.
        PCR comparison of Mycobacterium avium isolates obtained from patients and foods.
        Appl. Environ. Microbiol. 1999; 65 (10347056): 2650-2653
        • Yoo J.W.
        • Jo K.W.
        • Kim M.N.
        • Lee S.D.
        • Kim W.S.
        • Kim D.S.
        • Shim T.S.
        Increasing trend of isolation of non-tuberculous mycobacteria in a tertiary university hospital in South Korea.
        Tuberc. Respir. Dis. (Seoul). 2012; 72 (23101005): 409-415