In cases of severe wildlife population decline, a key question is whether recovery efforts will be impeded by genetic factors, such as inbreeding depression. Decades of excess mortality from gillnet fishing have driven Mexico’s vaquita porpoise ( Phocoena sinus ) to ~10 remaining individuals. We analyzed whole-genome sequences from 20 vaquitas and integrated genomic and demographic information into stochastic, individual-based simulations to quantify the species’ recovery potential. Our analysis suggests that the vaquita’s historical rarity has resulted in a low burden of segregating deleterious variation, reducing the risk of inbreeding depression. Similarly, genome-informed simulations suggest that the vaquita can recover if bycatch mortality is immediately halted. This study provides hope for vaquitas and other naturally rare endangered species and highlights the utility of genomics in predicting extinction risk.

The genetic consequences of species-wide declines are rarely quantified because the timing and extent of the decline varies across the species' range. The sea otter (Enhydra lutris) is a unique model in this regard. Their dramatic decline from thousands to fewer than 100 individuals per population occurred range-wide and nearly simultaneously due to the 18th-19th century fur trade. Consequently, each sea otter population represents an independent natural experiment of recovery after extreme population decline. We designed sequence capture probes for 50 Mb of sea otter exonic and neutral genomic regions. We sequenced 107 sea otters from five populations that span the species range to high coverage (18-76×) and three historical Californian samples from ~1500 and ~200 years ago to low coverage (1.5-3.5×). We observe distinct population structure and find that sea otters in California are the last survivors of a divergent lineage isolated for thousands of years and therefore warrant special conservation concern. We detect signals of extreme population decline

Phenotypic variation along environmental gradients can provide evidence suggesting local adaptation has shaped observed morphological disparities. These differences, in traits such as body and extremity size, as well as skin and coat pigmentation, may affect the overall fitness of individuals in their environments. The Virginia opossum (Didelphis virginiana) is a marsupial that shows phenotypic variation across its range, one that has recently expanded into temperate environments. It is unknown, however, whether the variation observed in the species fits adaptive ecogeographic patterns, or if phenotypic change is associated with any environmental factors. Using phenotypic measurements of over 300 museum specimens of Virginia opossum, collected throughout its distribution range, we applied regression analysis to determine if phenotypes change along a latitudinal gradient. Then, using predictors from remote-sensing databases and a random forest algorithm, we tested environmental models to find the most important variables driving the phenotypic variation. We found that despite the recent expansion into temperate environments, the phenotypic variation in the Virginia opossum follows a latitudinal gradient fitting three adaptive ecogeographic patterns codified under Bergmann’s, Allen’s and Gloger’s rules. Temperature seasonality was an important predictor of body size variation, with larger opossums occurring at high latitudes with more seasonal environments. Annual mean temperature predicted important variation in extremity size, with smaller extremities found in northern populations. Finally, we found that precipitation and temperature seasonality as well as low temperatures were strong environmental predictors of skin and coat pigmentation variation; darker opossums are distributed at low latitudes in warmer environments with higher precipitation seasonality. These results indicate that the adaptive mechanisms underlying the variation in body size, extremity size and pigmentation are related to the resource seasonality, heat conservation, and pathogen-resistance hypotheses, respectively. Our findings suggest that marsupials may be highly susceptible to environmental changes, and in the case of the Virginia opossum, the drastic phenotypic evolution in northern populations may have arisen rapidly, facilitating the colonization of seasonal and colder habitats of temperate North America.

Four fin whale sub-species are currently considered valid: Balaenoptera physalus physalus in the North Atlantic, B. p. velifera in the North Pacific, B. p. quoyi and B. p. patachonica in the Southern Hemisphere. The last, not genetically validated, was described as a pygmy-type sub-species, found in low to mid latitudes of the Southern Hemisphere. Genetic analyses across hemispheres show strong phylogeographic structure, yet low geographic coverage in middle latitudes of the Southern Hemisphere impeded an assessment within the area, as well as evaluating the validity of B. p. patachonica. New mtDNA sequences from the Southeastern Pacific allowed an improved coverage of the species’ distribution. Our phylogenetic analyses showed three main lineages and contrasting phylogeographic patterns between Northern and Southern Hemispheres. Absence of recurrent female mediated gene flow between hemispheres was found; however, rare dispersal events revealing old migrations were noted. The absence of genetic structure suggests the existence of one single taxa within the Southern Hemisphere. Thus, until further evidence supporting this subspecies can be produced, such as genetic, ecological, behavioral, or morphological data, we propose that all fin whales from the Southern Hemisphere, including those from middle latitudes of the Southeastern Pacific belong to B. p. quoyi subspecies. This information is important for the current assessment of fin whales, contributing to the evaluation of the taxonomic classification and the conservation of the species.

One of the most isolated populations of fin whales occurs in the Gulf of California (GOC) with 400-800 individuals. This population shows reduced neutral genetic variation in comparison to the North Pacific population and thus might also display limited adaptive polymorphism. We sampled 36 fin whales from the GOC and assessed genetic variation at exon 2 of the major histocompatibility complex class II DQB-1 genes responsible for eliciting immune responses. Three divergent alleles were found with higher nonsynonymous than synonymous substitution rates within the peptide-binding region positions as well as the likely retention of ancient alleles, indicating that positive selection has shaped diversity in this species. Limited levels of nonneutral polymorphism, in addition to previously described low levels of neutral polymorphism, are consistent with the results of previous studies on vertebrate populations that have remained small and demographically stable for a very long time. Such low genetic variation in the GOC fin whales could be explained by 2 demographic scenarios: an ancient isolated population with limited gene flow or a more recent founder event after the last glacial maximum with very restricted gene flow.

No abstract

Skin pigmentation and coat pigmentation are two of the best-studied examples of traits under natural selection given their quantifiable fitness interactions with the environment (e.g., camouflage) and signalling with other organisms (e.g., warning coloration). Previous morphological studies have found that skin pigmentation variation in the Virginia opossum (Didelphis virginiana) is associated with variation in precipitation and temperatures across its distribution range following Gloger's rule (lighter pigmentation in temperate environments). To investigate the molecular mechanism associated with skin pigmentation variation, we used RNA-Seq and quantified gene expression of wild opossums from tropical and temperate populations. Using differential expression analysis and a co-expression network approach, we found that expression variation in genes with melanocytic and immune functions is significantly associated with the degree of skin pigmentation variation and may be underlying this phenotypic difference. We also found evidence suggesting that the Wnt/β-catenin signalling pathway might be regulating the depigmentation observed in temperate populations. Based on our study results, we present several alternative hypotheses that may explain Gloger's rule pattern of skin pigmentation variation in opossum, including changes in pathogen diversity supporting a pathogen-resistant hypothesis, thermal stress associated with temperate environments, and pleiotropic and epistatic interactions between melanocytic and immune genes.

Little is known about how the spectrum and etiology of germline mutagenesis might vary among mammalian species. To shed light on this mystery, we quantify variation in mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k-mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clocklike mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these clocklike signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the mutation spectrum's phylogenetic signal when fit to non-context-dependent mutation spectrum data in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.