Two are better than one: combining landscape genomics and common gardens for detecting local adaptation in forest trees

Predicting likely species responses to an alteration of their local environment is key to decision‐making in resource management, ecosystem restoration and biodiversity conservation practice in the face of global human‐induced habitat disturbance. This is especially true for forest trees which are a dominant life form on Earth and play a central role in supporting diverse communities and structuring a wide range of ecosystems. In Europe, it is expected that most forest tree species will not be able to migrate North fast enough to follow the estimated temperature isocline shift given current predictions for rapid climate warming. In this context, a topical question for forest genetics research is to quantify the ability for tree species to adapt locally to strongly altered environmental conditions (Kremer et al. 2012 ). Identifying environmental factors driving local adaptation is, however, a major challenge for evolutionary biology and ecology in general but is particularly difficult in trees given their large individual and population size and long generation time. Empirical evaluation of local adaptation in trees has traditionally relied on fastidious long‐term common garden experiments (provenance trials) now supplemented by reference genome sequence analysis for a handful of economically valuable species. However, such resources have been lacking for most tree species despite their ecological importance in supporting whole ecosystems. In this issue of Molecular Ecology, De Kort et al. ( 2014 ) provide original and convincing empirical evidence of local adaptation to temperature in black alder, Alnus glutinosa L. Gaertn, a surprisingly understudied keystone species supporting riparian ecosystems. Here, De Kort et al. ( 2014 ) use an innovative empirical approach complementing state‐of‐the‐art landscape genomics analysis of A. glutinosa populations sampled in natura across a regional climate gradient with phenotypic trait assessment in a common garden experiment (Fig. 1 ). By combining the two methods, De Kort et al. ( 2014 ) were able to detect unequivocal association between temperature and phenotypic traits such as leaf size as well as with genetic loci putatively under divergent selection for temperature. The research by De Kort et al. ( 2014 ) provides valuable insight into adaptive response to temperature variation for an ecologically important species and demonstrates the usefulness of an integrated approach for empirical evaluation of local adaptation in nonmodel species (Sork et al. 2013 ).     

Pattern,prediction and parsimony in continental‐scale synthesis of pyromes: a reply to Gosper et al.

We (Murphy et al., 2013; Clarke et al., 2015) have recently developed a framework to understand the spatial distribution of fire regimes and plant fire‐response traits at large spatial scales. We integrated a range of data sources to create a continental‐scale overview of Australian pyromes from which to infer pyrogeographic drivers. Gosper et al. (in press) have criticized our approach, based on our misclassification of a vegetation type (eucalypt woodland), with distinct fire regime, in the Coolgardie bioregion of Western Australia. We argue that the intention of our integrative approach was to develop and refine conceptual models of Australian pyrogeography, not to produce a predictive map of fire regimes, and certainly not to guide local‐scale fire management. Like all models, continental‐scale syntheses of pyromes are imperfect, yet they still represent powerful tools for understanding the drivers of the spatial distribution of fire regimes.     

Tropical soils degraded by slash‐and‐burn cultivation can be recultivated when amended with ashes and compost

In many tropical regions, slash‐and‐burn agriculture is considered as a driver of deforestation; the forest is converted into agricultural land by cutting and burning the trees. However, the fields are abandoned after few years because of yield decrease and weed invasion. Consequently, new surfaces are regularly cleared from the primary forest. We propose a reclamation strategy for abandoned fields allowing and sustaining re‐cultivation. In the dry region of south‐western Madagascar, we tested, according to a split‐plot design, an alternative selective slash‐and‐burn cultivation technique coupled with compost amendment on 30–year‐old abandoned fields. Corn plants (Zea mays L.) were grown on four different types of soil amendments: no amendment (control), compost, ashes (as in traditional slash‐and‐burn cultivation), and compost + ashes additions. Furthermore, two tree cover treatments were applied: 0% tree cover (as in traditional slash‐and‐burn cultivation) and 50% tree cover (selective slash‐and‐burn). Both corn growth and soil fertility parameters were monitored during the growing season 2015 up to final harvest. The amendment compost + ashes strongly increased corn yield, which was multiplied by 4–5 in comparison with ashes or compost alone, reaching 1.5 t/ha compared to 0.25 and 0.35 t/ha for ashes and compost, respectively. On control plots, yield was negligible as expected on these degraded soils. Structural equation modeling evidenced that compost and ashes were complementary fertilizing pathways promoting soil fertility through positive effects on soil moisture, pH, organic matter, and microbial activity. Concerning the tree cover treatment, yield was reduced on shaded plots (50% tree cover) compared to sunny plots (0% tree cover) for all soil amendments, except ashes. To conclude, our results provide empirical evidence on the potential of recultivating tropical degraded soils with compost and ashes. This would help mitigating deforestation of the primary forest by increasing lifespan of agricultural lands.     

