RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site‐associated DNA sequencing (RAD‐Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F_ST 0.008; range 0.00–0.018) were concordant with previous microsatellite‐based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age‐1 and age‐2 families of full and half‐siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non‐native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.     

Effects of captivity,diet, and relocation on the gut bacterial communities of white‐footed mice

Microbes can have important impacts on their host's survival. Captive breeding programs for endangered species include periods of captivity that can ultimately have an impact on reintroduction success. No study to date has investigated the impacts of captive diet on the gut microbiota during the relocation process of generalist species. This study simulated a captive breeding program with white‐footed mice (Peromyscus leucopus) to describe the variability in gut microbial community structure and composition during captivity and relocation in their natural habitat, and compared it to wild individuals. Mice born in captivity were fed two different diets, a control with dry standardized pellets and a treatment with nonprocessed components that reflect a version of their wild diet that could be provided in captivity. The mice from the two groups were then relocated to their natural habitat. Relocated mice that had the treatment diet had more phylotypes in common with the wild‐host microbiota than mice under the control diet or mice kept in captivity. These results have broad implications for our understanding of microbial community dynamics and the effects of captivity on reintroduced animals, including the potential impact on the survival of endangered species. This study demonstrates that ex situ conservation actions should consider a more holistic perspective of an animal's biology including its microbes.     

Native drivers of fish life history traits are lost during the invasion process

Rapid adaptation to global change can counter vulnerability of species to population declines and extinction. Theoretically, under such circumstances both genetic variation and phenotypic plasticity can maintain population fitness, but empirical support for this is currently limited. Here, we aim to characterize the role of environmental and genetic diversity, and their prior evolutionary history (via haplogroup profiles) in shaping patterns of life history traits during biological invasion. Data were derived from both genetic and life history traits including a morphological analysis of 29 native and invasive populations of topmouth gudgeon Pseudorasbora parva coupled with climatic variables from each location. General additive models were constructed to explain distribution of somatic growth rate (SGR) data across native and invasive ranges, with model selection performed using Akaike's information criteria. Genetic and environmental drivers that structured the life history of populations in their native range were less influential in their invasive populations. For some vertebrates at least, fitness‐related trait shifts do not seem to be dependent on the level of genetic diversity or haplogroup makeup of the initial introduced propagule, nor of the availability of local environmental conditions being similar to those experienced in their native range. As long as local conditions are not beyond the species physiological threshold, its local establishment and invasive potential are likely to be determined by local drivers, such as density‐dependent effects linked to resource availability or to local biotic resistance.     

Is population subdivision different from speciation? From phylogeography to species delimitation

Species‐level diversity and the underlying mechanisms that lead to the formation of new species, that is, speciation, have often been confounded with intraspecific diversity and population subdivision. The delineation between intraspecific and interspecific divergence processes has received much less attention than species delimitation. The ramifications of confounding speciation and population subdivision are that the term speciation has been used to describe many different biological divergence processes, rendering the results, or inferences, between studies incomparable. Phylogeographic studies have advanced our understanding of how spatial variation in the pattern of biodiversity can begin, become structured, and persist through time. Studies of species delimitation have further provided statistical and model‐based approaches to determine the phylogeographic entities that merit species status. However, without a proper understanding and delineation between the processes that generate and maintain intraspecific and interspecific diversity in a study system, the delimitation of species may still not be biologically and evolutionarily relevant. I argue that variation in the continuity of the divergence process among biological systems could be a key factor leading to the enduring contention in delineating divergence patterns, or species delimitation, meriting future comparative studies to help us better understand the nature of biological species.     

Age‐related variation in the trophic characteristics of a marsupial carnivore,the Tasmanian devil Sarcophilus harrisii

Age‐related changes in diet have implications for competitive interactions and for predator–prey dynamics, affecting individuals and groups at different life stages. To quantify patterns of variation and ontogenetic change in the diets of Tasmanian devils Sarcophilus harrisii, a threatened marsupial carnivore, we analyzed variation in the stable isotope composition of whisker tissue samples taken from 91 individual devils from Wilmot, Tasmania from December 2014 to February 2017. Both δ¹³C and δ¹⁵N decreased with increasing age in weaned Tasmanian devils, indicating that as they age devils rely less on small mammals and birds, and more on large herbivores. Devils <12 months old had broader group isotopic niches, as estimated by Bayesian standard ellipses (SEA_B mode = 1.042) than devils from 12 to 23 months old (mode = 0.541) and devils ≥24 months old (mode = 0.532). Devils <24 months old had broader individual isotopic niches (SEA_B mode range 0.492–1.083) than devils ≥24 months old (mode range 0.092–0.240). A decrease in δ¹⁵N from the older whisker sections to the more recently grown sections in devils <24 months old likely reflects the period of weaning in this species, as this pattern was not observed in devils ≥24 months old. Our data reveal changes in the isotopic composition of devil whiskers with increasing age, accompanied by a reduction in isotopic variation both among population age classes and within individuals, reflecting the effect of weaning in early life, and a likely shift from an initially diverse diet of small mammals, birds, and invertebrates towards increasing consumption of larger herbivores in adulthood.     

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.     

