Who infects whom?-Reconstructing infection chains of Mycobacterium avium ssp. paratuberculosis in an endemically infected dairy herd by use of genomic data

High clonality of Mycobacterium avium subsp. paratuberculosis field isolates from red deer revealed by two different methodological approaches of comparative genomic analysis

Mycobacterium avium subsp. paratuberculosis (MAP) is the aetiological agent of paratuberculosis (Johne's disease) in both domestic and wild ruminants. In the present study, using a whole-genome sequence (WGS) approach, we investigated the genetic diversity of 15 Mycobacterium avium field strains isolated in the last 10 years from red deer inhabiting the Stelvio National Park and affected by paratuberculosis. Combining de novo assembly and a reference-based method, followed by a pangenome analysis, we highlight a very close relationship among 13 MAP field isolates, suggesting that a single infecting event occurred in this population. Moreover, two isolates have been classified as Mycobacterium avium subsp. hominissuis, distinct from the other MAPs under comparison but close to each other. This is the first time that this subspecies has been found in Italy in samples without evident epidemiological correlations, having been isolated in two different locations of the Stelvio National Park and in different years. Our study highlights the importance of a multidisciplinary approach incorporating molecular epidemiology and ecology into traditional infectious disease knowledge in order to investigate the nature of infectious disease in wildlife populations.

Environmental Sources and Transmission of Nontuberculous Mycobacteria

The field of environmental nontuberculous mycobacteria (NTM) is benefiting from a new era of genomics that has catapulted our understanding of preferred niches, transmission, and outbreak investigations. The ability to forecast environmental features that promote or reduce environmental NTM prevalence will greatly improve with coordinated environmental sampling and by elevating the necessity for uniform disease notifications. Studies that synergize environmental biology, isolate notifications, and comparative genomics in prospective, longitudinal studies, particularly during climate changes and weather events, will be useful to solve longstanding NTM public health quandaries.

Genetic Features of Mycobacterium avium subsp. paratuberculosis Strains Circulating in the West of France Deciphered by Whole-Genome Sequencing

Paratuberculosis is a chronic infection of the intestine, mainly the ileum, caused by Mycobacterium avium subsp. paratuberculosis in cattle and other ruminants. This enzootic disease is present worldwide and has a negative impact on the dairy cattle industry. For this subspecies, the current genotyping tools do not provide the needed resolution to investigate the genetic diversity of closely related strains. These limitations can be overcome by the application of whole-genome sequencing (WGS), particularly for clonal populations such as M. avium subsp. paratuberculosis. The purpose of the present study was to undertake a WGS analysis with a panel of 200 animal field M. avium subsp. paratuberculosis strains selected based on a previous large-scale longitudinal study of Prim'Holstein and Normande dairy breeds naturally infected with M. avium subsp. paratuberculosis in the West of France. The pangenome analysis revealed that M. avium subsp. paratuberculosis has a closed pangenome. The phylogeny, based on alignment of 2,786 nonhomoplasic single nucleotide polymorphisms (SNPs), showed that the strain population is structured into three clades independently of the cattle breed or geographic distribution. The increased resolution of phylogeny obtained by WGS confirmed the homoplasic nature of the markers variable-number tandem repeat (VNTR) and short sequence repeat (SSR) used for M. avium subsp. paratuberculosis genotyping. These phylogenetic data also revealed independent introductions of the different genotypes in two main waves since at least 2003. WGS applied to this sampling demonstrated the presence of mixed infections in herds and at the individual animal level. Collectively, the phylogeny results inferred with French isolates compared to M. avium subsp. paratuberculosis isolates from around the world suggest introductions of M. avium subsp. paratuberculosis genotypes through the animal trade. Relationships between genetic traits and epidemiological data can now be investigated to better understand transmission dynamics of the disease. IMPORTANCE Mycobacterium avium subsp. paratuberculosis causes Johne's disease in ruminants, which is present worldwide and has significant negative impacts on the dairy cattle industry and animal welfare. Prevention and control of M. avium subsp. paratuberculosis infection are hampered by knowledge gaps in strain virulence, genotype distribution, and transmission dynamics. This work has revealed new insights into M. avium subsp. paratuberculosis strains currently circulating in western France and how they are related to strains circulating globally. We applied whole-genome sequencing (WGS) to obtain comprehensive information on genome evolution and discrimination of closely related strains. This approach revealed the history of M. avium subsp. paratuberculosis infection in France, refined the pangenomic characteristics of M. avium subsp. paratuberculosis, and demonstrated the existence of mixed infection in animals. Finally, this study identified predominant genotypes, which allow a better understanding of disease transmission dynamics. This information will facilitate tracking of this pathogen on farms and across agricultural regions, thus informing transmission pathways and disease control points.

Can breeders prevent pathogen adaptation when selecting for increased resistance to infectious diseases?

Background

Recent research shows that genetic selection has high potential to reduce the prevalence of infectious diseases in livestock. However, like all interventions that target infectious diseases, genetic selection of livestock can exert selection pressure on pathogen populations. Such selection on the pathogen may lead to escape strategies and reduce the effect of selection of livestock for disease resistance. Thus, to successfully breed livestock for lower disease prevalence, it is essential to develop strategies that prevent the invasion of pathogen mutants that escape host resistance. Here we investigate the conditions under which such “escape mutants” can replace wild-type pathogens in a closed livestock population using a mathematical model of disease transmission.

Results

Assuming a single gene that confers sufficient resistance, results show that genetic selection for resistance in livestock typically leads to an “invasion window” within which an escape mutant of the pathogen can invade. The bounds of the invasion window are determined by the frequency of resistant hosts in the population. The lower bound occurs when the escape mutant has an advantage over the wild-type pathogen in the population. The upper bound occurs when local eradication of the pathogen is expected. The invasion window is smallest when host resistance is strong and when infection with the wild-type pathogen provides cross immunity to infection with the escape mutant.

Conclusions

To minimise opportunities for pathogens to adapt, under the assumptions of our model, the aim of disease control through genetic selection should be to achieve herd-level eradication of the infection faster than the rate of emergence of escape mutants of the pathogen. Especially for microparasitic infections, this could be achieved by placing animals into herds according to their genetic resistance, such that these herds stay completely out of the invasion window. In contrast to classical breeding theory, our model suggests that multi-trait selection with gradual improvement of each trait of the breeding goal might not be the best strategy when resistance to infectious disease is part of the breeding goal. Temporally, combining genetic selection with other interventions helps to make the invasion window smaller, and thereby reduces the risk of invasion of escape mutants.

Quantitative Risk Assessment of Exposure to Mycobacterium avium subsp. paratuberculosis (MAP) via Different Types of Milk for the Slovenian Consumer

This study aimed to assess the risk of exposure to Mycobacterium avium subsp. paratuberculosis (MAP) via milk for the Slovenian consumer. MAP is suspected to be associated with several diseases in humans, therefore the risk of exposure should be better understood. The primary source of MAP for humans is thought to be cattle, in which MAP causes paratuberculosis or Johne’s disease. We developed a stochastic quantitative risk assessment model using Monte Carlo simulations. Considering the assumptions and uncertainties, we estimated the overall risk of exposure to MAP via milk to be low. For people consuming raw milk from MAP positive farms, the risk was high. On-farm pasteurisation reduced the risk considerably, but not completely. The risk of exposure via pasteurised retail milk was most likely insignificant. However, with a higher paratuberculosis prevalence the risk would also increase. Given the popularity of raw milk vending machines and homemade dairy products, this risk should not be ignored. To reduce the risk, consumers should heat raw milk before consumption. To prevent a potential public health scare and safeguard farmers’ livelihoods, a reduction in paratuberculosis prevalence should be sought. Our results show that culling clinically infected cows was insufficient to reduce milk contamination with MAP.

Modeling Paratuberculosis Transmission in a Small Dairy Herd Typical of Slovenia Suggests That Different Models Should Be Used to Study Disease Spread in Herds of Different Sizes

This study aimed to investigate the possible dynamics of paratuberculosis or Johne’s disease in a typical Slovenian dairy herd of about 17 cows. Paratuberculosis is a worldwide endemic disease of cattle caused by Mycobacterium avium subsp. paratuberculosis (MAP) and is associated with significant economic losses. We developed a stochastic compartmental model with two pathways of disease progression, infections of adult cows and infections of young animals through horizontal and vertical transmission, and transmission through animal movements. The average proportions of subclinically and clinically infected cows were 4% and 0.47%, respectively. The prevalence within the herd, which included latently infected animals, averaged 7.13% and ranged from 0% to 70.59%. Under the given circumstances, the results showed a relatively high rate of spontaneous elimination (0.22 per herd per year) of the disease and a high rate of reinfection (0.18 per herd per year) facilitated by active animal trade. To our knowledge, this stochastic compartmental model is the first to be developed specifically to represent a small dairy herd and could apply to other countries with a similar structure of dairy farms. The results suggest that different models should be used to study MAP spread in herds of various sizes.

Mycobacterium avium subsp. paratuberculosis Proteome Changes Profoundly in Milk

Mycobacterium avium subspecies paratuberculosis (MAP) are detectable viable in milk and other dairy products. The molecular mechanisms allowing the adaptation of MAP in these products are still poorly understood. To obtain information about respective adaptation of MAP in milk, we differentially analyzed the proteomes of MAP cultivated for 48 h in either milk at 37 °C or 4 °C or Middlebrook 7H9 broth as a control. From a total of 2197 MAP proteins identified, 242 proteins were at least fivefold higher in abundance in milk. MAP responded to the nutritional shortage in milk with upregulation of 32% of proteins with function in metabolism and 17% in fatty acid metabolism/synthesis. Additionally, MAP upregulated clusters of 19% proteins with roles in stress responses and immune evasion, 19% in transcription/translation, and 13% in bacterial cell wall synthesis. Dut, MmpL4_1, and RecA were only detected in MAP incubated in milk, pointing to very important roles of these proteins for MAP coping with a stressful environment. Dut is essential and plays an exclusive role for growth, MmpL4_1 for virulence through secretion of specific lipids, and RecA for SOS response of mycobacteria. Further, 35 candidates with stable expression in all conditions were detected, which could serve as targets for detection. Data are available via ProteomeXchange with identifier PXD027444.