Reanalysis suggests that genomic islands of speciation are due to reduced diversity,not reduced gene flow

The metaphor of ‘genomic islands of speciation’ was first used to describe heterogeneous differentiation among loci between the genomes of closely related species. The biological model proposed to explain these differences was that the regions showing high levels of differentiation were resistant to gene flow between species, while the remainder of the genome was being homogenized by gene flow and consequently showed lower levels of differentiation. However, the conditions under which such differentiation can occur at multiple unlinked loci are restrictive; additionally, essentially, all previous analyses have been carried out using relative measures of divergence, which can be misleading when regions with different levels of recombination are compared. Here, we test the model of differential gene flow by asking whether absolute divergence is also higher in the previously identified ‘islands’. Using five species pairs for which full sequence data are available, we find that absolute measures of divergence are not higher in genomic islands. Instead, in all cases examined, we find reduced diversity in these regions, a consequence of which is that relative measures of divergence are abnormally high. These data therefore do not support a model of differential gene flow among loci, although islands of relative divergence may represent loci involved in local adaptation. Simulations using the program IMa2 further suggest that inferences of any gene flow may be incorrect in many comparisons. We instead present an alternative explanation for heterogeneous patterns of differentiation, one in which postspeciation selection generates patterns consistent with multiple aspects of the data.     

Evaluating noninvasive genetic sampling techniques to estimate large carnivore abundance

Monitoring large carnivores is difficult because of intrinsically low densities and can be dangerous if physical capture is required. Noninvasive genetic sampling (NGS) is a safe and cost‐effective alternative to physical capture. We evaluated the utility of two NGS methods (scat detection dogs and hair sampling) to obtain genetic samples for abundance estimation of coyotes, black bears and Canada lynx in three areas of Newfoundland, Canada. We calculated abundance estimates using program capwire , compared sampling costs, and the cost/sample for each method relative to species and study site, and performed simulations to determine the sampling intensity necessary to achieve abundance estimates with coefficients of variation (CV) of <10%. Scat sampling was effective for both coyotes and bears and hair snags effectively sampled bears in two of three study sites. Rub pads were ineffective in sampling coyotes and lynx. The precision of abundance estimates was dependent upon the number of captures/individual. Our simulations suggested that ~3.4 captures/individual will result in a < 10% CV for abundance estimates when populations are small (23–39), but fewer captures/individual may be sufficient for larger populations. We found scat sampling was more cost‐effective for sampling multiple species, but suggest that hair sampling may be less expensive at study sites with limited road access for bears. Given the dependence of sampling scheme on species and study site, the optimal sampling scheme is likely to be study‐specific warranting pilot studies in most circumstances.     

In situ genetic differentiation in a Hispaniolan lizard (Ameiva chrysolaema): a multilocus perspective

A previous phylogeographic study of mitochondrial haplotypes for the Hispaniolan lizard Ameiva chrysolaema revealed deep genetic structure associated with seawater inundation during the late Pliocene/early Pleistocene and evidence of subsequent population expansion into formerly inundated areas. We revisit hypotheses generated by our previous study using increased geographic sampling of populations and analysis of three nuclear markers (alpha-enolase intron 8, alpha-cardiac-actin intron 4, and beta-actin intron 3) in addition to mitochondrial haplotypes (ND2). Large genetic discontinuities correspond spatially and temporally with historical barriers to gene flow (sea inundations). NCPA cross-validation analysis and Bayesian multilocus analyses of divergence times (IMa and MCMCcoal) reveal two separate episodes of fragmentation associated with Pliocene and Pleistocene sea inundations, separating the species into historically separate Northern, East-Central, West-Central, and Southern population lineages. Multilocus Bayesian analysis using IMa indicates asymmetrical migration from the East-Central to the West-Central populations following secondary contact, consistent with expectations from the more pervasive sea inundation in the western region. The West-Central lineage has a genetic signature of population growth consistent with the expectation of geographic expansion into formerly inundated areas. Within each lineage, significant spatial genetic structure indicates isolation by distance at comparable temporal scales. This study adds to the growing body of evidence that vicariant speciation may be the prevailing source of lineage accumulation on oceanic islands. Thus, prior theories of island biogeography generally underestimate the role and temporal scale of intra-island vicariant processes.     

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species‐specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri‐specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure ), bidirectional estimates of gene flow (migrate ‐n and IMa2) and phylogeny‐based methods (*beast , buck y) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species‐specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri‐like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.     

Sensitivity of human germ-cell-tumor cell lines to human interferons

