Showing papers in "Molecular Ecology in 2014"
TL;DR: Testing the model of differential gene flow among loci by asking whether absolute divergence is also higher in the previously identified ‘islands’ finds that absolute measures of divergence are not higher in genomic islands, and simulations using the program IMa2 suggest that inferences of any gene flow may be incorrect in many comparisons.
Abstract: The metaphor of ‘genomic islands of speciation’ was first used to describe heterogeneous differentiation among loci between the genomes of closely related species. The biological model proposed to explain these differences was that the regions showing high levels of differentiation were resistant to gene flow between species, while the remainder of the genome was being homogenized by gene flow and consequently showed lower levels of differentiation. However, the conditions under which such differentiation can occur at multiple unlinked loci are restrictive; additionally, essentially, all previous analyses have been carried out using relative measures of divergence, which can be misleading when regions with different levels of recombination are compared. Here, we test the model of differential gene flow by asking whether absolute divergence is also higher in the previously identified ‘islands’. Using five species pairs for which full sequence data are available, we find that absolute measures of divergence are not higher in genomic islands. Instead, in all cases examined, we find reduced diversity in these regions, a consequence of which is that relative measures of divergence are abnormally high. These data therefore do not support a model of differential gene flow among loci, although islands of relative divergence may represent loci involved in local adaptation. Simulations using the program IMa2 further suggest that inferences of any gene flow may be incorrect in many comparisons. We instead present an alternative explanation for heterogeneous patterns of differentiation, one in which postspeciation selection generates patterns consistent with multiple aspects of the data.
948 citations
TL;DR: The underlying ecological processes that can generate patterns of IBE are described, its implications for a wide variety of disciplines are examined and several areas of future research are outlined that can answer pressing questions about the ecological basis of genetic diversity.
Abstract: The interactions between organisms and their environments can shape distributions of spatial genetic variation, resulting in patterns of isolation by environment (IBE) in which genetic and environmental distances are positively correlated, independent of geographic distance. IBE represents one of the most important patterns that results from the ways in which landscape heterogeneity influences gene flow and population connectivity, but it has only recently been examined in studies of ecological and landscape genetics. Nevertheless, the study of IBE presents valuable opportunities to investigate how spatial heterogeneity in ecological processes, agents of selection and environmental variables contributes to genetic divergence in nature. New and increasingly sophisticated studies of IBE in natural systems are poised to make significant contributions to our understanding of the role of ecology in genetic divergence and of modes of differentiation both within and between species. Here, we describe the underlying ecological processes that can generate patterns of IBE, examine its implications for a wide variety of disciplines and outline several areas of future research that can answer pressing questions about the ecological basis of genetic diversity.
619 citations
TL;DR: It is concluded that in species that exhibit IBD or have undergone range expansion, many of the published FST outliers based on FDIST2 and BayeScan are probably false positives, but FLK and Bayenv2 show great promise for accurately identifying loci under spatially divergent selection.
Abstract: FST outlier tests are a potentially powerful way to detect genetic loci under spatially divergent selection. Unfortunately, the extent to which these tests are robust to nonequilibrium demographic histories has been understudied. We developed a landscape genetics simulator to test the effects of isolation by distance (IBD) and range expansion on FST outlier methods. We evaluated the two most commonly used methods for the identification of FST outliers (FDIST2 and BayeScan, which assume samples are evolutionarily independent) and two recent methods (FLK and Bayenv2, which estimate and account for evolutionary nonindependence). Parameterization with a set of neutral loci (‘neutral parameterization’) always improved the performance of FLK and Bayenv2, while neutral parameterization caused FDIST2 to actually perform worse in the cases of IBD or range expansion. BayeScan was improved when the prior odds on neutrality was increased, regardless of the true odds in the data. On their best performance, however, the widely used methods had high false-positive rates for IBD and range expansion and were outperformed by methods that accounted for evolutionary nonindependence. In addition, default settings in FDIST2 and BayeScan resulted in many false positives suggesting balancing selection. However, all methods did very well if a large set of neutral loci is available to create empirical P-values. We conclude that in species that exhibit IBD or have undergone range expansion, many of the published FST outliers based on FDIST2 and BayeScan are probably false positives, but FLK and Bayenv2 show great promise for accurately identifying loci under spatially divergent selection.
493 citations
TL;DR: It is suggested that three mosquito species spanning the breadth of the Culicidae depend on their gut microbiome for development, and pyrosequencing of 16S rRNA provides new insights about the acquisition and structure of bacterial communities in mosquitoes.
Abstract: Field studies indicate adult mosquitoes (Culicidae) host low diversity communities of bacteria that vary greatly among individuals and species. In contrast, it remains unclear how adult mosquitoes acquire their microbiome, what influences community structure, and whether the microbiome is important for survival. Here, we used pyrosequencing of 16S rRNA to characterize the bacterial communities of three mosquito species reared under identical conditions. Two of these species, Aedes aegypti and Anopheles gambiae, are anautogenous and must blood-feed to produce eggs, while one, Georgecraigius atropalpus, is autogenous and produces eggs without blood feeding. Each mosquito species contained a low diversity community comprised primarily of aerobic bacteria acquired from the aquatic habitat in which larvae developed. Our results suggested that the communities in Ae. aegypti and An. gambiae larvae share more similarities with one another than with G. atropalpus. Studies with Ae. aegypti also strongly suggested that adults transstadially acquired several members of the larval bacterial community, but only four genera of bacteria present in blood fed females were detected on eggs. Functional assays showed that axenic larvae of each species failed to develop beyond the first instar. Experiments with Ae. aegypti indicated several members of the microbial community and Escherichia coli successfully colonized axenic larvae and rescued development. Overall, our results provide new insights about the acquisition and structure of bacterial communities in mosquitoes. They also indicate that three mosquito species spanning the breadth of the Culicidae depend on their gut microbiome for development.
386 citations
TL;DR: It is proposed that the potential functions of microbial symbionts in facilitating or restricting the use of host plants are constrained by their location (intracellular, gut or environmental), and by the fidelity of their associations with insect host lineages.
Abstract: Herbivory, defined as feeding on live plant tissues, is characteristic of highly successful and diverse groups of insects and represents an evolutionarily derived mode of feeding. Plants present various nutritional and defensive barriers against herbivory; nevertheless, insects have evolved a diverse array of mechanisms that enable them to feed and develop on live plant tissues. For decades, it has been suggested that insect-associated microbes may facilitate host plant use, and new molecular methodologies offer the possibility to elucidate such roles. Based on genomic data, specialized feeding on phloem and xylem sap is highly dependent on nutrient provisioning by intracellular symbionts, as exemplified by Buchnera in aphids, although it is unclear whether such symbionts play a substantive role in host plant specificity of their hosts. Microorganisms present in the gut or outside the insect body could provide more functions including digestion of plant polymers and detoxification of plant-produced toxins. However, the extent of such contributions to insect herbivory remains unclear. We propose that the potential functions of microbial symbionts in facilitating or restricting the use of host plants are constrained by their location (intracellular, gut or environmental), and by the fidelity of their associations with insect host lineages. Studies in the next decade, using molecular methods from environmental microbiology and genomics, will provide a more comprehensive picture of the role of microbial symbionts in insect herbivory.
354 citations
TL;DR: For advancing the field of polyploid population genetics, most priority should be given to development of new molecular approaches that allow efficient dosage determination, and to further development of analytical approaches to circumvent dosage uncertainty and to accommodate ‘flexible’ modes of inheritance.
Abstract: Despite the importance of polyploidy and the increasing availability of new genomic data, there remain important gaps in our knowledge of polyploid population genetics. These gaps arise from the complex nature of polyploid data (e.g. multiple alleles and loci, mixed inheritance patterns, association between ploidy and mating system variation). Furthermore, many of the standard tools for population genetics that have been developed for diploids are often not feasible for polyploids. This review aims to provide an overview of the state-of-the-art in polyploid population genetics and to identify the main areas where further development of molecular techniques and statistical theory is required. We review commonly used molecular tools (amplified fragment length polymorphism, microsatellites, Sanger sequencing, next-generation sequencing and derived technologies) and their challenges associated with their use in polyploid populations: that is, allele dosage determination, null alleles, difficulty of distinguishing orthologues from paralogues and copy number variation. In addition, we review the approaches that have been used for population genetic analysis in polyploids and their specific problems. These problems are in most cases directly associated with dosage uncertainty and the problem of inferring allele frequencies and assumptions regarding inheritance. This leads us to conclude that for advancing the field of polyploid population genetics, most priority should be given to development of new molecular approaches that allow efficient dosage determination, and to further development of analytical approaches to circumvent dosage uncertainty and to accommodate ‘flexible’ modes of inheritance. In addition, there is a need for more simulation-based studies that test what kinds of biases could result from both existing and novel approaches.
300 citations
TL;DR: High‐throughput 16S rRNA gene sequencing is used to examine the host and environmental influences on the skin microbiota of the cohabiting amphibian species Anaxyrus boreas, Pseudacris regilla, Taricha torosa and Lithobates catesbeianus from the Central Valley in California, and populations of Rana cascadae tadpoles.
Abstract: Skin-associated bacteria of amphibians are increasingly recognized for their role in defence against pathogens, yet we have little understanding of their basic ecology. Here, we use high-throughput 16S rRNA gene sequencing to examine the host and environmental influences on the skin microbiota of the cohabiting amphibian species Anaxyrus boreas, Pseudacris regilla, Taricha torosa and Lithobates catesbeianus from the Central Valley in California. We also studied populations of Rana cascadae over a large geographic range in the Klamath Mountain range of Northern California, and across developmental stages within a single site. Dominant bacterial phylotypes on amphibian skin included taxa from Bacteroidetes, Gammaproteobacteria, Alphaproteobacteria, Firmicutes, Sphingobacteria and Actinobacteria. Amphibian species identity was the strongest predictor of microbial community composition. Secondarily, within a given amphibian species, wetland site explained significant variation. Amphibian-associated microbiota differed systematically from microbial assemblages in their environments. Rana cascadae tadpoles have skin bacterial communities distinct from postmetamorphic conspecifics, indicating a strong developmental shift in the skin microbes following metamorphosis. Establishing patterns observed in the skin microbiota of wild amphibians and environmental factors that underlie them is necessary to understand skin symbiont community assembly, and ultimately, the role skin microbiota play in the extended host phenotype including disease resistance.
288 citations
TL;DR: It is found that it is possible to greatly reduce error rates by considering the results of all three methods when identifying outlier loci, and the relative ranking between the methods is impacted by the consideration of polygenic selection.
Abstract: The recent availability of next-generation sequencing (NGS) has made possible the use of dense genetic markers to identify regions of the genome that may be under the influence of selection. Several statistical methods have been developed recently for this purpose. Here, we present the results of an individual-based simulation study investigating the power and error rate of popular or recent genome scan methods: linear regression, Bayescan, BayEnv and LFMM. Contrary to previous studies, we focus on complex, hierarchical population structure and on polygenic selection. Additionally, we use a false discovery rate (FDR)-based framework, which provides an unified testing framework across frequentist and Bayesian methods. Finally, we investigate the influence of population allele frequencies versus individual genotype data specification for LFMM and the linear regression. The relative ranking between the methods is impacted by the consideration of polygenic selection, compared to a monogenic scenario. For strongly hierarchical scenarios with confounding effects between demography and environmental variables, the power of the methods can be very low. Except for one scenario, Bayescan exhibited moderate power and error rate. BayEnv performance was good under nonhierarchical scenarios, while LFMM provided the best compromise between power and error rate across scenarios. We found that it is possible to greatly reduce error rates by considering the results of all three methods when identifying outlier loci.
288 citations
TL;DR: The different forms of evolutionary rate heterogeneity are described and explained how they can be accommodated in molecular‐clock analyses, and an outline of the various clock methods and models that are available is provided.
Abstract: The molecular clock presents a means of estimating evolutionary rates and timescales using genetic data. These estimates can lead to important insights into evolutionary processes and mechanisms, as well as providing a framework for further biological analyses. To deal with rate variation among genes and among lineages, a diverse range of molecular-clock methods have been developed. These methods have been implemented in various software packages and differ in their statistical properties, ability to handle different models of rate variation, capacity to incorporate various forms of calibrating information and tractability for analysing large data sets. Choosing a suitable molecular-clock model can be a challenging exercise, but a number of model-selection techniques are available. In this review, we describe the different forms of evolutionary rate heterogeneity and explain how they can be accommodated in molecular-clock analyses. We provide an outline of the various clock methods and models that are available, including the strict clock, local clocks, discrete clocks and relaxed clocks. Techniques for calibration and clock-model selection are also described, along with methods for handling multilocus data sets. We conclude our review with some comments about the future of molecular clocks.
267 citations
TL;DR: Evidence suggests that overharvest drives the decay of genetic diversity across a wide range of marine fishes, which may lead to a long-term impact of fishing on their evolutionary potential, particularly if abundance remains low and diversity continues to decay.
Abstract: While population declines can drive the loss of genetic diversity under some circumstances, it has been unclear whether this loss is a general consequence of overharvest in highly abundant marine fishes We compiled data from 11 049 loci across 140 species and found that allelic richness was lower in overfished populations within 9 of 12 genera and families A multiple linear regression showed that allelic richness was on average 12% lower (P < 00001) in overharvested populations after accounting for the effects of body size, latitude and other factors Heterozygosity was on average 2% lower (P = 0030) Simulations confirmed that these patterns are consistent with a recent bottleneck in abundant species and also showed that our analysis likely underestimates the loss of rare alleles by a factor of two or three This evidence suggests that overharvest drives the decay of genetic diversity across a wide range of marine fishes Such reductions of genetic diversity in some of the world's most abundant species may lead to a long-term impact of fishing on their evolutionary potential, particularly if abundance remains low and diversity continues to decay
263 citations
TL;DR: In this article, the authors explored patterns of phylogenetic correlation in the distribution of gut bacteria among species of turtle ants (genus Cephalotes), which host a dense gut microbial community.
Abstract: Correlation between gut microbiota and host phylogeny could reflect codiversification over shared evolutionary history or a selective environment that is more similar in related hosts. These alternatives imply substantial differences in the relationship between host and symbiont, but can they be distinguished based on patterns in the community data themselves? We explored patterns of phylogenetic correlation in the distribution of gut bacteria among species of turtle ants (genus Cephalotes), which host a dense gut microbial community. We used 16S rRNA pyrosequencing from 25 Cephalotes species to show that their gut community is remarkably stable, from the colony to the genus level. Despite this overall similarity, the existing differences among species' microbiota significantly correlated with host phylogeny. We introduced a novel analytical technique to test whether these phylogenetic correlations are derived from recent bacterial evolution, as would be expected in the case of codiversification, or from broader shifts more likely to reflect environmental filters imposed by factors such as diet or habitat. We also tested this technique on a published data set of ape microbiota, confirming earlier results while revealing previously undescribed patterns of phylogenetic correlation. Our results indicated a high degree of partner fidelity in the Cephalotes microbiota, suggesting that vertical transmission of the entire community could play an important role in the evolution and maintenance of the association. As additional comparative microbiota data become available, the techniques presented here can be used to explore trends in the evolution of host-associated microbial communities.
TL;DR: It is concluded that microbial community composition must be viewed within an informed context of host ecology and physiology, and that this is of particular importance with respect to research planning and sampling design.
Abstract: High-throughput sequencing approaches have enabled characterizations of the community composition of numerous gut microbial communities, which in turn has enhanced interest in their diversity and functional relationships in different groups of vertebrates. Although fishes represent the greatest taxonomic and ecological diversity of vertebrates, our understanding of their gut microbiota and its functional significance has lagged well behind that of terrestrial vertebrates. In order to highlight emerging issues, we provide an overview of research on fish gut microbiotas and the biology of their hosts. We conclude that microbial community composition must be viewed within an informed context of host ecology and physiology, and that this is of particular importance with respect to research planning and sampling design.
TL;DR: It is confirmed that both diet and phylogeny drive the evolution of mammalian gut microbiota, with cases of convergence in global composition, but also examples of phylogenetic inertia.
Abstract: Mammals have diversified into many dietary niches. Specialized myrmecophagous (ant- and termite-eating) placental mammals represent a textbook example of evolutionary convergence driven by extreme diet specialization. Armadillos, anteaters, aardvarks, pangolins and aardwolves thus provide a model system for understanding the potential role of gut microbiota in the convergent adaptation to myrmecophagy. Here, we expand upon previous mammalian gut microbiome studies by using high-throughput barcoded Illumina sequencing of the 16S rRNA gene to characterize the composition of gut microbiota in 15 species representing all placental myrmecophagous lineages and their close relatives from zoo- and field-collected samples. We confirm that both diet and phylogeny drive the evolution of mammalian gut microbiota, with cases of convergence in global composition, but also examples of phylogenetic inertia. Our results reveal specialized placental myrmecophages as a spectacular case of large-scale convergence in gut microbiome composition. Indeed, neighbour-net networks and beta-diversity plots based on UniFrac distances show significant clustering of myrmecophagous species (anteaters, aardvarks and aardwolves), even though they belong to phylogenetically distant lineages representing different orders. The aardwolf, which diverged from carnivorous hyenas only in the last 10 million years, experienced a convergent shift in the composition of its gut microbiome to become more similar to other myrmecophages. These results confirm diet adaptation to be a major driving factor of convergence in gut microbiome composition over evolutionary timescales. This study sets the scene for future metagenomic studies aiming at evaluating potential convergence in functional gene content in the microbiomes of specialized mammalian myrmecophages.
TL;DR: It is found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event per chromosome have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution.
Abstract: Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high-density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best-order map contained 4215 markers, with a total distance of 3132 cM and a mean genetic distance between markers of 0.12 cM. Facilitated by the array being designed to include markers from most scaffolds, we obtained a second-generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super-scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cM/Mb and correlated closely with chromosome size, from 2 cM/Mb for chromosomes >100 Mb to >10 cM/Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cM was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.
TL;DR: This work challenges the assertion that the original RAD protocol minimizes the impact of PCR artefacts relative to that of other RAD protocols, and highlights the strengths and weaknesses of four different approaches to RADseq which are a representative sample of all RAD variants.
Abstract: 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. AL 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.
TL;DR: Overall, fungal root symbionts appear to be more constrained by dispersal and biogeography than by host availability, which is consistent with previous work on rhizospheric fungi and bacteria.
Abstract: Microbial communities in plant roots provide critical links between above- and belowground processes in terrestrial ecosystems. Variation in root communities has been attributed to plant host effects and microbial host preferences, as well as to factors pertaining to soil conditions, microbial biogeography and the presence of viable microbial propagules. To address hypotheses regarding the influence of plant host and soil biogeography on root fungal and bacterial communities, we designed a trap-plant bioassay experiment. Replicate Populus, Quercus and Pinus plants were grown in three soils originating from alternate field sites. Fungal and bacterial community profiles in the root of each replicate were assessed through multiplex 454 amplicon sequencing of four loci (i.e., 16S, SSU, ITS, LSU rDNA). Soil origin had a larger effect on fungal community composition than did host species, but the opposite was true for bacterial communities. Populus hosted the highest diversity of rhizospheric fungi and bacteria. Root communities on Quercus and Pinus were more similar to each other than to Populus. Overall, fungal root symbionts appear to be more constrained by dispersal and biogeography than by host availability.
TL;DR: The results complement and extend existing theory on alternative stable states during population divergence, distinct phases of speciation and the rapid emergence of multilocus barriers to gene flow, and are a step towards aligning population genomic theory with modern empirical studies.
Abstract: A long-standing problem in evolutionary biology has been determining whether and how gradual, incremental changes at the gene level can account for rapid speciation and bursts of adaptive radiation. Using genome-scale computer simulations, we extend previous theory showing how gradual adaptive change can generate nonlinear population transitions, resulting in the rapid formation of new, reproductively isolated species. We show that these transitions occur via a mechanism rooted in a basic property of biological heredity: the organization of genes in genomes. Genomic organization of genes facilitates two processes: (i) the build-up of statistical associations among large numbers of genes and (ii) the action of divergent selection on persistent combinations of alleles. When a population has accumulated a critical amount of standing, divergently selected variation, the combination of these two processes allows many mutations of small effect to act synergistically and precipitously split one population into two discontinuous, reproductively isolated groups. Periods of allopatry, chromosomal linkage among loci, and large-effect alleles can facilitate this process under some conditions, but are not required for it. Our results complement and extend existing theory on alternative stable states during population divergence, distinct phases of speciation and the rapid emergence of multilocus barriers to gene flow. The results are thus a step towards aligning population genomic theory with modern empirical studies.
TL;DR: Pyrosequencing is applied to the elucidation of Symbiodinium diversity via analysis of the internal transcribed spacer 2 (ITS2) region, a multicopy genetic marker commonly used to analyse Symbiod inium diversity, and an operational taxonomic unit (OTU)‐based pipeline is developed to provisionally resolve ecologically discrete entities from intragenomic variation.
Abstract: The persistence of coral reef ecosystems relies on the symbiotic relationship between scleractinian corals and intracellular, photosynthetic dinoflagellates in the genus Symbiodinium. Genetic evidence indicates that these symbionts are biologically diverse and exhibit discrete patterns of environmental and host distribution. This makes the assessment of Symbiodinium diversity critical to understanding the symbiosis ecology of corals. Here, we applied pyrosequencing to the elucidation of Symbiodinium diversity via analysis of the internal transcribed spacer 2 (ITS2) region, a multicopy genetic marker commonly used to analyse Symbiodinium diversity. Replicated data generated from isoclonal Symbiodinium cultures showed that all genomes contained numerous, yet mostly rare, ITS2 sequence variants. Pyrosequencing data were consistent with more traditional denaturing gradient gel electrophoresis (DGGE) approaches to the screening of ITS2 PCR amplifications, where the most common sequences appeared as the most intense bands. Further, we developed an operational taxonomic unit (OTU)-based pipeline for Symbiodinium ITS2 diversity typing to provisionally resolve ecologically discrete entities from intragenomic variation. A genetic distance cut-off of 0.03 collapsed intragenomic ITS2 variants of isoclonal cultures into single OTUs. When applied to the analysis of field-collected coral samples, our analyses confirm that much of the commonly observed Symbiodinium ITS2 diversity can be attributed to intragenomic variation. We conclude that by analysing Symbiodinium populations in an OTU-based framework, we can improve objectivity, comparability and simplicity when assessing ITS2 diversity in field-based studies.
TL;DR: The authors' results indicate a number of well‐supported dependencies between abundances of certain AMF taxa and soil properties such as pH, soil fertility and texture, and a surprising lack of effect of available soil phosphorus on the AMF community profiles.
Abstract: Arbuscular mycorrhizal fungi (AMF) are ubiquitous soil fungi, forming mutualistic symbiosis with a majority of terrestrial plant species. They are abundant in nearly all soils, less diverse than soil prokaryotes and other intensively studied soil organisms and thus are promising candidates for universal indicators of land management legacies and soil quality degradation. However, insufficient data on how the composition of indigenous AMF varies along soil and landscape gradients have hampered the definition of baselines and effect thresholds to date. Here, indigenous AMF communities in 154 agricultural soils collected across Switzerland were profiled by quantitative real-time PCR with taxon-specific markers for six widespread AMF species. To identify the key determinants of AMF community composition, the profiles were related to soil properties, land management and site geography. Our results indicate a number of well-supported dependencies between abundances of certain AMF taxa and soil properties such as pH, soil fertility and texture, and a surprising lack of effect of available soil phosphorus on the AMF community profiles. Site geography, especially the altitude and large geographical distance, strongly affected AMF communities. Unexpected was the apparent lack of a strong land management effect on the AMF communities as compared to the other predictors, which could be due to the rarity of highly intensive and unsustainable land management in Swiss agriculture. In spite of the extensive coverage of large geographical and soil gradients, we did not identify any taxon suitable as an indicator of land use among the six taxa we studied.
TL;DR: It was found that the heritability captured by SNP markers asymptoted at about half the SNPs available, suggesting that denser marker panels are not necessarily required for precise and unbiased heritability estimates.
Abstract: The estimation of quantitative genetic parameters in wild populations is generally limited by the accuracy and completeness of the available pedigree information. Using relatedness at genomewide markers can potentially remove this limitation and lead to less biased and more precise estimates. We estimated heritability, maternal genetic effects and genetic correlations for body size traits in an unmanaged long-term study population of Soay sheep on St Kilda using three increasingly complete and accurate estimates of relatedness: (i) Pedigree 1, using observation-derived maternal links and microsatellite-derived paternal links; (ii) Pedigree 2, using SNP-derived assignment of both maternity and paternity; and (iii) whole-genome relatedness at 37 037 autosomal SNPs. In initial analyses, heritability estimates were strikingly similar for all three methods, while standard errors were systematically lower in analyses based on Pedigree 2 and genomic relatedness. Genetic correlations were generally strong, differed little between the three estimates of relatedness and the standard errors declined only very slightly with improved relatedness information. When partitioning maternal effects into separate genetic and environmental components, maternal genetic effects found in juvenile traits increased substantially across the three relatedness estimates. Heritability declined compared to parallel models where only a maternal environment effect was fitted, suggesting that maternal genetic effects are confounded with direct genetic effects and that more accurate estimates of relatedness were better able to separate maternal genetic effects from direct genetic effects. We found that the heritability captured by SNP markers asymptoted at about half the SNPs available, suggesting that denser marker panels are not necessarily required for precise and unbiased heritability estimates. Finally, we present guidelines for the use of genomic relatedness in future quantitative genetics studies in natural populations.
TL;DR: The various ecological sources of divergent selection and genomic changes are outlined, showing that gene loss and changes in gene expression and in genomic architecture are important adaptation processes, in addition to the more widely recognized processes of amino acid substitution and gene duplication.
Abstract: Fungi are ideal model organisms for dissecting the genomic bases of adaptive divergence in eukaryotes. They have simple morphologies and small genomes, occupy contrasting, well-identified ecological niches and tend to have short generation times, and many are amenable to experimental approaches. Fungi also display diverse lifestyles, from saprotrophs to pathogens or mutualists, and they play extremely important roles in both ecosystems and human activities, as wood decayers, mycorrhizal fungi, lichens, endophytes, plant and animal pathogens, and in fermentation or drug production. We review here recent insights into the patterns and mechanisms of adaptive divergence in fungi, including sources of divergence, genomic variation and, ultimately, speciation. We outline the various ecological sources of divergent selection and genomic changes, showing that gene loss and changes in gene expression and in genomic architecture are important adaptation processes, in addition to the more widely recognized processes of amino acid substitution and gene duplication. We also review recent findings regarding the interspecific acquisition of genomic variation and suggesting an important role for introgression, hybridization and horizontal gene transfers (HGTs). We show that transposable elements can mediate several of these genomic changes, thus constituting important factors for adaptation. Finally, we review the consequences of divergent selection in terms of speciation, arguing that genetic incompatibilities may not be as widespread as generally thought and that pleiotropy between adaptation and reproductive isolation is an important route of speciation in fungal pathogens.
TL;DR: It is suggested that bacterial communities may converge during assembly supporting determinism, whereas fungal communities show no such patterns, and differing drivers for fungal and bacterial trajectories during early primary succession in recently deglaciated soils are highlighted.
Abstract: Early community assembly of soil microbial communities is essential for pedogenesis and development of organic legacies. We examined fungal and bacterial successions along a well-established temperate glacier forefront chronosequence representing ~70 years of deglaciation to determine community assembly. As microbial communities may be heavily structured by establishing vegetation, we included nonvegetated soils as well as soils from underneath four plant species with differing mycorrhizal ecologies (Abies lasiocarpa, ectomycorrhizal; Luetkea pectinata, arbuscular mycorrhizal; Phyllodoce empetriformis, ericoid mycorrhizal; Saxifraga ferruginea, nonmycorrhizal). Our main objectives were to contrast fungal and bacterial successional dynamics and community assembly as well as to decouple the effects of plant establishment and time since deglaciation on microbial trajectories using high-throughput sequencing. Our data indicate that distance from glacier terminus has large effects on biomass accumulation, community membership, and distribution for both fungi and bacteria. Surprisingly, presence of plants rather than their identity was more important in structuring bacterial communities along the chronosequence and played only a very minor role in structuring the fungal communities. Further, our analyses suggest that bacterial communities may converge during assembly supporting determinism, whereas fungal communities show no such patterns. Although fungal communities provided little evidence of convergence in community structure, many taxa were nonrandomly distributed across the glacier foreland; similar taxon-level responses were observed in bacterial communities. Overall, our data highlight differing drivers for fungal and bacterial trajectories during early primary succession in recently deglaciated soils.
TL;DR: In this paper, the effects of admixture and spatial isolation on how biological diversity is organized in a group of Lycaeides butterflies are quantified using DNA sequences and genetic ancestry.
Abstract: Detailed information about the geographic distribution of genetic and genomic variation is necessary to better understand the organization and structure of biological diversity. In particular, spatial isolation within species and hybridization between them can blur species boundaries and create evolutionary relationships that are inconsistent with a strictly bifurcating tree model. Here, we analyse genome-wide DNA sequence and genetic ancestry variation in Lycaeides butterflies to quantify the effects of admixture and spatial isolation on how biological diversity is organized in this group. We document geographically widespread and pervasive historical admixture, with more restricted recent hybridization. This includes evidence supporting previously known and unknown instances of admixture. The genome composition of admixed individuals varies much more among than within populations, and tree- and genetic ancestry-based analyses indicate that multiple distinct admixed lineages or populations exist. We find that most genetic variants in Lycaeides are rare (minor allele frequency <0.5%). Because the spatial and taxonomic distributions of alleles reflect demographic and selective processes since mutation, rare alleles, which are presumably younger than common alleles, were spatially and taxonomically restricted compared with common variants. Thus, we show patterns of genetic variation in this group are multifaceted, and we argue that this complexity challenges simplistic notions concerning the organization of biological diversity into discrete, easily delineated and hierarchically structured entities.
TL;DR: This study determines the landscape and ecological factors affecting gene flow in the western slimy salamander and demonstrates that composite representations of ecologically explicit processes can provide novel insight and can better explain genetic differentiation than ecologically implicit landscape resistance surfaces.
Abstract: Landscape genetics has seen tremendous advances since its introduction, but parameterization and optimization of resistance surfaces still poses significant challenges. Despite increased availability and resolution of spatial data, few studies have integrated empirical data to directly represent ecological processes as genetic resistance surfaces. In our study, we determine the landscape and ecological factors affecting gene flow in the western slimy salamander (Plethodon albagula). We used field data to derive resistance surfaces representing salamander abundance and rate of water loss through combinations of canopy cover, topographic wetness, topographic position, solar exposure and distance from ravine. These ecologically explicit composite surfaces directly represent an ecological process or physiological limitation of our organism. Using generalized linear mixed-effects models, we optimized resistance surfaces using a nonlinear optimization algorithm to minimize model AIC. We found clear support for the resistance surface representing the rate of water loss experienced by adult salamanders in the summer. Resistance was lowest at intermediate levels of water loss and higher when the rate of water loss was predicted to be low or high. This pattern may arise from the compensatory movement behaviour of salamanders through suboptimal habitat, but also reflects the physiological limitations of salamanders and their sensitivity to extreme environmental conditions. Our study demonstrates that composite representations of ecologically explicit processes can provide novel insight and can better explain genetic differentiation than ecologically implicit landscape resistance surfaces. Additionally, our study underscores the fact that spatial estimates of habitat suitability or abundance may not serve as adequate proxies for describing gene flow, as predicted abundance was a poor predictor of genetic differentiation.
TL;DR: These findings offer a new explanation for heterogeneous genomic divergence and thus challenge the standard view that peaks in population divergence harbour divergently selected loci and that low-divergence regions result from balancing selection or localized introgression.
Abstract: Parallel adaptation is common and may often occur from shared genetic variation, but the genomic consequences of this process remain poorly understood. We first use individual-based simulations to demonstrate that comparisons between populations adapted in parallel to similar environments from shared variation reveal a characteristic genomic signature around a selected locus: a low-divergence valley centred at the locus and flanked by twin peaks of high divergence. This signature is initiated by the hitchhiking of haplotype tracts differing between derived populations in the broader neighbourhood of the selected locus (driving the high-divergence twin peaks) and shared haplotype tracts in the tight neighbourhood of the locus (driving the low-divergence valley). This initial hitchhiking signature is reinforced over time because the selected locus acts as a barrier to gene flow from the source to the derived populations, thus promoting divergence by drift in its close neighbourhood. We next empirically confirm the peak-valley-peak signature by combining targeted and RAD sequence data at three candidate adaptation genes in multiple marine (source) and freshwater (derived) populations of threespine stickleback. Finally, we use a genome-wide screen for the peak-valley-peak signature to discover additional genome regions involved in parallel marine-freshwater divergence. Our findings offer a new explanation for heterogeneous genomic divergence and thus challenge the standard view that peaks in population divergence harbour divergently selected loci and that low-divergence regions result from balancing selection or localized introgression. We anticipate that genome scans for peak-valley-peak divergence signatures will promote the discovery of adaptation genes in other organisms.
TL;DR: Changes in pH, moisture and organic nutrients availability caused shifts in the composition of eukaryotic microbial communities at multiple trophic levels, and a co‐occurrence based network analysis revealed complex potential interaction patterns between OTUs from different troPHic levels.
Abstract: To understand the fine-scale effects of changes in nutrient availability on eukaryotic soil microorganisms communities, a multiple barcoding approach was used to analyse soil samples from four different treatments in a long-term fertilization experiment. We performed PCR amplification on soil DNA with primer pairs specifically targeting the 18S rRNA genes of all eukaryotes and three protist groups (Cercozoa, Chrysophyceae-Synurophyceae and Kinetoplastida) as well as the ITS gene of fungi and the 23S plastid rRNA gene of photoautotrophic microorganisms. Amplicons were pyrosequenced, and a total of 88 706 quality filtered reads were clustered into 1232 operational taxonomic units (OTU) across the six data sets. Comparisons of the taxonomic coverage achieved based on overlapping assignment of OTUs revealed that half of the eukaryotic taxa identified were missed by the universal eukaryotic barcoding marker. There were only little differences in OTU richness observed between organic- (farmyard manure), mineral- and nonfertilized soils. However, the community compositions appeared to be strongly structured by organic fertilization in all data sets other than that generated using the universal eukaryotic 18S rRNA gene primers, whereas mineral fertilization had only a minor effect. In addition, a co-occurrence based network analysis revealed complex potential interaction patterns between OTUs from different trophic levels, for example between fungivorous flagellates and fungi. Our results demonstrate that changes in pH, moisture and organic nutrients availability caused shifts in the composition of eukaryotic microbial communities at multiple trophic levels.
TL;DR: In this article, the authors examined the bacterial community composition along two transects, one from the inner Pearl River estuary to the open waters of the South China Sea (SCS) and the other from the Luzon Strait to the SCS basin, using 454 pyrosequencing of the 16S rRNA and 16S RRNA gene (V1-3 regions) and thereby characterizing the active and total bacterial community, respectively.
Abstract: To test the hypothesis that different drivers shape the diversity and biogeography of the total and active bacterial community, we examined the bacterial community composition along two transects, one from the inner Pearl River estuary to the open waters of the South China Sea (SCS) and the other from the Luzon Strait to the SCS basin, using 454 pyrosequencing of the 16S rRNA and 16S rRNA gene (V1-3 regions) and thereby characterizing the active and total bacterial community, respectively. The diversity and biogeographic patterns differed substantially between the active and total bacterial communities. Although the composition of both the total and active bacterial community was strongly correlated with environmental factors and weakly correlated with geographic distance, the active bacterial community displayed higher environmental sensitivity than the total community and particularly a greater distance effect largely caused by the active assemblage from deep waters. The 16S rRNA vs. rDNA relationships indicated that the active bacteria were low in relative abundance in the SCS. This might be due to a high competition between active bacterial taxa as indicated by our community network models. Based on these analyses, we speculate that high competition could cause some dispersal limitation of the active bacterial community resulting in a distinct distance-decay relationship. Altogether, our results indicated that the biogeographic distribution of bacteria in the SCS is the result of both environmental control and distance decay.
TL;DR: In this paper, 381 SNPs located in transcribed regions were used to assess large and fine-scale population structure in the European hake (Merluccius merluccii), a widely distributed demersal species of high priority for the European fishery.
Abstract: Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome-wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large- and fine-scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (FCT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (FCT range 0.275�0.705) and fine-scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.
TL;DR: It is shown that MHC IIb polymorphism is associated with among‐individual variation in gut microbiota within a single wild vertebrate population of a small fish, the threespine stickleback, and speculated that macroparasite‐driven selection on MHC has the potential to indirectly alter the host gut microbiota, and vice versa.
Abstract: Animals harbour diverse communities of symbiotic bacteria, which differ dramatically among host individuals. This heterogeneity poses an immunological challenge: distinguishing between mutualistic and pathogenic members of diverse and host-specific microbial communities. We propose that Major Histocompatibility class II (MHC) genotypes contribute to recognition and regulation of gut microbes, and thus, MHC polymorphism contributes to microbial variation among hosts. Here, we show that MHC IIb polymorphism is associated with among-individual variation in gut microbiota within a single wild vertebrate population of a small fish, the threespine stickleback. We sampled stickleback from Cedar Lake, on Vancouver Island, and used next-generation sequencing to genotype the sticklebacks' gut microbiota (16S sequencing) and their MHC class IIb exon 2 sequences. The presence of certain MHC motifs was associated with altered relative abundance (increase or decrease) of some microbial Families. The effect sizes are modest and entail a minority of microbial taxa, but these results represent the first indication that MHC genotype may affect gut microbiota composition in natural populations (MHC-microbe associations have also been found in a few studies of lab mice). Surprisingly, these MHC effects were frequently sex-dependent. Finally, hosts with more diverse MHC motifs had less diverse gut microbiota. One implication is that MHC might influence the efficacy of therapeutic strategies to treat dysbiosis-associated disease, including the outcome of microbial transplants between healthy and diseased patients. We also speculate that macroparasite-driven selection on MHC has the potential to indirectly alter the host gut microbiota, and vice versa.
TL;DR: It is concluded that stored pollen is not evolved for microbially mediated nutrient conversion, but is a preservative environment due primarily to added honey, nectar, bee secretions and properties of pollen itself.
Abstract: Honey bee hives are filled with stored pollen, honey, plant resins and wax, all antimicrobial to differing degrees. Stored pollen is the nutritionally rich currency used for colony growth and consists of 40–50% simple sugars. Many studies speculate that prior to consumption by bees, stored pollen undergoes long-term nutrient conversion, becoming more nutritious ‘bee bread’ as microbes predigest the pollen. We quantified both structural and functional aspects associated with this hypothesis using behavioural assays, bacterial plate counts, microscopy and 454 amplicon sequencing of the 16S rRNA gene from both newly collected and hive-stored pollen. We found that bees preferentially consume fresh pollen stored for 96 h. The estimated microbe to pollen grain surface area ratio was 1:1 000 000 indicating a negligible effect of microbial metabolism on hive-stored pollen. Consistent with these findings, hive-stored pollen grains did not appear compromised according to microscopy. Based on year round 454 amplicon sequencing, bacterial communities of newly collected and hive-stored pollen did not differ, indicating the lack of an emergent microbial community co-evolved to digest stored pollen. In accord with previous culturing and 16S cloning, acid resistant and osmotolerant bacteria like Lactobacillus kunkeei were found in greatest abundance in stored pollen, consistent with the harsh character of this microenvironment. We conclude that stored pollen is not evolved for microbially mediated nutrient conversion, but is a preservative environment due primarily to added honey, nectar, bee secretions and properties of pollen itself.