key: cord-348841-qxkmngyk authors: Kozakiewicz, Christopher P.; Burridge, Christopher P.; Funk, W. Chris; VandeWoude, Sue; Craft, Meggan E.; Crooks, Kevin R.; Ernest, Holly B.; Fountain‐Jones, Nicholas M.; Carver, Scott title: Pathogens in space: Advancing understanding of pathogen dynamics and disease ecology through landscape genetics date: 2018-07-28 journal: Evol Appl DOI: 10.1111/eva.12678 sha: doc_id: 348841 cord_uid: qxkmngyk Landscape genetics has provided many insights into how heterogeneous landscape features drive processes influencing spatial genetic variation in free‐living organisms. This rapidly developing field has focused heavily on vertebrates, and expansion of this scope to the study of infectious diseases holds great potential for landscape geneticists and disease ecologists alike. The potential application of landscape genetics to infectious agents has garnered attention at formative stages in the development of landscape genetics, but systematic examination is lacking. We comprehensively review how landscape genetics is being used to better understand pathogen dynamics. We characterize the field and evaluate the types of questions addressed, approaches used and systems studied. We also review the now established landscape genetic methods and their realized and potential applications to disease ecology. Lastly, we identify emerging frontiers in the landscape genetic study of infectious agents, including recent phylogeographic approaches and frameworks for studying complex multihost and host‐vector systems. Our review emphasizes the expanding utility of landscape genetic methods available for elucidating key pathogen dynamics (particularly transmission and spread) and also how landscape genetic studies of pathogens can provide insight into host population dynamics. Through this review, we convey how increasing awareness of the complementarity of landscape genetics and disease ecology among practitioners of each field promises to drive important cross‐disciplinary advances. its formal inception in 2003, facilitated by technological advances that have increased the availability of molecular and landscape data in conjunction with more powerful computational and analytical approaches. Landscape genetics is fuelled by a steady stream of new ideas and methodologies, which, while exciting, can contribute to a lack of consensus or consistency in some key aspects. These aspects include the formulation of research questions, sampling strategies, analytical methods (Balkenhol, Waits, & Dezzani, 2009; Richardson, Brady, Wang, & Spear, 2016; Wagner & Fortin, 2013) and even the identity of the field itself (Dyer, 2015; Storfer et al., 2007) . In fact, landscape genetics has yet to develop its own comprehensive, unifying theory for linking spatial and temporal landscape heterogeneity to genetic variation . While these issues are expected to be remedied as the field matures, many suggestions have been made to facilitate this progress. These have included calls for an increase in cross-disciplinary collaboration (Balkenhol, Gugerli et al., 2009) and an expansion of the scope of landscape genetic research beyond its current emphasis on vertebrates Dyer, 2015) and, particularly, mammals (Kozakiewicz, Carver, & Burridge, 2018) . One logical avenue for cross-disciplinary expansion of landscape genetics is in disease ecology (Biek & Real, 2010) . Elucidating the specific influences of landscape features on pathogen transmission can provide key insights into the processes that affect disease risk and incidence. However, accomplishing this has been a challenge for disease ecologists (Ostfeld, Glass, & Keesing, 2005) . Indeed, the field of spatial epidemiology has only recently begun to emphasize the use of explicit landscape approaches in studies of spatial heterogeneity in infectious disease (i.e., "landscape epidemiology"; Ostfeld et al., 2005; Meentemeyer, Haas, & Václavík, 2012) . A major challenge for the study of landscape epidemiology, a field which does not traditionally implement genetic approaches, is that it is typically dependent on the ability to identify the location and timing of transmission events such that they can be compared to landscape features of interest. Transmission events are essentially impossible to observe, so disease ecologists often assume that contacts between infected and susceptible individuals are a reasonable proxy for transmission. Such contacts generally must be inferred indirectly using methods such as proximity collars, mark-recapture or telemetry, often using spatial overlap as a proxy for contact (Craft & Caillaud, 2011) . These methods are logistically challenging to employ, and whether an inferred contact resulted in transmission is uncertain (Craft, 2015) . Further, much landscape epidemiological research uses infection or exposure data to indicate past transmission, but these methods provide static snapshots of pathogen prevalence and may be inappropriate for inferring how transmission or spread has occurred (or is occurring) over time (Meentemeyer et al., 2012) . The spatial distribution and movement of hosts are major factors affecting the likelihood, timing and spatial patterns of pathogen transmission and spread (Dougherty, Seidel, Carlson, Spiegel, & Getz, 2018) . Landscape genetics can identify landscape factors that are important drivers of host population structure. These landscape factors can determine the spatial configuration of a population, its density, its connectivity with other populations, its demographic structure and its genetic health-all of which have implications for the dynamics of microorganisms infecting the host species (Ellis, Václavík, & Meentemeyer, 2010; Prentice, Marion, White, Davidson, & Hutchings, 2014; Spielman, Brook, Briscoe, & Frankham, 2004) . Further, pathogen dynamics can be inferred directly using pathogen genetic data (Archie, Luikart, & Ezenwa, 2009; DeCandia, Dobson, & vonHoldt, 2018) and incorporated into landscape genetic analyses. Understanding specifically how infectious agents respond to the influence of landscape factors on hosts enables us to predict how such agents might spread based on present landscape configurations, as well as under potential future landscape scenarios (Real & Biek, 2007) . This knowledge can subsequently inform management efforts at the population level (such as vaccination targeted at key regions, culling), as well as broader decisions relating to the management of the landscape itself, which is a key aim of landscape genetics generally (Manel & Holderegger, 2013; Segelbacher et al., 2010) . Landscape genetics is being applied by managers at relatively low rates compared to related ecological fields such as landscape ecology, conservation biology and telemetry research (Bowman et al., 2016) . Therefore, studies that contribute to the management of disease agents within populations could increase the practical impacts of landscape genetics significantly. However, the conceptual underpinnings of pathogen landscape genetics are not fully developed, and the methodologies employed are diverse and potentially confusing for new practitioners. Here, we investigate how landscape genetic techniques are being used to better understand dynamics of microorganisms infecting host species. In conducting this review, we aim to both advocate and facilitate landscape genetic research involving disease-causing organisms. We first evaluate the use of landscape genetics in disease ecology, including the types of questions addressed, the approaches used and the infectious agents studied. We then review established landscape genetic methods and their realized and potential applications to disease ecology. At last, we identify emerging frontiers in the landscape genetic study of pathogens that hold significant potential for advancing research in this field. Landscape genetics was first implemented in the study of rabies virus by Real et al. (2005) , offering an approach to overcome many feasibility issues associated with understanding landscape influences on pathogen transmission. The landscape genetic approach to studying disease was later reviewed by Biek and Real (2010) , who were optimistic about its growth and future use. In particular, they noted that microparasites, such as viruses, are well-suited to landscape genetic study due to their rapid mutation rate and potential spatial genetic structure that can be compared to heterogeneous landscape features at fine temporal and spatial scales. Analyses could be conducted using both pathogenic organisms and agents that do not cause significant diseases in their hosts (Biek, Drummond, & Poss, 2006) . They also identified that methodologies such as GIS, which are commonly employed both in the wider landscape genetics literature and in spatial studies of infectious disease, had not been widely implemented in molecular epidemiology (Archie et al., 2009) . Further, other popular landscape genetic tools, such as those focused on differential landscape permeability (e.g., least-cost paths), were greatly underused despite compatibility with pathogen spatial genetic data. Similar to landscape genetics, landscape epidemiology is an interdisciplinary field undergoing rapid development driven by technological advancements, and arguably still working to develop clear directions for future research (Meentemeyer et al., 2012) . It is therefore likely that the interface of these two fields (i.e., where landscape genetics is used in epidemiology) is similarly challenged, perhaps to the extent that its potential is remaining unrealized. We thus believe it is timely to revisit the body of research that combines landscape genetics and landscape epidemiology, leveraging the work done both prior and subsequent to Biek and Real's (2010) earlier review into clear directions for future research. We conducted a literature search in February 2018 using the ISI Web of Science database with the following terms: TS=(("landscape genetic*" OR "landscape genom*") AND (disease* OR pathogen* OR parasit* OR virus* OR virol* OR epidem* OR infect* OR transmi*)) The search returned 133 results. We read each article and retained the 51 empirical papers that used landscape genetic methods to address questions related to pathogens (see Supporting Information Appendix S1). We excluded reviews (n = 15), meeting abstracts (n = 1), purely methods-based papers (n = 6) and articles that identified as or mentioned landscape genetics but did not sufficiently incorporate landscape factors or genetic data into the study (n = 32), studies that referred to any of our pathogen-related search terms without it being a primary motivation for the study (n = 21), and studies that used words like "transmit" or "parasite" outside of the context of infectious agents (such as the transmission of behaviours) (n = 6). One paper was excluded due to a lack of access at our institutions. Studies that qualitatively discussed landscape with respect to genetic variation were kept, although one might argue that landscape genetics requires quantitative testing of landscape effects. We classified each paper according to the type of host system studied (plant, wild animal, domestic animal and human), the type of pathogen studied (bacterium, protozoan, virus, prion, fungus, macroparasite and transmissible cancer) and the source of genetic data (host, pathogen and vector), and we estimated the severity of disease that each studied pathogen causes in its sampled host or vector. We also categorized each article according to its general conceptual approach. Most examples described in this study were found in our literature search, while several other examples were cited by papers from our search and subsequently also discussed here. Following publication of the first study using landscape genetics to investigate disease in 2005, there was little further research in this area until 2009, which saw a rapid increase in the number of publications (Figure 1a ). This increase coincided with two prominent review articles (Archie et al., 2009; Biek & Real, 2010 ) that were strong proponents of a landscape genetics approach to disease ecology and expressed optimism about its future use. The rate of publication has remained relatively steady (and arguably low) since then, with none of the subsequent 7 years recording more publications than in 2009, when six papers were published. However, 10 articles using landscape genetics to investigate disease were published in 2017, potentially indicating increasing interest in this area of research. A majority of studies (27 of 51) used genetic data from the host for comparison with landscape features (Figure 1b ). This is likely because DNA is easier to obtain from larger, free-living hosts than for pathogens, and methods for genotyping and characterizing host spatial genetic variation are more familiar to landscape geneticists, who predominantly study free-living organisms (Storfer, Murphy, Spear, Holderegger, & Waits, 2010) . Among pathogens that are associated with a particular animal vector, the vector is often genotyped (9 of 14 studies of vector-borne diseases), as vectors such as ticks or mosquitos are also easily sampled, and vector gene flow can be used as a proxy for pathogen spread. Vectors can be targeted for population control as a means of limiting pathogen spread, which makes their study of immediate relevance to wildlife and livestock managers (Townson et al., 2005) . Pathogen genetic data are used in only 16 of 51 pathogen landscape genetic studies, which was somewhat surprising considering that the pathogen is the primary motivation behind many of the reviewed studies. One study included both host and pathogen genetic data (Talbot, Vonhof, Broders, Fenton, & Keyghobadi, 2017) . Viruses were the most frequently studied type of infectious agent (14 of 51 studies; Figure 1c ). In general, viruses evolve more rapidly than other microparasites, which makes them well-suited to study of genetic variation for inference of transmission history (Archie et al., 2009; Grenfell et al., 2004) . However, a majority of landscape genetic studies involving viruses used host genetic data, potentially reflecting the relative difficulty of obtaining viral data, which we discuss later in this section. Instead, the high representation of viruses is largely due to the considerable effort devoted to studying rabies, which comprised half of all landscape genetic studies on viral systems. Rabies is one of the most well-known wildlife pathogens globally, due to its negative impacts on wildlife, domestic animal and human health (Gordon et al., 2004) . Large outbreaks have occurred in North American and European wildlife in recent years, where considerable resources have been devoted to its management (Holmala & Kauhala, 2006; Slate et al., 2009) . Animals infected with rabies also often exhibit behavioural changes that may make them easier to identify (Lefèvre et al., 2009) , potentially aiding sampling of infected individuals. F I G U R E 1 Papers using landscape genetic approaches for the study of infectious agents. (a) Number of publications per year that met our search criteria. (b) Number of publications using genetic data from each of the host, agent or vector species. (c) Number of publications studying pathogens by type, with genetic data source indicated for each type ("unspecified" typically involves studies of a hypothetical agent or estimates of overall pathogen exposure, such as inferred by immune-linked loci). (d) Number of publications adopting each of our broadly identified conceptual approaches for applying landscape genetics to the study of pathogens/ infectious agents-using host/vector genetics to predict agent spread, using host/vector genetics to explain agent spread/distribution and using pathogen genetics to directly study agent spread We broadly define three distinct conceptual approaches by which landscape genetics has been used to study infectious agents ( Figure 1d ). These are the prediction of agent spread using genetic information from the host or vector; the use of host or vector genetic information to explain existing spatial variation in infection risk or prevalence; and the use of genetic information from the infectious agent to directly study transmission and spread. The remainder of this section will address each of these approaches in turn. Because the spread of many microparasites (particularly directly was unrelated to landscape features tested, determining that current rabies oral vaccination plans should be expanded given the high potential for long-distance host movement. In another rabies study, landscape genetics was used to characterize striped skunk dispersal across riverine and highway barriers to assess their utility as barriers to pathogen spread (Talbot, Garant, Paquette, Mainguy, & Pelletier, 2012) . Using host or vector genetic data to predict pathogen spread is attractive as it avoids sampling of the agent itself, which may be substantially more difficult, especially in wildlife populations. Identification of infected hosts often requires laboratory testing and may require specific, potentially invasive sampling approaches (e.g., necropsy) for accurate diagnosis. In addition, extensive sampling may be required to obtain adequate sample sizes when prevalence is low and must be conducted strategically to capture spatial heterogene- Therefore, studies using host or vector data alone have limitations for inferring or predicting pathogen spread, or lack thereof, directly. However, host landscape genetic studies can provide indications of the potential risk of spread of infectious agents, and the understanding gained about host movements can inform subsequent studies of pathogen dynamics. Spatial variation in pathogen prevalence or infection risk can be represented in much the same way as any landscape variable , making spatial data relating to presence of an infectious agent well-suited for incorporation into host landscape genetic models. While spatial heterogeneity in pathogen prevalence could also be considered a component of the landscape that may influence spatial genetic variation in the host, typically only adaptive loci are investigated in this context. More commonly, host neutral genetic variation is used to explain spatial patterns of infection risk or prevalence. A prominent example is a study of chronic wasting disease (CWD) in white-tailed deer. Blanchong et al. (2008) found that populations with lower CWD prevalence showed higher genetic differentiation from those that had high CWD prevalence. This genetic differentiation was found to be associated with roads and rivers, which were likely barriers to both host gene flow and CWD spread. These inferences have subsequently informed and been verified by additional landscape epidemiological research (Robinson, Samuel, Rolley, & Shelton, 2013) . Spatial heterogeneity in pathogen infection risk can also drive microevolutionary responses in the host (Epstein et al., 2016; Monello et al., 2017) . Host species are constantly being challenged by parasitic organisms, which, if not overcome, cause disease and can have fitness consequences. This can create strong selection that acts on various genes, and geographic variation in selection at loci that are known to be associated with adaptive immune genes may reflect variation in pathogen pressure, and individual infection or disease risk (Fumagalli et al., 2011) . This variation may be tested for association with environmental features such as temperature, humidity or urbanization (Tonteri, Vasemägi, Lumme, & Primmer, 2010) , enabling insights into how future changes in climate or land use might influence overall pathogen prevalence. Using the sampled disease agent as the source of genetic data is the most direct way to infer pathogen spread across landscapes, but can be challenging to accomplish. Genetic material may be absent from, or uninformative in some infectious agents, such as prions or clonally transmissible cancers, necessitating genetic analysis of the host (Kelly et al., 2014; Storfer et al., 2017) . In addition to the aforementioned difficulties with pathogen diagnosis, pathogen nucleic acid can be difficult to isolate from samples taken from the host or vector and would ideally be present in the blood, saliva or other easily collected sample. Samples may also require enrichment to obtain sufficient quantities of genetic material for analysis, which can be difficult to accomplish for many pathogens, particularly viruses. However, genetic information from viruses may be particularly useful for molecular epidemiologic analyses due to their rapid mutation rate that can closely infer transmission history (Archie et al., 2009; Brunker, Hampson, Horton, & Biek, 2012) . Further, viruses are prominent emerging pathogens and have relatively small genomes, aiding whole genome-analysis. There are a variety of methods available for implementing landscape genetics, some designed specifically for landscape genetics, while others have been adapted from other fields. The rapid development of landscape genetics means that new methods are regularly emerging, and it is difficult to comprehensively review all of them. However, there are some well-established methodological approaches that have either seen wide use for some time or are becoming increasingly popular at the cutting edge of the field . We describe the approaches (Table 1) and discuss their implementation in the study of pathogen transmission and spread. In landscape genetics, simulation models are usually agent-based and spatially explicit (Landguth, Cushman, & Balkenhol, 2016) . Genetic data are modelled for individuals which have discrete spatial locations with respect to one another and with respect to environmental heterogeneity. Individuals move, behave and reproduce according to their own attributes in response to other individuals and in response to the simulated environment, and the model simulates changes in allele frequencies in response to these parameters. Landscape genetic simulation modelling has been used to test and validate methodological approaches (Cushman, Wasserman, Landguth, & Shirk, 2013; Zeller et al., 2016) , address theoretical questions about how and why landscape heterogeneity influences genetics (Landguth et al., 2010) , and evaluate and explain empirical observations (Shirk, Cushman, & Landguth, 2012) . Further, simulation modelling can predict how a system might respond to certain changes, such as habitat fragmentation or future management activities. Simulation modelling has been widely implemented in the study of pathogenic and nonpathogenic disease, beginning with medical research in the 1960s (Elveback & Varma, 1965 (Rees et al., 2008) . The spread of particular host genes relevant to disease can also be simulated to inform management efforts. For instance, Landguth, Holden, Mahalovich, and Cushman (2017) used landscape genetic simulations to determine optimal planting regimes to maximize the spread of blister rust resistant genes among whitebark pine populations. Such simulations could undoubtedly be applied to vector species in particular, such as predicting the spread of pesticide resistance genes in mosquitos (Chang et al., 2016) and selecting appropriate sites for introduction of genetically modified vectors (Lavery, Harrington, & Scott, 2008) . In addition, with the need to develop further landscape genetic frameworks for the study of pathogens, simulation modelling can prove useful in testing and validating these techniques, as it has done in the broader landscape genetics field (Cushman et al., 2013; Zeller et al., 2016) . For example, Leo, Gonzalez, Millien, and Cristescu (2016) used landscape genetic simulations to validate their multitaxa integrated landscape genetic framework, which appears to be a promising solution to the challenge of studying pathogens with multiple hosts and/or vectors. Landscape genetic simulations may also include epidemiological parameters such as mortality or activity responses to infection, or limited infectious periods, which may otherwise confound conventional (i.e., nonsimulation) landscape genetic approaches. Edge detection methods, such as Monmonier's maximum difference algorithm, (Monmonier, 1973) have also been used to detect landscape barriers to transmission in pathogen studies (Carrel et al., 2015; Joannon et al., 2010) . Ancestry estimates from model-based clustering algorithms can assign individuals to their populations of origin, enabling inference of landscape barrier permeability through the identification of migrants and thus estimation of the risk of pathogen spread across the barrier. Most of the studies implementing clustering and assignment methods did not use approaches that incorporate environmental data. Instead, spatially or nonspatially explicit methods were typically used to identify genetic discontinuities and relationships with landscape barriers were inferred ad hoc, or analyses proceeded to entirely different methods that explicitly include environmental data. Associations between genetic discontinuities and landscape be applied to pathogens directly without these potential constraints. Resistance surfaces are commonly used in landscape genetics for modelling hypotheses concerning the influence of landscape features (from GIS landscape variables) on functional connectivity using techniques such as least-cost paths (Adriaensen et al., 2003) or circuit theory (McRae, Dickson, Keitt, & Shah, 2008) . These techniques produce measures of landscape or "effective" distance among populations or individuals for each hypothesis, which can be tested against observed genetic variation. The primary applications of resistance surface modelling in landscape genetics have been the identification of dispersal corridors and predicting the impacts of landscape and environmental change, such as habitat fragmentation or climate change, on connectivity. Similar to that, landscape genetic resistance surfaces can identify transmission corridors or future patterns of spread (e.g., Streicker et al., 2016) , and such tools have been identified previously as having great utility for pathogen landscape genetic studies (Biek & Real, 2010) . However, resistance surface modelling remains infrequently applied among pathogen studies. Careful consideration is required for identifying the most relevant landscape variables to be tested and correctly parameterizing (assigning costs to) the resistance surface(s) so that these variables are represented in a biologically meaningful way. Developing landscape resistance hypotheses for transmitted agents may be more difficult as their interactions with the landscape are often indirect, mediated by the ecology of hosts and vectors. Pathogen ecological niche models offer an empirical approach for constructing resistance surfaces based on ecological factors influencing pathogen prevalence Fountain-Jones, Pearse et al., 2017) , but these also may not adequately represent host/vector movements. Our literature search returned only one study that explicitly modelled landscape resistance based on pathogen-specific biology, testing elevation (as a proxy for temperature) as a predictor of Plasmodium spread, in addition to resistance surfaces that modelled human movements and mosquito vector ecology (Lo et al., 2017 Graph theoretical approaches, which describe connections (edges) between discrete objects (nodes) (Newman, 2003) , are a flexible yet powerful tool for use in landscape genetics (Dyer, Nason, & Garrick, 2010; Garroway, Bowman, Carr, & Wilson, 2008) . In landscape genetics, nodes can represent individuals, populations or habitat patches, possessing genetic parameters such as diversity measures (Dyer et al., 2010) , or landscape parameters such as percentage habitat or habitat quality (Murphy, Dezzani, Pilliod, & Storfer, 2010) . Similar to that, edges can represent genetic relationships between nodes such as genetic distances, gene flow or dispersal (Decout, Manel, Miaud, & Luque, 2012) , or spatial/landscape relationships such as geographic distance or landscape resistance (Dyer et al., 2010) . Distinct from other landscape genetic analytical approaches, graphs allow inferences based on the overall shape, or topology, of the network, which can provide unique insights into systemwide processes, such as hierarchical population structure (Dyer & Nason, 2004) . Network topology may be used to identify populations or habitat patches that form important "stepping stones" for maintaining genetic connectivity across an entire system. Such an approach enables experimental simulation whereby nodes may be selectively removed and the overall effect on the system's topology (e.g., overall connectivity, population structure) assessed. Metrics pertaining to the importance of individual nodes to network topology can be correlated with variables such as landscape to identify important drivers of network processes. Despite their unique applications, graph theory and network approaches are relatively underutilized in landscape genetics compared to methods specifically derived from population genetics and landscape ecology. However, among studies of infectious agents, network approaches in wildlife are becoming increasingly popular (Craft, 2015; Craft & Caillaud, 2011) . Epidemiological network models are typically based on host contact networks, which are usually constructed using direct observations or indirect techniques such as mark-recapture, telemetry or proximity loggers, and pathogens are simulated on these contact networks. Such approaches have already incorporated landscape and other environmental features. In addition, the potential for inferring host contacts in network models using pathogen genetic markers (see below) has been acknowledged in recent reviews (Craft, 2015; Gilbertson, Fountain-Jones, & Craft, 2018; White, Forester, & Craft, 2017) , and some studies have directly compared host contact network parameters to parasite genotypes (Bull, Godfrey, & Gordon, 2012) . Despite this, to our knowledge, no published studies have used network models to investigate pathogen movement within a landscape genetic framework. While landscape genetics initially was used to investigate spatial genetic patterns using relatively few neutral markers, the more modern advent of landscape genomics allows the study of variation across the entire genome and effectively expands the scope of landscape genetics to include the study of functional, adaptive genetic variation. Next-generation sequencing (NGS) techniques such as restriction-site-associated DNA sequencing (RADseq) require minimal prior knowledge of the genome under study and can genotype thousands of SNPs randomly distributed across the genome. Some of these SNPs will by chance be located within or near (and thus linked to) genes or regulatory regions that are under selection. Genomewide association studies (GWAS) can make use of this information to identify loci linked to phenotypic variation such as disease susceptibility. Genotyping of candidate loci identified using quantitative trait locus mapping and GWAS can be expanded across a large number of individuals using methods such as targeted sequence capture (Grover, Salmon, & Wendel, 2012) , and these data can be tested in a landscape genomic framework for associations with environmental variables. Loci exhibiting a signature of selection can be identified using outlier tests (Excoffier, Hofer, & Foll, 2009; Luu, Bazin, & Blum, 2017) , which search for loci with allelic frequencies that are outliers relative to the majority. Such loci are considered potentially under selection and may then be tested a posteriori for correlations with environmental variables. Newer methods have focused on explicitly incorporating environmental variables into landscape genomic analyses, known as genetic-environment association (GEA) tests (Lotterhos & Whitlock, 2015; Rellstab, Gugerli, Eckert, Hancock, & Holderegger, 2015) . GEA analyses test for correlations between environmental variables and individual genotypes, which eliminates problems due to underlying population structure that must be controlled when using outlier tests. NGS approaches also generate thousands of neutral loci, which provide greater power to detect fine-scale neutral genetic structure than conventional studies based on relatively few loci (Allendorf, Hohenlohe, & Luikart, 2010) . However, for studies with a particular focus on functional genetic variation, NGS approaches can also be adapted specifically for this purpose through targeted sequencing of the exome (Roffler et al., 2016) (Roffler et al., 2016) . The spread of functional alleles has also been incorporated into landscape genetic simulations , enhancing predictions of future pathogen spread and its effects on host populations. This small body of research is promising for expansion of landscape genomic studies designed to couple pathogen-related functional genetic variation with landscape variables. While we believe that there remains much unexplored utility in established landscape genetic methods for the study of pathogen dynamics as we have described above, we also note new frontiers with significant potential for expanding research in this area. We complete this review by discussing three particularly promising frontiers. Studies relating pathogen genetic data directly to the landscape using resistance surfaces are challenged by the mediating influence of distinct host and vector traits, as well as relative differences in the contributions of multiple host and/or vector species to microparasite gene flow. This necessitates frameworks that more holistically incorporate multiple host and vector factors into studies of pathogen gene flow, which can expand the potential insights provided by landscape genetic studies of infectious agents (Figure 2 ). Single or multiple host or vector species can be added as "landscape variables" (e.g., as resistance surfaces) in addition to physical landscape and environmental variables to test as factors shaping spatial pathogen genetic structure. Resistance surfaces for tests of microparasite gene flow can represent host/vector distributions or abundance, ideally inferred from empirically derived ecological niche or species distribution models. Optimally, host/vector movement would be represented (Dougherty et al., 2018) , using outputs from agent-based movement models informed by telemetry or mark-recapture data, or host/vector landscape genetic data representing spatial patterns of gene flow. We note that the common issue in conventional landscape genetics of spatio-temporal mismatches between landscape processes and genetic change (Anderson et al., 2010; Landguth et al., 2010) would apply even more strongly here. Researchers must simultaneously consider the potentially different spatial and temporal scales over which host and pathogen genetic changes (and poten- Explain observed spatial patterns of prevalence least-cost path models of water bird movement estimated from ecological niche models, and road networks representing human movement, as potential predictors of avian influenza spread (Young et al., 2017 The rapid mutation of microparasites relative to their hosts has potential to provide greater power to detect subtle variation in host movement patterns in response to the landscape, as well as earlier detectability of changes in host movements (such as in response to a new barrier) that are yet to be reflected in host genetic structure (Landguth et al., 2010) . In addition, movements of nonreproducing hosts are difficult to detect using host genetic markers, but instead might be inferred using markers from directly transmitted microorganisms. Such an approach has demonstrated the utility of a chronic, relatively apathogenic infection of felids (feline immunodeficiency virus) for identifying demographic structure of mountain lions and recent population history (Biek et al., 2006) , and has identified movement of bobcats across a highway barrier that was not detectable using host markers (Lee et al., 2012) . However, these approaches have not been broadly applied, particularly in the study of landscape effects. Phylogenetic approaches can reconstruct very recent epidemic histories, providing insights into particular transmission events and pathways that may be contextualized temporally and spatially (Corman et al., 2014; Faria et al., 2014; Carroll et al., 2015; Magee, Beard, Suchard, Lemey, & Scotch, 2015; Fountain-Jones, Packer, et al., 2017; Fountain-Jones, Pearse et al., 2017) . The majority of such work has been conducted on RNA viruses owing to their small, rapidly mutating genomes, requiring relatively little sequencing effort to detect contemporary phylogenetic signals. Other pathogens that evolve more slowly, such as bacteria or fungal pathogens, require the sequencing of larger portions of their genomes to capture equivalent phylogenetic signals (Biek, Pybus, Lloyd-Smith, & Didelot, 2015) . While this is becoming increasingly feasible (Kao, Haydon, Lycett, & Murcia, 2014) , more complex computational analysis is required to make meaningful conclusions. Several approaches may be used for relating phylogenetic information with landscape variables. Neighbour joining trees can identify clusters for quantifying population-level landscape genetic relationships (Joannon et al., 2010) . The calculation of genetic distances based on maximum likelihood trees (Carrel, Emch, Tung, Jobe, & Wan, 2012; Real et al., 2005; Young et al., 2017) results in distance matrices that can be correlated with landscape resistance matrices using conventional landscape genetic approaches. Relaxed random walk phylogeographic approaches (Lemey, Rambaut, Welch, & Suchard, 2010 ) that can reconstruct pathogen dispersal have been linked to landscape predictors using a "phylogeographic GLM" method (Faria, Suchard, Rambaut, Streicker, & Lemey, 2013; Jacquot, Nomikou, Palmarini, Mertens, & Biek, 2017) . The phylogeographic GLM approach has enabled a better understanding of how landscape and hosts can constrain pathogen spread. For example, using the phylogeographic GLM approach on viral genomic data, roads and rivers, coupled with dog distribution, were found to impact rabies spread in Tanzania (Brunker et al., 2018) . However, this approach is limited to discrete sampling locations and is computationally intensive (Dellicour, Rose, & Pybus, 2016) . A recent framework by Dellicour et al. (2016) modifies the phylogeographic GLM approach to use resistance surfaces to efficiently quantify landscape resistance along transmission pathways inferred by continuous phylogeographic analyses. These landscape resistances are then correlated with temporal estimates of transmission along these routes to estimate how the landscape has shaped rates and directions of pathogen spread. Such approaches are yet to be broadly applied, but appear to be important developments that should see increasing application in the future. Overall, landscape genetics has been relatively underutilized in disease ecology research. We believe this is partly due to a lack of cross-disciplinary awareness between the two fields, but also a lack of a clear landscape genetic framework specifically designed for tackling pathogen systems, which are often complex and do not facilitate easy translation of existing landscape genetic tools. However, we note there has been a recent effort to develop new frameworks for such research, expanding the utility of the landscape genetic toolset. These tools will increase our capacity to study complex multihost and host-vector systems, improving the integration of multiple genetic datasets and accounting for interspecific interactions. Improved understanding of host-parasite associations will facilitate the use of microparasite genetic markers to provide insights into host processes that may be difficult to detect using conventional host landscape genetics. Identification of idealized systems that are designed to target specific ecological questions will also facilitate progress in this field. Recent methods that enable the incorporation of quantitative landscape data into spatio-temporal phylogenetic reconstructions of recent transmission events, coupled with advances in high-throughput sequencing, hold great promise for studying how the landscape shapes transmission processes. We believe that these recent developments represent a renewed interest in advancing landscape genetic research in pathogen systems, which we expect will translate to continued growth of research in this area. We thank Louis Bernatchez for inviting us to contribute this paper, and Nicolas Bierne and two anonymous reviewers for their helpful suggestions to improve the manuscript. This work was supported by a National Research Foundation Ecology of Infectious Diseases research programme grant (DEB 1413925) awarded to S.V., S.C., W.C.F., K.R.C., M.E.C. and H.B.E. C.P.K. was supported by an Australian Government Research Training Program Scholarship. The application of "least-cost" modelling as a functional landscape model Genomics and the future of conservation genetics Considering spatial and temporal scale in landscape-genetic studies of gene flow Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology Infecting epidemiology with genetics: A new frontier in disease ecology Landscape genetics: Concepts, methods, applications Current status, future opportunities, and remaining challenges in landscape genetics Identifying future research needs in landscape genetics: Where to from here? Statistical approaches in landscape genetics: An evaluation of methods for linking landscape and genetic data Contrasting landscape epidemiology of two sympatric rabies virus strains A virus reveals population structure and recent demographic history of its carnivore host Measurably evolving pathogens in the genomic era The landscape genetics of infectious disease emergence and spread Landscape genetics and the spatial distribution of chronic wasting disease On applications of landscape genetics Colonization history, host distribution, anthropogenic influence and landscape features shape populations of white pine blister rust, an invasive alien tree pathogen Integrating the landscape epidemiology and genetics of RNA viruses: Rabies in domestic dogs as a model Landscape attributes governing local transmission of an endemic zoonosis: Rabies virus in domestic dogs Social networks and the spread of Salmonella in a sleepy lizard population A host-pathogen simulation model: Powdery mildew of grapevine Population-environment drivers of H5N1 avian influenza molecular change in Vietnam The geography of malaria genetics in the Democratic Republic of Congo: A complex and fragmented landscape Temporal and spatial analysis of the 2014-2015 Ebola virus outbreak in West Africa Landscape genetic structure and evolutionary genetics of insecticide resistance gene mutations in Anopheles sinensis Rooting the phylogenetic tree of Middle East Respiratory Syndrome coronavirus by characterization of a conspecific virus from an African bat Genetic structure and rabies spread potential in raccoons: The role of landscape barriers and sex-biased dispersal Infectious disease transmission and contact networks in wildlife and livestock Network models: An underutilized tool in wildlife epidemiology? Interdisciplinary Perspectives on Infectious Diseases Differential permeability of rivers to raccoon gene flow corresponds to rabies incidence in Ontario Re-evaluating causal modeling with Mantel tests in landscape genetics Toward an integrative molecular approach to wildlife disease Integrative approach for landscape-based graph connectivity analysis: A case study with the common frog (Rana temporaria) in human-dominated landscapes Explaining the geographic spread of emerging epidemics: A framework for comparing viral phylogenies and environmental landscape data Landscape-genetic analysis of population structure in the Texas gray fox oral rabies vaccination zone Going through the motions: Incorporating movement analyses into disease research Nonstationary patterns of isolation-by-distance: inferring measures of local genetic differentiation with Bayesian kriging Is there such a thing as landscape genetics? Molecular Ecology Population graphs: The graph theoretic shape of genetic structure Landscape modelling of gene flow: Improved power using conditional genetic distance derived from the topology of population networks When is connectivity important? A case study of the spatial pattern of sudden oak death Simulation of mathematical models for public health problems Rapid evolutionary response to a transmissible cancer in Tasmanian devils Inferring the ecological niche of Toxoplasma gondii and Bartonella spp. in wild felids Detecting loci under selection in a hierarchically structured population The early spread and epidemic ignition of HIV-1 in human populations Simultaneously reconstructing viral cross-species transmission history and identifying the underlying constraints Urban landscapes can change virus gene flow and evolution in a fragmentationsensitive carnivore Linking social and spatial networks to viral community phylogenetics reveals subtype-specific transmission dynamics in African lions Towards an eco-phylogenetic framework for infectious disease ecology Clustering and assignment methods in landscape genetics Using spatial Bayesian methods to determine the genetic structure of a continuously distributed population: Clusters or isolation by distance? Signatures of environmental genetic adaptation pinpoint pathogens as the main selective pressure through human evolution Applications of graph theory to landscape genetics Fine-scale genetic structure in a wild bird population: The role of limited dispersal and environmentally based selection as causal factors Incorporating genomic methods into contact networks to reveal new insights into animal behaviour and infectious disease dynamics Temporal dynamics of rabies in a wildlife host and the risk of cross-species transmission Unifying the epidemiological and evolutionary dynamics of pathogens Targeted sequence capture as a powerful tool for evolutionary analysis GENELAND: A computer package for landscape genetics Ecology of wildlife rabies in Bluetongue virus spread in Europe is a consequence of climatic, landscape and vertebrate host factors as revealed by phylogeographic inference POPS: Prediction of population genetic structure. France: Program documentation Forecasting changes in population genetic structure of alpine plants in response to global warming Barriers to gene flow between emerging populations of watermelon mosaic virus in southeastern France Why infectious disease research needs community ecology Revealing cryptic spatial patterns in genetic variability by a new multivariate method Supersize me: How whole-genome sequencing and big data are transforming epidemiology Genetic assessment of environmental features that influence deer dispersal: Implications for prion-infected populations Underrepresentation of avian studies in landscape genetics Simulation modeling in landscape genetics Quantifying the lag time to detect barriers in landscape genetics Using landscape genetics simulations for planting blister rust resistant whitebark pine in the US Northern Rocky Mountains Signals of heterogeneous selection at an MHC locus in geographically proximate ecotypes of sockeye salmon Ethical, social, and cultural considerations for site selection for research with genetically modified mosquitoes Targeted enrichment for pathogen detection and characterization in three felid species Gene flow and pathogen transmission among bobcats (Lynx rufus) in a fragmented urban landscape Chapter 3. Invasion of the body snatchers: The diversity and evolution of manipulative strategies in host-parasite interactions Emergence of new virulent populations of apple scab from nonagricultural disease reservoirs Phylogeography takes a relaxed random walk in continuous space and time Multitaxa integrated landscape genetics for zoonotic infectious diseases: Deciphering variables influencing disease emergence The genetic structure of a Venturia inaequalis population in a heterogeneous host population composed of different Malus species Wetlands explain most in the genetic divergence pattern of Oncomelania hupensis Transmission dynamics of co-endemic Plasmodium vivax and P. falciparum in Ethiopia and prevalence of antimalarial resistant genotypes The relative power of genome scans to detect local adaptation depends on sampling design and statistical method pcadapt: An R package to perform genome scans for selection based on principal component analysis Combining phylogeography and spatial epidemiology to uncover predictors of H5N1 influenza A virus diffusion Ten years of landscape genetics Landscape genetics: Combining landscape ecology and population genetics Host dispersal as the driver of parasite genetic structure: A paradigm lost? Using circuit theory to model connectivity in ecology, evolution, and conservation Landscape epidemiology of emerging infectious diseases in natural and humanaltered ecosystems Pathogenmediated selection in free-ranging elk populations infected by chronic wasting disease Maximum-difference barriers: An alternative numerical regionalization method Landscape genetics of high mountain frog metapopulations Spatial epidemiology: An emerging (or re-emerging) discipline. Trends in Ecology and Evolution Landscape genetics of California mule deer (Odocoileus hemionus): The roles of ecological and historical factors in generating differentiation Visualizing spatial population structure with estimated effective migration surfaces Mapping Averaged Pairwise Information (MAPI): A new exploratory tool to uncover spatial structure Demographic processes drive increases in wildlife disease following population reduction Inference of population structure using multilocus genotype data Spatial dynamics and genetics of infectious diseases on heterogeneous landscapes Unifying the spatial population dynamics and molecular evolution of epidemic rabies virus Assessing a landscape barrier using genetic simulation modelling: Implications for raccoon rabies management A practical guide to environmental association analysis in landscape genomics Navigating the pitfalls and promise of landscape genetics Recent range expansion and agricultural landscape heterogeneity have only minimal effect on the spatial genetic structure of the plant pathogenic fungus Mycosphaerella fijiensis Inferences on pathogenic fungus population structures from microsatellite data: New insights from spatial genetics approaches Modelling the dispersal of the two main hosts of the raccoon rabies variant in heterogeneous environments with landscape genetics Using landscape epidemiological models to understand the distribution of chronic wasting disease in the Midwestern USA SNP discovery in candidate adaptive genes using exon capture in a free-ranging alpine ungulate Prediction and prevention of parasitic diseases using a landscape genomics framework Applications of landscape genetics in conservation biology: Concepts and challenges Simulating pattern-process relationships to validate landscape genetic models Oral rabies vaccination in North America: Opportunities, complexities, and challenges Does inbreeding and loss of genetic diversity decrease disease resistance? Conservation Genetics Landscape genetics of the Tasmanian devil: Implications for spread of an infectious cancer Putting the "landscape" in landscape genetics Landscape genetics: Where are we now? Host-pathogen evolutionary signatures reveal dynamics and future invasions of vampire bat rabies Lack of genetic structure and female-specific effect of dispersal barriers in a rabies vector, the striped skunk (Mephitis mephitis) Comparative analysis of landscape effects on spatial genetic structure of the big brown bat and one of its cimicid ectoparasites Integrating microorganism and macroorganism dispersal: Modes, techniques and challenges with particular focus on co-dispersal Beyond MHC: Signals of elevated selection pressure on Atlantic salmon (Salmo salar) immune-relevant loci Exploiting the potential of vector control for disease prevention Juxtaposition between host population structures: Implications for disease transmission in a sympatric cervid community A conceptual framework for the spatial analysis of landscape genetic data The role of parasite-driven selection in shaping landscape genomic structure in red grouse (Lagopus lagopus scotica) Using contact networks to explore mechanisms of parasite transmission in wildlife How's the flu getting through? Landscape genetics suggests both humans and birds spread H5N1 in Egypt Using simulations to evaluate Mantel-based methods for assessing landscape resistance to gene flow Additional supporting information may be found online in the Supporting Information section at the end of the article.