Hannah Trewby

Whole Genome Sequencing Reveals Local Transmission Patterns of Mycobacterium bovis in Sympatric Cattle and Badger Populations

Whole genome sequencing (WGS) technology holds great promise as a tool for the forensic epidemiology of bacterial pathogens. It is likely to be particularly useful for studying the transmission dynamics of an observed epidemic involving a largely unsampled ‘reservoir’ host, as for bovine tuberculosis (bTB) in British and Irish cattle and badgers. BTB is caused by Mycobacterium bovis, a member of the M. tuberculosis complex that also includes the aetiological agent for human TB. In this study, we identified a spatio-temporally linked group of 26 cattle and 4 badgers infected with the same Variable Number Tandem Repeat (VNTR) type of M. bovis. Single-nucleotide polymorphisms (SNPs) between sequences identified differences that were consistent with bacterial lineages being persistent on or near farms for several years, despite multiple clear whole herd tests in the interim. Comparing WGS data to mathematical models showed good correlations between genetic divergence and spatial distance, but poor correspondence to the network of cattle movements or within-herd contacts. Badger isolates showed between zero and four SNP differences from the nearest cattle isolate, providing evidence for recent transmissions between the two hosts. This is the first direct genetic evidence of M. bovis persistence on farms over multiple outbreaks with a continued, ongoing interaction with local badgers. However, despite unprecedented resolution, directionality of transmission cannot be inferred at this stage. Despite the often notoriously long timescales between time of infection and time of sampling for TB, our results suggest that WGS data alone can provide insights into TB epidemiology even where detailed contact data are not available, and that more extensive sampling and analysis will allow for quantification of the extent and direction of transmission between cattle and badgers.

Limitations of variable number of tandem repeat typing identified through whole genome sequencing of Mycobacterium avium subsp. paratuberculosis on a national and herd level.

BackgroundMycobacterium avium subsp. paratuberculosis (MAP), the causative bacterium of Johne’s disease in dairy cattle, is widespread in the Canadian dairy industry and has significant economic and animal welfare implications. An understanding of the population dynamics of MAP can be used to identify introduction events, improve control efforts and target transmission pathways, although this requires an adequate understanding of MAP diversity and distribution between herds and across the country. Whole genome sequencing (WGS) offers a detailed assessment of the SNP-level diversity and genetic relationship of isolates, whereas several molecular typing techniques used to investigate the molecular epidemiology of MAP, such as variable number of tandem repeat (VNTR) typing, target relatively unstable repetitive elements in the genome that may be too unpredictable to draw accurate conclusions. The objective of this study was to evaluate the diversity of bovine MAP isolates in Canadian dairy herds using WGS and then determine if VNTR typing can distinguish truly related and unrelated isolates.ResultsPhylogenetic analysis based on 3,039 SNPs identified through WGS of 124 MAP isolates identified eight genetically distinct subtypes in dairy herds from seven Canadian provinces, with the dominant type including over 80% of MAP isolates. VNTR typing of 527 MAP isolates identified 12 types, including “bison type” isolates, from seven different herds. At a national level, MAP isolates differed from each other by 1–2 to 239–240 SNPs, regardless of whether they belonged to the same or different VNTR types. A herd-level analysis of MAP isolates demonstrated that VNTR typing may both over-estimate and under-estimate the relatedness of MAP isolates found within a single herd.ConclusionsThe presence of multiple MAP subtypes in Canada suggests multiple introductions into the country including what has now become one dominant type, an important finding for Johne’s disease control. VNTR typing often failed to identify closely and distantly related isolates, limiting the applicability of using this typing scheme to study the molecular epidemiology of MAP at a national and herd-level.

Analysis of the long control region of bovine papillomavirus type 1 associated with sarcoids in equine hosts indicates multiple cross-species transmission events and phylogeographical structure.

