Showing papers in "Molecular Ecology in 2018"

Cryptic species as a window into the paradigm shift of the species concept

Abstract: The species concept is the cornerstone of biodiversity science, and any paradigm shift in the delimitation of species affects many research fields. Many biologists now are embracing a new "species" paradigm as separately evolving populations using different delimitation criteria. Individual criteria can emerge during different periods of speciation; some may never evolve. As such, a paradigm shift in the species concept relates to this inherent heterogeneity in the speciation process and species category-which is fundamentally overlooked in biodiversity research. Cryptic species fall within this paradigm shift: they are continuously being reported from diverse animal phyla but are poorly considered in current tests of ecological and evolutionary theory. The aim of this review is to integrate cryptic species in biodiversity science. In the first section, we address that the absence of morphological diversification is an evolutionary phenomenon, a "process" counterpart to the long-studied mechanisms of morphological diversification. In the next section regarding taxonomy, we show that molecular delimitation of cryptic species is heavily biased towards distance-based methods. We also stress the importance of formally naming of cryptic species for better integration into research fields that use species as units of analysis. Finally, we show that incorporating cryptic species leads to novel insights regarding biodiversity patterns and processes, including large-scale biodiversity assessments, geographic variation in species distribution and species coexistence. It is time for incorporating multicriteria species approaches aiming to understand speciation across space and taxa, thus allowing integration into biodiversity conservation while accommodating for species uncertainty.

Comparing methods for detecting multilocus adaptation with multivariate genotype-environment associations.

Abstract: Identifying adaptive loci can provide insight into the mechanisms underlying local adaptation. Genotype-environment association (GEA) methods, which identify these loci based on correlations between genetic and environmental data, are particularly promising. Univariate methods have dominated GEA, despite the high dimensional nature of genotype and environment. Multivariate methods, which analyse many loci simultaneously, may be better suited to these data as they consider how sets of markers covary in response to environment. These methods may also be more effective at detecting adaptive processes that result in weak, multilocus signatures. Here, we evaluate four multivariate methods and five univariate and differentiation-based approaches, using published simulations of multilocus selection. We found that Random Forest performed poorly for GEA. Univariate GEAs performed better, but had low detection rates for loci under weak selection. Constrained ordinations, particularly redundancy analysis (RDA), showed a superior combination of low false-positive and high true-positive rates across all levels of selection. These results were robust across the demographic histories, sampling designs, sample sizes and weak population structure tested here. The value of combining detections from different methods was variable and depended on the study goals and knowledge of the drivers of selection. Re-analysis of genomic data from grey wolves highlighted the unique, covarying sets of adaptive loci that could be identified using RDA. Although additional testing is needed, this study indicates that RDA is an effective means of detecting adaptation, including signatures of weak, multilocus selection, providing a powerful tool for investigating the genetic basis of local adaptation.

These aren't the loci you'e looking for: Principles of effective SNP filtering for molecular ecologists.

Abstract: Sequencing reduced-representation libraries of restriction site-associated DNA (RADseq) to identify single nucleotide polymorphisms (SNPs) is quickly becoming a standard methodology for molecular ecologists. Because of the scale of RADseq data sets, putative loci cannot be assessed individually, making the process of filtering noise and correctly identifying biologically meaningful signal more difficult. Artefacts introduced during library preparation and/or bioinformatic processing of SNP data can create patterns that are incorrectly interpreted as indicative of population structure or natural selection. Therefore, it is crucial to carefully consider types of errors that may be introduced during laboratory work and data processing, and how to minimize, detect and remove these errors. Here, we discuss issues inherent to RADseq methodologies that can result in artefacts during library preparation and locus reconstruction resulting in erroneous SNP calls and, ultimately, genotyping error. Further, we describe steps that can be implemented to create a rigorously filtered data set consisting of markers accurately representing independent loci and compare the effect of different combinations of filters on four RAD data sets. At last, we stress the importance of publishing raw sequence data along with final filtered data sets in addition to detailed documentation of filtering steps and quality control measures.

Scaling up: A guide to high-throughput genomic approaches for biodiversity analysis.