We examined the sensitivity of four human germ-cell-tumor cell lines exhibiting different stages of differentiation to human interferons (IFNs) in vitro. The cell lines were derived from two embryonal carcinomas (NEC 8 and NEC 14), a choriocarcinoma (IMa), and a yolk-sac tumor (HUOT). Treatment with poly I:C induced IFN production in IMa and HUOT cells, but not in NEC-8 and NEC-14 cells. In the two embryonal-carcinoma cell lines, the addition of IFN-alpha, IFN-beta, and IFN-gamma did not prevent infection by vesicular stomatitis virus and encephalomyocarditis virus. Also, in these two lines, 2-5A synthetase was not induced by the addition of IFN-alpha. In contrast, both IMa and HUOT showed sensitivity to the antiviral action of IFN-alpha and IFN-beta against the two viruses, and 2-5A synthetase was induced by IFN-alpha. IFNs added at doses of up to 1000 IU/ml had no antiproliferative effect on NEC 8, NEC 14, and HUOT, whereas colony formation by IMa cells was greatly suppressed by all three forms of IFN. These results indicate that the production of and sensitivity to IFN are developmentally regulated and are related to the level of differentiation of human germ-cell stem cells.     

Hierarchical distance sampling to estimate population sizes of common lizards across a desert ecoregion

Multispecies wildlife monitoring across large geographical regions is important for effective conservation planning in response to expected impacts from climate change and land use. Unlike many species of birds, mammals, and amphibians which can be efficiently sampled using automated sensors including cameras and sound recorders, reptiles are often much more challenging to detect, in part because of their typically cryptic behavior and generally small body sizes. Although many lizard species are more active during the day which makes them easier to detect using visual encounter surveys, they may be unavailable for sampling during certain periods of the day or year due to their sensitivity to temperature. In recognition of these sampling challenges, we demonstrate application of a recent innovation in distance sampling that adjusts for temporary emigration between repeat survey visits. We used transect surveys to survey lizards at 229 sites throughout the Mojave Desert in California, USA, 2016. We estimated a total population size of 82 million (90% CI: 65–99 million) for the three most common species of lizards across this 66,830 km² ecoregion. We mapped how density at the 1‐km² scale was predicted to vary with vegetation cover and human development. We validated these results against independent surveys from the southern portion of our study area. Our methods and results demonstrate how multispecies monitoring programs spanning arid ecoregions can better incorporate information about reptiles.     

Estimating occupancy dynamics for large‐scale monitoring networks: amphibian breeding occupancy across protected areas in the northeast United States

Regional monitoring strategies frequently employ a nested sampling design where a finite set of study areas from throughout a region are selected and intensive sampling occurs within a subset of sites within the individual study areas. This sampling protocol naturally lends itself to a hierarchical analysis to account for dependence among subsamples. Implementing such an analysis using a classic likelihood framework is computationally challenging when accounting for detection errors in species occurrence models. Bayesian methods offer an alternative approach for fitting models that readily allows for spatial structure to be incorporated. We demonstrate a general approach for estimating occupancy when data come from a nested sampling design. We analyzed data from a regional monitoring program of wood frogs (Lithobates sylvaticus) and spotted salamanders (Ambystoma maculatum) in vernal pools using static and dynamic occupancy models. We analyzed observations from 2004 to 2013 that were collected within 14 protected areas located throughout the northeast United States. We use the data set to estimate trends in occupancy at both the regional and individual protected area levels. We show that occupancy at the regional level was relatively stable for both species. However, substantial variation occurred among study areas, with some populations declining and some increasing for both species. In addition, When the hierarchical study design is not accounted for, one would conclude stronger support for latitudinal gradient in trends than when using our approach that accounts for the nested design. In contrast to the model that does not account for nesting, the nested model did not include an effect of latitude in the 95% credible interval. These results shed light on the range‐level population status of these pond‐breeding amphibians, and our approach provides a framework that can be used to examine drivers of local and regional occurrence dynamics.     

Interpreting the estimated timing of migration events between hybridizing species

The question of whether speciation can occur in the presence of gene flow has long been a contentious one. However, measuring the amount and timing of gene flow remains challenging. The computer program IMa2 allows researchers to estimate the timing of migration events for each locus during analyses, and these estimates have been used to infer the timing of introgression and mode of speciation. We use simulated data sets to examine the degree to which gene-flow timing estimates can be used for these purposes, and what demographic conditions and data sets may be most amenable to gene-flow timing estimation. We find that the 90% highest posterior density (HPD) interval of gene-flow timing is almost always substantially wider than the actual window of gene flow, and increasing the information content of the data set in terms of number of loci, number of sequences sampled or locus length (and thus number of variable sites) has little impact on the posterior distribution over the range of values we tested. Even when simulated gene flow only occurred over the most recent 0.01% of the species' history, the HPD interval usually encompasses the inferred divergence time. Our results indicate that gene-flow timing estimates made using the method currently implemented in IMa2 cannot reliably be used to make inferences about the timing of introgression between diverged species or to distinguish between speciation with gene flow and allopatric speciation followed by one or more episodes of gene flow.     