Challenges and prospects for interpreting long‐term phytoplankton diversity changes in Lake Zurich (Switzerland)

相似文献   

Demystifying the RAD fad

We are writing in response to the population and phylogenomics meeting review by Andrews & Luikart ( 2014 ) entitled ‘Recent novel approaches for population genomics data analysis’. Restriction‐site‐associated DNA (RAD) sequencing has become a powerful and useful approach in molecular ecology, with several different published methods now available to molecular ecologists, none of which can be considered the best option in all situations. A&L report that the original RAD protocol of Miller et al. ( 2007 ) and Baird et al. ( 2008 ) is superior to all other RAD variants because putative PCR duplicates can be identified (see Baxter et al. 2011 ), thereby reducing the impact of PCR artefacts on allele frequency estimates (Andrews & Luikart 2014 ). In response, we (i) challenge the assertion that the original RAD protocol minimizes the impact of PCR artefacts relative to that of other RAD protocols, (ii) present additional biases in RADseq that are at least as important as PCR artefacts in selecting a RAD protocol and (iii) highlight the strengths and weaknesses of four different approaches to RADseq which are a representative sample of all RAD variants: the original RAD protocol (mbRAD, Miller et al. 2007 ; Baird et al. 2008 ), double digest RAD (ddRAD, Peterson et al. 2012 ), ezRAD (Toonen et al. 2013 ) and 2bRAD (Wang et al. 2012 ). With an understanding of the strengths and weaknesses of different RAD protocols, researchers can make a more informed decision when selecting a RAD protocol.     

The discovery of Antarctic RNA viruses: a new game changer

Antarctic ecosystems are dominated by micro‐organisms, and viruses play particularly important roles in the food webs. Since the first report in 2009 (López‐Bueno et al. 2009 ), ‘omic’‐based studies have greatly enlightened our understanding of Antarctic aquatic microbial diversity and ecosystem function (Wilkins et al. 2013 ; Cavicchioli 2015 ). This has included the discovery of many new eukaryotic viruses (López‐Bueno et al. 2009 ), virophage predators of algal viruses (Yau et al. 2011 ), bacteria with resistance to phage (Lauro et al. 2011 ) and mechanisms of haloarchaeal evasion, defence and adaptation to viruses (Tschitschko et al. 2015 ). In this issue of Molecular Ecology, López‐Bueno et al. ( 2015 ) report the first discovery of RNA viruses from an Antarctic aquatic environment. High sequence coverage enabled genome variation to be assessed for four positive‐sense single‐stranded RNA viruses from the order Picornavirales. By examining the populations present in the water column and in the lake's catchment area, populations of ‘quasispecies’ were able to be linked to local environmental factors. In view of the importance of viruses in Antarctic ecosystems but lack of data describing them, this study represents a significant advance in the field.     

Ecologically flexible endemics dominate Indo‐Pacific bird communities

Reeve et al. (2016, Ecography, 39 , 990‐997) found that ecologically flexible endemics dominate Indo‐Pacific bird communities. This negative relationship between local abundance and global range size contrasts strongly with the positive range size‐abundance relationship “rule,” which would predict community dominance by globally widespread species. Theuerkauf et al. (2017, Journal of Biogeography, 44 , 2161–2163) provide new data from New Caledonia which they claim invalidate our study. They find positive relationships between local abundance and local range size, which they attribute to endemic species having narrower habitat niches than globally widespread species. We reanalysed their data using global range sizes, corroborating the pattern we originally reported: negative relationships between local abundance and global range size, driven by a subset of adaptable endemic species. We stress the importance of being explicit about the scale of ecological mechanisms, and ensuring that the scale of analysis matches the scale of interpretation.     

