Boundary effects cause false signals of range expansions in population genomic data.

Studying range expansions (REs) is central for understanding genetic variation through space and time as well as for identifying refugia and biological invasions. Range expansions are characterized by serial founder events causing clines of decreasing diversity away from the center of origin and asymmetries in the two-dimensional allele frequency spectra. These asymmetries, summarized by the directionality index (ψ), are sensitive to REs and persist for longer than clines in genetic diversity. In continuous and finite meta-populations, genetic drift tends to be stronger at the edges of the species distribution in equilibrium populations and populations undergoing REs alike. Such boundary effects (BEs) are expected to affect geographic patterns in genetic diversity and ψ. Here we demonstrate that BEs cause high false positive rates in equilibrium meta-populations when testing for REs. In the simulations, the absolute value of ψ (|ψ|) in equilibrium data sets was proportional to the fixation index (FST). By fitting signatures of REs as a function of ɛ=|ψ|/FST and geographic clines in ψ, strong evidence for REs could be detected in data from a recent rapid invasion of the cane toad, Rhinella marina, in Australia, but not in 28 previously published empirical data sets from Australian scincid lizards that were significant for the standard RE tests. Thus, while clinal variation in ψ is still the most sensitive statistic to REs, to detect true signatures of REs in natural populations, its magnitude needs to be considered in relation to the overall levels of genetic structuring in the data.

Genetic insight into a polygenic trait using a novel genome-wide association approach in a wild amphibian population.

Body size variation is central in the evolution of life-history traits in amphibians, but the underlying genetic architecture of this complex trait is still largely unknown. Herein, we studied the genetic basis of body size and fecundity of the alternative morphotypes in a wild population of the Greek smooth newt (Lissotriton graecus). By combining a genome-wide association approach with linkage disequilibrium network analysis, we were able to identify clusters of highly correlated loci thus maximizing sequence data for downstream analysis. The putatively associated variants explained 12.8% to 44.5% of the total phenotypic variation in body size and were mapped to genes with functional roles in the regulation of gene expression and cell cycle processes. Our study is the first to provide insights into the genetic basis of complex traits in newts and provides a useful tool to identify loci potentially involved in fitness-related traits in small data sets from natural populations in non-model species.

Parallel evolution despite low genetic diversity in three-spined sticklebacks.

When populations repeatedly adapt to similar environments they can evolve similar phenotypes based on shared genetic mechanisms (parallel evolution). The likelihood of parallel evolution is affected by demographic history, as it depends on the standing genetic variation of the source population. The three-spined stickleback (Gasterosteus aculeatus) repeatedly colonized and adapted to brackish and freshwater. Most parallel evolution studies in G. aculeatus were conducted at high latitudes, where freshwater populations maintain connectivity to the source marine populations. Here, we analysed southern and northern European marine and freshwater populations to test two hypotheses. First, that southern European freshwater populations (which currently lack connection to marine populations) lost genetic diversity due to bottlenecks and inbreeding compared to their northern counterparts. Second, that the degree of genetic parallelism is higher among northern than southern European freshwater populations, as the latter have been subjected to strong drift due to isolation. The results show that southern populations exhibit lower genetic diversity but a higher degree of genetic parallelism than northern populations. Hence, they confirm the hypothesis that southern populations have lost genetic diversity, but this loss probably happened after they had already adapted to freshwater conditions, explaining the high degree of genetic parallelism in the south.

Cast away in the Adriatic: Low degree of parallel genetic differentiation in three-spined sticklebacks.

The three-spined stickleback (Gasterosteus aculeatus) has repeatedly and independently adapted to freshwater habitats from standing genetic variation (SGV) following colonization from the sea. However, in the Mediterranean Sea G. aculeatus is believed to have gone extinct, and thus the spread of locally adapted alleles between different freshwater populations via the sea since then has been highly unlikely. This is expected to limit parallel evolution, that is the extent to which phylogenetically related alleles can be shared among independently colonized freshwater populations. Using whole genome and 2b-RAD sequencing data, we compared levels of genetic differentiation and genetic parallelism of 15 Adriatic stickleback populations to 19 Pacific, Atlantic and Caspian populations, where gene flow between freshwater populations across extant marine populations is still possible. Our findings support previous studies suggesting that Adriatic populations are highly differentiated (average FST  ≈ 0.45), of low genetic diversity and connectivity, and likely to stem from multiple independent colonizations during the Pleistocene. Linkage disequilibrium network analyses in combination with linear mixed models nevertheless revealed several parallel marine-freshwater differentiated genomic regions, although still not to the extent observed elsewhere in the world. We hypothesize that current levels of genetic parallelism in the Adriatic lineages are a relic of freshwater adaptation from SGV prior to the extinction of marine sticklebacks in the Mediterranean that has persisted despite substantial genetic drift experienced by the Adriatic stickleback isolates.

