Search

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil's geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex.

Emerging infectious diseases (EIDs) not only cause catastrophic declines in wildlife populations but also generate selective pressures that may result in rapid evolutionary responses. One such EID is devil facial tumour disease (DFTD) in the Tasmanian devil. DFTD is almost always fatal and has reduced the average lifespan of individuals by around 2 years, likely causing strong selection for traits that reduce susceptibility to the disease, but population decline has also left Tasmanian devils vulnerable to inbreeding depression. We analysed 22 years of data from an ongoing study of a population of Tasmanian devils on Freycinet Peninsula, Tasmania, to (1) identify whether DFTD may be causing selection on body size, by estimating phenotypic and genetic correlations between DFTD and size traits, (2) estimate the additive genetic variance of susceptibility to DFTD, and (3) investigate whether size traits or susceptibility to DFTD were under inbreeding depression. We found a positive phenotypic relationship between head width and susceptibility to DFTD, but this was not underpinned by a genetic correlation. Conversely, we found a negative phenotypic relationship between body weight and susceptibility to DFTD, and there was evidence for a negative genetic correlation between susceptibility to DFTD and body weight. There was additive genetic variance in susceptibility to DFTD, head width and body weight, but there was no evidence for inbreeding depression in any of these traits. These results suggest that Tasmanian devils have the potential to respond adaptively to DFTD, although the realised evolutionary response will critically further depend on the evolution of DFTD itself. [ABSTRACT FROM AUTHOR]

Identifying the genetic architecture of complex phenotypes is a central goal of modern biology, particularly for disease-related traits. Genome-wide association methods are a classical approach for identifying the genomic basis of variation in disease phenotypes, but such analyses are particularly challenging in natural populations due to sample size difficulties. Extensive mark-recapture data, strong linkage disequilibrium and a lethal transmissible cancer make the Tasmanian devil (Sarcophilus harrisii) an ideal model for such an association study. We used a RAD-capture approach to genotype 624 devils at ~16,000 loci and then used association analyses to assess the heritability of three cancer-related phenotypes: infection case-control (where cases were infected devils and controls were devils that were never infected), age of first infection and survival following infection. The SNP array explained much of the phenotypic variance for female survival (>80%) and female case-control (>61%). We found that a few large-effect SNPs explained much of the variance for female survival (~5 SNPs explained >61% of the total variance), whereas more SNPs (~56) of smaller effect explained less of the variance for female case-control (~23% of the total variance). By contrast, these same SNPs did not account for a significant proportion of phenotypic variance in males, suggesting that the genetic bases of these traits and/or selection differ across sexes. Loci involved with cell adhesion and cell-cycle regulation underlay trait variation, suggesting that the devil immune system is rapidly evolving to recognize and potentially suppress cancer growth through these pathways. Overall, our study provided necessary data for genomics-based conservation and management in Tasmanian devils.

The threatened eastern wolf is found predominantly in protected areas of central Ontario and has an evolutionary history obscured by interbreeding with coyotes and gray wolves, which challenges its conservation status and subsequent management. Here, we used a population genomics approach to uncover spatial patterns of variation in 281 canids in central Ontario and the Great Lakes region. This represents the first genome-wide single nucleotide polymorphism (SNP) dataset with substantial sample sizes of representative populations. Although they comprise their own genetic cluster, we found evidence of eastern wolf dispersal outside of the boundaries of protected areas, in that the frequency of eastern wolf genetic variation decreases with increasing distance from provincial parks. We detected eastern wolf alleles in admixed coyotes along the northeastern regions of Lake Huron and Lake Ontario. Our analyses confirm the unique genomic composition of eastern wolves, which are mostly restricted to small fragmented patches of protected habitat in central Ontario. We hope this work will encourage an innovative discussion regarding a plan for managed introgression, which could conserve eastern wolf genetic material in any genome regardless of their potential mosaic ancestry composition and the habitats that promote them.