Papillomaviruses are a family of slowly evolving DNA viruses and their evolution is commonly linked to that of their host species. However, whilst bovine papillomavirus-1 (BPV-1) primarily causes warts in its natural host, the cow, it can also cause locally aggressive and invasive skin tumours in equids, known as sarcoids, and thus provides a rare contemporary example of cross-species transmission of a papillomavirus. Here, we describe the first phylogenetic analysis of BPV-1 in equine sarcoids to our knowledge, allowing us to explore the evolutionary history of BPV-1 and investigate its cross-species association with equids. A phylogenetic analysis of the BPV-1 transcriptional promoter region (the long control region or LCR) was conducted on 15 bovine and 116 equine samples from four continents. Incorporating previous estimates for evolutionary rates in papillomavirus implied that the genetic diversity in the LCR variants was ancient and predated domestication of both equids and cattle. The phylogeny demonstrated geographical segregation into an ancestral group (African, South American and Australian samples), and a more recently derived, largely European clade. Whilst our data are consistent with BPV-1 originating in cattle, we found evidence of multiple, probably relatively recent, cross-species transmission events into horses. We also demonstrated the high prevalence of one particular sequence variant (variant 20), and suggest this may indicate that this variant shows a fitness advantage in equids. Although strong host specificity remains the norm in papillomaviruses, our results demonstrate that exceptions to this rule exist and can become epidemiologically relevant.

Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis

Highlights • We performed whole genome sequencing (WGS) of Mycobacterium bovis for a single molecular (VNTR) type.• Under-sampling of one lineage was caused by switching between VNTR-types.• Pairwise SNP distances showed a weak genetic isolation by distance pattern.• Bayesian phylogeographic inference was feasible despite a low substitution rate.• WGS studies of M. bovis need to account for slow evolution and molecular type switching.

What has molecular epidemiology ever done for wildlife disease research? Past contributions and future directions

The increasing availability of novel molecular techniques has transformed the study of human health and disease epidemiology. However, uptake of such approaches has been more conservative in the field of wildlife disease epidemiology. We consider the reasons for this and discuss current and potential applications of molecular techniques in a variety of relevant areas within the field of wildlife disease research. These include conducting wildlife disease surveillance, identifying sources of pathogen emergence, uncovering host-pathogen dynamics and managing current outbreaks, including the development and monitoring of wildlife vaccines. We highlight key examples of applications of molecular epidemiological approaches to wildlife disease scenarios and draw parallels from human disease research to suggest potential future directions. The potential value of next generation sequencing technologies to the field of wildlife disease research is discussed, and initial applications are highlighted, balanced against consideration of the challenges involved. Using a wide range of examples drawn from research into human, livestock and wildlife diseases, we demonstrate the value of using molecular epidemiological approaches at all scales of wildlife disease research, from pathogen strains circulating at a global scale to intra-individual host-pathogen dynamics. The potential future contribution of these technologies to the field of wildlife disease epidemiology is substantial. In particular, they are likely to play an increasingly important role in helping us to address a principal challenge in the management of wildlife diseases which is how to tease apart the transmission dynamics of complex multi-host systems in order to develop effective and sustainable interventions.

Genome-Wide Diversity and Phylogeography of Mycobacterium avium subsp. paratuberculosis in Canadian Dairy Cattle.

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative bacterium of Johne’s disease (JD) in ruminants. The control of JD in the dairy industry is challenging, but can be improved with a better understanding of the diversity and distribution of MAP subtypes. Previously established molecular typing techniques used to differentiate MAP have not been sufficiently discriminatory and/or reliable to accurately assess the population structure. In this study, the genetic diversity of 182 MAP isolates representing all Canadian provinces was compared to the known global diversity, using single nucleotide polymorphisms identified through whole genome sequencing. MAP isolates from Canada represented a subset of the known global diversity, as there were global isolates intermingled with Canadian isolates, as well as multiple global subtypes that were not found in Canada. One Type III and six “Bison type” isolates were found in Canada as well as one Type II subtype that represented 86% of all Canadian isolates. Rarefaction estimated larger subtype richness in Québec than in other Canadian provinces using a strict definition of MAP subtypes and lower subtype richness in the Atlantic region using a relaxed definition. Significant phylogeographic clustering was observed at the inter-provincial but not at the intra-provincial level, although most major clades were found in all provinces. The large number of shared subtypes among provinces suggests that cattle movement is a major driver of MAP transmission at the herd level, which is further supported by the lack of spatial clustering on an intra-provincial scale.

Application of high-throughput sequencing to whole rabies viral genome characterisation and its use for phylogenetic re-evaluation of a raccoon strain incursion into the province of Ontario.