Stable inversion clines in a grasshopper species group despite complex geographical history.

Chromosomal inversions are known to play roles in adaptation and differentiation in many species. They involve clusters of correlated genes (i.e., loci in linkage disequilibrium, LD) possibly associated with environmental variables. The grasshopper "species complex" Trimerotropis pallidipennis comprises several genetic lineages distributed from North to South America in arid and semi-arid high-altitude environments. The southernmost lineage, Trimerotropis sp., segregates for four to seven putative inversions that display clinal variation, possibly through adaptation to temperate environments. We analysed chromosomal, mitochondrial and genome-wide single nucleotide polymorphism data in 19 Trimerotropis sp. populations mainly distributed along two altitudinal gradients (MS and Ju). Populations across Argentina comprise two main chromosomally and genetically differentiated lineages: one distributed across the southernmost border of the "Andes Centrales," adding evidence for a differentiation hotspot in this area; and the other widely distributed in Argentina. Within the latter, network analytical approaches to LD found three clusters of correlated loci (LD-clusters), with inversion karyotypes explaining >79% of the genetic variation. Outlier loci associated with environmental variables mapped to two of these LD-clusters. Furthermore, despite the complex geographical history indicated by population genetic analyses, the clines in inversion karyotypes have remained stable for more than 20 generations, implicating their role in adaptation and differentiation within this lineage. We hypothesize that these clines could be the consequence of a coupling between extrinsic postzygotic barriers and spatially varying selection along environmental gradients resulting in a hybrid zone. These results provide a framework for future investigations about candidate genes implicated in rapid adaptation to new environments.

Population Structure Limits Parallel Evolution in Sticklebacks.

Population genetic theory predicts that small effective population sizes (Ne) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e., parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically codistributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbors less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius: 2,996 versus 379 single nucleotide polymorphisms located within 26 versus 9 genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius, respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution.

Genetic population structure constrains local adaptation in sticklebacks.

Repeated and independent adaptation to specific environmental conditions from standing genetic variation is common. However, if genetic variation is limited, the evolution of similar locally adapted traits may be restricted to genetically different and potentially less optimal solutions or prevented from happening altogether. Using a quantitative trait locus (QTL) mapping approach, we identified the genomic regions responsible for the repeated pelvic reduction (PR) in three crosses between nine-spined stickleback populations expressing full and reduced pelvic structures. In one cross, PR mapped to linkage group 7 (LG7) containing the gene Pitx1, known to control pelvic reduction also in the three-spined stickleback. In the two other crosses, PR was polygenic and attributed to 10 novel QTL, of which 90% were unique to specific crosses. When screening the genomes from 27 different populations for deletions in the Pitx1 regulatory element, these were only found in the population in which PR mapped to LG7, even though the morphological data indicated large-effect QTL for PR in several other populations as well. Consistent with the available theory and simulations parameterized on empirical data, we hypothesize that the observed variability in genetic architecture of PR is due to heterogeneity in the spatial distribution of standing genetic variation caused by >2× stronger population structuring among freshwater populations and >10× stronger genetic isolation by distance in the sea in nine-spined sticklebacks as compared to three-spined sticklebacks.

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe.

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat-related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model-based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.

On the causes of geographically heterogeneous parallel evolution in sticklebacks.