Abstract: The purpose of this review is to present the most common and emerging DNA-based methods used to generate data for biodiversity and biomonitoring studies. As environmental assessment and monitoring programmes may require biodiversity information at multiple levels, we pay particular attention to the DNA metabarcoding method and discuss a number of bioinformatic tools and considerations for producing DNA-based indicators using operational taxonomic units (OTUs), taxa at a variety of ranks and community composition. By developing the capacity to harness the advantages provided by the newest technologies, investigators can "scale up" by increasing the number of samples and replicates processed, the frequency of sampling over time and space, and even the depth of sampling such as by sequencing more reads per sample or more markers per sample. The ability to scale up is made possible by the reduced hands-on time and cost per sample provided by the newest kits, platforms and software tools. Results gleaned from broad-scale monitoring will provide opportunities to address key scientific questions linked to biodiversity and its dynamics across time and space as well as being more relevant for policymakers, enabling science-based decision-making, and provide a greater socio-economic impact. As genomic approaches are continually evolving, we provide this guide to methods used in biodiversity genomics.

The search for loci under selection: trends, biases and progress

Abstract: Detecting genetic variants under selection using FST outlier analysis (OA) and environmental association analyses (EAAs) are popular approaches that provide insight into the genetic basis of local adaptation. Despite the frequent use of OA and EAA approaches and their increasing attractiveness for detecting signatures of selection, their application to field-based empirical data have not been synthesized. Here, we review 66 empirical studies that use Single Nucleotide Polymorphisms (SNPs) in OA and EAA. We report trends and biases across biological systems, sequencing methods, approaches, parameters, environmental variables and their influence on detecting signatures of selection. We found striking variability in both the use and reporting of environmental data and statistical parameters. For example, linkage disequilibrium among SNPs and numbers of unique SNP associations identified with EAA were rarely reported. The proportion of putatively adaptive SNPs detected varied widely among studies, and decreased with the number of SNPs analysed. We found that genomic sampling effort had a greater impact than biological sampling effort on the proportion of identified SNPs under selection. OA identified a higher proportion of outliers when more individuals were sampled, but this was not the case for EAA. To facilitate repeatability, interpretation and synthesis of studies detecting selection, we recommend that future studies consistently report geographical coordinates, environmental data, model parameters, linkage disequilibrium, and measures of genetic structure. Identifying standards for how OA and EAA studies are designed and reported will aid future transparency and comparability of SNP-based selection studies and help to progress landscape and evolutionary genomics.

Rates of gut microbiome divergence in mammals.

Abstract: The variation and taxonomic diversity among mammalian gut microbiomes raises several questions about the factors that contribute to the rates and patterns of change in these microbial communities. By comparing the microbiome compositions of 112 species representing 14 mammalian orders, we assessed how host and ecological factors contribute to microbiome diversification. Except in rare cases, the same bacterial phyla predominate in mammalian gut microbiomes, and there has been some convergence of microbiome compositions according to dietary category across all mammalians lineages except Chiropterans (bats), which possess high proportions of Proteobacteria and tend to be most similar to one another regardless of diet. At lower taxonomic ranks (families, genera, 97% OTUs), bacteria are more likely to be associated with a particular mammalian lineage than with a particular dietary category, resulting in a strong phylogenetic signal in the degree to which microbiomes diverge. Despite different physiologies, the gut microbiomes of several mammalian lineages have diverged at roughly the same rate over the past 75 million years; however, the gut microbiomes of Cetartiodactyla (ruminants, whales, hippopotami) have evolved much faster and those of Chiropterans much slower. Contrary to expectations, the number of dietary transitions within a lineage does not influence rates of microbiome divergence, but instead, some of the most dramatic changes are associated with the loss of bacterial taxa, such as those accompanying the transition from terrestrial to marine lifestyles and the evolution of hominids.

The role of gene flow in rapid and repeated evolution of cave-related traits in Mexican tetra, Astyanax mexicanus

Abstract: Understanding the molecular basis of repeatedly evolved phenotypes can yield key insights into the evolutionary process. Quantifying gene flow between populations is especially important in interpreting mechanisms of repeated phenotypic evolution, and genomic analyses have revealed that admixture occurs more frequently between diverging lineages than previously thought. In this study, we resequenced 47 whole genomes of the Mexican tetra from three cave populations, two surface populations and outgroup samples. We confirmed that cave populations are polyphyletic and two Astyanax mexicanus lineages are present in our data set. The two lineages likely diverged much more recently than previous mitochondrial estimates of 5-7 mya. Divergence of cave populations from their phylogenetically closest surface population likely occurred between ~161 and 191 k generations ago. The favoured demographic model for most population pairs accounts for divergence with secondary contact and heterogeneous gene flow across the genome, and we rigorously identified gene flow among all lineages sampled. Therefore, the evolution of cave-related traits occurred more rapidly than previously thought, and trogolomorphic traits are maintained despite gene flow with surface populations. The recency of these estimated divergence events suggests that selection may drive the evolution of cave-derived traits, as opposed to disuse and drift. Finally, we show that a key trogolomorphic phenotype QTL is enriched for genomic regions with low divergence between caves, suggesting that regions important for cave phenotypes may be transferred between caves via gene flow. Our study shows that gene flow must be considered in studies of independent, repeated trait evolution.