How many marker loci are necessary? Analysis of dominant marker data sets using two popular population genetic algorithms

The number of marker loci required to answer a given research question satisfactorily is especially important for dominant markers since they have a lower information content than co‐dominant marker systems. In this study, we used simulated dominant marker data sets to determine the number of dominant marker loci needed to obtain satisfactory results from two popular population genetic analyses: STRUCTURE and AMOVA (analysis of molecular variance). Factors such as migration, level of population differentiation, and unequal sampling were varied in the data sets to mirror a range of realistic research scenarios. AMOVA performed well under all scenarios with a modest quantity of markers while STRUCTURE required a greater number, especially when populations were closely related. The popular ΔK method of determining the number of genetically distinct groups worked well when sampling was balanced, but underestimated the true number of groups with unbalanced sampling. These results provide a window through which to interpret previous work with dominant markers and we provide a protocol for determining the number of markers needed for future dominant marker studies.     

Efficient distinction of invasive aquatic plant species from non‐invasive related species using DNA barcoding

Biological invasions are regarded as threats to global biodiversity. Among invasive aliens, a number of plant species belonging to the genera Myriophyllum, Ludwigia and Cabomba, and to the Hydrocharitaceae family pose a particular ecological threat to water bodies. Therefore, one would try to prevent them from entering a country. However, many related species are commercially traded, and distinguishing invasive from non‐invasive species based on morphology alone is often difficult for plants in a vegetative stage. In this regard, DNA barcoding could become a good alternative. In this study, 242 samples belonging to 26 species from 10 genera of aquatic plants were assessed using the chloroplast loci trnH‐psbA, matK and rbcL. Despite testing a large number of primer sets and several PCR protocols, the matK locus could not be amplified or sequenced reliably and therefore was left out of the analysis. Using the other two loci, eight invasive species could be distinguished from their respective related species, a ninth one failed to produce sequences of sufficient quality. Based on the criteria of universal application, high sequence divergence and level of species discrimination, the trnH‐psbA noncoding spacer was the best performing barcode in the aquatic plant species studied. Thus, DNA barcoding may be helpful with enforcing a ban on trade of such invasive species, such as is already in place in the Netherlands. This will become even more so once DNA barcoding would be turned into machinery routinely operable by a nonspecialist in botany and molecular genetics.     

Genetic structure among remnant populations of a migratory passerine,the Northern Wheatear Oenanthe oenanthe

Continuous animal populations often become fragmented due to anthropogenic habitat alterations. These small, fragmented populations are fragile due to demographic and genetic factors, whereas immigration can enhance their long‐term viability. Previously, we showed that high philopatry affected the local dynamics of three small and remnant subpopulations of Northern Wheatears in The Netherlands. Here, we show that these three populations together with an additional larger population in the European lowlands are highly genetically differentiated based on 22 microsatellite markers. In contrast, we found no evidence for differentiation using two mitochondrial DNA markers. An IMa2 analysis indicates that gene flow has occurred regularly among our sampled populations. As immigration of colour‐ringed birds among our sampled populations is rare at best, our results suggest that the populations have recently become isolated from one another. Low dispersal rates in highly mobile birds may occur when suitable habitat becomes highly fragmented, and will accentuate stochastic demographic processes and inbreeding, both reducing population viability. As dispersal rates are low among populations of Northern Wheatears in The Netherlands, there is only a small probability of recolonization of habitat patches where populations have become locally extinct.     

Simulations inform design of regional occupancy‐based monitoring for a sparsely distributed,territorial species

Sparsely distributed species attract conservation concern, but insufficient information on population trends challenges conservation and funding prioritization. Occupancy‐based monitoring is attractive for these species, but appropriate sampling design and inference depend on particulars of the study system. We employed spatially explicit simulations to identify minimum levels of sampling effort for a regional occupancy monitoring study design, using white‐headed woodpeckers (Picoides albolvartus), a sparsely distributed, territorial species threatened by habitat decline and degradation, as a case study. We compared the original design with commonly proposed alternatives with varying targets of inference (i.e., species range, space use, or abundance) and spatial extent of sampling. Sampling effort needed to achieve adequate power to observe a long‐term population trend (≥80% chance to observe a 2% yearly decline over 20 years) with the previously used study design consisted of annually monitoring ≥120 transects using a single‐survey approach or ≥90 transects surveyed twice per year using a repeat‐survey approach. Designs that shifted inference toward finer‐resolution trends in abundance and extended the spatial extent of sampling by shortening transects, employing a single‐survey approach to monitoring, and incorporating a panel design (33% of units surveyed per year) improved power and reduced error in estimating abundance trends. In contrast, efforts to monitor coarse‐scale trends in species range or space use with repeat surveys provided extremely limited statistical power. Synthesis and applications. Sampling resolutions that approximate home range size, spatially extensive sampling, and designs that target inference of abundance trends rather than range dynamics are probably best suited and most feasible for broad‐scale occupancy‐based monitoring of sparsely distributed territorial animal species.     

