BUSCADOR RECOLECTA

Vector transmission from DSF grapevines data set Data from vector transmission experiment for Xylella fastidiosa from DSF transgenic grapevines. File is an .rds file for loading directly into R. The file is a list of length 2. The first item in the list is metadata describing each variable of the data set, including total sample size. The second item is the data set itself. zeilinger_dsf_transmission_dataset.rds, Effective management of vector-borne plant pathogens often relies on disease-resistant cultivars. While heterogeneity in host resistance and in pathogen population density at the host population level play important and well-recognized roles in epidemiology, the effects of resistance traits on pathogen distribution at the individual host level, and the epidemiological consequences in turn, are poorly understood. Transgenic disease-resistant plants that produce bacterial Diffusible Signaling Factor (DSF) could provide resistance to the vector-borne bacterium Xylella fastidiosa by impeding plant colonization and reducing virulence. However, the effects of constitutive in planta production of DSF on insect vector transmission has remained unresolved. We investigated the transmission biology of X. fastidiosa in DSF and wild-type (WT) grapevines with the efficient vector Graphocephala atropunctata. We also developed a novel Bayesian hierarchical model to improve statistical inference on the multiple components of the vector transmission process. We found that insect vectors had a greater colonization efficiency on DSF plants—meaning they acquired a greater population size of X. fastidiosa—than on WT plants. However, DSF plants also maintained much lower X. fastidiosa populations. These apparently conflicting processes resulted in a lower but highly variable probability of transmission from DSF plants compared to WT plants. Our Bayesian model improved statistical inference compared to widely used frequentist statistics in part by estimating and correcting for imperfect detection of X. fastidiosa in plant and insect tissues. Overall, our results support current models on the roles that DSF plays in vector transmission of X. fastidiosa. In line with our hypothesis, DSF production reduced mean X. fastidiosa population density but increased heterogeneity within host plants. While DSF-producing plants could potentially improve disease management, our results suggest that they could, under some conditions, facilitate X. fastidiosa spread., Peer reviewed

data for dryads Raw data used for comparative analyses on eye colour in owls, Birds, due to their multiple colourful displays, constitute a classic paradigm for the study of colour evolution. Although avian eyes are remarkably coloured, the functional basis behind inter‐specific variability in iris colouration remains poorly understood. Owls are an ideal system to shed light on the role of ecology in promoting iris colour evolution as they show inter‐specific variation in iris colour and in niche specialization with some species being strictly nocturnal and others active during the day. Owls perching for hunting at night might be unnoticed by both predators and their prey if they had dark irises, which would predict that dark irises were more likely to evolve in strictly nocturnal species than in diurnal ones. Using phylogenetic comparative models, we tested the camouflage hypothesis for eye colour. Ancestral state reconstruction revealed that the owl ancestor of the family Strigidae was more likely bright‐irided whereas the ancestor of the family Tytonidae was more likely dark‐irided. We found that iris colour and activity rhythm have more likely evolved in concert than independently, and a non‐significant trend of dark eyes to evolve more easily in owl species presenting strictly nocturnal habits than in diurnal species. The transition from diurnality to nocturnality was a previous requisite for the evolution of dark irises in owls. Taken together our results are only partly consistent with the camouflage hypothesis suggesting that dark irises in owls have primarily evolved to enhance concealment in nocturnal conditions., Peer reviewed

genolike_map20_q20_snp16_maf002_ind32_430K.beagle Input file of genotype likelihoods in Beagle format (.beagle.gz) for PCANGSD analysis. ORJUSTRUx8_biall_dp450_q40_perpopmiss075_hwe00001_maf002_LD02_neutral01 Dataset of 16,858 independent neutral SNPs from 64 samples in vcf format for STRUCTURE analysis. ORJU06x12_biall_dp450_q40_maf002_nomissing_hwe00001_noZ_neutral01 Dataset of 15,209 neutral SNPs used in the PCA intended for population structure correction in the redundancy analysis (RDA). ORJU06x12_biall_dp450_q40_maf002_nomissing_hwe00001_noZ Dataset of 15,252 SNPs used in the redundancy analysis. ORJU06x12_biall_dp450_q40_maf002_nomissing_hwe00001_sel01_BAYESCENV Subset of 49 SNP outliers identified with BayScEnv for the second redundancy analysis in the manuscript. ORJUGEO_coords_ALLSAMPLES Table of sample location and geographic coordinates for all sequenced samples of the study. ORJUGEO_coords_ClimVariables Sample data and WorldClim environmental data (Fick, S.E. and R.J. Hijmans, 2017. Worldclim 2: New 1-km spatial resolution climate surfaces for global land areas. International Journal of Climatology.) for each one of the samples included in the genotype-environment association analyses., The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (i) drift in geographic isolation, (ii) differentiation along continuous selective gradients, and (iii) isolation by adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci., Peer reviewed