Evolution of heat-shock protein expression underlying adaptive responses to environmental stress

Abstract: Heat-shock proteins (Hsps) and their cognates are primary mitigators of cell stress. With increasingly severe impacts of climate change and other human modifications of the biosphere, the ability of the heat-shock system to affect evolutionary fitness in environments outside the laboratory and to evolve in response is topic of growing importance. Since the last major reviews, several advances have occurred. First, demonstrations of the heat-shock response outside the laboratory now include many additional taxa and environments. Many of these demonstrations are only correlative, however. More importantly, technical advances in "omic" quantification of nucleic acids and proteins, genomewide association analysis, and manipulation of genes and their expression have enabled the field to move beyond correlation. Several consequent advances are already evident: The pathway from heat-shock gene expression to stress tolerance in nature can be extremely complex, mediated through multiple biological processes and systems, and even multiple species. The underlying genes are more numerous, diverse and variable than previously appreciated, especially with respect to their regulatory variation and epigenetic changes. The impacts and limitations (e.g., due to trade-offs) of natural selection on these genes have become more obvious and better established. At last, as evolutionary capacitors, Hsps may have distinctive impacts on the evolution of other genes and ecological consequences.

Transitions in sex determination and sex chromosomes across vertebrate species

Abstract: Despite the prevalence of sexual reproduction across eukaryotes, there is a remarkable diversity of sex-determination mechanisms. The underlying causes of this diversity remain unclear, and it is unknown whether there are convergent trends in the directionality of turnover in sex-determination mechanisms. We used the recently assembled Tree of Sex database to assess patterns in the evolution of sex-determination systems in the remarkably diverse vertebrate clades of teleost fish, squamate reptiles and amphibians. Contrary to theoretical predictions, we find no evidence that the evolution of separate sexes is irreversible, as transitions from separate sexes to hermaphroditism occur at higher rates than the reverse in fish. We also find that transitions from environmental sex determination to genetic sex determination occur at higher rates than the reverse in both squamates and fish, suggesting that genetic sex determination is more stable. However, our data are not consistent with the hypothesis that heteromorphic sex chromosomes are an "evolutionary trap." Rather, we find similar transition rates between homomorphic and heteromorphic sex chromosomes in both fish and amphibians, and to environmental sex determination from heteromorphic vs. homomorphic sex chromosome systems in fish. Finally, we find that transitions between male and female heterogamety occur at similar rates in amphibians and squamates, while transitions to male heterogamety occur at higher rates in fish. Together, these results provide the most comprehensive view to date of the evolution of vertebrate sex determination in a phylogenetic context, providing new insight into long-standing questions about the evolution of sexual reproduction.

Why the COI barcode should be the community DNA metabarcode for the metazoa.

Abstract: Metabarcoding of complex metazoan communities is increasingly being used to measure biodiversity in terrestrial, freshwater and marine ecosystems, revolutionizing our ability to observe patterns and infer processes regarding the origin and conservation of biodiversity. A fundamentally important question is which genetic marker to amplify, and although the mitochondrial cytochrome oxidase subunit I (COI) gene is one of the more widely used markers in metabarcoding for the Metazoa, doubts have recently been raised about its suitability. We argue that (a) the extensive coverage of reference sequence databases for COI; (b) the variation it presents; (c) the comparative advantages for denoising protein-coding genes; and (d) recent advances in DNA sequencing protocols argue in favour of standardizing for the use of COI for metazoan community samples. We also highlight where research efforts should focus to maximize the utility of metabarcoding.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics.

Abstract: Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.

Two key features influencing community assembly processes at regional scale: Initial state and degree of change in environmental conditions.