Mock communities highlight the diversity of host‐associated eukaryotes

Host‐associated microbes are ubiquitous. Every multicellular eukaryote, and even many unicellular eukaryotes (protists), hosts a diverse community of microbes. High‐throughput sequencing (HTS) tools have illuminated the vast diversity of host‐associated microbes and shown that they have widespread influence on host biology, ecology and evolution (McFall‐Ngai et al. 2013 ). Bacteria receive most of the attention, but protists are also important components of microbial communities associated with humans (Parfrey et al. 2011 ) and other hosts. As HTS tools are increasingly used to study eukaryotes, the presence of numerous and diverse host‐associated eukaryotes is emerging as a common theme across ecosystems. Indeed, HTS studies demonstrate that host‐associated lineages account for between 2 and 12% of overall eukaryotic sequences detected in soil, marine and freshwater data sets, with much higher relative abundances observed in some samples (Ramirez et al. 2014 ; Simon et al. 2015 ; de Vargas et al. 2015 ). Previous studies in soil detected large numbers of predominantly parasitic lineages such as Apicomplexa, but did not delve into their origin [e.g. (Ramirez et al. 2014 )]. In this issue of Molecular Ecology, Geisen et al. ( 2015 ) use mock communities to show that many of the eukaryotic organisms detected by environmental sequencing in soils are potentially associated with animal hosts rather than free‐living. By isolating the host‐associated fraction of soil microbial communities, Geisen and colleagues help explain the surprisingly high diversity of parasitic eukaryotic lineages often detected in soil/terrestrial studies using high‐throughput sequencing (HTS) and reinforce the ubiquity of these host‐associated microbes. It is clear that we can no longer assume that organisms detected in bulk environmental sequencing are free‐living, but instead need to design studies that specifically enumerate the diversity and function of host‐associated eukaryotes. Doing so will allow the field to determine the role host‐associated eukaryotes play in soils and other environments and to evaluate hypotheses on assembly of host‐associated communities, disease ecology and more.     

Accelerating forest loss in Southeast Asian Massif in the 21st century: A case study in Nan Province,Thailand

Farmers are carving a new agricultural frontier from the forests in the Southeast Asian Massif (SAM) in the 21st century, triggering significant environment degradation at the local scale; however, this frontier has been missed by existing global land use and forest loss analyses. In this paper, we chose Thailand's Nan Province, which is located in the geometric center of SAM, as a case study, and combined high resolution forest cover change product with a fine‐scale land cover map to investigate land use dynamics with respect to topography in this region. Our results show that total forest loss in Nan Province during 2001–2016 was 66,072 ha (9.1% of the forest cover in 2000), and that the majority of this lost forest (92%) had been converted into crop (mainly corn) fields by 2017. Annual forest loss is significantly correlated with global corn price (p < 0.01), re‐confirming agricultural expansion as a key driver of forest loss in Nan Province. Along with the increasing global corn price, forest loss in Nan Province has accelerated at a rate of 2,616 ± 730 ha per decade (p < 0.01). Global corn price peaked in 2012, in which year annual forest loss also reached its peak (7,523 ha); since then, the location of forest loss has moved to steeper land at higher elevations. Spatially, forest loss driven by this smallholder agricultural expansion emerges as many small patches that are not recognizable even at a moderate spatial resolution (e.g. 300 m). It explains how existing global land use/cover change products have missed the widespread and rapid forest loss in SAM. It also highlights the importance of high‐resolution observations to evaluate the environmental impacts of agricultural expansion and forest loss in SAM, including, but not limited to, the impacts on the global carbon cycle, regional hydrology, and local environmental degradation.     

Positive range–abundance relationships in Indo‐Pacific bird communities