The three-spined stickleback (Gasterosteus aculeatus) is an important model system for the study of parallel evolution in the wild, having repeatedly colonized and adapted to freshwater from the sea throughout the northern hemisphere. Previous studies identified numerous genomic regions showing consistent genetic differentiation between freshwater and marine ecotypes but these had typically limited geographic sampling and mostly focused on the Eastern Pacific region. We analysed population genomic data from global samples of the three-spined stickleback marine and freshwater ecotypes to detect loci involved in parallel evolution at different geographic scales. Most signatures of parallel evolution were unique to the Eastern Pacific and trans-oceanic marine-freshwater differentiation was restricted to a limited number of shared genomic regions, including three chromosomal inversions. On the basis of simulations and empirical data, we demonstrate that this could result from the stochastic loss of freshwater-adapted alleles during the invasion of the Atlantic basin and selection against freshwater-adapted variants in the sea, both of which can reduce standing genetic variation available for freshwater adaptation outside the Eastern Pacific region. Moreover, the elevated linkage disequilibrium associated with marine-freshwater differentiation in the Eastern Pacific is consistent with secondary contact between marine and freshwater populations that evolved in isolation from each other during past glacial periods. Thus, contrary to what earlier studies from the Eastern Pacific region have led us to believe, parallel marine-freshwater differentiation in sticklebacks is far less prevalent and pronounced in all other parts of the species global distribution range.

Accounting for heteroscedasticity and censoring in chromosome partitioning analyses.

A fundamental assumption in quantitative genetics is that traits are controlled by many loci of small effect. Using genomic data, this assumption can be tested using chromosome partitioning analyses, where the proportion of genetic variance for a trait explained by each chromosome (h2 c ), is regressed on its size. However, as h2 c -estimates are necessarily positive (censoring) and the variance increases with chromosome size (heteroscedasticity), two fundamental assumptions of ordinary least squares (OLS) regression are violated. Using simulated and empirical data we demonstrate that these violations lead to incorrect inference of genetic architecture. The degree of bias depends mainly on the number of chromosomes and their size distribution and is therefore specific to the species; using published data across many different species we estimate that not accounting for this effect overall resulted in 28% false positives. We introduce a new and computationally efficient resampling method that corrects for inflation caused by heteroscedasticity and censoring and that works under a large range of dataset sizes and genetic architectures in empirical datasets. Our new method substantially improves the robustness of inferences from chromosome partitioning analyses.

Inferences of genetic architecture of bill morphology in house sparrow using a high-density SNP array point to a polygenic basis.

Understanding the genetic architecture of quantitative traits can provide insights into the mechanisms driving phenotypic evolution. Bill morphology is an ecologically important and phenotypically variable trait, which is highly heritable and closely linked to individual fitness. Thus, bill morphology traits are suitable candidates for gene mapping analyses. Previous studies have revealed several genes that may influence bill morphology, but the similarity of gene and allele effects between species and populations is unknown. Here, we develop a custom 200K SNP array and use it to examine the genetic basis of bill morphology in 1857 house sparrow individuals from a large-scale, island metapopulation off the coast of Northern Norway. We found high genomic heritabilities for bill depth and length, which were comparable with previous pedigree estimates. Candidate gene and genomewide association analyses yielded six significant loci, four of which have previously been associated with craniofacial development. Three of these loci are involved in bone morphogenic protein (BMP) signalling, suggesting a role for BMP genes in regulating bill morphology. However, these loci individually explain a small amount of variance. In combination with results from genome partitioning analyses, this indicates that bill morphology is a polygenic trait. Any studies of eco-evolutionary processes in bill morphology are therefore dependent on methods that can accommodate polygenic inheritance of the phenotype and molecular-scale evolution of genetic architecture.

Linkage disequilibrium clustering-based approach for association mapping with tightly linked genomewide data.

Genomewide association studies (GWAS) aim to identify genetic markers strongly associated with quantitative traits by utilizing linkage disequilibrium (LD) between candidate genes and markers. However, because of LD between nearby genetic markers, the standard GWAS approaches typically detect a number of correlated SNPs covering long genomic regions, making corrections for multiple testing overly conservative. Additionally, the high dimensionality of modern GWAS data poses considerable challenges for GWAS procedures such as permutation tests, which are computationally intensive. We propose a cluster-based GWAS approach that first divides the genome into many large nonoverlapping windows and uses linkage disequilibrium network analysis in combination with principal component (PC) analysis as dimensional reduction tools to summarize the SNP data to independent PCs within clusters of loci connected by high LD. We then introduce single- and multilocus models that can efficiently conduct the association tests on such high-dimensional data. The methods can be adapted to different model structures and used to analyse samples collected from the wild or from biparental F2 populations, which are commonly used in ecological genetics mapping studies. We demonstrate the performance of our approaches with two publicly available data sets from a plant (Arabidopsis thaliana) and a fish (Pungitius pungitius), as well as with simulated data.