Abstract: Belowground microbial communities strongly influence ecosystem function such that predicting function may rely on understanding ecological processes that assemble communities. Uncertainty remains, however, in what governs the relative contributions of different ecological processes. To help fill this knowledge gap, we test the general hypothesis that both initial state and degree of change in environmental conditions govern the relative contributions of different ecological assembly processes. To do so, we leveraged regional-scale nutrient and organic matter addition experiments and used soil organic matter (SOM) as a proxy of integrated soil environmental conditions. Consistent with our hypothesis, we found that both the initial amount of SOM and the degree of change in SOM-in response to nutrient addition-influenced the relative contributions of different ecological assembly processes. These influences were most clearly observed at the regional scale, suggesting potential scale dependence. More specifically, nutrient additions homogenized bacterial community composition due to enhanced influences of homogenizing dispersal when SOM content was initially high. In contrast, nutrient additions led to divergence in community composition due to variable selection when initial SOM was low and/or when SOM increased significantly in response to nutrient additions. Our findings indicate important connections among initial conditions, degree of change in environmental variables and microbial community assembly processes that may influence ecosystem processes. These conceptual inferences highlight a need to strengthen connections between ecological theory and biogeochemical modelling.

Demographic and genetic approaches to study dispersal in wild animal populations: A methodological review

Abstract: Dispersal is a central process in ecology and evolution At the individual level, the three stages of the dispersal process (ie, emigration, transience and immigration) are affected by complex interactions between phenotypes and environmental factors Condition- and context-dependent dispersal have far-reaching consequences, both for the demography and the genetic structuring of natural populations and for adaptive processes From an applied point of view, dispersal also deeply affects the spatial dynamics of populations and their ability to respond to land-use changes, habitat degradation and climate change For these reasons, dispersal has received considerable attention from ecologists and evolutionary biologists Demographic and genetic methods allow quantifying non-effective (ie, followed or not by a successful reproduction) and effective (ie, with a successful reproduction) dispersal and to investigate how individual and environmental factors affect the different stages of the dispersal process Over the past decade, demographic and genetic methods designed to quantify dispersal have rapidly evolved but interactions between researchers from the two fields are limited We here review recent developments in both demographic and genetic methods to study dispersal in wild animal populations We present their strengths and limits, as well as their applicability depending on study objectives and population characteristics We propose a unified framework allowing researchers to combine methods and select the more suitable tools to address a broad range of important topics about the ecology and evolution of dispersal and its consequences on animal population dynamics and genetics

Fungal community reveals less dispersal limitation and potentially more connected network than that of bacteria in bamboo forest soils

Abstract: A central aim of this microbial ecology research was to investigate the mechanisms shaping the assembly of soil microbial communities. Despite the importance of bacterial and fungal mediation of carbon cycling in forest ecosystems, knowledge concerning their distribution patterns and underlying mechanisms remains insufficient. Here, soils were sampled from six bamboo forests across the main planting area of Moso bamboo in southern China. The bacterial and fungal diversities were assessed by sequencing 16S rRNA and ITS gene amplicons, respectively, with an Illumina MiSeq. Based on structural equation modelling, dispersal limitation had strongest impact on bacterial beta diversity, while the mean annual precipitation had a smaller impact by directly or indirectly mediating the soil organic carbon density. However, only the mean annual temperature and precipitation played direct roles in fungal beta diversity. Moreover, the co-occurrence network analyses revealed a possibly much higher network connectivity in the fungal network than in the bacteria. With less dispersal limitation, stronger environmental selection and a potentially more connected network, the fungal community had more important roles in the soil carbon metabolisms in bamboo forests. Fungal beta diversity and the clustering coefficient explained approximately 14.4% and 6.1% of the variation in the carbon metabolic profiles among sites, respectively, but that of bacteria only explained approximately 1.7% and 1.8%, respectively. This study explored soil microbial spatial patterns along with the underlying mechanisms of dispersal limitation, selection and connectivity of ecological networks, thus providing novel insights into the study of the distinct functional traits of different microbial taxa.

Agriculture shapes the trophic niche of a bat preying on multiple pest arthropods across Europe: Evidence from DNA metabarcoding

Abstract: The interaction between agricultural production and wildlife can shape, and even condition, the functioning of both systems. In this study, we i) explored the degree to which a widespread European bat, namely the common bent-wing bat Miniopterus schreibersii, consumes crop-damaging insects at a continental scale, and ii) tested whether its dietary niche is shaped by the extension and type of agricultural fields. We employed a dual-primer DNA metabarcoding approach to characterize arthropod 16S and COI DNA sequences within bat faecal pellets collected across 16 Southern European localities, to first characterize the bat species' dietary niche, second measure the incidence of agricultural pests across their ranges and third assess whether geographical dietary variation responds to climatic, landscape diversity, agriculture type and vegetation productivity factors. We detected 12 arthropod orders, among which lepidopterans were predominant. We identified >200 species, 44 of which are known to cause agricultural damage. Pest species were detected at all but one sampling site and in 94% of the analysed samples. Furthermore, the dietary diversity of M. schreibersii exhibited a negative linear relation with the area of intensive agricultural fields, thus suggesting crops restrict the dietary niche of bats to prey taxa associated with agricultural production within their foraging range. Overall, our results imply that M. schreibersii might be a valuable asset for biological pest suppression in a variety of agricultural productions and highlight the dynamic interplay between wildlife and agricultural systems.