Reeve et al. (2016, Ecography, 39 , 990–997) recently reported negative range–abundance relationships in Indo‐Pacific bird communities and speculated that geographical isolation facilitates the evolution of broad‐niched, small‐ranged and abundant species. We tested this relationship using a large independent data set on range and abundance of birds across New Caledonia (over 4,000 bird census points for 17,300 km²). In contradiction to Reeve et al. (2016, Ecography, 39 , 990–997), we found clear evidence that range–abundance relationships are positive and endemic species have narrower habitat niches than wide‐range species. Our findings are likely valid also for other islands in the Indo‐Pacific.     

Pollen‐based biomes for Beringia 18,000, 6000 and 0 14C yr bp†

The objective biomization method developed by Prentice et al. (1996) for Europe was extended using modern pollen samples from Beringia and then applied to fossil pollen data to reconstruct palaeovegetation patterns at 6000 and 18,000 ¹⁴C yr bp . The predicted modern distribution of tundra, taiga and cool conifer forests in Alaska and north‐western Canada generally corresponds well to actual vegetation patterns, although sites in regions characterized today by a mosaic of forest and tundra vegetation tend to be preferentially assigned to tundra. Siberian larch forests are delimited less well, probably due to the extreme under‐representation of Larix in pollen spectra. The biome distribution across Beringia at 6000 ¹⁴C yr bp was broadly similar to today, with little change in the northern forest limit, except for a possible northward advance in the Mackenzie delta region. The western forest limit in Alaska was probably east of its modern position. At 18,000 ¹⁴C yr bp the whole of Beringia was covered by tundra. However, the importance of the various plant functional types varied from site to site, supporting the idea that the vegetation cover was a mosaic of different tundra types.     

Host plant richness explains diversity of ectomycorrhizal fungi: Response to the comment of Tedersoo et al. (2014)

Exploring the relationships between the biodiversity of groups of interacting organisms yields insight into ecosystem stability and function (Hooper et al. 2000 ; Wardle 2006 ). We demonstrated positive relationships between host plant richness and ectomycorrhizal (EM) fungal diversity both in a field study in subtropical China (Gutianshan) and in a meta‐analysis of temperate and tropical studies (Gao et al. 2013 ). However, based on re‐evaluation of our data sets, Tedersoo et al. ( 2014 ) argue that the observed positive correlation between EM fungal richness and EM plant richness at Gutianshan and also in our metastudies was based mainly from (i) a sampling design with inconsistent species pool and (ii) poor data compilation for the meta‐analysis. Accordingly, we checked our data sets and repeated the analysis performed by Tedersoo et al. ( 2014 ). In contrast to Tedersoo et al. ( 2014 ), our re‐analysis still confirms a positive effect of plant richness on EM fungal diversity in Gutianshan, temperate and tropical ecosystems, respectively.     

Will reduced host connectivity curb the spread of a devastating epidemic?