Inference of genetic architecture from chromosome partitioning analyses is sensitive to genome variation, sample size, heritability and effect size distribution.

Genomewide association studies have contributed immensely to our understanding of the genetic basis of complex traits. One major conclusion arising from these studies is that most traits are controlled by many loci of small effect, confirming the infinitesimal model of quantitative genetics. A popular approach to test for polygenic architecture involves so-called "chromosome partitioning" where phenotypic variance explained by each chromosome is regressed on the size of the chromosome. First developed for humans, this has now been repeatedly used in other species, but there has been no evaluation of the suitability of this method in species that can differ in their genome characteristics such as number and size of chromosomes. Nor has the influence of sample size, heritability of the trait, effect size distribution of loci controlling the trait or the physical distribution of the causal loci in the genome been examined. Using simulated data, we show that these characteristics have major influence on the inferences of the genetic architecture of traits we can infer using chromosome partitioning analyses. In particular, small variation in chromosome size, small sample size, low heritability, a skewed effect size distribution and clustering of loci can lead to a loss of power and consequently altered inference from chromosome partitioning analyses. Future studies employing this approach need to consider and derive an appropriate null model for their study system, taking these parameters into consideration. Our simulation results can provide some guidelines on these matters, but further studies examining a broader parameter space are needed.

Divergent evolution and niche differentiation within the common peatmoss Sphagnum magellanicum.

PREMISE OF THE STUDY: Populations with phenotypic polymorphism in discrete characters may be good models for investigating genome evolution and speciation. Sphagnum magellanicum Brid. is found throughout the northern hemisphere, and despite considerable variation in morphological characters, it is considered one of the least taxonomically controversial peatmoss species. We have observed two main morphs of the species associated with different microhabitats. Here we investigated the genomic and environmental basis of this intraspecific morphological variation. METHODS: We conducted transplant and common garden experiments to test whether the two morphs are genetically differentiated. We then used RAD-sequencing to quantify the genomic divergence between the morphs and approximate Bayesian computation (ABC) to infer the most likely demographic scenario explaining the genome-wide differentiation of the two morphs. KEY RESULTS: We found that genomic differentiation between the two morphs is unexpectedly high and that several of the differentiated morphological characters have a genetic basis. Using simulation approaches, we found support for a scenario of ancient divergence followed by recent secondary contact. CONCLUSIONS: We show that the two morphs represent the two main genetic clusters previously found worldwide. Our results demonstrate that relatively minor morphological differentiation in a presumed phenotypically plastic peatmoss may be associated with massive divergence across the genome.

Controlling for P-value inflation in allele frequency change in experimental evolution and artificial selection experiments.

Experimental evolution studies can be used to explore genomic response to artificial and natural selection. In such studies, loci that display larger allele frequency change than expected by genetic drift alone are assumed to be directly or indirectly associated with traits under selection. However, such studies report surprisingly many loci under selection, suggesting that current tests for allele frequency change may be subject to P-value inflation and hence be anticonservative. One factor known from genomewide association (GWA) studies to cause P-value inflation is population stratification, such as relatedness among individuals. Here, we suggest that by treating presence of an individual in a population after selection as a binary response variable, existing GWA methods can be used to account for relatedness when estimating allele frequency change. We show that accounting for relatedness like this effectively reduces false-positives in tests for allele frequency change in simulated data with varying levels of population structure. However, once relatedness has been accounted for, the power to detect causal loci under selection is low. Finally, we demonstrate the presence of P-value inflation in allele frequency change in empirical data spanning multiple generations from an artificial selection experiment on tarsus length in two free-living populations of house sparrow and correct for this using genomic control. Our results indicate that since allele frequencies in large parts of the genome may change when selection acts on a heritable trait, such selection is likely to have considerable and immediate consequences for the eco-evolutionary dynamics of the affected populations.

Linkage disequilibrium network analysis (LDna) gives a global view of chromosomal inversions, local adaptation and geographic structure.