Evaluating genomic signatures of “the large X‐effect” during complex speciation

Abstract: The ubiquity of the “two rules of speciation”— Haldane’s rule and the large X-effect— imply a general, special role for sex chromosomes in the evolution of intrinsic postzygotic reproductive isolation. The recent proliferation of genome-scale analyses has revealed two further general observations: (1) complex speciation involving some form of gene flow is not uncommon; and (2) sex chromosomes in male- and in female-heterogametic taxa tend to show elevated differentiation relative to autosomes. Together these observations are consistent with speciation histories in which population genetic differentiation at autosomal loci is reduced by gene flow while natural selection against hybrid incompatibilities renders sex chromosomes relatively refractory to gene flow. Here I summarize multi-locus population genetic and population genomic evidence for greater differentiation on the X (or Z) versus the autosomes and consider the possible causes. I review common population genetic circumstances involving no selection and/or no interspecific gene flow that are nevertheless expected to elevate differentiation on sex chromosomes relative to autosomes. I then review theory for why large X-effects exist for hybrid incompatibilities and, more generally, for loci mediating local adaptation. The observed levels of sex chromosome versus autosomal differentiation, in many cases, appear consistent with simple explanations requiring neither large X-effects nor gene flow. Discerning signatures of large X-effects during complex speciation will therefore require analyses that go beyond chromosome-scale summaries of population genetic differentiation, explicitly test for differential introgression, and/or integrate experimental genetic data.

Sex chromosomes and speciation in birds and other ZW systems.

Abstract: Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.

Seascape genomics reveals adaptive divergence in a connected and commercially important mollusc, the greenlip abalone (Haliotis laevigata), along a longitudinal environmental gradient

Abstract: Populations of broadcast spawning marine organisms often have large sizes and are exposed to reduced genetic drift. Under such scenarios, strong selection associated with spatial environmental heterogeneity is expected to drive localized adaptive divergence, even in the face of connectivity. We tested this hypothesis using a seascape genomics approach in the commercially important greenlip abalone (Haliotis laevigata). We assessed how its population structure has been influenced by environmental heterogeneity along a zonal coastal boundary in southern Australia linked by strong oceanographic connectivity. Our data sets include 9,109 filtered SNPs for 371 abalones from 13 localities and environmental mapping across ~800 km. Genotype-environment association analyses and outlier tests defined 8,786 putatively neutral and 323 candidate adaptive loci. From a neutral perspective, the species is better represented by a metapopulation with very low differentiation (global FST = 0.0081) and weak isolation by distance following a stepping-stone model. For the candidate adaptive loci, however, model-based and model-free approaches indicated five divergent population clusters. After controlling for spatial distance, the distribution of putatively adaptive variation was strongly correlated to selection linked to minimum sea surface temperature and oxygen concentration. Around 80 candidates were annotated to genes with functions related to high temperature and/or low oxygen tolerance, including genes that influence the resilience of abalone species found in other biogeographic regions. Our study includes a documented example about the uptake of genomic information in fisheries management and supports the hypothesis of adaptive divergence due to coastal environmental heterogeneity in a connected metapopulation of a broadcast spawner.

The joint evolution of the Myxozoa and their alternate hosts: a cnidarian recipe for success and vast biodiversity.

Abstract: The relationships between parasites and their hosts are intimate, dynamic and complex; the evolution of one is inevitably linked to the other. Despite multiple origins of parasitism in the Cnidaria, only parasites belonging to the Myxozoa are characterized by a complex life cycle, alternating between fish and invertebrate hosts, as well as by high species diversity. This inspired us to examine the history of adaptive radiations in myxozoans and their hosts by determining the degree of congruence between their phylogenies and by timing the emergence of myxozoan lineages in relation to their hosts. Recent genomic analyses suggested a common origin of Polypodium hydriforme, a cnidarian parasite of acipenseriform fishes, and the Myxozoa, and proposed fish as original hosts for both sister lineages. We demonstrate that the Myxozoa emerged long before fish populated Earth and that phylogenetic congruence with their invertebrate hosts is evident down to the most basal branches of the tree, indicating bryozoans and annelids as original hosts and challenging previous evolutionary hypotheses. We provide evidence that, following invertebrate invasion, fish hosts were acquired multiple times, leading to parallel cospeciation patterns in all major phylogenetic lineages. We identify the acquisition of vertebrate hosts that facilitate alternative transmission and dispersion strategies as reason for the distinct success of the Myxozoa, and identify massive host specification-linked parasite diversification events. The results of this study transform our understanding of the origins and evolution of parasitism in the most basal metazoan parasites known.

