Spatially structured statistical network models for landscape genetics

Multi-generational dispersal and dynamic patch occupancy reveals spatial and temporal stability of seascapes

The success of non-native species (NNS) invasions depends on patterns of dispersal and connectivity, which underpin genetic diversity, population establishment and growth. In the marine environment, both global environmental change and increasing anthropogenic activity can alter hydrodynamic patterns, leading to significant inter-annual variability in dispersal pathways. Despite this, multi-generational dispersal is rarely explicitly considered in attempts to understand NNS spread or in the design of management interventions. Here, we present a novel approach to quantifying species spread that considers range expansion and network formation across time using the non-native Pacific oyster, Magallana gigas (Thunberg 1793), as a model. We combined biophysical modelling, dynamic patch occupancy models, consideration of environmental factors, and graph network theory to model multi-generational dispersal in northwest Europe over 13 generations. Results revealed that M. gigas has a capacity for rapid range expansion through the creation of an ecological network of dispersal pathways that remains stable through time. Maximum network size was achieved in four generations, after which connectivity patterns remained temporally stable. Multi-generational connectivity could therefore be divided into two periods: network growth (2000-2003) and network stability (2004-2012). Our study is the first to examine how dispersal trajectories affect the temporal stability of ecological networks across biogeographic scales, and provides an approach for the assignment of site-based prioritisation of non-native species management at different stages of the invasion timeline. More broadly, the framework we present can be applied to other fields (e.g. Marine Protected Area design, management of threatened species and species range expansion due to climate change) as a means of characterising and defining ecological network structure, functioning and stability.

Why More Biologists Must Embrace Quantitative Modeling

Biology as a field has transformed since the time of its foundation from an organized enterprise cataloging the diversity of the natural world to a quantitatively rigorous science seeking to answer complex questions about the functions of organisms and their interactions with each other and their environments. As the mathematical rigor of biological analyses has improved, quantitative models have been developed to describe multi-mechanistic systems and to test complex hypotheses. However, applications of quantitative models have been uneven across fields, and many biologists lack the foundational training necessary to apply them in their research or to interpret their results to inform biological problem-solving efforts. This gap in scientific training has created a false dichotomy of "biologists" and "modelers" that only exacerbates the barriers to working biologists seeking additional training in quantitative modeling. Here, we make the argument that all biologists are modelers and are capable of using sophisticated quantitative modeling in their work. We highlight four benefits of conducting biological research within the framework of quantitative models, identify the potential producers and consumers of information produced by such models, and make recommendations for strategies to overcome barriers to their widespread implementation. Improved understanding of quantitative modeling could guide the producers of biological information to better apply biological measurements through analyses that evaluate mechanisms, and allow consumers of biological information to better judge the quality and applications of the information they receive. As our explanations of biological phenomena increase in complexity, so too must we embrace modeling as a foundational skill.

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease

Wildlife diseases are a major global threat to biodiversity. Boreal toads (Anaxyrus [Bufo] boreas) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity (π = 0.00034-0.00040) and effective population sizes (Ne = 8.9-38.4) were low, likely due to post-Pleistocene founder effects and Bd-related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd.

Fencing Can Alter Gene Flow of Asian Elephant Populations within Protected Areas

The Asian elephant is mostly confined to mountainous ranges and therefore risks population fragmentation if hard protected area (PA) boundaries near steep slopes prevent movement. We tested whether elephant gene flow is (i) controlled by slope and (ii) affected by the interplay between barriers and slope. We used 176 unique genotypes obtained non-invasively from fresh elephant dung to assess individual-by-individual genetic distance across the Western Ghats of India, a biodiversity hotspot. To assess landscape distance, 36 resistance models were produced by transforming a slope raster. Core areas and corridors were calculated from the raster that provided the best correlation between the genetic and distance matrices. The influence of the closure of PAs on gene flow was examined for one region, the Nilgiri Biosphere Reserve. The best resistance raster obtained by transforming the slope occupancy model was better than Euclidean distance for explaining genetic distance, indicating that slope partially controls gene flow. Fencing elephant PAs on hilly terrain reduces core areas and disrupts corridors. Consequently, hard PA boundaries abutting slopes can fragment elephant populations, but this can be ameliorated by protecting the adjacent flatter terrain.

Scale-dependent influence of the sagebrush community on genetic connectivity of the sagebrush obligate Gunnison sage-grouse

Habitat fragmentation and degradation impacts an organism's ability to navigate the landscape, ultimately resulting in decreased gene flow and increased extinction risk. Understanding how landscape composition impacts gene flow (i.e., connectivity) and interacts with scale is essential to conservation decision‐making. We used a landscape genetics approach implementing a recently developed statistical model based on the generalized Wishart probability distribution to identify the primary landscape features affecting gene flow and estimate the degree to which each component influences connectivity for Gunnison sage‐grouse (Centrocercus minimus). We were interested in two spatial scales: among distinct populations rangewide and among leks (i.e., breeding grounds) within the largest population, Gunnison Basin. Populations and leks are nested within a landscape fragmented by rough terrain and anthropogenic features, although requisite sagebrush habitat is more contiguous within populations. Our best fit models for each scale confirm the importance of sagebrush habitat in connectivity, although the important sagebrush characteristics differ. For Gunnison Basin, taller shrubs and higher quality nesting habitat were the primary drivers of connectivity, while more sagebrush cover and less conifer cover facilitated connectivity rangewide. Our findings support previous assumptions that Gunnison sage‐grouse range contraction is largely the result of habitat loss and degradation. Importantly, we report direct estimates of resistance for landscape components that can be used to create resistance surfaces for prioritization of specific locations for conservation or management (i.e., habitat preservation, restoration, or development) or as we demonstrated, can be combined with simulation techniques to predict impacts to connectivity from potential management actions.