Adaptive evolution can facilitate species' range expansions across environmentally heterogeneous landscapes. However, serial founder effects can limit the efficacy of selection, and the evolution of increased dispersal during range expansions may result in gene flow swamping local adaptation. Here, we study how genetic drift, gene flow and selection interact during the cane toad's (Rhinella marina) invasion across the heterogeneous landscape of Australia. Following its introduction in 1935, the cane toad colonised eastern Australia and established several stable range edges. The ongoing, more rapid range expansion in north‐central Australia has occurred concomitant with an evolved increase in dispersal capacity. Using reduced representation genomic data of Australian cane toads from the expansion front and from two areas of their established range, we test the hypothesis that high gene flow constrains local adaptation at the expansion front relative to established areas. Genetic analyses indicate the three study areas are genetically distinct but show similar levels of allelic richness, heterozygosity and inbreeding. Markedly higher gene flow or recency of colonisation at the expansion front have likely hindered local adaptation at the time of sampling, as indicated by reduced slopes of genetic‐environment associations (GEAs) estimated using a novel application of geographically weighted regression that accounts for allele surfing; GEA slopes are significantly steeper in established parts of the range. Our work bolsters evidence supporting adaptation of invasive species post‐introduction and adds novel evidence for differing strengths of evolutionary forces among geographic areas with different invasion histories. [ABSTRACT FROM AUTHOR]

High-throughput techniques based on restriction site-associated DNA sequencing (RADseq) are enabling the low-cost discovery and genotyping of thousands of genetic markers for any species, including non-model organisms, which is revolutionizing ecological, evolutionary and conservation genetics. Technical differences among these methods lead to important considerations for all steps of genomics studies, from the specific scientific questions that can be addressed, and the costs of library preparation and sequencing, to the types of bias and error inherent in the resulting data. In this Review, we provide a comprehensive discussion of RADseq methods to aid researchers in choosing among the many different approaches and avoiding erroneous scientific conclusions from RADseq data, a problem that has plagued other genetic marker types in the past.

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Minimally invasive samples are often the best option for collecting genetic material from species of conservation concern, but they perform poorly in many genomic sequencing methods due to their tendency to yield low DNA quality and quantity. Genotyping‐in‐thousands by sequencing (GT‐seq) is a powerful amplicon sequencing method that can genotype large numbers of variable‐quality samples at a standardized set of single nucleotide polymorphism (SNP) loci. Here, we develop, optimize, and validate a GT‐seq panel for the federally threatened northern Idaho ground squirrel (Urocitellus brunneus) to provide a standardized approach for future genetic monitoring and assessment of recovery goals using minimally invasive samples. The optimized panel consists of 224 neutral and 81 putatively adaptive SNPs. DNA collected from buccal swabs from 2016 to 2020 had 73% genotyping success, while samples collected from hair from 2002 to 2006 had little to no DNA remaining and did not genotype successfully. We evaluated our GT‐seq panel by measuring genotype discordance rates compared to RADseq and whole‐genome sequencing. GT‐seq and other sequencing methods had similar population diversity and FST estimates, but GT‐seq consistently called more heterozygotes than expected, resulting in negative FIS values at the population level. Genetic ancestry assignment was consistent when estimated with different sequencing methods and numbers of loci. Our GT‐seq panel is an effective and efficient genotyping tool that will aid in the monitoring and recovery of this threatened species, and our results provide insights for applying GT‐seq for minimally invasive DNA sampling techniques in other rare animals. [ABSTRACT FROM AUTHOR]

Striking examples of local adaptation at fine geographic scales are increasingly being documented in natural populations. However, the relative contributions made by natural selection, phenotype-dependent dispersal (when individuals disperse with respect to a habitat preference), and mate preference in generating and maintaining microgeographic adaptation and divergence are not well studied. Here, we develop quantitative genetics models and individual-based simulations (IBSs) to uncover the evolutionary forces that possibly drive microgeographic divergence. We also perform Bayesian estimation of the parameters in our IBS using empirical data on habitat-specific variation in bill morphology in the island scrub-jay (Aphelocoma insularis) to apply our models to a natural system. We find that natural selection and phenotype-dependent dispersal can generate the patterns of divergence we observe in the island scrub-jay. However, mate preference for a mate with similar bill morphology, even though observed in the species, does not play a significant role in driving divergence. Our modeling approach provides insights into phenotypic evolution occurring over small spatial scales relative to dispersal ranges, suggesting that adaptive divergence at microgeographic scales may be common across a wider range of taxa than previously thought. Our quantitative genetic models help to inform future theoretical and empirical work to determine how selection, habitat preference, and mate preference contribute to local adaptation and microgeographic divergence. [ABSTRACT FROM AUTHOR]