Faster-X evolution: Theory and evidence from Drosophila.

Abstract: A faster rate of adaptive evolution of X-linked genes compared with autosomal genes can be caused by the fixation of recessive or partially recessive advantageous mutations, due to the full expression of X-linked mutations in hemizygous males. Other processes, including recombination rate and mutation rate differences between X chromosomes and autosomes, may also cause faster evolution of X-linked genes. We review population genetics theory concerning the expected relative values of variability and rates of evolution of X-linked and autosomal DNA sequences. The theoretical predictions are compared with data from population genomic studies of several species of Drosophila. We conclude that there is evidence for adaptive faster-X evolution of several classes of functionally significant nucleotides. We also find evidence for potential differences in mutation rates between X-linked and autosomal genes, due to differences in mutational bias towards GC to AT mutations. Many aspects of the data are consistent with the male hemizygosity model, although not all possible confounding factors can be excluded.

Genomewide SNP markers breathe new life into phylogeography and species delimitation for the problematic short-necked turtles (Chelidae: Emydura) of eastern Australia.

Abstract: Understanding the evolutionary history of diversifying lineages and the delineation of evolutionarily significant units and species remains major challenges for evolutionary biology. Low-cost representational sampling of the genome for single nucleotide polymorphisms shows great potential at the temporal scales that are typically the focus of species delimitation and phylogeography. We apply these markers to a case study of a freshwater turtle, Emydura macquarii, whose systematics has so far defied resolution, to bring to light a dynamic system of substantive allopatric lineages diverging on independent evolutionary trajectories, but held back in the process of speciation by low level and episodic exchange of alleles across drainage divides on various timescales. In the context of low-level episodic gene flow, speciation is often reticulate, rather than a bifurcating process. We argue that species delimitation needs to take into account the pattern of ancestry and descent of diverging lineages in allopatry together with the recent and contemporary processes of dispersal and gene flow that retard and obscure that divergence. Underpinned by a strong focus on lineage diagnosability, this combined approach provides a means for addressing the challenges of incompletely isolated populations with uncommon, but recurrent gene flow in studies of species delimitation, a situation likely to be frequently encountered. Taxonomic decisions in cases of allopatry often require subjective judgements. Our strategy, which adds an additional level of objectivity before that subjectivity is applied, reduces the risk of taxonomic inflation that can accompany lineage approaches to species delimitation.

Asymmetric oceanographic processes mediate connectivity and population genetic structure, as revealed by RADseq, in a highly dispersive marine invertebrate (Parastichopus californicus)

Abstract: Marine populations are typically characterized by weak genetic differentiation due to the potential for long-distance dispersal favouring high levels of gene flow. However, strong directional advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (FST = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R2 = .49) than geographic distance (IBD) (R2 = .18), and directional processes played an important role in shaping fine-scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine species.

Antibiotics reduce genetic diversity of core species in the honeybee gut microbiome.

Abstract: The gut microbiome plays a key role in animal health, and perturbing it can have detrimental effects. One major source of perturbation to microbiomes, in humans and human-associated animals, is exposure to antibiotics. Most studies of how antibiotics affect the microbiome have used amplicon sequencing of highly conserved 16S rRNA sequences, as in a recent study showing that antibiotic treatment severely alters the species-level composition of the honeybee gut microbiome. But because the standard 16S rRNA-based methods cannot resolve closely related strains, strain-level changes could not be evaluated. To address this gap, we used amplicon sequencing of protein-coding genes to assess effects of antibiotics on fine-scale genetic diversity of the honeybee gut microbiota. We followed the population dynamics of alleles within two dominant core species of the bee gut community, Gilliamella apicola and Snodgrassella alvi, following antibiotic perturbation. Whereas we observed a large reduction of genetic diversity in G. apicola, S. alvi diversity was mostly unaffected. The reduction of G. apicola diversity accompanied an increase in the frequency of several alleles, suggesting resistance to antibiotic treatment. We find that antibiotic perturbation can cause major shifts in diversity, and that the extent of these shifts can vary substantially across species. Thus, antibiotics impact not only species composition, but also allelic diversity within species, potentially affecting hosts if variants with particular functions are reduced or eliminated. Overall, we show that amplicon sequencing of protein-coding genes, without clustering into operational taxonomic units (OTUs), provides an accurate picture of the fine-scale dynamics of microbial communities over time.

