Geography is more important than host plant use for the population genetic structure of a generalist insect herbivore

Population divergence can occur due to mechanisms associated with geographic isolation and/or due to selection associated with different ecological niches. Much of the evidence for selection‐driven speciation has come from studies of specialist insect herbivores that use different host plant species; however, the influence of host plant use on population divergence of generalist herbivores remains poorly understood. We tested how diet breadth, host plant species and geographic distance influence population divergence of the fall webworm (Hyphantria cunea; FW). FW is a broadly distributed, extreme generalist herbivore consisting of two morphotypes that have been argued to represent two different species: black‐headed and red‐headed. We characterized the differentiation of FW populations at two geographic scales. We first analysed the influence of host plant and geographic distance on genetic divergence across a broad continental scale for both colour types. We further analysed the influence of host plant, diet breadth and geographic distance on divergence at a finer geographic scale focusing on red‐headed FW in Colorado. We found clear genetic and morphological distinction between red‐ and black‐headed FW, and Colorado FW formed a genetic cluster distinct from other locations. Although both geographic distance and host plant use were correlated with genetic distance, geographic distance accounted for up to 3× more variation in genetic distance than did host plant use. As a rare study investigating the genetic structure of a widespread generalist herbivore over a broad geographic range (up to 3,000 km), our study supports a strong role for geographic isolation in divergence in this system.     

Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (F_ST ~0.1), over a short coastline (<60 km) in the North Sea‐Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post‐glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.     

Hybrid zones and the speciation continuum in Heliconius butterflies

Tropical butterflies in the genus Heliconius have long been models in the study of the stages of speciation. Heliconius are unpalatable to predators, and many species are notable for multiple geographic populations with striking warning colour pattern differences associated with Müllerian mimicry. A speciation continuum is evident in Heliconius hybrid zones. Examples range from hybrid zones across which (a) there is little genetic differentiation other than at mimicry loci, but where hybrids are common, (b) to ‘bimodal‘ hybrid zones with strong genetic divergence and few hybrids, (c) through to ‘good’ sympatric species, with hybridization extremely rare or absent. Now, in this issue of Molecular Ecology, Arias et al. ( 2012 ) have found an intermediate case in Colombian Heliconius cydno showing evidence for assortative mating and molecular differences, but where hybrids are abundant.     

Strong population structure but no equilibrium yet: Genetic connectivity and phylogeography in the kelp Saccharina latissima (Laminariales,Phaeophyta)

