Genomic approaches to the study of population demography rely on accurate SNP calling and by-proxy the site frequency spectrum (SFS). Two main questions for the design of such studies remain poorly investigated: do reduced genomic sequencing summary statistics reflect that of whole genome, and how do sequencing strategies and derived summary statistics impact demographic inferences? To address those questions, we applied the ddRAD sequencing approach to 254 individuals and whole genome resequencing approach to 35 mountain goat (Oreamnos americanus) individuals across the species range with a known demographic history. We identified SNPs with 5 different variant callers and used ANGSD to estimate the genotype likelihoods (GLs). We tested combinations of SNP filtering by linkage disequilibrium (LD), minor allele frequency (MAF) and the genomic region. We compared the resulting suite of summary statistics reflective of the SFS and quantified the relationship to demographic inferences by estimating the contemporary effective population size (Ne), isolation-by-distance and population structure, FST, and explicit modelling of the demographic history with δaδi. Filtering had a larger effect than sequencing strategy, with the former strongly influencing summary statistics. Estimates of contemporary Ne and isolation-by-distance patterns were largely robust to the choice of sequencing, pipeline, and filtering. Despite the high variance in summary statistics, whole genome and reduced representation approaches were overall similar in supporting a glacial induced vicariance and low Ne in mountain goats. We discuss why whole genome resequencing data is preferable, and reiterate support the use of GLs, in part because it limits user-determined filters.

Nutrition is at the center of interactions between organisms and their environment. During this time of unprecedented anthropogenic environmental change, biological indicators of nutritional changes and their effects on ecosystems are crucial if we are to understand and mitigate these changes. In this study, we developed a qPCR assay to examine nutritional responsiveness of ten genes identified as potential indicators of nutritional state in the freshwater zooplankton Daphnia pulex. We grew animals in six ecologically relevant treatments: nutrient replete, low carbon (food), low phosphorus, low nitrogen, low calcium, and high Cyanobacteria. We measured the growth rate and elemental composition of these animals and extracted their RNA. We then selected ten nutrient sensitive genes, two per limiting nutrient, and two reference genes from an RNA sequencing dataset of Daphnia pulex grown under the same six nutritional treatments. We then designed and validated qPCR primers for the ten indicator and two reference genes. We ran qPCR using these primers on cDNA from our experimental animals and found that the differential expression patterns of these genes could discriminate between our six nutritional states with high levels of accuracy. These results represent a compelling proof of concept for the use of gene-based nutritional biomarkers, paving the way for the use of genetically based nutritional state biomarkers in the study of ecology

The loss of genetic diversity is a challenge many species are facing, and genomics is a potential tool that can inform and prioritize decision making. Caribou populations have experienced significant recent declines throughout Québec, Canada, and some are considered threatened or endangered. We calculated the ancestral and contemporary patterns of genomic diversity of five caribou populations and applied a comparative framework to assess the interplay between demography and genomic diversity. We calculated a caribou specific mutation rate, μ, by extracting orthologous genes from related ungulates. Whole genome re-sequencing was completed on 67 caribou and genotype likelihoods were estimated. We calculated nucleotide diversity, θπ and estimated the coalescent or ancestral Ne, which ranged from 12,030 to 15,513. When compared to the census size, NC, the endangered Gaspésie Mountain caribou population had the highest Ne:NC ratio which is consistent with recent work suggesting high ancestral Ne:NC is of conservation concern. These ratios were highly correlated with genomic signatures (i.e. Tajima’s D) and explicit demographic model parameters. Values of contemporary Ne, estimated from linkage-disequilibrium, ranged from 11 to 162, with Gaspésie having among the highest contemporary Ne:NC ratio. Importantly, conservation genetics theory would predict this population to be of less concern based on this ratio. Of note, F varied only slightly between populations, and runs of homozygosity were not abundant in the genome. Our study highlights how genomic patterns are nuanced and misleading if viewed only through a contemporary lens; a holistic view should integrate ancestral Ne and Tajima’s D into conservation decisions.