Natural variation in stomata size contributes to the local adaptation of water-use efficiency in Arabidopsis thaliana.

Abstract: Stomata control gas exchanges between the plant and the atmosphere. How natural variation in stomata size and density contributes to resolve trade-offs between carbon uptake and water loss in response to local climatic variation is not yet understood. We developed an automated confocal microscopy approach to characterize natural genetic variation in stomatal patterning in 330 fully sequenced Arabidopsis thaliana accessions collected throughout the European range of the species. We compared this to variation in water-use efficiency, measured as carbon isotope discrimination (δ13 C). We detect substantial genetic variation for stomata size and density segregating within Arabidopsis thaliana. A positive correlation between stomata size and δ13 C further suggests that this variation has consequences on water-use efficiency. Genome wide association analyses indicate a complex genetic architecture underlying not only variation in stomatal patterning but also to its covariation with carbon uptake parameters. Yet, we report two novel QTL affecting δ13 C independently of stomatal patterning. This suggests that, in A. thaliana, both morphological and physiological variants contribute to genetic variance in water-use efficiency. Patterns of regional differentiation and covariation with climatic parameters indicate that natural selection has contributed to shape some of this variation, especially in Southern Sweden, where water availability is more limited in spring relative to summer. These conditions are expected to favour the evolution of drought avoidance mechanisms over drought escape strategies.

Macroalgal-bacterial interactions: Role of dimethylsulfoniopropionate in microbial gardening by Ulva (Chlorophyta).

Abstract: The marine macroalga Ulva mutabilis (Chlorophyta) develops into callus-like colonies consisting of undifferentiated cells and abnormal cell walls under axenic conditions. Ulva mutabilis is routinely cultured with two bacteria, the Roseovarius sp. MS2 strain and the Maribacter sp. MS6 strain, which release morphogenetic compounds and ensure proper algal morphogenesis. Using this tripartite community as an emerging model system, we tested the hypothesis that the bacterial-algal interactions evolved as a result of mutually taking advantage of signals in the environment. Our study aimed to determine whether cross-kingdom crosstalk is mediated by the attraction of bacteria through algal chemotactic signals. Roseovarius sp. MS2 senses the known osmolyte dimethylsulfoniopropionate (DMSP) released by Ulva into the growth medium. Roseovarius sp. is attracted by DMSP and takes it up rapidly such that DMSP can only be determined in axenic growth media. As DMSP did not promote bacterial growth under the tested conditions, Roseovarius benefited solely from glycerol as the carbon source provided by Ulva. Roseovarius quickly catabolized DMSP into methanethiol (MeSH) and dimethylsulphide (DMS). We conclude that many bacteria can use DMSP as a reliable signal indicating a food source and promote the subsequent development and morphogenesis in Ulva.

Environmental DNA detection of rare and invasive fish species in two Great Lakes tributaries.

Abstract: The extraction and characterization of DNA from aquatic environmental samples offers an alternative, noninvasive approach for the detection of rare species. Environmental DNA, coupled with PCR and next-generation sequencing ("metabarcoding"), has proven to be very sensitive for the detection of rare aquatic species. Our study used a custom-designed group-specific primer set and next-generation sequencing for the detection of three species at risk (Eastern Sand Darter, Ammocrypta pellucida; Northern Madtom, Noturus stigmosus; and Silver Shiner, Notropis photogenis), one invasive species (Round Goby, Neogobius melanostomus) and an additional 78 native species from two large Great Lakes tributary rivers in southern Ontario, Canada: the Grand River and the Sydenham River. Of 82 fish species detected in both rivers using capture-based and eDNA methods, our eDNA method detected 86.2% and 72.0% of the fish species in the Grand River and the Sydenham River, respectively, which included our four target species. Our analyses also identified significant positive and negative species co-occurrence patterns between our target species and other identified species. Our results demonstrate that eDNA metabarcoding that targets the fish community as well as individual species of interest provides a better understanding of factors affecting the target species spatial distribution in an ecosystem than possible with only target species data. Additionally, eDNA is easily implemented as an initial survey tool, or alongside capture-based methods, for improved mapping of species distribution patterns.

