Confidently identifying the correct K value using the ΔK method: When does K = 2?

Populations delineated based on genetic data are commonly used for wildlife conservation and management. Many studies use the program structure combined with the ΔK method to identify the most probable number of populations (K). We recently found K = 2 was identified more often when studies used ΔK compared to studies that did not. We suggested two reasons for this: hierarchical population structure leads to underestimation, or the ΔK method does not evaluate K = 1 causing an overestimation. The present contribution aims to develop a better understanding of the limits of the method using one, two and three population simulations across migration scenarios. From these simulations we identified the “best K” using model likelihood and ΔK. Our findings show that mean probability plots and ΔK are unable to resolve the correct number of populations once migration rate exceeds 0.005. We also found a strong bias towards selecting K = 2 using the ΔK method. We used these data to identify the range of values where the ΔK statistic identifies a value of K that is not well supported. Finally, using the simulations and a review of empirical data, we found that the magnitude of ΔK corresponds to the level of divergence between populations. Based on our findings, we suggest researchers should use the ΔK method cautiously; they need to report all relevant data, including the magnitude of ΔK, and an estimate of connectivity for the research community to assess whether meaningful genetic structure exists within the context of management and conservation.     

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.     

Assessing changes in genomic divergence following a century of human‐mediated secondary contact among wild and captive‐bred ducks

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human‐mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD‐seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game‐farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game‐farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait‐capture array targeting thousands of loci in century‐old (1842–1915) and contemporary (2009–2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game‐farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game‐farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens.     

Genetic variation in Plethodon cinereus and Plethodon hubrichti from in and around a contact zone

Climate change poses several challenges to biological communities including changes in the frequency of encounters between closely related congeners as a result of range shifts. When climate change leads to increased hybridization, hybrid dysfunction or genetic swamping may increase extinction risk—particularly in range‐restricted species with low vagility. The Peaks of Otter Salamander, Plethodon hubrichti, is a fully terrestrial woodland salamander that is restricted to ~18 km of ridgeline in the mountains of southwestern Virginia, and its range is surrounded by the abundant and widespread Eastern Red‐backed Salamander, Plethodon cinereus. In order to determine whether these two species are hybridizing and how their range limits may be shifting, we assessed variation at eight microsatellite loci and a 1,008 bp region of Cytochrome B in both species at allopatric reference sites and within a contact zone. Our results show that hybridization between P. hubrichti and P. cinereus either does not occur or is very rare. However, we find that diversity and differentiation are substantially higher in the mountaintop endemic P. hubrichti than in the widespread P. cinereus, despite similar movement ability for the two species as assessed by a homing experiment. Furthermore, estimation of divergence times between reference and contact zone populations via approximate Bayesian computation is consistent with the idea that P. cinereus has expanded into the range of P. hubrichti. Given the apparent recent colonization of the contact zone by P. cinereus, future monitoring of P. cinereus range limits should be a priority for the management of P. hubrichti populations.     

Timing,frequency, and duration of incubation recesses in dabbling ducks

Nest attendance is an important determinant of avian reproductive success, and identifying factors that influence the frequency and duration of incubation recesses furthers our understanding of how incubating birds balance their needs with those of their offspring. We characterized the frequency and timing (start time, end time, and duration) of incubation recesses for mallard (Anas platyrhynchos) and gadwall (Mareca strepera) hens breeding in Suisun Marsh, California, USA, and examined the influences of day of year, ambient temperature at the nest, incubation day, and clutch size on recess frequency and timing using linear mixed models. Mallard, on average, took more recesses per day (1.69 ± 0.80, mean ± standard deviation) than did gadwall (1.39 ± 0.69), and 45% of mallard nest‐days were characterized by two recesses, while only 27% of gadwall nest‐days were characterized by two recesses. Mallard morning recesses started at 06:14 ± 02:46 and lasted 106.11 ± 2.01 min, whereas mallard afternoon recesses started at 16:39 ± 02:11 and lasted 155.39 ± 1.99 min. Gadwall morning recesses started at 06:30 ± 02:46 and lasted 91.28 ± 2.32 min, and gadwall afternoon recesses started at 16:31 ± 01:57 and lasted 192.69 ± 1.89 min. Mallard and gadwall started recesses earlier in the day with increasing ambient temperature, but later in the day as the season progressed. Recess duration decreased as the season progressed and as clutch size increased, and increased with ambient temperature at the nest. The impending darkness of sunset appeared to be a strong cue for ending a recess and returning to the nest, because hens returned to their nests earlier than expected when recesses were expected to end after sunset. Within hens, the timing of incubation recesses was repeatable across incubation days and was most repeatable for mallard afternoon recesses and on days in which hens took only one recess. Hens were most likely to be away from nests between 04:00 and 07:00 and between 16:00 and 19:00; therefore, investigators should search for nests between 07:00 and 16:00. Our analyses identified important factors influencing incubation recess timing in dabbling ducks and have important implications for nest monitoring programs.     

Contrasting the seasonal and elevational prevalence of generalist avian haemosporidia in co‐occurring host species

Understanding the ecology and evolution of parasites is contingent on identifying the selection pressures they face across their infection landscape. Such a task is made challenging by the fact that these pressures will likely vary across time and space, as a result of seasonal and geographical differences in host susceptibility or transmission opportunities. Avian haemosporidian blood parasites are capable of infecting multiple co‐occurring hosts within their ranges, yet whether their distribution across time and space varies similarly in their different host species remains unclear. Here, we applied a new PCR method to detect avian haemosporidia (genera Haemoproteus, Leucocytozoon, and Plasmodium) and to determine parasite prevalence in two closely related and co‐occurring host species, blue tits (Cyanistes caeruleus, N = 529) and great tits (Parus major, N = 443). Our samples were collected between autumn and spring, along an elevational gradient in the French Pyrenees and over a three‐year period. Most parasites were found to infect both host species, and while these generalist parasites displayed similar elevational patterns of prevalence in the two host species, this was not always the case for seasonal prevalence patterns. For example, Leucocytozoon group A parasites showed inverse seasonal prevalence when comparing between the two host species, being highest in winter and spring in blue tits but higher in autumn in great tits. While Plasmodium relictum prevalence was overall lower in spring relative to winter or autumn in both species, spring prevalence was also lower in blue tits than in great tits. Together, these results reveal how generalist parasites can exhibit host‐specific epidemiology, which is likely to complicate predictions of host–parasite co‐evolution.     

Recognizing pulses of extinction from clusters of last occurrences

The distribution of last occurrences of fossil taxa in a stratigraphic column are used to infer the pattern, timing and tempo of extinction from the fossil record. Clusters of last occurrences are commonly interpreted as an abrupt pulse of extinction. However, stratigraphic architecture alone can produce clusters of last occurrences that can be misinterpreted as an extinction pulse. These clusters will typically occur in strata that immediately underlie facies changes and sequence‐stratigraphic surfaces. It has been proposed that a basin‐wide analysis of the fossil record within a sequence‐stratigraphic framework can be used to distinguish between clusters of last occurrences caused solely by extinction pulses from those generated by sequence‐stratigraphic architecture. A basin‐wide approach makes it possible to observe lateral facies shifts in response to sea‐level change, mitigating the effects of stratigraphic architecture. Using computer simulations of plausible Late Ordovician mass‐extinction scenarios tuned to an inferred Late Ordovician sea‐level curve, we show that stratigraphically‐generated clusters of last occurrences are observed even in basin‐wide analyses of the simulated fossil records due to the basin‐wide loss of preferred facies for many taxa. Nonetheless, we demonstrate that by coarsening the stratigraphic resolution to the systems‐tract level and identifying facies preferences of simulated taxa, we can filter out taxa whose last occurrences coincide with the basin‐wide loss of their preferred facies. This enables consistent identification of the underlying extinction pattern for a wide variety of extinction scenarios. Applying this approach to empirical field data can help to constrain underlying extinction patterns from the fossil record.     

Insufficient potassium and sulfur supply threaten the productivity of perennial forage grasses in smallholder farms on tropical sandy soils

Plant and Soil - Perennial forage grass production has the potential to improve smallholder livelihoods in the tropics. However, nutrient management is often challenging, especially on infertile...