Raw data on intact seeds in shorebirds This file contains data of the faecal samples collected from the field in different locations with date, faeces sample mass, plant species name, and the number of seeds per taxon. This file was made with Microsoft Excel 2016 shorebird.xlsx, Shorebirds (Charadriiformes) undergo rapid migrations with potential for long-distance dispersal (LDD) of plants. We studied the frequency of endozoochory by shorebirds in different parts of Europe covering a broad latitudinal range and different seasons. We assessed whether plants dispersed conformed to morphological dispersal syndromes. A total of 409 excreta samples (271 faeces and 138 pellets) were collected from redshank (Tringa totanus), black-winged stilt (Himantopus himantopus), pied avocet (Recurvirostra avosetta), northern lapwing (Vanellus vanellus), Eurasian curlew (Numenius arquata) and black-tailed godwit (Limosa limosa) in south-west Spain, north-west England, southern Ireland and Iceland in 2005 and 2016, and intact seeds were extracted and identified. Godwits were sampled just before or after migratory movements between England and Iceland. The germinability of seeds was tested. Intact diaspores were recovered from all bird species and study areas, and were present in 13% of samples overall. Thirteen plant families were represented, including Charophyceae and 26 angiosperm taxa. Only four species had an "endozoochory syndrome". Four alien species were recorded. Ellenberg values classified three species as aquatic and 20 as terrestrial. Overall, 89% of seeds were from terrestrial plants, and 11% from aquatic plants. Average seed length was higher in redshank pellets than in their faeces. Six species were germinated, none of which had an endozoochory syndrome. Seeds were recorded during spring and autumn migration. Plant species recorded have broad latitudinal ranges consistent with LDD via shorebirds. Crucially, morphological syndromes do not adequately predict LDD potential, and more empirical work is required to identify which plants are dispersed by shorebirds. Incorporating endozoochory by shorebirds and other migratory waterbirds into plant distribution models would allow us to better understand the natural processes that facilitated colonization of oceanic islands, or to improve predictions of how plants will respond to climate change, or how alien species spread., Peer reviewed

Adult wood frog traits Jumping performance, adult size, and other morphometrics AdultWoodFrogTraits2014.xls, Human-modified habitats rarely yield outcomes that are aligned with conservation ideals. Landscapes that are subdivided by roads are no exception, precipitating negative impacts on populations due to fragmentation, pollution, and road kill. Although many populations in human modified habitats show evidence for local adaptation, rarely does environmental change yield outright benefits for populations of conservation interest. Contrary to expectations, we report surprising benefits experienced by amphibian populations breeding and dwelling in proximity to roads. We show that roadside populations of the wood frog, Rana sylvatica, exhibit better locomotor performance and higher measures of traits related to fitness compared with frogs from less disturbed environments located further away from roads. These results contrast previous evidence for maladaptation in roadside populations of wood frogs studied elsewhere. Our results indicate that altered habitats might not be unequivocally detrimental, and at times might contribute to metapopulation success. While the frequency of such beneficial outcomes remains unknown, their occurrence underscores the complexity of inferring consequences of environmental change., National Science Foundation, Award: DEB-1655092, Peer reviewed

