Seeing the whole picture: What molecular ecology is gaining from whole genomes

The genome of Anoplarchus purpurescens (Stichaeidae) reflects its carnivorous diet

Digestion is driven by digestive enzymes and digestive enzyme gene copy number can provide insights on the genomic underpinnings of dietary specialization. The "Adaptive Modulation Hypothesis" (AMH) proposes that digestive enzyme activity, which increases with increased gene copy number, should correlate with substrate quantity in the diet. To test the AMH and reveal some of the genetics of herbivory vs carnivory, we sequenced, assembled, and annotated the genome of Anoplarchus purpurescens, a carnivorous prickleback fish in the family Stichaeidae, and compared the gene copy number for key digestive enzymes to that of Cebidichthys violaceus, a herbivorous fish from the same family. A highly contiguous genome assembly of high quality (N50 = 10.6 Mb) was produced for A. purpurescens, using combined long-read and short-read technology, with an estimated 33,842 protein-coding genes. The digestive enzymes that we examined include pancreatic α-amylase, carboxyl ester lipase, alanyl aminopeptidase, trypsin, and chymotrypsin. Anoplarchus purpurescens had fewer copies of pancreatic α-amylase (carbohydrate digestion) than C. violaceus (1 vs. 3 copies). Moreover, A. purpurescens had one fewer copy of carboxyl ester lipase (plant lipid digestion) than C. violaceus (4 vs. 5). We observed an expansion in copy number for several protein digestion genes in A. purpurescens compared to C. violaceus, including trypsin (5 vs. 3) and total aminopeptidases (6 vs. 5). Collectively, these genomic differences coincide with measured digestive enzyme activities (phenotypes) in the two species and they support the AMH. Moreover, this genomic resource is now available to better understand fish biology and dietary specialization.

Species-wide genomics of kākāpō provides tools to accelerate recovery

The kākāpō is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analysed whole-genome sequence data for nearly all individuals living in early 2018 (169 individuals) to generate a high-quality species-wide genetic variant callset. We leverage extensive long-term metadata to quantify genome-wide diversity of the species over time and present new approaches using probabilistic programming, combined with a phenotype dataset spanning five decades, to disentangle phenotypic variance into environmental and genetic effects while quantifying uncertainty in small populations. We find associations for growth, disease susceptibility, clutch size and egg fertility within genic regions previously shown to influence these traits in other species. Finally, we generate breeding values to predict phenotype and illustrate that active management over the past 45 years has maintained both genome-wide diversity and diversity in breeding values and, hence, evolutionary potential. We provide new pathways for informing future conservation management decisions for kākāpō, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of the small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

A comparison of genomic diversity and demographic history of the North Atlantic and Southwest Atlantic southern right whales

Right whales (genus Eubalaena) were among the first, and most extensively pursued, targets of commercial whaling. However, understanding the impacts of this persecution requires knowledge of the demographic histories of these species prior to exploitation. We used deep whole genome sequencing (~40×) of 12 North Atlantic (E. glacialis) and 10 Southwest Atlantic southern (E. australis) right whales to quantify contemporary levels of genetic diversity and infer their demographic histories over time. Using coalescent- and identity-by-descent-based modelling to estimate ancestral effective population sizes from genomic data, we demonstrate that North Atlantic right whales have lived with smaller effective population sizes (Ne ) than southern right whales in the Southwest Atlantic since their divergence and describe the decline in both populations around the time of whaling. North Atlantic right whales exhibit reduced genetic diversity and longer runs of homozygosity leading to higher inbreeding coefficients compared to the sampled population of southern right whales. This study represents the first comprehensive assessment of genome-wide diversity of right whales in the western Atlantic and underscores the benefits of high coverage, genome-wide datasets to help resolve long-standing questions about how historical changes in effective population size over different time scales shape contemporary diversity estimates. This knowledge is crucial to improve our understanding of the right whales' history and inform our approaches to address contemporary conservation issues. Understanding and quantifying the cumulative impact of long-term small Ne , low levels of diversity and recent inbreeding on North Atlantic right whale recovery will be important next steps.

Museum specimens of a landlocked pinniped reveal recent loss of genetic diversity and unexpected population connections

The Saimaa ringed seal (Pusa hispida saimensis) is endemic to Lake Saimaa in Finland. The subspecies is thought to have originated when parts of the ringed seal population of the Baltic region were trapped in lakes emerging due to postglacial bedrock rebound around 9000 years ago. During the 20th century, the population experienced a drastic human-induced bottleneck. Today encompassing a little over 400 seals with extremely low genetic diversity, it is classified as endangered. We sequenced sections of the mitochondrial control region from 60 up to 125-years-old museum specimens of the Saimaa ringed seal. The generated dataset was combined with publicly available sequences. We studied how genetic variation has changed through time in this subspecies and how it is phylogenetically related to other ringed seal populations from the Baltic Sea, Lake Ladoga, North America, Svalbard, and the White Sea. We observed temporal fluctuations in haplotype frequencies and loss of haplotypes accompanied by a recent reduction in female effective population size. In apparent contrast with the traditionally held view of the Baltic origin of the population, the Saimaa ringed seal mtDNA variation also shows affinities to North American ringed seals. Our results suggest that the Saimaa ringed seal has experienced recent genetic drift associated with small population size. The results further suggest that extant Baltic ringed seal is not representative of the ancestral population of the Saimaa ringed seal, which calls for re-evaluation of the deep history of this subspecies.

The genomic basis of reproductive and migratory behaviour in a polymorphic salmonid

Recent ecotypic differentiation provides unique opportunities to investigate the genomic basis and architecture of local adaptation, while offering insights into how species form and persist. Sockeye salmon (Oncorhynchus nerka) exhibit migratory and resident ("kokanee") ecotypes, which are further distinguished into shore-spawning and stream-spawning reproductive ecotypes. Here, we analysed 36 sockeye (stream-spawning) and kokanee (stream- and shore-spawning) genomes from a system where they co-occur and have recent common ancestry (Okanagan Lake/River in British Columbia, Canada) to investigate the genomic basis of reproductive and migratory behaviour. Examination of the genomic landscape of differentiation, differences in allele frequencies and genotype-phenotype associations revealed three main blocks of sequence differentiation on chromosomes 7, 12 and 20, associated with migratory behaviour, spawning location and spawning timing. Structural variants identified in these same areas suggest they could contribute to ecotypic differentiation directly as causal variants or via maintenance of their genomic architecture through recombination suppression mechanisms. Genes in these regions were related to spatial memory and swimming endurance (SYNGAP, TPM3), as well as eye and brain development (including SIX6), potentially associated with differences in migratory behaviour and visual habitats across spawning locations, respectively. Additional genes (GREB1L, ROCK1) identified here have been associated with timing of migration in other salmonids and could explain variation in timing of O. nerka spawning. Together, these results based on the joint analysis of sequence and structural variation represent a significant advance in our understanding of the genomic landscape of ecotypic differentiation at different stages in the speciation continuum.