Using multilevel models to identify drivers of landscape‐genetic structure among management areas

Landscape genetics offers a powerful approach to understanding species' dispersal patterns. However, a central obstacle is to account for ecological processes operating at multiple spatial scales, while keeping research outcomes applicable to conservation management. We address this challenge by applying a novel multilevel regression approach to model landscape drivers of genetic structure at both the resolution of individuals and at a spatial resolution relevant to management (i.e. local government management areas: LGAs) for the koala (Phascolartos cinereus) in Australia. Our approach allows for the simultaneous incorporation of drivers of landscape‐genetic relationships operating at multiple spatial resolutions. Using microsatellite data for 1106 koalas, we show that, at the individual resolution, foliage projective cover (FPC) facilitates high gene flow (i.e. low resistance) until it falls below approximately 30%. Out of six additional land‐cover variables, only highways and freeways further explained genetic distance after accounting for the effect of FPC. At the LGA resolution, there was significant variation in isolation‐by‐resistance (IBR) relationships in terms of their slopes and intercepts. This was predominantly explained by the average resistance distance among LGAs, with a weaker effect of historical forest cover. Rates of recent landscape change did not further explain variation in IBR relationships among LGAs. By using a novel multilevel model, we disentangle the effect of landscape resistance on gene flow at the fine resolution (i.e. among individuals) from effects occurring at coarser resolutions (i.e. among LGAs). This has important implications for our ability to identify appropriate scale‐dependent management actions.     

Current approaches using genetic distances produce poor estimates of landscape resistance to interindividual dispersal

Landscape resistance reflects how difficult it is for genes to move across an area with particular attributes (e.g. land cover, slope). An increasingly popular approach to estimate resistance uses Mantel and partial Mantel tests or causal modelling to relate observed genetic distances to effective distances under alternative sets of resistance parameters. Relatively few alternative sets of resistance parameters are tested, leading to relatively poor coverage of the parameter space. Although this approach does not explicitly model key stochastic processes of gene flow, including mating, dispersal, drift and inheritance, bias and precision of the resulting resistance parameters have not been assessed. We formally describe the most commonly used model as a set of equations and provide a formal approach for estimating resistance parameters. Our optimization finds the maximum Mantel r when an optimum exists and identifies the same resistance values as current approaches when the alternatives evaluated are near the optimum. Unfortunately, even where an optimum existed, estimates from the most commonly used model were imprecise and were typically much smaller than the simulated true resistance to dispersal. Causal modelling using Mantel significance tests also typically failed to support the true resistance to dispersal values. For a large range of scenarios, current approaches using a simple correlational model between genetic and effective distances do not yield accurate estimates of resistance to dispersal. We suggest that analysts consider the processes important to gene flow for their study species, model those processes explicitly and evaluate the quality of estimates resulting from their model.     

Urban land use limits regional bumble bee gene flow

Potential declines in native pollinator communities and increased reliance on pollinator‐dependent crops have raised concerns about native pollinator conservation and dispersal across human‐altered landscapes. Bumble bees are one of the most effective native pollinators and are often the first to be extirpated in human‐altered habitats, yet little is known about how bumble bees move across fine spatial scales and what landscapes promote or limit their gene flow. In this study, we examine regional genetic differentiation and fine‐scale relatedness patterns of the yellow‐faced bumble bee, Bombus vosnesenskii, to investigate how current and historic habitat composition impact gene flow. We conducted our study across a landscape mosaic of natural, agricultural and urban/suburban habitats, and we show that B. vosnesenskii exhibits low but significant levels of differentiation across the study system (F_ST = 0.019, D_est = 0.049). Most importantly, we reveal significant relationships between pairwise F_ST and resistance models created from contemporary land use maps. Specifically, B. vosnesenskii gene flow is most limited by commercial, industrial and transportation‐related impervious cover. Finally, our fine‐scale analysis reveals significant but declining relatedness between individuals at the 1–9 km spatial scale, most likely due to local queen dispersal. Overall, our results indicate that B. vosnesenskii exhibits considerable local dispersal and that regional gene flow is significantly limited by impervious cover associated with urbanization.     

Geographical variation in the spatial synchrony of a forest-defoliating insect: isolation of environmental and spatial drivers