Kelp aquaculture is globally developing steadily as human food source, along with other applications. One of the newer crop species is Saccharina latissima, a northern hemisphere kelp inhabiting temperate to arctic rocky shores. To protect and document its natural genetic variation at the onset of this novel aquaculture, as well as increase knowledge on its taxonomy and phylogeography, we collected new genetic data, both nuclear and mitochondrial, and combined it with previous knowledge to estimate genetic connectivity and infer colonization history. Isolation‐with‐migration coalescent analyses demonstrate that gene flow among the sampled locations is virtually nonexistent. An updated scenario for the origin and colonization history of S. latissima is developed as follows: We propose that the species (or species complex) originated in the northwest Pacific, crossed to the northeast Pacific in the Miocene, and then crossed the Bering Strait after its opening ~5.5 Ma into the Arctic and northeast Atlantic. It subsequently crossed the Atlantic from east to west. During the Pleistocene, it was compressed in the south with evidence for northern refugia in Europe. Postglacial recolonization led to secondary contact in the Canadian Arctic. Saccharina cichorioides is shown to probably belong to the S. latissima species complex and to derive from ancestral populations in the Asian North Pacific. Our novel approach of comparing inferred gene flow based on coalescent analysis versus Wright's island model suggests that equilibrium levels of differentiation have not yet been reached in Europe and, hence, that genetic differentiation is expected to increase further if populations are left undisturbed.     

Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation‐with‐gene‐flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence‐with‐gene‐flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F_ST outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between‐population recombination and gene exchange are reduced over large ‘divergence hitchhiking’ (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F_ST outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation‐with‐gene‐flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.     

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

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.     

Genetic evidence for the origin of Aedes aegypti,the yellow fever mosquito,in the southwestern Indian Ocean

Aedes aegypti is among the best‐studied mosquitoes due to its critical role as a vector of human pathogens and ease of laboratory rearing. Until now, this species was thought to have originated in continental Africa, and subsequently colonized much of the world following the establishment of global trade routes. However, populations of this mosquito on the islands in the southwestern Indian Ocean (SWIO), where the species occurs with its nearest relatives referred to as the Aegypti Group, have received little study. We re‐evaluated the evolutionary history of Ae. aegypti and these relatives, using three data sets: nucleotide sequence data, 18,489 SNPs and 12 microsatellites. We found that: (a) the Aegypti Group diverged 16 MYA (95% HPD: 7–28 MYA) from its nearest African/Asian ancestor; (b) SWIO populations of Ae. aegypti are basal to continental African populations; (c) after diverging 7 MYA (95% HPD: 4–15 MYA) from its nearest formally described relative (Ae. mascarensis), Ae. aegypti moved to continental Africa less than 85,000 years ago, where it recently (<1,000 years ago) split into two recognized subspecies Ae. aegypti formosus and a human commensal, Ae. aegypti aegypti; (d) the Madagascar samples form a clade more distant from all other Ae. aegypti than the named species Ae. mascarensis, implying that Madagascar may harbour a new cryptic species; and (e) there is evidence of introgression between Ae. mascarensis and Ae. aegypti on Réunion, and between the two subspecies elsewhere in the SWIO, a likely consequence of recent introductions of domestic Ae. aegypti aegypti from Asia.     

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within‐population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in‐depth examples of how molecular ecology has been successfully applied to specific amphibian systems.     

Climate change and the spread of vector‐borne diseases: using approximate Bayesian computation to compare invasion scenarios for the bluetongue virus vector Culicoides imicola in Italy

Bluetongue (BT) is a commonly cited example of a disease with a distribution believed to have recently expanded in response to global warming. The BT virus is transmitted to ruminants by biting midges of the genus Culicoides, and it has been hypothesized that the emergence of BT in Mediterranean Europe during the last two decades is a consequence of the recent colonization of the region by Culicoides imicola and linked to climate change. To better understand the mechanism responsible for the northward spread of BT, we tested the hypothesis of a recent colonization of Italy by C. imicola, by obtaining samples from more than 60 localities across Italy, Corsica, Southern France, and Northern Africa (the hypothesized source point for the recent invasion of C. imicola), and by genotyping them with 10 newly identified microsatellite loci. The patterns of genetic variation within and among the sampled populations were characterized and used in a rigorous approximate Bayesian computation framework to compare three competing historical hypotheses related to the arrival and establishment of C. imicola in Italy. The hypothesis of an ancient presence of the insect vector was strongly favoured by this analysis, with an associated P ≥ 99%, suggesting that causes other than the northward range expansion of C. imicola may have supported the emergence of BT in southern Europe. Overall, this study illustrates the potential of molecular genetic markers for exploring the assumed link between climate change and the spread of diseases.     

Genome‐wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.     

Repetitive sequences in phylogeographic inference: a reply to Saltonstall and Lambertini (2012)

In 2011, Vachon and Freeland presented the results of a study that illustrated the importance of considering mutation patterns, including the potential for homoplasy, when deciding whether to include or exclude repetitive sequences from phylogeographic inferences. Saltonstall and Lambertini (2012) criticized this study by suggesting that some of the analyses and interpretations were flawed. In this reply, we explain why we disagree with most of their criticisms and identify some inconsistencies in their analyses. Most importantly, we reiterate the need to examine underlying assumptions, particularly with respect to mutation patterns, when using molecular genetic data to untangle evolutionary relationships.     

Phylogeography of ninespine sticklebacks (Pungitius pungitius) in North America: glacial refugia and the origins of adaptive traits

The current geographical distribution of the ninespine stickleback (Pungitius pungitius) was shaped in large part by the glaciation events of the Pleistocene epoch (2.6 Mya–10 Kya). Previous efforts to elucidate the phylogeographical history of the ninespine stickleback in North America have focused on a limited set of morphological traits, some of which are likely subject to widespread convergent evolution, thereby potentially obscuring relationships among populations. In this study, we used genetic information from both mitochondrial DNA (mtDNA) sequences and nuclear microsatellite markers to determine the phylogenetic relationships among ninespine stickleback populations. We found that ninespine sticklebacks in North America probably dispersed from at least three glacial refugia—the Mississippi, Bering, and Atlantic refugia—not two as previously thought. However, by applying a molecular clock to our mtDNA data, we found that these three groups diverged long before the most recent glacial period. Our new phylogeny serves as a critical framework for examining the evolution of derived traits in this species, including adaptive phenotypes that evolved multiple times in different lineages. In particular, we inferred that loss of the pelvic (hind fin) skeleton probably evolved independently in populations descended from each of the three putative North American refugia.     

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad‐scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.     

Molecular population genetics of the melanic plumage polymorphism in Arctic skuas (Stercorarius parasiticus): evidence for divergent selection on plumage colour

The Arctic skua (Stercorarius parasiticus) is a classic example of an avian plumage polymorphism, with variation in melanin‐based ventral plumage coloration defining pale, intermediate and dark morphs in adults of both sexes. However, despite several decades of field research, there is an incomplete understanding of how the polymorphism in ventral plumage colour is maintained and the selective forces involved. Here, we investigate selection on a locus (MC1R) that is strongly associated with plumage colour variation in Arctic skuas using patterns of nucleotide variation and comparison to neutral loci (nuclear introns and mtDNA). We find that three linked nonsynonymous mutations in MC1R, including the single mutation described previously, are associated with plumage colour in the Arctic skua. The position of nonsynonymous mutations on a MC1R haplotype network implies that divergent selection drove the initial evolution of the colour morphs. Comparisons of F_STs of MC1R vs. nuclear introns among five skua populations differing in proportion of dark morphs along an approximate north–south cline reveal a signature of divergent selection on MC1R. In contrast, we find limited evidence for balancing selection on MC1R within populations, although the power is low. Our results provide strong evidence for both past and ongoing selection on MC1R, and, by implication, plumage colour in Arctic skuas. The results suggest that a fruitful avenue for future ecological studies will be analysis of selection on morphs in colonies at the extremes along the morph ratio cline.     

Asymmetric hybridization and gene flow between Joshua trees (Agavaceae: Yucca) reflect differences in pollinator host specificity

The angiosperms are by far the largest group of terrestrial plants. Their spectacular diversity is often attributed to specialized pollination. Obligate pollination mutualisms where both a plant and its pollinator are dependent upon one another for reproduction are thought to be prone to rapid diversification through co‐evolution and pollinator isolation. However, few studies have evaluated the degree to which pollinators actually mediate reproductive isolation in these systems. Here, we examine evidence for hybridization and gene flow between two subspecies of Joshua tree (Yucca brevifolia brevifolia and Yucca brevifolia jaegeriana) pollinated by two sister species of yucca moth. Previous work indicated that the pollinators differ in host specificity, and DNA sequence data suggested asymmetric introgression between the tree subspecies. Through intensive sampling in a zone of sympatry, a large number of morphologically intermediate trees were identified. These included trees with floral characters typical of Y. b. jaegeriana, but vegetative features typical of Y. b. brevifolia. The opposite combination—Y. b. brevifolia flowers with Y. b. jaegeriana vegetative morphology—never occurred. Microsatellite genotyping revealed a high frequency of genetically admixed, hybrid trees. Coalescent‐based estimates of migration indicated significant gene flow between the subspecies and that the direction of gene flow matches differences in pollinator host fidelity. The data suggest that pollinator behaviour determines the magnitude and direction of gene flow between the two subspecies, but that specialized pollination alone is not sufficient to maintain species boundaries. Natural selection may be required to maintain phenotypic differences in the face of ongoing gene flow.     

landgenreport: a new r function to simplify landscape genetic analysis using resistance surface layers

We describe functions recently added to the r package popgenreport that can be used to perform a landscape genetic analysis (LGA) based on landscape resistance surfaces, which aims to detect the effect of landscape features on gene flow. These functions for the first time implement a LGA in a single framework. Although the approach has been shown to be a valuable tool to study gene flow in landscapes, it has not been widely used to date, despite the type of data being widely available. In part, this is likely due to the necessity to use several software packages to perform landscape genetic analyses. To apply LGA functions, two types of data sets are required: a data set with spatially referenced and genotyped individuals, and a resistance layer representing the effect of the landscape. The function outputs three pairwise distance matrices from these data: a genetic distance matrix, a cost distance matrix and a Euclidean distance matrix. Statistical tests are performed to test whether the cost matrix contributes to the understanding of the observed population structure. A full report on the analysis and outputs in the form of plots and tables of all intermediate steps of the LGA is produced. It is possible to customize the LGA to allow for different cost path approaches and measures of genetic distances. The package is written in the r language and is available through the Comprehensive r Archive. Comprehensive tutorials and information on how to install and use the package are provided at the authors’ website ( www.popgenreport.org ).     

Four cleaved amplified polymorphic sequence (CAPS) markers for the detection of the Juglans ailantifolia chloroplast in putatively native J. cinerea populations

Hybridization between butternut (Juglans cinerea), a forest tree native to eastern North America, and Japanese walnut (J. ailantifolia), a tree tolerant to the lethal fungal disease butternut canker, casts doubt on the genetic identity of the remaining butternuts. We report a diagnostic test to distinguish the J. cinerea chloroplast from the J. ailantifolia chloroplast using cleaved amplified polymorphic sequences resolvable in 1.5% agarose gels. J. ailantifolia maternal ancestry in naturally regenerated stands provides a site selection criterion for studies of introgression dynamics when the non-native parent and the hybrids tolerate a disease to which the native species is susceptible.     

Mature male parr contribution to the effective size of an anadromous Atlantic salmon (Salmo salar) population over 30 years

We describe temporal changes in the genetic composition of a small anadromous Atlantic salmon (Salmo salar) population from South Newfoundland, an area where salmon populations are considered threatened (COSEWIC 2010). We examined the genetic variability (13 microsatellite loci) in 869 out‐migrating smolt and post‐spawning kelt samples, collected from 1985 to 2011 for a total of 22 annual collections and a 30 year span of assigned cohorts. We estimated the annual effective number of breeders (N_b) and the generational effective population size (N_e) through genetic methods and demographically using the adult sex ratio. Comparisons between genetic and demographic estimates show that the adult spawners inadequately explain the observed N_e estimates, suggesting that mature male parr are significantly increasing N_b and N_e over the study period. Spawning as parr appears to be a viable and important strategy in the near absence of adult males.     

Gene flow and selection interact to promote adaptive divergence in regions of low recombination

Adaptation to new environments often occurs in the face of gene flow. Under these conditions, gene flow and recombination can impede adaptation by breaking down linkage disequilibrium between locally adapted alleles. Theory predicts that this decay can be halted or slowed if adaptive alleles are tightly linked in regions of low recombination, potentially favouring divergence and adaptive evolution in these regions over others. Here, we compiled a global genomic data set of over 1,300 individual threespine stickleback from 52 populations and compared the tendency for adaptive alleles to occur in regions of low recombination between populations that diverged with or without gene flow. In support of theory, we found that putatively adaptive alleles (F_ST and d_XY outliers) tend to occur more often in regions of low recombination in populations where divergent selection and gene flow have jointly occurred. This result remained significant when we employed different genomic window sizes, controlled for the effects of mutation rate and gene density, controlled for overall genetic differentiation, varied the genetic map used to estimate recombination and used a continuous (rather than discrete) measure of geographic distance as proxy for gene flow/shared ancestry. We argue that our study provides the first statistical evidence that the interaction of gene flow and selection biases divergence toward regions of low recombination.