The white‐nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans, is threatening the cave‐dwelling bat fauna of North America by killing individuals by the thousands in hibernacula each winter since its appearance in New York State less than ten years ago. Epidemiological models predict that WNS will reach the western coast of the USA by 2035, potentially eliminating most populations of susceptible bat species in its path (Frick et al. 2015; O'Regan et al. 2015). These models were built and validated using distributional data from the early years of the epidemic, which spread throughout eastern North America following a route driven by cave density and winter severity (Maher et al. 2012). In this issue of Molecular Ecology, Wilder et al. (2015) refine these findings by showing that connectivity among host populations, as assessed by population genetic markers, is crucial in determining the spread of the pathogen. Because host connectivity is much reduced in the hitherto disease free western half of North America, Wilder et al. make the reassuring prediction that the disease will spread more slowly west of the Great Plains.     

Forest floor vegetation response to nitrogen deposition in Europe

Chronic nitrogen (N) deposition is a threat to biodiversity that results from the eutrophication of ecosystems. We studied long‐term monitoring data from 28 forest sites with a total of 1,335 permanent forest floor vegetation plots from northern Fennoscandia to southern Italy to analyse temporal trends in vascular plant species cover and diversity. We found that the cover of plant species which prefer nutrient‐poor soils (oligotrophic species) decreased the more the measured N deposition exceeded the empirical critical load (CL) for eutrophication effects (P = 0.002). Although species preferring nutrient‐rich sites (eutrophic species) did not experience a significantly increase in cover (P = 0.440), in comparison to oligotrophic species they had a marginally higher proportion among new occurring species (P = 0.091). The observed gradual replacement of oligotrophic species by eutrophic species as a response to N deposition seems to be a general pattern, as it was consistent on the European scale. Contrary to species cover changes, neither the decrease in species richness nor of homogeneity correlated with nitrogen CL exceedance (ExCL_empN). We assume that the lack of diversity changes resulted from the restricted time period of our observations. Although existing habitat‐specific empirical CL still hold some uncertainty, we exemplify that they are useful indicators for the sensitivity of forest floor vegetation to N deposition.     

Using RAD‐seq to recognize sex‐specific markers and sex chromosome systems

Next‐generation sequencing methods have initiated a revolution in molecular ecology and evolution (Tautz et al. 2010 ). Among the most impressive of these sequencing innovations is restriction site‐associated DNA sequencing or RAD‐seq (Baird et al. 2008 ; Andrews et al. 2016 ). RAD‐seq uses the Illumina sequencing platform to sequence fragments of DNA cut by a specific restriction enzyme and can generate tens of thousands of molecular genetic markers for analysis. One of the many uses of RAD‐seq data has been to identify sex‐specific genetic markers, markers found in one sex but not the other (Baxter et al. 2011 ; Gamble & Zarkower 2014 ). Sex‐specific markers are a powerful tool for biologists. At their most basic, they can be used to identify the sex of an individual via PCR. This is useful in cases where a species lacks obvious sexual dimorphism at some or all life history stages. For example, such tests have been important for studying sex differences in life history (Sheldon 1998 ; Mossman & Waser 1999 ), the management and breeding of endangered species (Taberlet et al. 1993 ; Griffiths & Tiwari 1995 ; Robertson et al. 2006 ) and sexing embryonic material (Hacker et al. 1995 ; Smith et al. 1999 ). Furthermore, sex‐specific markers allow recognition of the sex chromosome system in cases where standard cytogenetic methods fail (Charlesworth & Mank 2010 ; Gamble & Zarkower 2014 ). Thus, species with male‐specific markers have male heterogamety (XY) while species with female‐specific markers have female heterogamety (ZW). In this issue, Fowler & Buonaccorsi ( 2016 ) illustrate the ease by which RAD‐seq data can generate sex‐specific genetic markers in rockfish (Sebastes). Moreover, by examining RAD‐seq data from two closely related rockfish species, Sebastes chrysomelas and Sebastes carnatus (Fig.  1 ), Fowler & Buonaccorsi ( 2016 ) uncover shared sex‐specific markers and a conserved sex chromosome system.     

Phylogenetic structural equation modelling reveals no need for an ‘origin’ of the leaf economics spectrum

The leaf economics spectrum (LES) is a prominent ecophysiological paradigm that describes global variation in leaf physiology across plant ecological strategies using a handful of key traits. Nearly a decade ago, Shipley et al. (2006) used structural equation modelling to explore the causal functional relationships among LES traits that give rise to their strong global covariation. They concluded that an unmeasured trait drives LES covariation, sparking efforts to identify the latent physiological trait underlying the ‘origin’ of the LES. Here, we use newly developed phylogenetic structural equation modelling approaches to reassess these conclusions using both global LES data as well as data collected across scales in the genus Helianthus. For global LES data, accounting for phylogenetic non‐independence indicates that no additional unmeasured traits are required to explain LES covariation. Across datasets in Helianthus, trait relationships are highly variable, indicating that global‐scale models may poorly describe LES covariation at non‐global scales.     

Why some parasites are widespread and abundant while others are local and rare?

Abundances and distributions of species are usually associated. This implies that as a species declines in abundance so does the number of sites it occupies. Conversely, when there is an increase in a species' range size, it is usually followed by an increase in population size (Gaston et al. 2000 ). This ecological phenomenon, also known as the abundance–occupancy relationship (AOR), is well documented in several species of animals and plants (Gaston et al. 2000 ) but has been little investigated in parasites. In this issue of Molecular Ecology, Drovetski et al. ( 2014 ) investigated the AOR in avian haemosporidians (vector‐borne blood parasites) using data from four well‐sampled bird communities. In support of the AOR, the research group found that the abundance of parasite cytochrome b lineages (a commonly used proxy for species identification within this group of parasites) was positively linked with the abundance of susceptible avian host species and that the most abundant haemospordian lineages were those with larger ranges. Drovetski et al. ( 2014 ) also found evidence for both hypotheses proposed to explain the AOR in parasites: the trade‐off hypothesis (TOH) and the niche‐breadth hypothesis (NBH). Interestingly, the main predictor of the AOR was the number of susceptible hosts (i.e. number of infected birds) and not the number of host species the parasites were able to exploit.     

High‐throughput SNPs for all: genotyping‐in‐thousands

Understanding the genetic structure of species is essential for conservation. It is only with this information that managers, academics, user groups and land‐use planners can understand the spatial scale of migration and local adaptation, source‐sink dynamics and effective population size. Such information is essential for a multitude of applications including delineating management units, balancing management priorities, discovering cryptic species and implementing captive breeding programmes. Species can range from locally adapted by hundreds of metres (Pavey et al. 2010 ) to complete species panmixia (Côté et al. 2013 ). Even more remarkable is that this essential information can be obtained without fully sequenced or annotated genomes, but from mere (putatively) nonfunctional variants. First with allozymes, then microsatellites and now SNPs, this neutral genetic variation carries a wealth of information about migration and drift. For many of us, it may be somewhat difficult to remember our understanding of species conservation before the widespread usage of these useful tools. However most species on earth have yet to give us that ‘peek under the curtain’. With the current diversity on earth estimated to be nearly 9 million species (Mora et al. 2011 ), we have a long way to go for a comprehensive meta‐phylogeographic understanding. A method presented in this issue by Campbell and colleagues (Campbell et al. 2015 ) is a tool that will accelerate the pace in this area. Genotyping‐in‐thousands (GT‐seq) leverages recent advancements in sequencing technology to save many hours and dollars over previous methods to generate this important neutral genetic information.     

Ensemble modeling to predict habitat suitability for a large‐scale disturbance specialist

To conserve habitat for disturbance specialist species, ecologists must identify where individuals will likely settle in newly disturbed areas. Habitat suitability models can predict which sites at new disturbances will most likely attract specialists. Without validation data from newly disturbed areas, however, the best approach for maximizing predictive accuracy can be unclear (Northwestern U.S.A.). We predicted habitat suitability for nesting Black‐backed Woodpeckers (Picoides arcticus; a burned‐forest specialist) at 20 recently (≤6 years postwildfire) burned locations in Montana using models calibrated with data from three locations in Washington, Oregon, and Idaho. We developed 8 models using three techniques (weighted logistic regression, Maxent, and Mahalanobis D² models) and various combinations of four environmental variables describing burn severity, the north–south orientation of topographic slope, and prefire canopy cover. After translating model predictions into binary classifications (0 = low suitability to unsuitable, 1 = high to moderate suitability), we compiled “ensemble predictions,” consisting of the number of models (0–8) predicting any given site as highly suitable. The suitability status for 40% of the area burned by eastside Montana wildfires was consistent across models and therefore robust to uncertainty in the relative accuracy of particular models and in alternative ecological hypotheses they described. Ensemble predictions exhibited two desirable properties: (1) a positive relationship with apparent rates of nest occurrence at calibration locations and (2) declining model agreement outside surveyed environments consistent with our reduced confidence in novel (i.e., “no‐analogue”) environments. Areas of disagreement among models suggested where future surveys could help validate and refine models for an improved understanding of Black‐backed Woodpecker nesting habitat relationships. Ensemble predictions presented here can help guide managers attempting to balance salvage logging with habitat conservation in burned‐forest landscapes where black‐backed woodpecker nest location data are not immediately available. Ensemble modeling represents a promising tool for guiding conservation of large‐scale disturbance specialists.