ddRAD reads Quality-filtered reads for each sample in FASTA format. ddRAD_reads.zip Genotype_counts.pl Perl script for calculating genotype counts from a PED file to be used in the VCF2LR program. Genotype_counts.zip Assortativity.pl Perl script for calculating the assortativity coefficient between two groups or sectors and example file. Assortativity.zip SNP dataset SNP dataset in PED and VCF formats. SNP_dataset.zip, Analyzing the impact of anthropogenic and natural river barriers on the dispersal of aquatic and semi-aquatic species may be critical for their conservation, but no adequate genetic methods have been developed for quantifying the effect of specific barriers on current connectivity. Knowledge of kinship relationships between individuals and reconstructions of pedigrees obtained using genomic data can be extremely useful, not only for studying the social organization of animals, but also inferring how the last few generations of offspring have dispersed. In this study, we used kinship data to analyze connectivity patterns in a small semi-aquatic mammal, the Pyrenean desman, in an area comprising two river systems with close headwaters and dams of various sizes. Using a large SNP dataset from 70 specimens, we obtained kinship categories and reconstructed pedigrees. To quantify the barrier effect of specific obstacles, we constructed kinship networks and devised a method based on the assortativity coefficient, which measures the proportion between observed and expected kinship relationships across a barrier. The estimation of this parameter enabled us to infer that the most important barrier in the area was the watershed divide between the rivers, followed by a dam on one of the rivers. Other barriers did not significantly reduce the expected number of kinship relationships across them. This strategy and the information obtained with it may be crucial in determining the most important connectivity problems in an area and help develop conservation plans aimed at improving genetic exchange between populations of threatened species., Peer reviewed

Arabidopsis_data_279accessions_20jan16, Genetic diversity provides insight into heterogeneous demographic and adaptive history across organisms’ distribution ranges. For this reason, decomposing single species into genetic units may represent a powerful tool to better understand biogeographical patterns as well as improve predictions of the effects of GCC (global climate change) on biodiversity loss. Using 279 georeferenced Iberian accessions, we used classes of three intraspecific genetic units of the annual plant Arabidopsis thaliana obtained from the genetic analyses of nuclear SNPs (single nucleotide polymorphisms), chloroplast SNPs, and the vernalization requirement for flowering. We used SDM (species distribution models), including climate, vegetation, and soil data, at the whole-species and genetic-unit levels. We compared model outputs for present environmental conditions and with a particularly severe GCC scenario. SDM accuracy was high for genetic units with smaller distribution ranges. Kernel density plots identified the environmental variables underpinning potential distribution ranges of genetic units. Combinations of environmental variables accounted for potential distribution ranges of genetic units, which shrank dramatically with GCC at almost all levels. Only two genetic clusters increased their potential distribution ranges with GCC. The application of SDM to intraspecific genetic units provides a detailed picture on the biogeographical patterns of distinct genetic groups based on different genetic criteria. Our approach also allowed us to pinpoint the genetic changes, in terms of genetic background and physiological requirements for flowering, that Iberian A. thaliana may experience with a GCC scenario applying SDM to intraspecific genetic units., Peer reviewed

AFLP profiles for the plants of Lavandula latifolia sampled for the study Markers are identified by primer combination and size (base pairs) AFLP.scores.xlsx MSAP profiles for the plants of Lavandula latifolia sampled for the study Markers are identified by primer combination and size (base pairs). MSAP.scores.xlsx Spatial coordinates of Lavandula latifolia plants sampled for the study Coordinates are expressed in meters and the origin (0,0) is the lower left corner of the plot as depicted in the figures. Plant.coordinates.xlsx, Little is known on the potential of ecological disturbance to cause genetic and epigenetic changes in plant populations. We take advantage of a long-term field experiment initiated in 1986 to study the demography of the shrub Lavandula latifolia, and compare genetic and epigenetic characteristics of plants in two adjacent subplots, one experimentally disturbed and one left undisturbed, 20 years after disturbance. Experimental setup was comparable to an unreplicated ‘Before-After-Control-Impact’ (BACI) design where a single pair of perturbed and control areas were compared. When sampled in 2005, plants in the two subplots had roughly similar ages, but they had established in contrasting environments: dense conspecific population (‘Undisturbed’ subpopulation) versus open area with all conspecifics removed (‘Disturbed’ subpopulation). Plants were characterized genetically and epigenetically using amplified fragment length polymorphism (AFLP) and two classes of methylation-sensitive AFLP (MSAP) markers. Subpopulations were similar in genetic diversity but differed in epigenetic diversity and multilocus genetic and epigenetic characteristics. Epigenetic divergence between subpopulations was statistically unrelated to genetic divergence. Bayesian clustering revealed an abrupt linear boundary between subpopulations closely coincident with the arbitrary demarcation line between subplots drawn 20 years back, which supports that genetic and epigenetic divergence between subpopulations was caused by artificial disturbance. There was significant fine-scale spatial structuring of MSAP markers in both subpopulations, which in the Undisturbed one was indistinguishable from that of AFLP markers. Genetic differences between subpopulations could be explained by divergent selection alone, while the concerted action of divergent selection and disturbance-driven appearance of new methylation variants in the Disturbed subpopulation is proposed to explain epigenetic differences. This study provides the first empirical evidence to date suggesting that relatively mild disturbances could leave genetic and epigenetic signatures on the next adult generation of long-lived plants., Peer reviewed