Algorithmic single‐locus species delimitation: effects of sampling effort,variation and nonmonophyly in four methods and 1870 species of beetles

The vast number of undescribed species and the fast rate of biodiversity loss call for new approaches to speed up alpha taxonomy. A plethora of methods for delimiting species or operational taxonomic units (OTUs) based on sequence data have been published in recent years. We test the ability of four delimitation methods (BIN, ABGD, GMYC, PTP) to reproduce established species boundaries on a carefully curated DNA barcode data set of 1870 North European beetle species. We also explore how sampling effort, intraspecific variation, nearest neighbour divergence and nonmonophyly affect the OTU delimitations. All methods produced approximately 90% identity between species and OTUs. The effects of variation and sampling differed between methods. ABGD was sensitive to singleton sequences, while GMYC showed tendencies for oversplitting. The best fit between species and OTUs was achieved using simple rules to find consensus between discordant OTU delimitations. Using several approaches simultaneously allows the methods to compensate for each other's weaknesses. Barcode‐based OTU‐picking is an efficient way to delimit putative species from large data sets where the use of more sophisticated methods based on multilocus or genomic data is not feasible.     

Increased divergence but reduced variation on the Z chromosome relative to autosomes in Ficedula flycatchers: differential introgression or the faster-Z effect?

Recent multilocus studies of congeneric birds have shown a pattern of elevated interspecific divergence on the Z chromosome compared to the autosomes. In contrast, intraspecifically, birds exhibit less polymorphism on the Z chromosome relative to the autosomes. We show that the four black-and-white Ficedula flycatcher species show greater genetic divergence on the Z chromosome than on the autosomes, and that the ratios of intraspecific polymorphism at Z-linked versus autosomal markers are below the neutral expectation of 75%. In all species pairs, we found more fixed substitutions and fewer shared polymorphisms on the Z chromosome than on the autosomes. Finally, using isolation with migration (IMa) models we estimated gene flow among the four closely related flycatcher species. The results suggest that different pattern of evolution of Z chromosomes and autosomes is best explained by the faster-Z hypothesis, since the estimated long-term gene flow parameters were close to zero in all comparisons.     

Vocal activity rate index: a useful method to infer terrestrial bird abundance with acoustic monitoring

Autonomous recording units have been widely used in a large number of bird studies in recent years, but challenges remain in estimating abundance based on acoustic monitoring. We tested whether vocal activity rate index (VAR; the number of songs per unit time for a species), recorded using autonomous recording units, was related to population abundance in two terrestrial bird species, the European Bee‐eater Merops apiaster and the Dupont's Lark Chersophilus duponti. We took recordings at sites where censuses were also carried out to estimate local populations around recorders. We found a positive and significant relationship for the two monitored species. Although our results are not conclusive, the strong and significant relationship found for both monitored species suggests that VAR may be used to infer bird abundance around recorders in terrestrial species. We describe five logical steps for using the VAR with autonomous recording units in other species to guide future studies.     

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations

Thanks to genome‐scale diversity data, present‐day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up‐to‐date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self‐fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype–environment correlations and five designed to detect adaptive differentiation (based on F_ST or similar measures). We show that genotype–environment correlation methods have substantially more power to detect selection than differentiation‐based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype–environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.     

A tangled tale of two teal: population history of the grey Anas gracilis and chestnut teal A. castanea of Australia