Acknowledging uncertainty in evolutionary reconstructions of ecological niches

Reconstructing ecological niche evolution can provide insight into the biogeography and diversification of evolving lineages. However, comparative phylogenetic methods may infer the history of ecological niche evolution inaccurately because (a) species' niches are often poorly characterized; and (b) phylogenetic comparative methods rely on niche summary statistics rather than full estimates of species' environmental tolerances. Here, we propose a new framework for coding ecological niches and reconstructing their evolution that explicitly acknowledges and incorporates the uncertainty introduced by incomplete niche characterization. Then, we modify existing ancestral state inference methods to leverage full estimates of environmental tolerances. We provide a worked empirical example of our method, investigating ecological niche evolution in the New World orioles (Aves: Passeriformes: Icterus spp.). Temperature and precipitation tolerances were generally broad and conserved among orioles, with niche reduction and specialization limited to a few terminal branches. Tools for performing these reconstructions are available in a new R package called nichevol.     

In silico and empirical evaluation of twelve metabarcoding primer sets for insectivorous diet analyses

During the most recent decade, environmental DNA metabarcoding approaches have been both developed and improved to minimize the biological and technical biases in these protocols. However, challenges remain, notably those relating to primer design. In the current study, we comprehensively assessed the performance of ten COI and two 16S primer pairs for eDNA metabarcoding, including novel and previously published primers. We used a combined approach of in silico, in vivo‐mock community (33 arthropod taxa from 16 orders), and guano‐based analyses to identify primer sets that would maximize arthropod detection and taxonomic identification, successfully identify the predator (bat) species, and minimize the time and financial costs of the experiment. We focused on two insectivorous bat species that live together in mixed colonies: the greater horseshoe bat (Rhinolophus ferrumequinum) and Geoffroy's bat (Myotis emarginatus). We found that primer degeneracy is the main factor that influences arthropod detection in silico and mock community analyses, while amplicon length is critical for the detection of arthropods from degraded DNA samples. Our guano‐based results highlight the importance of detecting and identifying both predator and prey, as guano samples can be contaminated by other insectivorous species. Moreover, we demonstrate that amplifying bat DNA does not reduce the primers' capacity to detect arthropods. We therefore recommend the simultaneous identification of predator and prey. Finally, our results suggest that up to one‐third of prey occurrences may be unreliable and are probably not of primary interest in diet studies, which may decrease the relevance of combining several primer sets instead of using a single efficient one. In conclusion, this study provides a pragmatic framework for eDNA primer selection with respect to scientific and methodological constraints.     

High predictability of direct competition between marine diatoms under different temperatures and nutrient states

The distribution of marine phytoplankton will shift alongside changes in marine environments, leading to altered species frequencies and community composition. An understanding of the response of mixed populations to abiotic changes is required to adequately predict how environmental change may affect the future composition of phytoplankton communities. This study investigated the growth and competitive ability of two marine diatoms, Phaeodactylum tricornutum and Thalassiosira pseudonana, along a temperature gradient (9–35°C) spanning the thermal niches of both species under both high‐nitrogen nutrient‐replete and low‐nitrogen nutrient‐limited conditions. Across this temperature gradient, the competitive outcome under both nutrient conditions at any assay temperature, and the critical temperature at which competitive advantage shifted from one species to the other, was well predicted by the temperature dependencies of the growth rates of the two species measured in monocultures. The temperature at which the competitive advantage switched from P. tricornutum to T. pseudonana increased from 18.8°C under replete conditions to 25.3°C under nutrient‐limited conditions. Thus, P. tricornutum was a better competitor over a wider temperature range in a low N environment. Being able to determine the competitive outcomes from physiological responses of single species to environmental changes has the potential to significantly improve the predictive power of phytoplankton spatial distribution and community composition models.