Bibliometric Analysis of Global Research on Ecological Networks in Nature Conservation from 1990 to 2020

As a nature-based solution to land-use sustainability, ecological networks (ENs) have received substantial attention from researchers, planners, and decision-makers worldwide. To portray the global research on ENs in nature conservation during the period of 1990–2020, 1371 papers in 53 subject categories were reviewed with bibliometric methods and CiteSpace. The results showed a successive growth of publications at an annually averaged rate of 18.9% during the past three decades. Co-citation analysis indicated that the most popular topic was connectivity, on which the studies concentrated on quantifying connectivity, identifying priority areas, and integrating conservation planning. A recent hotspot is to study the landscape fragmentation effects on natural habitats or biodiversity under land-use changes in urbanized areas. Multidisciplinary approaches have been increasingly used to tackle more complex interplays among economic, social, ecological, and cultural factors, with the aim of alleviating ecological service losses attributed to human activities. Spatiotemporal dynamics and participatory design of ENs at different scales have become an emerging trend. In order to address increasing pressures on biodiversity or landscape connectivity brought about by land use and climate change, it is suggested to develop more research on the evaluation and management of the resilience of ENs.

What processes must we understand to forecast regional-scale population dynamics?

An urgent challenge facing biologists is predicting the regional-scale population dynamics of species facing environmental change. Biologists suggest that we must move beyond predictions based on phenomenological models and instead base predictions on underlying processes. For example, population biologists, evolutionary biologists, community ecologists and ecophysiologists all argue that the respective processes they study are essential. Must our models include processes from all of these fields? We argue that answering this critical question is ultimately an empirical exercise requiring a substantial amount of data that have not been integrated for any system to date. To motivate and facilitate the necessary data collection and integration, we first review the potential importance of each mechanism for skilful prediction. We then develop a conceptual framework based on reaction norms, and propose a hierarchical Bayesian statistical framework to integrate processes affecting reaction norms at different scales. The ambitious research programme we advocate is rapidly becoming feasible due to novel collaborations, datasets and analytical tools.

The use and misuse of regression models in landscape genetic analyses

The field of landscape genetics has been rapidly evolving, adopting and adapting analytical frameworks to address research questions. Current studies are increasingly using regression-based frameworks to infer the individual contributions of landscape and habitat variables on genetic differentiation. This paper outlines appropriate and inappropriate uses of multiple regression for these purposes, and demonstrates through simulation the limitations of different analytical frameworks for making correct inference. Of particular concern are recent studies seeking to explain genetic differences by fitting regression models with effective distance variables calculated independently on separate landscape resistance surfaces. When moving across the landscape, organisms cannot respond independently and uniquely to habitat and landscape features. Analyses seeking to understand how landscape features affect gene flow should model a single conductance or resistance surface as a parameterized function of relevant spatial covariates, and estimate the values of these parameters by linking a single set of resistance distances to observed genetic dissimilarity via a loss function. While this loss function may involve a regression-like step, the associated nuisance parameters are not interpretable in terms of organismal movement and should not be conflated with what is actually of interest: the mapping between spatial covariates and conductance/resistance. The growth and evolution of landscape genetics as a field has been rapid and exciting. It is the goal of this paper to highlight past missteps and demonstrate limitations of current approaches to ensure that future use of regression models will appropriately consider the process being modeled, which will provide clarity to model interpretation.

A novel quantitative framework for riverscape genetics

Riverscape genetics, which applies concepts in landscape genetics to riverine ecosystems, lack appropriate quantitative methods that address the spatial autocorrelation structure of linear stream networks and account for bidirectional geneflow. To address these challenges, we present a general framework for the design and analysis of riverscape genetic studies. Our framework starts with the estimation of pairwise genetic distance at sample sites and the development of a spatially structured ecological network (SSEN) on which riverscape covariates are measured. We then introduce the novel bidirectional geneflow in riverscapes (BGR) model that uses principles of isolation-by-resistance to quantify the effects of environmental covariates on genetic connectivity, with spatial covariance defined using simultaneous autoregressive models on the SSEN and the generalized Wishart distribution to model pairwise distance matrices arising through a random walk model of geneflow. We highlight the utility of this framework in an analysis of riverscape genetics for brook trout (Salvelinus fontinalis) in north central Pennsylvania, USA. Using the fixation index (F_ST ) as the measure of genetic distance, we estimated the effects of 12 riverscape covariates on geneflow by evaluating the relative support of eight competing BGR models. We then compared the performance of the top-ranked BGR model to results obtained from comparable analyses using multiple regression on distance matrices (MRM) and the program STRUCTURE. We found that the BGR model had more power to detect covariate effects, particularly for variables that were only partial barriers to geneflow and/or uncommon in the riverscape, making it more informative for assessing patterns of population connectivity and identifying threats to species conservation. This case study highlights the utility of our modeling framework over other quantitative methods in riverscape genetics, particularly the ability to rigorously test hypotheses about factors that influence geneflow and probabilistically estimate the effect of riverscape covariates, including stream flow direction. This framework is flexible across taxa and riverine networks, is easily executable, and provides intuitive results that can be used to investigate the likely outcomes of current and future management scenarios.