Signatures of local adaptation along environmental gradients in a range-expanding damselfly (Ischnura elegans)

Abstract: Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range-expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site-associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V-ATPase) and visual processes (long-wavelength-sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate-mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients.

Oral microbiomes from hunter-gatherers and traditional farmers reveal shifts in commensal balance and pathogen load linked to diet.

Abstract: Maladaptation to modern diets has been implicated in several chronic disorders. Given the higher prevalence of disease such as dental caries and chronic gum diseases in industrialized societies, we sought to investigate the impact of different subsistence strategies on oral health and physiology, as documented by the oral microbiome. To control for confounding variables such as environment and host genetics, we sampled saliva from three pairs of populations of hunter-gatherers and traditional farmers living in close proximity in the Philippines. Deep shotgun sequencing of salivary DNA generated high-coverage microbiomes along with human genomes. Comparing these microbiomes with publicly available data from individuals living on a Western diet revealed that abundance ratios of core species were significantly correlated with subsistence strategy, with hunter-gatherers and Westerners occupying either end of a gradient of Neisseria against Haemophilus, and traditional farmers falling in between. Species found preferentially in hunter-gatherers included microbes often considered as oral pathogens, despite their hosts' apparent good oral health. Discriminant analysis of gene functions revealed vitamin B5 autotrophy and urease-mediated pH regulation as candidate adaptations of the microbiome to the hunter-gatherer and Western diets, respectively. These results suggest that major transitions in diet selected for different communities of commensals and likely played a role in the emergence of modern oral pathogens.

Spatial population genomics of the brown rat (Rattus norvegicus) in New York City.

Abstract: Human commensal species such as rodent pests are often widely distributed across cities and threaten both infrastructure and public health. Spatially explicit population genomic methods provide insights into movements for cryptic pests that drive evolutionary connectivity across multiple spatial scales. We examined spatial patterns of neutral genomewide variation in brown rats (Rattus norvegicus) across Manhattan, New York City (NYC), using 262 samples and 61,401 SNPs to understand (i) relatedness among nearby individuals and the extent of spatial genetic structure in a discrete urban landscape; (ii) the geographic origin of NYC rats, using a large, previously published data set of global rat genotypes; and (iii) heterogeneity in gene flow across the city, particularly deviations from isolation by distance. We found that rats separated by ≤200 m exhibit strong spatial autocorrelation (r = .3, p = .001) and the effects of localized genetic drift extend to a range of 1,400 m. Across Manhattan, rats exhibited a homogeneous population origin from rats that likely invaded from Great Britain. While traditional approaches identified a single evolutionary cluster with clinal structure across Manhattan, recently developed methods (e.g., fineSTRUCTURE, sPCA, EEMS) provided evidence of reduced dispersal across the island's less residential Midtown region resulting in fine-scale genetic structuring (FST = 0.01) and two evolutionary clusters (Uptown and Downtown Manhattan). Thus, while some urban populations of human commensals may appear to be continuously distributed, landscape heterogeneity within cities can drive differences in habitat quality and dispersal, with implications for the spatial distribution of genomic variation, population management and the study of widely distributed pests.

A comparison of genomic islands of differentiation across three young avian species pairs

Abstract: Detailed evaluations of genomic variation between sister species often reveal distinct chromosomal regions of high relative differentiation (ie, "islands of differentiation" in FST ), but there is much debate regarding the causes of this pattern We briefly review the prominent models of genomic islands of differentiation and compare patterns of genomic differentiation in three closely related pairs of New World warblers with the goal of evaluating support for the four models Each pair (MacGillivray's/mourning warblers; Townsend's/black-throated green warblers; and Audubon's/myrtle warblers) consists of forms that were likely separated in western and eastern North American refugia during cycles of Pleistocene glaciations and have now come into contact in western Canada, where each forms a narrow hybrid zone We show strong differences between pairs in their patterns of genomic heterogeneity in FST , suggesting differing selective forces and/or differing genomic responses to similar selective forces among the three pairs Across most of the genome, levels of within-group nucleotide diversity (πWithin ) are almost as large as levels of between-group nucleotide distance (πBetween ) within each pair, suggesting recent common ancestry and/or gene flow In two pairs, a pattern of the FST peaks having low πBetween suggests that selective sweeps spread between geographically differentiated groups, followed by local differentiation This "sweep-before-differentiation" model is consistent with signatures of gene flow within the yellow-rumped warbler species complex These findings add to our growing understanding of speciation as a complex process that can involve phases of adaptive introgression among partially differentiated populations