Despite the pervasiveness of spatial synchrony of population fluctuations in virtually every taxon, it remains difficult to disentangle its underlying mechanisms, such as environmental perturbations and dispersal. We used multiple regression of distance matrices (MRMs) to statistically partition the importance of several factors potentially synchronizing the dynamics of the gypsy moth, an invasive species in North America, exhibiting outbreaks that are partially synchronized over long distances (approx. 900 km). The factors considered in the MRM were synchrony in weather conditions, spatial proximity and forest-type similarity. We found that the most likely driver of outbreak synchrony is synchronous precipitation. Proximity played no apparent role in influencing outbreak synchrony after accounting for precipitation, suggesting dispersal does not drive outbreak synchrony. Because a previous modelling study indicated weather might indirectly synchronize outbreaks through synchronization of oak masting and generalist predators that feed upon acorns, we also examined the influence of weather and proximity on synchrony of acorn production. As we found for outbreak synchrony, synchrony in oak masting increased with synchrony in precipitation, though it also increased with proximity. We conclude that precipitation could synchronize gypsy moth populations directly, as in a Moran effect, or indirectly, through effects on oak masting, generalist predators or diseases.     

Bacterial farming by the fungus Morchella crassipes

The interactions between bacteria and fungi, the main actors of the soil microbiome, remain poorly studied. Here, we show that the saprotrophic and ectomycorrhizal soil fungus Morchella crassipes acts as a bacterial farmer of Pseudomonas putida, which serves as a model soil bacterium. Farming by M. crassipes consists of bacterial dispersal, bacterial rearing with fungal exudates, as well as harvesting and translocation of bacterial carbon. The different phases were confirmed experimentally using cell counting and ¹³C probing. Common criteria met by other non-human farming systems are also valid for M. crassipes farming, including habitual planting, cultivation and harvesting. Specific traits include delocalization of food production and consumption and separation of roles in the colony (source versus sink areas), which are also found in human agriculture. Our study evidences a hitherto unknown mutualistic association in which bacteria gain through dispersal and rearing, while the fungus gains through the harvesting of an additional carbon source and increased stress resistance of the mycelium. This type of interaction between fungi and bacteria may play a key role in soils.     

稃的有无对大赖草种子扩散与萌发特性的影响

对大赖草[Leymus racemosus（Lam．）Tzvel．]带稃种子和去稃种子的扩散性、萌发率、活力及吸水性和失水性的差异进行了比较。结果表明：在静止空气中带稃种子的降落速度显著小于去稃种子,分别为3．96和4．81m·s-1;而在风速1和4m·s-1的水平气流中带稃种子的扩散距离（15．26和11．86cm）均显著大于去稃种子（27．80和21．60cm）。培养20d后,带稃种子和去稃种子的萌发率分别为99％和97％,差异不显著;在自然条件及60℃高温条件下,带稃种子和去稃种子的活力无显著差异。带稃种子的吸水饱和时间和吸水量均显著大于去稃种子,而其失水速率小于后者但二者差异不显著。研究结果显示：稃对大赖草种子萌发和活力均无显著影响,但可增加种子的风媒扩散能力、有效保持种子含水量并能降低种子失水速率,对大赖草适应干旱的沙漠环境具有重要作用。     

Turkey as a crossroad for Greater Flamingos Phoenicopterus roseus: evidence from population trends and ring-resightings (Aves: Phoenicopteridae)

The Greater Flamingo Phoenicopterus roseus is a waterbird commonly found in saline and brackish lagoons throughout the Mediterranean Region. We have gathered existing data on Greater Flamingos in Turkey and carried out field surveys to present the most up to date information on wintering (1999–2014) and breeding (1969– 2014). The wintering population of flamingos shows an increasing trend with 54,947±20,794 individuals mainly concentrated in the Gediz, Büyük Menderes and Çukurova deltas, respectively. Breeding attempts were recorded in at least seven wetlands in Turkey in the past, yet after 1999 most of the colonies were abandoned due to basin scale intensive water management practices in Central Anatolia. Currently, only Tuz Lake and Gediz Delta are used as regular breeding sites, while breeding has been recorded sporadically in Ac?göl and Ak?ehir Lakes. The breeding colony of Tuz Lake is of prime importance at the Mediterranean scale, with the number of young chicks in 2011, 2012 and 2013 accounting for the highest number of fledglings in the Mediterranean Region and West Africa (18,418, 20,274 and 20,292 respectively). Finally, building upon the previous findings about Turkey and the western Mediterranean metapopulation links, recent resightings of Turkish flamingos (despite the limited numbers) confirm post-fledging and natal dispersal reaching the western Mediterranean Basin and West Africa. Flamingos from Turkey were also found to disperse to Israel and to a region outside the known flyways of the western Mediterranean and West African flamingos (i.e. to Israel and UAE). Thus, Turkey, due to its geographic position, appears to be a crossroad between the western and eastern Mediterranean Region and southwest Asia.     