13C values of amino acids Sample 13C values (‰) of amino acids. Each sample was analyzed in triplicate (mean ± SD). NA indicates missing values or replicates. The values are compiled from this study and Larsen et al. (Plos One 2013) (marked with asterisk under sample code). Larsen_13CAA.csv, Supplementation of nutrients by symbionts enables consumers to thrive on resources that might otherwise be insufficient to meet nutritional demands. Such nutritional subsidies by intracellular symbionts has been well studied; however, supplementation of de novo synthesized nutrients to hosts by extracellular gut symbionts is poorly documented, especially for generalists with relatively undifferentiated intestinal tracts. Although gut symbionts facilitate degradation of resources that would otherwise remain inaccessible to the host, such digestive actions alone cannot make up for dietary insufficiencies of macronutrients such as essential amino acids (EAA). Documenting whether gut symbionts also function as partners for symbiotic EAA supplementation is important because the question of how some detritivores are able to subsist on nutritionally insufficient diets has remained unresolved. To answer this poorly-understood nutritional aspect of symbiont-host interactions, we studied the enchytraeid worm, a bulk soil feeder that thrives in arctic peatlands. In a combined field and laboratory study, we employed stable isotope fingerprinting of amino acids to identify the biosynthetic origins of amino acids to bacteria, fungi and plants in enchytraeids. Enchytraeids collected from arctic peatlands derived more than 80% of their EAA from bacteria. In a controlled feeding study with the enchytraeid Enchytraeus crypticus, EAA derived almost exclusively from gut bacteria when the worms fed on higher fiber diets, whereas most of the enchytraeids' EAA derived from dietary sources when fed on lower fiber diets. Our gene sequencing results of gut microbiota showed that the worms harbor several taxa in their gut lumen absent from their diets and substrates. Almost all gut taxa are candidates for EAA supplementation because almost all belong to clades capable of biosynthesizing EAA. Our study provides the first evidence of extensive symbiotic supplementation of EAA by microbial gut symbionts, and demonstrate that symbiotic bacteria in the gut lumen appear to function as partners for both symbiotic EAA supplementation as well as for digestion of insoluble plant fibers., Peer reviewed

De novo assembly of Rhabdomys pumilio transcriptome and detailed differential gene expression results This dataset contains a fasta file, including the de novo assembly of Rhabdomys pumilio transcriptome (built from RNA-Seq data using the Trinity suite of de novo transcriptome assembly tools - https://github.com/trinityrnaseq/trinityrnaseq/wiki), as well as a tab separated text file indicating the gene annotations corresponding to each assembly contig, obtained by mapping the de novo assembly to the human reference transcriptome, as described in the manuscript. A third Excel file contains the detailed results of RNA-Seq differential gene expression study in skin obtained from three dorsal regions, as described in the manuscript. Rhabdomys_pumilio.zip, Mammalian colour patterns are among the most recognizable characteristics found in nature and can have a profound impact on fitness. However, little is known about the mechanisms underlying the formation and subsequent evolution of these patterns. Here we show that, in the African striped mouse (Rhabdomys pumilio), periodic dorsal stripes result from underlying differences in melanocyte maturation, which give rise to spatial variation in hair colour. We identify the transcription factor ALX3 as a regulator of this process. In embryonic dorsal skin, patterned expression of Alx3 precedes pigment stripes and acts to directly repress Mitf, a master regulator of melanocyte differentiation, thereby giving rise to light-coloured hair. Moreover, Alx3 is upregulated in the light stripes of chipmunks, which have independently evolved a similar dorsal pattern. Our results show a previously undescribed mechanism for modulating spatial variation in hair colour and provide insights into how phenotypic novelty evolves., Peer reviewed