Raccoon rabies remains a serious public health problem throughout much of the eastern seaboard of North America due to the urban nature of the reservoir host and the many challenges inherent in multi-jurisdictional efforts to administer co-ordinated and comprehensive wildlife rabies control programmes. Better understanding of the mechanisms of spread of rabies virus can play a significant role in guiding such control efforts. To facilitate a detailed molecular epidemiological study of raccoon rabies virus movements across eastern North America, we developed a methodology to efficiently determine whole genome sequences of hundreds of viral samples. The workflow combines the generation of a limited number of overlapping amplicons covering the complete viral genome and use of high throughput sequencing technology. The value of this approach is demonstrated through a retrospective phylogenetic analysis of an outbreak of raccoon rabies which occurred in the province of Ontario between 1999 and 2005. As demonstrated by the number of single nucleotide polymorphisms detected, whole genome sequence data were far more effective than single gene sequences in discriminating between samples and this facilitated the generation of more robust and informative phylogenies that yielded insights into the spatio-temporal pattern of viral spread. With minor modification this approach could be applied to other rabies virus variants thereby facilitating greatly improved phylogenetic inference and thus better understanding of the spread of this serious zoonotic disease. Such information will inform the most appropriate strategies for rabies control in wildlife reservoirs.

The impact of sheep breed on the risk of classical scrapie

The risk of classical scrapie in sheep is associated with polymorphisms in the prion protein (PrP) gene. In recent years, large-scale selective breeding programmes for sheep at lower risk of disease have been undertaken across the European Union. We analysed large-scale datasets on scrapie and sheep demography to investigate additional effects of sheep breed on scrapie risk. There was evidence for variation between certain breeds in the scrapie risk of some PrP genotypes, which could be caused by innate breed differences or distinct scrapie strains circulating within them. While the PrP genotypes of cases are generally consistent across breeds, some exceptions provide evidence that scrapie strain may influence affected PrP genotypes to a greater extent than innate breed differences. There was a significant association between the breed-level incidence of scrapie and the frequency of susceptible PrP genotypes in breeds. Our results lend support to selective breeding programmes which aim to reduce the frequency of high-risk PrP genotypes with measures not varying by sheep breed.

Processes underlying rabies virus incursions across US–Canada Border as revealed by whole-genome phylogeography

Disease control programs aim to constrain and reduce the spread of infection. Human disease interventions such as wildlife vaccination play a major role in determining the limits of a pathogen’s spatial distribution. Over the past few decades, a raccoon-specific variant of rabies virus (RRV) has invaded large areas of eastern North America. Although expansion into Canada has been largely prevented through vaccination along the US border, several outbreaks have occurred in Canada. Applying phylogeographic approaches to 289 RRV whole-genome sequences derived from isolates collected in Canada and adjacent US states, we examined the processes underlying these outbreaks. RRV incursions were attributable predominantly to systematic virus leakage of local strains across areas along the border where vaccination has been conducted but also to single stochastic events such as long-distance translocations. These results demonstrate the utility of phylogeographic analysis of pathogen genomes for understanding transboundary outbreaks.

Relative abundance of Mycobacterium bovis molecular types in cattle: a simulation study of potential epidemiological drivers

BackgroundThe patterns of relative species abundance are commonly studied in ecology and epidemiology to provide insights into underlying dynamical processes. Molecular types (MVLA-types) of Mycobacterium bovis, the causal agent of bovine tuberculosis, are now routinely recorded in culture-confirmed bovine tuberculosis cases in Northern Ireland. In this study, we use ecological approaches and simulation modelling to investigate the distribution of relative abundances of MVLA-types and its potential drivers. We explore four biologically plausible hypotheses regarding the processes driving molecular type relative abundances: sampling and speciation; structuring of the pathogen population; historical changes in population size; and transmission heterogeneity (superspreading).ResultsNorthern Irish herd-level MVLA-type surveillance shows a right-skewed distribution of MVLA-types, with a small number of types present at very high frequencies and the majority of types very rare. We demonstrate that this skew is too extreme to be accounted for by simple neutral ecological processes. Simulation results indicate that the process of MVLA-type speciation and the manner in which the MVLA-typing loci were chosen in Northern Ireland cannot account for the observed skew. Similarly, we find that pathogen population structure, assuming for example a reservoir of infection in a separate host, would drive the relative abundance distribution in the opposite direction to that observed, generating more even abundances of molecular types. However, we find that historical increases in bovine tuberculosis prevalence and/or transmission heterogeneity (superspreading) are both capable of generating the skewed MVLA-type distribution, consistent with findings of previous work examining the distribution of molecular types in human tuberculosis.ConclusionAlthough the distribution of MVLA-type abundances does not fit classical neutral predictions, our simulations show that increases in pathogen population size and/or superspreading are consistent with the pattern observed, even in the absence of selective pressures acting on the system.