Role of evolutionary and ecological factors in the reproductive success and the spatial genetic structure of the temperate gorgonian Paramuricea clavata

Dispersal and mating features strongly influence the evolutionary dynamics and the spatial genetic structure (SGS) of marine populations. For the first time in a marine invertebrate, we examined individual reproductive success, by conducting larval paternity assignments after a natural spawning event, combined with a small‐scale SGS analysis within a population of the gorgonian Paramuricea clavata. Thirty four percent of the larvae were sired by male colonies surrounding the brooding female colonies, revealing that the bulk of the mating was accomplished by males from outside the studied area. Male success increased with male height and decreased with increasing male to female distance. The parentage analyses, with a strong level of self‐recruitment (25%), unveiled the occurrence of a complex family structure at a small spatial scale, consistent with the limited larval dispersal of this species. However, no evidence of small scale SGS was revealed despite this family structure. Furthermore, temporal genetic structure was not observed, which appears to be related to the rather large effective population size. The low level of inbreeding found suggests a pattern of random mating in this species, which disagrees with expectations that limited larval dispersal should lead to biparental inbreeding. Surface brooding and investment in sexual reproduction in P. clavata contribute to multiple paternity (on average 6.4 fathers were assigned per brood), which enhance genetic diversity of the brood. Several factors may have contributed to the lack of biparental inbreeding in our study such as (i) the lack of sperm limitation at a small scale, (ii) multiple paternity, and (iii) the large effective population size. Thus, our results indicate that limited larval dispersal and complex family structure do not necessarily lead to biparental inbreeding and SGS. In the framework of conservation purposes, our results suggested that colony size, proximity among colonies and the population size should be taken into consideration for restoration projects.     

The magnitude of local adaptation under genotype‐dependent dispersal

Dispersal moves individuals from patches where their immediate ancestors were successful to sites where their genotypes are untested. As a result, dispersal generally reduces fitness, a phenomenon known as “migration load.” The strength of migration load depends on the pattern of dispersal and can be dramatically lessened or reversed when individuals move preferentially toward patches conferring higher fitness. Evolutionary ecologists have long modeled nonrandom dispersal, focusing primarily on its effects on population density over space, the maintenance of genetic variation, and reproductive isolation. Here, we build upon previous work by calculating how the extent of local adaptation and the migration load are affected when individuals differ in their dispersal rate in a genotype‐dependent manner that alters their match to their environment. Examining a one‐locus, two‐patch model, we show that local adaptation occurs through a combination of natural selection and adaptive dispersal. For a substantial portion of parameter space, adaptive dispersal can be the predominant force generating local adaptation. Furthermore, genetic load may be largely averted with adaptive dispersal whenever individuals move before selective deaths occur. Thus, to understand the mechanisms driving local adaptation, biologists must account for the extent and nature of nonrandom, genotype‐dependent dispersal, and the potential for adaptation via spatial sorting of genotypes.     

Limited contemporary gene flow and high self‐replenishment drives peripheral isolation in an endemic coral reef fish

Extensive ongoing degradation of coral reef habitats worldwide has lead to declines in abundance of coral reef fishes and local extinction of some species. Those most vulnerable are ecological specialists and endemic species. Determining connectivity between locations is vital to understanding recovery and long‐term persistence of these species following local extinction. This study explored population connectivity in the ecologically‐specialized endemic three‐striped butterflyfish (Chaetodon tricinctus) using mt and msatDNA (nuclear microsatellites) to distinguish evolutionary versus contemporary gene flow, estimate self‐replenishment and measure genetic diversity among locations at the remote Australian offshore coral reefs of Middleton Reef (MR), Elizabeth Reef (ER), Lord Howe Island (LHI), and Norfolk Island (NI). Mt and msatDNA suggested genetic differentiation of the most peripheral location (NI) from the remaining three locations (MR, ER, LHI). Despite high levels of mtDNA gene flow, there is limited msatDNA gene flow with evidence of high levels of self‐replenishment (≥76%) at all four locations. Taken together, this suggests prolonged population recovery times following population declines. The peripheral population (NI) is most vulnerable to local extinction due to its relative isolation, extreme levels of self‐replenishment (95%), and low contemporary abundance.