Two Australian species of teal (Anseriformes: Anatidae: Anas ), the grey teal Anas gracilis and the chestnut teal A. castanea , are remarkable for the zero or near-zero divergence recorded between them in earlier surveys of mitochondrial DNA (mtDNA) diversity. We confirmed this result through wider geographical and population sampling as well as nucleotide sampling in the more rapidly evolving mtDNA control region. Any data set where two species share polymorphism as is the case here can be explained by a model of gene flow through hybridization on one hand or by incomplete lineage sorting on the other hand. Ideally, analysis of such shared polymorphism would simultaneously estimate the likelihood of both phenomena. To do this, we used the underlying principle of the IMa package to explore ramifications to understanding population histories of A. gracilis and A. castanea . We cannot reject that hybridization occurs between the two species but an equally or more plausible finding for their nearly zero divergence is incomplete sorting following very recent divergence between the two, probably in the mid-late Pleistocene. Our data add to studies that explore intermediate stages in the evolution of reciprocal monophyly and paraphyletic or polyphyletic relationships in mtDNA diversity among widespread Australian birds.     

Fine‐scale habitat heterogeneity explains the local distribution of two Amazonian frog species of concern for conservation

To test the hypothesis that subtle differences in abiotic requirements can result in almost total spatial segregation, we sampled two species of diurnal frogs, Atelopus spumarius and Allobates sumtuosus, in a primary forest reserve in central Brazilian Amazonia. We conducted visual and acoustic surveys on three occasions over 2 months, in 40 streamside (riparian) plots distributed throughout the reserve's two major drainage basins, using a grid system that covers 64 km². On average, drainages differed in the pH of stream water and the number of connected and isolated streamside pools. Differences in abiotic characteristics of drainages were associated with the spatial distribution of frog species. The occurrence and density of Allobates sumtuosus was negatively related to stream pH and discharge and positively related to the number of isolated pools in plots. The occurrence and density of Atelopus spumarius was associated with streams with high discharge and pH near neutral. These results indicate that although very large reserves will probably contain sufficient landscape heterogeneity to accommodate different species of diurnal frogs, due to strongly patchy distributions, in situ studies using fine‐scale species‐distribution models will be necessary to assess the adequacy of small reserves in Amazonia that cover hundreds of square kilometers or less for the conservation of some anuran species.     

Five species of passerine bird differ in their ability to detect Batesian mimics

Multiple predators affect the evolution of aposematic signals in nature and these predators may substantially differ in terms of ecological and cognitive parameters. However, most experimental studies testing the evolution of Batesian mimics use only a single species of predator (usually the great tit or a domestic chick). Therefore, in the present study, we experimentally tested the responses of five passerine predators to an artificially made Batesian mimic (a cockroach equipped with the warning pattern of the red firebug) with respect to their dietary ecology. Half of the individuals of each species were fed on unmodified roaches before the experiment, whereas the other half were fed with mealworms and thus had no previous experience with roaches. We found that Batesian mimics were better protected than inconspicuous prey against inexperienced great tits and robins alone. The other three bird species showed high level of neophobia; therefore, the effect of warning coloration could not be assessed. We also found that experienced birds attacked a greater number of Batesian mimics compared to inexperienced individuals of all tested species, with the exception of blackcaps. In the great tits, robins, and blue tits, a significant number of experienced birds attacked the Batesian mimic, which was possibly the result of a learned search image for a roach. Our results suggest that using a limited array of predators to describe evolutionary processes forming the diversity of antipredatory strategies of the prey may be biased and need not describe the situation occurring in nature.     

Evaluating the interaction of faecal pellet deposition rates and DNA degradation rates to optimize sampling design for DNA‐based mark–recapture analysis of Sonoran pronghorn

Knowledge of population demographics is important for species management but can be challenging in low‐density, wide‐ranging species. Population monitoring of the endangered Sonoran pronghorn (Antilocapra americana sonoriensis) is critical for assessing the success of recovery efforts, and noninvasive DNA sampling (NDS) could be more cost‐effective and less intrusive than traditional methods. We evaluated faecal pellet deposition rates and faecal DNA degradation rates to maximize sampling efficiency for DNA‐based mark–recapture analyses. Deposition data were collected at five watering holes using sampling intervals of 1–7 days and averaged one pellet pile per pronghorn per day. To evaluate nuclear DNA (nDNA) degradation, 20 faecal samples were exposed to local environmental conditions and sampled at eight time points from one to 124 days. Average amplification success rates for six nDNA microsatellite loci were 81% for samples on day one, 63% by day seven, 2% by day 14 and 0% by day 60. We evaluated the efficiency of different sampling intervals (1–10 days) by estimating the number of successful samples, success rate of individual identification and laboratory costs per successful sample. Cost per successful sample increased and success and efficiency declined as the sampling interval increased. Results indicate NDS of faecal pellets is a feasible method for individual identification, population estimation and demographic monitoring of Sonoran pronghorn. We recommend collecting samples >7 days old and estimate that a sampling interval of 4–7 days in summer conditions (i.e. extreme heat and exposure to UV light) will achieve desired sample sizes for mark–recapture analysis while also maximizing efficiency.