Recent advances in sequencing allow population-genomic data to be generated for virtually any species. However, approaches to analyse such data lag behind the ability to generate it, particularly in nonmodel species. Linkage disequilibrium (LD, the nonrandom association of alleles from different loci) is a highly sensitive indicator of many evolutionary phenomena including chromosomal inversions, local adaptation and geographical structure. Here, we present linkage disequilibrium network analysis (LDna), which accesses information on LD shared between multiple loci genomewide. In LD networks, vertices represent loci, and connections between vertices represent the LD between them. We analysed such networks in two test cases: a new restriction-site-associated DNA sequence (RAD-seq) data set for Anopheles baimaii, a Southeast Asian malaria vector; and a well-characterized single nucleotide polymorphism (SNP) data set from 21 three-spined stickleback individuals. In each case, we readily identified five distinct LD network clusters (single-outlier clusters, SOCs), each comprising many loci connected by high LD. In A. baimaii, further population-genetic analyses supported the inference that each SOC corresponds to a large inversion, consistent with previous cytological studies. For sticklebacks, we inferred that each SOC was associated with a distinct evolutionary phenomenon: two chromosomal inversions, local adaptation, population-demographic history and geographic structure. LDna is thus a useful exploratory tool, able to give a global overview of LD associated with diverse evolutionary phenomena and identify loci potentially involved. LDna does not require a linkage map or reference genome, so it is applicable to any population-genomic data set, making it especially valuable for nonmodel species.

Mosquito genomics. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes.

Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts.

Parallel evolution of local adaptation and reproductive isolation in the face of gene flow.

Parallel evolution of similar phenotypes provides strong evidence for the operation of natural selection. Where these phenotypes contribute to reproductive isolation, they further support a role for divergent, habitat-associated selection in speciation. However, the observation of pairs of divergent ecotypes currently occupying contrasting habitats in distinct geographical regions is not sufficient to infer parallel origins. Here we show striking parallel phenotypic divergence between populations of the rocky-shore gastropod, Littorina saxatilis, occupying contrasting habitats exposed to either wave action or crab predation. This divergence is associated with barriers to gene exchange but, nevertheless, genetic variation is more strongly structured by geography than by ecotype. Using approximate Bayesian analysis of sequence data and amplified fragment length polymorphism markers, we show that the ecotypes are likely to have arisen in the face of continuous gene flow and that the demographic separation of ecotypes has occurred in parallel at both regional and local scales. Parameter estimates suggest a long delay between colonization of a locality and ecotype formation, perhaps because the postglacial spread of crab populations was slower than the spread of snails. Adaptive differentiation may not be fully genetically independent despite being demographically parallel. These results provide new insight into a major model of ecologically driven speciation.

Molecular characterization and identification of members of the Anopheles subpictus complex in Sri Lanka.

BACKGROUND: Anopheles subpictus sensu lato is a major malaria vector in South and Southeast Asia. Based initially on polytene chromosome inversion polymorphism, and subsequently on morphological characterization, four sibling species A-D were reported from India. The present study uses molecular methods to further characterize and identify sibling species in Sri Lanka. METHODS: Mosquitoes from Sri Lanka were morphologically identified to species and sequenced for the ribosomal internal transcribed spacer-2 (ITS2) and the mitochondrial cytochrome c oxidase subunit-I (COI) genes. These sequences, together with others from GenBank, were used to construct phylogenetic trees and parsimony haplotype networks and to test for genetic population structure. RESULTS: Both ITS2 and COI sequences revealed two divergent clades indicating that the Subpictus complex in Sri Lanka is composed of two genetically distinct species that correspond to species A and species B from India. Phylogenetic analysis showed that species A and species B do not form a monophyletic clade but instead share genetic similarity with Anopheles vagus and Anopheles sundaicus s.l., respectively. An allele specific identification method based on ITS2 variation was developed for the reliable identification of species A and B in Sri Lanka. CONCLUSION: Further multidisciplinary studies are needed to establish the species status of all chromosomal forms in the Subpictus complex. This study emphasizes the difficulties in using morphological characters for species identification in An. subpictus s.l. in Sri Lanka and demonstrates the utility of an allele specific identification method that can be used to characterize the differential bio-ecological traits of species A and B in Sri Lanka.