Ancient mitochondrial genomes unveil the origins and evolutionary history of New Zealand's enigmatic takahē and moho

Many avian species endemic to Aotearoa New Zealand were driven to extinction or reduced to relict populations following successive waves of human arrival, due to hunting, habitat destruction and the introduction of mammalian predators. Among the affected species were the large flightless South Island takahē (Porphyrio hochstetteri) and the moho (North Island takahē; P. mantelli), with the latter rendered extinct and the former reduced to a single relictual population. Little is known about the evolutionary history of these species prior to their decline and/or extinction. Here we sequenced mitochondrial genomes from takahē and moho subfossils (12 takahē and 4 moho) and retrieved comparable sequence data from takahē museum skins (n = 5) and contemporary individuals (n = 17) to examine the phylogeny and recent evolutionary history of these species. Our analyses suggest that prehistoric takahē populations lacked deep phylogeographic structure, in contrast to moho, which exhibited significant spatial genetic structure, albeit based on limited sample sizes (n = 4). Temporal genetic comparisons show that takahē have lost much of their mitochondrial genetic diversity, likely due to a sudden demographic decline soon after human arrival (~750 years ago). Time-calibrated phylogenetic analyses strongly support a sister species relationship between takahē and moho, suggesting these flightless taxa diverged around 1.5 million years ago, following a single colonisation of New Zealand by a flighted Porphyrio ancestor approximately 4 million years ago. This study highlights the utility of palaeogenetic approaches for informing the conservation and systematic understanding of endangered species whose ranges have been severely restricted by anthropogenic impacts.

A framework for identifying fertility gene targets for mammalian pest control

Fertility-targeted gene drives have been proposed as an ethical genetic approach for managing wild populations of vertebrate pests for public health and conservation benefit. This manuscript introduces a framework to identify and evaluate target gene suitability based on biological gene function, gene expression and results from mouse knockout models. This framework identified 16 genes essential for male fertility and 12 genes important for female fertility that may be feasible targets for mammalian gene drives and other non-drive genetic pest control technology. Further, a comparative genomics analysis demonstrates the conservation of the identified genes across several globally significant invasive mammals. In addition to providing important considerations for identifying candidate genes, our framework and the genes identified in this study may have utility in developing additional pest control tools such as wildlife contraceptives.

Genomics reveals repeated landlocking of diadromous fish on an isolated island

Landlocking of diadromous fish in freshwater systems can have significant genomic consequences. For instance, the loss of the migratory life stage can dramatically reduce gene flow across populations, leading to increased genetic structuring and stronger effects of local adaptation. These genomic consequences have been well-studied in some mainland systems, but the evolutionary impacts of landlocking in island ecosystems are largely unknown. In this study, we used a genotyping-by-sequencing (GBS) approach to examine the evolutionary history of landlocking in common smelt (Retropinna retropinna) on Chatham Island, a small isolated oceanic island 800 kilometres east of mainland New Zealand. We examined the relationship between Chatham Island and mainland smelt and used coalescent analyses to test the number and timing of landlocking events on Chatham Island. Our genomic analysis, based on 21,135 SNPs across 169 individuals, revealed that the Chatham Island smelt was genomically distinct from the mainland New Zealand fish, consistent with a single ancestral colonisation event of Chatham Island in the Pleistocene. Significant genetic structure was also evident within the Chatham Island smelt, with a diadromous Chatham Island smelt group, along with three geographically structured landlocked groups. Coalescent demographic analysis supported three independent landlocking events, with this loss of diadromy significantly pre-dating human colonisation. Our results illustrate how landlocking of diadromous fish can occur repeatedly across a narrow spatial scale, and highlight a unique system to study the genomic basis of repeated adaptation.

Fish biogeography and hybridization: do contemporary distributions predict introgression history?

Freshwater ecosystems frequently house diverse assemblages of closely related fish taxa, which can be particularly prone to hybridization and introgression. While extensive introgression may be expected among biogeographically proximate lineages, recent analyses imply that contemporary distributions do not always accurately predict hybridization history. Here, we use the ABBA-BABA approach to test biogeographic hypotheses regarding the extent of hybridization in the recent evolution of New Zealand's species-rich freshwater Galaxias vulgaris fish complex. Genome-wide comparisons reveal significant increases in introgression associated with increasing geographic overlap of taxa. The estimator DP, which assesses the net proportion of a genome originating from introgression, shows a particularly strong relationship with biogeographic overlap (R2 = .43; p = .005). Our analyses nevertheless reveal surprisingly substantial signatures of introgression among taxa that currently have disjunct distributions within drainages (e.g., separate subcatchments). These "anomalies" imply that current biogeography is not always an accurate predictor of introgression history. Our study suggests that both modern and ancient biogeographic shifts, including recent anthropogenic range fragmentation and tectonically driven riven capture events, have influenced introgression histories in this dynamic freshwater fish radiation.

A framework for identifying fertility gene targets for mammalian pest control

Fertility-targeted gene drives have been proposed as an ethical genetic approach for managing wild populations of vertebrate pests for public health and conservation benefit.This manuscript introduces a framework to identify and evaluate target gene suitability based on biological gene function, gene expression, and results from mouse knockout models.This framework identified 16 genes essential for male fertility and 12 genes important for female fertility that may be feasible targets for mammalian gene drives and other non-drive genetic pest control technology. Further, a comparative genomics analysis demonstrates the conservation of the identified genes across several globally significant invasive mammals.In addition to providing important considerations for identifying candidate genes, our framework and the genes identified in this study may have utility in developing additional pest control tools such as wildlife contraceptives.

Correlated evolution in an ectomycorrhizal host-symbiont system

Mechanisms of diversification in fungi are relatively poorly known. Many ectomycorrhizal symbionts show preference for particular host genera or families, so host-symbiont selection may be an important driver of fungal diversification in ectomycorrhizal systems. However, whether ectomycorrhizal hosts and symbionts show correlated evolutionary patterns remains untested, and it is unknown whether fungal specialisation also occurs in systems dominated by hosts from the same genus. We use metabarcoding of ectomycorrhizal fungi collected with hyphal ingrowth bags from Nothofagus forests across southern New Zealand to investigate host-symbiont specialisation and correlated evolution. We examine how ectomycorrhizal communities differ between host species and look for patterns of host-symbiont cophylogeny. We found substantial differences in ectomycorrhizal communities associated with different host taxa, particularly between hosts from different subgenera (Lophozonia and Fuscospora), but also between more closely related hosts. Twenty-four per cent of fungal taxa tested showed affiliations to particular hosts, and tests for cophylogeny revealed significant correlations between host relatedness and the fungal phylogeny that extended to substantial evolutionary depth. These results provide new evidence of correlated evolution in ectomycorrhizal systems, indicating that preferences among closely related host species may represent an important evolutionary driver for local lineage diversification in ectomycorrhizal fungi.

Fun in facets: A flexible new tool for population genomics in R

The landscape of analytical tools for population genomics continues to evolve. However, these tools are scattered across programming languages, making them largely inaccessible for many biologists. In this issue of Molecular Ecology Resources, Hemstrom and Jones, 2022 (Mol Ecol Resour; 962) introduce a new R package, snpR. This package combines a large number of existing analyses, to provide a one-stop shop for population genomics. F-statistics, admixture analyses, effective population size inferences, genome-wide association studies (GWAS), and parentage analyses are all implemented natively within the package. A variety of third-party software can also be run without leaving the R environment. The authors pay particular attention to data structure - avoiding redundancy - and allowing analyses to be run across multiple sample or single-nucleotide polymorphism (SNP) groupings. Because of its great accessibility and wide range of analyses, snpR has the potential to become a favourite within the Molecular Ecology community.

Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe

Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.

Paternal hypoxia exposure primes offspring for increased hypoxia resistance

Background

In a time of rapid environmental change, understanding how the challenges experienced by one generation can influence the fitness of future generations is critically needed. Using tolerance assays and transcriptomic and methylome approaches, we use zebrafish as a model to investigate cross-generational acclimation to hypoxia.

Results

We show that short-term paternal exposure to hypoxia endows offspring with greater tolerance to acute hypoxia. We detected two hemoglobin genes that are significantly upregulated by more than 6-fold in the offspring of hypoxia exposed males. Moreover, the offspring which maintained equilibrium the longest showed greatest upregulation in hemoglobin expression. We did not detect differential methylation at any of the differentially expressed genes, suggesting that other epigenetic mechanisms are responsible for alterations in gene expression.

Conclusions

Overall, our findings suggest that an epigenetic memory of past hypoxia exposure is maintained and that this environmentally induced information is transferred to subsequent generations, pre-acclimating progeny to cope with hypoxic conditions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01389-x.

Southern Hemisphere coasts are biologically connected by frequent, long-distance rafting events

Globally, species distributions are shifting in response to environmental change,1 and those that cannot disperse risk extinction.2 Many taxa, including marine species, are showing poleward range shifts as the climate warms.3 In the Southern Hemisphere, however, circumpolar oceanic fronts can present barriers to dispersal.4 Although passive, southward movement of species across this barrier has been considered unlikely,5,6 the recent discovery of buoyant kelp rafts on beaches in Antarctica7,8 demonstrates that such journeys are possible. Rafting is a key process by which diverse taxa-including terrestrial, e.g., Lindo,9 Godinot,10 and Censky et al.,11 and marine, e.g., Carlton et al.12 and Gillespie et al.13 species-can cross oceans.14 Kelp rafts can carry passengers7,15-17 and thus can act as vectors for long-distance dispersal of coastal organisms. The small numbers of kelp rafts previously found in Antarctica7,8 do not, however, shed much light on the frequency of such dispersal events.18 We use a combination of high-resolution phylogenomic analyses (>220,000 SNPs) and oceanographic modeling to show that long-distance biological dispersal events in Southern Ocean are not rare. We document tens of kelp (Durvillaea antarctica) rafting events of thousands of kilometers each, over several decades (1950-2019), with many kelp rafts apparently still reproductively viable. Modeling of dispersal trajectories from genomically inferred source locations shows that distant landmasses are well connected, for example South Georgia and New Zealand, and the Kerguelen Islands and Tasmania. Our findings illustrate the power of genomic approaches to track, and modeling to show frequencies of, long-distance dispersal events.

Are Cell Junctions Implicated in the Regulation of Vitellogenin Uptake? Insights from an RNAseq-Based Study in Eel, Anguilla australis

At the onset of puberty, ovarian follicles become competent to incorporate large amounts of vitellogenin (Vtg). Using an RNAseq-based approach, transcriptomes from pre-vitellogenic (PV) and early vitellogenic (EV) ovaries from wild-caught eel, Anguilla australis, were compared to investigate the expression of specific genes encoding cell junction proteins that could be involved in regulating Vtg uptake. Partial support was found for the mechanical barrier hypothesis proposing that the access of Vtg to the oolemma is restricted by a tight junction (TJ) network within the granulosa cell layer, which changes between the PV and EV stage. Among 25 genes encoding TJ-constituting proteins, five were down-regulated and two were up-regulated. A chemical barrier hypothesis stating that gap junctions (GJs) are involved in modulating Vtg uptake was not supported, as only five GJs were found to be expressed in the ovary with no significant changes in expression between stages. Furthermore, the endocytic pathway was found to be up-regulated during the PV-EV transition. Finally, the study showed that gene expression patterns may help identify suitable candidates involved in the regulation of Vtg uptake, and provided novel sequence data for A. australis, including putative Vtg receptors corresponding to Lr8 and Lrp13 members of the low-density lipoprotein receptor family.

Viromes of Freshwater Fish with Lacustrine and Diadromous Life Histories Differ in Composition

Viruses that infect fish are understudied, yet they provide important evolutionary context to the viruses that infect terrestrial vertebrates. We surveyed gill tissue meta-transcriptomes collected from two species of native freshwater fish from Aotearoa New Zealand—Retropinna retropinna and Gobiomorphus cotidianus. A total of 64 fish were used for gill tissue meta-transcriptomic sequencing, from populations with contrasting life histories—landlocked (i.e., lacustrine) and diadromous—on the South Island and Chatham Islands. We observed that both viral richness and taxonomic diversity were significantly associated with life history and host species, with lacustrine R. retropinna characterised by higher viral alpha diversity than diadromous R. retropinna. Additionally, we observed transcripts of fish viruses from 12 vertebrate host-associated virus families, and phylogenetically placed eight novel RNA viruses and three novel DNA viruses in the Astroviridae, Paramyxoviridae, Orthomyxoviridae, Rhabdoviridae, Totiviridae, Poxviridae, Alloherpesviridae, and Adintoviridae in their evolutionary contexts. These results represent an important survey of the viruses that infect two widespread native fish species in New Zealand, and provide insight useful for future fish virus surveys.

Concordant phylogeographic responses to large-scale coastal disturbance in intertidal macroalgae and their epibiota

Major ecological disturbance events can provide opportunities to assess multispecies responses to upheaval. In particular, catastrophic disturbances that regionally extirpate habitat-forming species can potentially influence the genetic diversity of large numbers of codistributed taxa. However, due to the rarity of such disturbance events over ecological timeframes, the genetic dynamics of multispecies recolonization processes have remained little understood. Here, we use single nucleotide polymorphism (SNP) data from multiple coastal species to track the dynamics of cocolonization events in response to ancient earthquake disturbance in southern New Zealand. Specifically, we use a comparative phylogeographic approach to understand the extent to which epifauna (with varying ecological associations with their macroalgal hosts) share comparable spatial and temporal recolonization patterns. Our study reveals concordant disturbance-related phylogeographic breaks in two intertidal macroalgal species along with two associated epibiotic species (a chiton and an isopod). By contrast, two codistributed species, one of which is an epibiotic amphipod and the other a subtidal macroalga, show few, if any, genetic effects of palaeoseismic coastal uplift. Phylogeographic model selection reveals similar post-uplift recolonization routes for the epibiotic chiton and isopod and their macroalgal hosts. Additionally, codemographic analyses support synchronous population expansions of these four phylogeographically similar taxa. Our findings indicate that coastal paleoseismic activity has driven concordant impacts on multiple codistributed species, with concerted recolonization events probably facilitated by macroalgal rafting. These results highlight that high-resolution comparative genomic data can help reconstruct concerted multispecies responses to recent ecological disturbance.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Genomic signatures of parallel alpine adaptation in recently evolved flightless insects

Natural selection along elevational gradients has potential to drive predictable adaptations across distinct lineages, but the extent of such repeated evolution remains poorly studied for many widespread alpine taxa. We present parallel genomic analyses of two recently evolved flightless alpine insect lineages to test for molecular signatures of repeated alpine adaptation. Specifically, we compare low-elevation vs. alpine stonefly ecotypes from parallel stream populations in which flightless upland ecotypes have been independently derived. We map 67,922 polymorphic genetic markers, generated across 176 Zelandoperla fenestrata specimens from two independent alpine stream populations in New Zealand's Rock and Pillar Range, to a newly developed plecopteran reference genome. Genome-wide scans revealed 31 regions with outlier single nucleotide polymorphisms (SNPs) differentiating lowland vs. alpine ecotypes in Lug Creek, and 37 regions with outliers differentiating ecotypes in Six Mile Creek. Of these regions, 13% (8/60) yielded outlier SNPs across both within-stream ecotype comparisons, implying comparable genomic shifts contribute to this repeated alpine adaptation. Candidate genes closely linked to repeated outlier regions include several with documented roles in insect wing-development (e.g., dishevelled), suggesting that they may contribute to repeated alpine wing reduction. Additional candidate genes have been shown to influence insect fecundity (e.g., ovo) and lifespan (e.g., Mrp4), implying that they might contribute to life history differentiation between upland and lowland ecotypes. Additional outlier genes have potential roles in the evolution of reproductive isolation among ecotypes (hedgehog and Desaturase 1). These results demonstrate how replicated outlier tests across independent lineages can potentially contribute to the discovery of genes underpinning repeated adaptation.

Genomic inference of contemporary effective population size in a large island population of collared flycatchers (Ficedula albicollis)

Due to its central importance to many aspects of evolutionary biology and population genetics, the long-term effective population size (N_e ) has been estimated for numerous species and populations. However, estimating contemporary N_e is difficult and in practice this parameter is often unknown. In principle, contemporary N_e can be estimated using either analyses of temporal changes in allele frequencies, or the extent of linkage disequilibrium (LD) between unlinked markers. We applied these approaches to estimate contemporary N_e of a relatively recently founded island population of collared flycatchers (Ficedula albicollis). We sequenced the genomes of 85 birds sampled in 1993 and 2015, and applied several temporal methods to estimate N_e at a few thousand (4000-7000). The approach based on LD provided higher estimates of N_e (20,000-32,000) and was associated with high variance, often resulting in infinite N_e . We conclude that whole-genome sequencing data offers new possibilities to estimate high (>1000) contemporary N_e , but also note that such estimates remain challenging, in particular for LD-based methods for contemporary N_e estimation.

Genomics Reveals Widespread Ecological Speciation in Flightless Insects

Recent genomic analyses have highlighted parallel divergence in response to ecological gradients, but the extent to which altitude can underpin such repeated speciation remains unclear. Wing reduction and flight loss have apparently evolved repeatedly in montane insect assemblages and have been suggested as important drivers of hexapod diversification. We test this hypothesis using genomic analyses of a widespread wing-polymorphic stonefly species complex in New Zealand. We identified over 50,000 polymorphic genetic markers generated across almost 200 Zelandoperla fenestrata stonefly specimens using a newly generated plecopteran reference genome, to reveal widespread parallel speciation between sympatric full-winged and wing-reduced ecotypes. Rather than the existence of a single, widespread, flightless taxon (Zelandoperla pennulata), evolutionary genomic data reveal that wing-reduced upland lineages have speciated repeatedly and independently from full-winged Z. fenestrata. This repeated evolution of reproductive isolation between local ecotype pairs that lack mitochondrial DNA differentiation suggests that ecological speciation has evolved recently. A cluster of outlier single-nucleotide polymorphisms detected in independently wing-reduced lineages, tightly linked in an approximately 85 kb genomic region that includes the developmental “supergene” doublesex, suggests that this “island of divergence” may play a key role in rapid ecological speciation. [Ecological speciation; genome assembly; genomic island of differentiation; genotyping-by-sequencing; incipient species; plecoptera; wing reduction.]

Tissue-specific patterns of regulatory changes underlying gene expression differences among Ficedula flycatchers and their naturally occurring F1 hybrids

Changes in interacting cis- and trans-regulatory elements are important candidates for Dobzhansky-Muller hybrid incompatibilities and may contribute to hybrid dysfunction by giving rise to misexpression in hybrids. To gain insight into the molecular mechanisms and determinants of gene expression evolution in natural populations, we analyzed the transcriptome from multiple tissues of two recently diverged Ficedula flycatcher species and their naturally occurring F₁ hybrids. Differential gene expression analysis revealed that the extent of differentiation between species and the set of differentially expressed genes varied across tissues. Common to all tissues, a higher proportion of Z-linked genes than autosomal genes showed differential expression, providing evidence for a fast-Z effect. We further found clear signatures of hybrid misexpression in brain, heart, kidney, and liver. However, while testis showed the highest divergence of gene expression among tissues, it showed no clear signature of misexpression in F₁ hybrids, even though these hybrids were found to be sterile. It is therefore unlikely that incompatibilities between cis-trans regulatory changes explain the observed sterility. Instead, we found evidence that cis-regulatory changes play a significant role in the evolution of gene expression in testis, which illustrates the tissue-specific nature of cis-regulatory evolution bypassing constraints associated with pleiotropic effects of genes.

The search for sexually antagonistic genes: Practical insights from studies of local adaptation and statistical genomics

Sexually antagonistic (SA) genetic variation-in which alleles favored in one sex are disfavored in the other-is predicted to be common and has been documented in several animal and plant populations, yet we currently know little about its pervasiveness among species or its population genetic basis. Recent applications of genomics in studies of SA genetic variation have highlighted considerable methodological challenges to the identification and characterization of SA genes, raising questions about the feasibility of genomic approaches for inferring SA selection. The related fields of local adaptation and statistical genomics have previously dealt with similar challenges, and lessons from these disciplines can therefore help overcome current difficulties in applying genomics to study SA genetic variation. Here, we integrate theoretical and analytical concepts from local adaptation and statistical genomics research-including F ST and F IS statistics, genome-wide association studies, pedigree analyses, reciprocal transplant studies, and evolve-and-resequence experiments-to evaluate methods for identifying SA genes and genome-wide signals of SA genetic variation. We begin by developing theoretical models for between-sex F ST and F IS, including explicit null distributions for each statistic, and using them to critically evaluate putative multilocus signals of sex-specific selection in previously published datasets. We then highlight new statistics that address some of the limitations of F ST and F IS, along with applications of more direct approaches for characterizing SA genetic variation, which incorporate explicit fitness measurements. We finish by presenting practical guidelines for the validation and evolutionary analysis of candidate SA genes and discussing promising empirical systems for future work.

Roe deer on ice: Selection despite limited effective population size through the Pleistocene

Roe deer (Capreolus spp.) are a little odd. They are one of only a few placental mammals-and the only genus among even-toed ungulates-capable of putting embryonic development "on ice", also known as embryonic diapause (Figure 1). It would seem such an unusual trait is probably the product of natural selection, but a big question is, how does selection for important traits, such as diapause, interact with the historical demography of a species? In a 'From the Cover' article in this issue of Molecular Ecology, de Jong et al. (2020) demonstrate that selection is acting on genes associated with reproductive biology in roe deer, despite heightened genetic drift due to reduced effective population size through the Pleistocene.

The genomic footprint of coastal earthquake uplift

Theory suggests that catastrophic earth-history events can drive rapid biological evolution, but empirical evidence for such processes is scarce. Destructive geological events such as earthquakes can represent large-scale natural experiments for inferring such evolutionary processes. We capitalized on a major prehistoric (800 yr BP) geological uplift event affecting a southern New Zealand coastline to test for the lasting genomic impacts of disturbance. Genome-wide analyses of three co-distributed keystone kelp taxa revealed that post-earthquake recolonization drove the evolution of novel, large-scale intertidal spatial genetic 'sectors' which are tightly linked to geological fault boundaries. Demographic simulations confirmed that, following widespread extirpation, parallel expansions into newly vacant habitats rapidly restructured genome-wide diversity. Interspecific differences in recolonization mode and tempo reflect differing ecological constraints relating to habitat choice and dispersal capacity among taxa. This study highlights the rapid and enduring evolutionary effects of catastrophic ecosystem disturbance and reveals the key role of range expansion in reshaping spatial genetic patterns.

Population structure of the New Zealand whelk, Cominella glandiformis (Gastropoda: Buccinidae), suggests sporadic dispersal of a direct developer

We examined phylogeographic structure in the direct-developing New Zealand endemic intertidal mud whelk, Cominella glandiformis. Two hundred and ninety-six whelks from 12 sites were collected from sheltered shores around New Zealand’s four largest islands (North Island, South Island, Stewart Island and Chatham Island), encompassing the geographical range of this species. Despite being direct developers, gene flow among C. glandiformis populations may occur over short distances by adult floating, and over larger distances by rafting of egg masses. Primers were developed to amplify variable microsatellite regions at six loci. All loci were variable, with 8–34 alleles/loci. Observed and expected heterozygosities were high across all alleles, with minimal evidence of null alleles. The average number of alleles varied from 3.5 (Chatham Island) to 7.5 (Waitemata Harbour). Strong genetic structure was evident, with distinct ‘eastern’ and ‘western’ groups. Each group extended over a large geographic area, including regions of unsuitable habitat, but were linked by oceanic currents. We suggest that the intraspecific geographic genetic structure in C. glandiformis has arisen due a combination of ocean currents (promoting gene flow between geographically distant regions) and upwelling areas (limiting gene flow between certain regions).

Species in the faeces: DNA metabarcoding as a method to determine the diet of the endangered yellow-eyed penguin

Abstract Context. Diet variability is a significant driver of seabird decline; however, data on seabird diet composition and trends have been affected by changes in precision and resolution owing to the evolution of different sampling methods over time. We investigated the effectiveness of applying a passive molecular diet method using faeces obtained from the endangered yellow-eyed penguin. Aims. To assess the feasibility of applying DNA metabarcoding methods to yellow-eyed penguin faeces to evaluate diet, and to compare the reliability of diet results derived from adults and chicks, and from latrine versus fresh faecal samples. Methods. We collected 313 faecal samples from yellow-eyed penguins resident on the Otago coast of New Zealand from October 2016 to August 2017. We used polymerase chain reaction (PCR) with mitochondrial 16S cephalopod and chordate primers to amplify prey DNA present in the faecal samples, and tested the completeness of our assembled reference databases based on previous diet research. Amplified prey DNA sequences were then assigned to taxa from our reference databases by using QIIME2. Key results. Mitochondrial 16S chordate PCR primers were effective at identifying 29 fish taxa, with 98.3% of amplified sequences being identified to species or genus level in 193 samples (61.7% collected). There was no significant difference in the number, occurrence or proportion of ray-finned fish prey DNA sequences derived from fresh samples or latrines. Mitochondrial 16S cephalopod PCR primers classified 1.98% of amplified DNA sequences as targets, with 96.5% of these target sequences being identified to species or genus level in 48 samples (15.3% collected), and five taxa identified. Conclusions. We recommend the collection of latrine samples to enable long-term monitoring of the diet of yellow-eyed penguins, which will optimise the trade-off between wildlife disturbance and dietary resolution. Further refinement is needed to identify cephalopod dietary components for yellow-eyed penguins, because our cephalopod primers were not as specific as those used for ray-finned fishes, amplifying a large number (>98%) of non-cephalopod species. Implications. DNA metabarcoding offers a robust and comprehensive alternative to other, more intrusive, seabird diet-assessment methods, but still requires parallel studies to provide critical information on prey size, true diet composition and diet quality.

SNP analyses reveal a diverse pool of potential colonists to earthquake-uplifted coastlines

In species that form dense populations, major disturbance events are expected to increase the chance of establishment for immigrant lineages. Real-time tests of the impact of disturbance on patterns of genetic structure are, however, scarce. Central to testing these concepts is determining the pool of potential immigrants dispersing into a disturbed area. In 2016, a 7.8 magnitude earthquake occurred on the South Island of New Zealand. Affecting approximately 100 km of coastline, this quake caused extensive uplift (several metres high), extirpating many intertidal populations, including keystone intertidal kelp species. Following the uplift, we set out to determine the geographic origins of detached kelp specimens which rafted into the disturbed zone. Specifically, we used genotyping-by-sequencing (GBS) approaches to compare beach-cast southern bull-kelp (Durvillaea antarctica and Durvillaea poha) samples to established populations throughout the species' ranges, and thus infer the geographic origins of potential colonists reaching the disturbed coast. Our findings revealed an ongoing supply of diverse lineages dispersing to the newly uplifted coastline, suggesting potential for establishment of "exotic" lineages following disturbance. Furthermore, we found that some drifting individuals of each species came from far-distant regions, some >1,200 km away. These results show that diverse lineages - in many cases from very distant sources - can compete for new space in the wake of an exceptional disturbance event, illustrating the potential of long-distance dispersal as a key mechanism for reassembly of coastal ecosystems. Furthermore, our findings demonstrate that high-resolution genomic baselines can be used to robustly assign the provenance of dispersing individuals.

Genomics detects population structure within and between ocean basins in a circumpolar seabird: The white-chinned petrel

The Southern Ocean represents a continuous stretch of circumpolar marine habitat, but the potential physical and ecological drivers of evolutionary genetic differentiation across this vast ecosystem remain unclear. We tested for genetic structure across the full circumpolar range of the white-chinned petrel (Procellaria aequinoctialis) to unravel the potential drivers of population differentiation and test alternative population differentiation hypotheses. Following range-wide comprehensive sampling, we applied genomic (genotyping-by-sequencing or GBS; 60,709 loci) and standard mitochondrial-marker approaches (cytochrome b and first domain of control region) to quantify genetic diversity within and among island populations, test for isolation by distance, and quantify the number of genetic clusters using neutral and outlier (non-neutral) loci. Our results supported the multi-region hypothesis, with a range of analyses showing clear three-region genetic population structure, split by ocean basin, within two evolutionary units. The most significant differentiation between these regions confirmed previous work distinguishing New Zealand and nominate subspecies. Although there was little evidence of structure within the island groups of the Indian or Atlantic oceans, a small set of highly-discriminatory outlier loci could assign petrels to ocean basin and potentially to island group, though the latter needs further verification. Genomic data hold the key to revealing substantial regional genetic structure within wide-ranging circumpolar species previously assumed to be panmictic.

Local adaptation and the evolution of inversions on sex chromosomes and autosomes

Spatially varying selection with gene flow can favour the evolution of inversions that bind locally adapted alleles together, facilitate local adaptation and ultimately drive genomic divergence between species. Several studies have shown that the rates of spread and establishment of new inversions capturing locally adaptive alleles depend on a suite of evolutionary factors, including the strength of selection for local adaptation, rates of gene flow and recombination, and the deleterious mutation load carried by inversions. Because the balance of these factors is expected to differ between X (or Z) chromosomes and autosomes, opportunities for inversion evolution are likely to systematically differ between these genomic regions, though such scenarios have not been formally modelled. Here, we consider the evolutionary dynamics of X-linked and autosomal inversions in populations evolving at a balance between migration and local selection. We identify three factors that lead to asymmetric rates of X-linked and autosome inversion establishment: (1) sex-biased migration, (2) dominance of locally adapted alleles and (3) chromosome-specific deleterious mutation loads. This theory predicts an elevated rate of fixation, and depressed opportunities for polymorphism, for X-linked inversions. Our survey of data on the genomic distribution of polymorphic and fixed inversions supports both theoretical predictions.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Biased Inference of Selection Due to GC-Biased Gene Conversion and the Rate of Protein Evolution in Flycatchers When Accounting for It

The rate of recombination impacts on rates of protein evolution for at least two reasons: it affects the efficacy of selection due to linkage and influences sequence evolution through the process of GC-biased gene conversion (gBGC). We studied how recombination, via gBGC, affects inferences of selection in gene sequences using comparative genomic and population genomic data from the collared flycatcher (Ficedula albicollis). We separately analyzed different mutation categories (“strong”-to-“weak,” “weak-to-strong,” and GC-conservative changes) and found that gBGC impacts on the distribution of fitness effects of new mutations, and leads to that the rate of adaptive evolution and the proportion of adaptive mutations among nonsynonymous substitutions are underestimated by 22–33%. It also biases inferences of demographic history based on the site frequency spectrum. In light of this impact, we suggest that inferences of selection (and demography) in lineages with pronounced gBGC should be based on GC-conservative changes only. Doing so, we estimate that 10% of nonsynonymous mutations are effectively neutral and that 27% of nonsynonymous substitutions have been fixed by positive selection in the flycatcher lineage. We also find that gene expression level, sex-bias in expression, and the number of protein–protein interactions, but not Hill–Robertson interference (HRI), are strong determinants of selective constraint and rate of adaptation of collared flycatcher genes. This study therefore illustrates the importance of disentangling the effects of different evolutionary forces and genetic factors in interpretation of sequence data, and from that infer the role of natural selection in DNA sequence evolution.

Covariation in levels of nucleotide diversity in homologous regions of the avian genome long after completion of lineage sorting

Closely related species may show similar levels of genetic diversity in homologous regions of the genome owing to shared ancestral variation still segregating in the extant species. However, after completion of lineage sorting, such covariation is not necessarily expected. On the other hand, if the processes that govern genetic diversity are conserved, diversity may potentially covary even among distantly related species. We mapped regions of conserved synteny between the genomes of two divergent bird species-collared flycatcher and hooded crow-and identified more than 600 Mb of homologous regions (66% of the genome). From analyses of whole-genome resequencing data in large population samples of both species we found nucleotide diversity in 200 kb windows to be well correlated (Spearman's ρ = 0.407). The correlation remained highly similar after excluding coding sequences. To explain this covariation, we suggest that a stable avian karyotype and a conserved landscape of recombination rate variation render the diversity-reducing effects of linked selection similar in divergent bird lineages. Principal component regression analysis of several potential explanatory variables driving heterogeneity in flycatcher diversity levels revealed the strongest effects from recombination rate variation and density of coding sequence targets for selection, consistent with linked selection. It is also possible that a stable karyotype is associated with a conserved genomic mutation environment contributing to covariation in diversity levels between lineages. Our observations imply that genetic diversity is to some extent predictable.

Divergence in gene expression within and between two closely related flycatcher species

Relatively little is known about the character of gene expression evolution as species diverge. It is for instance unclear if gene expression generally evolves in a clock‐like manner (by stabilizing selection or neutral evolution) or if there are frequent episodes of directional selection. To gain insights into the evolutionary divergence of gene expression, we sequenced and compared the transcriptomes of multiple organs from population samples of collared (Ficedula albicollis) and pied flycatchers (F. hypoleuca), two species which diverged less than one million years ago. Ordination analysis separated samples by organ rather than by species. Organs differed in their degrees of expression variance within species and expression divergence between species. Variance was negatively correlated with expression breadth and protein interactivity, suggesting that pleiotropic constraints reduce gene expression variance within species. Variance was correlated with between‐species divergence, consistent with a pattern expected from stabilizing selection and neutral evolution. Using an expression P_ST approach, we identified genes differentially expressed between species and found 16 genes uniquely expressed in one of the species. For one of these, DPP7, uniquely expressed in collared flycatcher, the absence of expression in pied flycatcher could be associated with a ≈20‐kb deletion including 11 of 13 exons. This study of a young vertebrate speciation model system expands our knowledge of how gene expression evolves as natural populations become reproductively isolated.

Linked selection and recombination rate variation drive the evolution of the genomic landscape of differentiation across the speciation continuum of Ficedula flycatchers

Speciation is a continuous process during which genetic changes gradually accumulate in the genomes of diverging species. Recent studies have documented highly heterogeneous differentiation landscapes, with distinct regions of elevated differentiation (“differentiation islands”) widespread across genomes. However, it remains unclear which processes drive the evolution of differentiation islands; how the differentiation landscape evolves as speciation advances; and ultimately, how differentiation islands are related to speciation. Here, we addressed these questions based on population genetic analyses of 200 resequenced genomes from 10 populations of four Ficedula flycatcher sister species. We show that a heterogeneous differentiation landscape starts emerging among populations within species, and differentiation islands evolve recurrently in the very same genomic regions among independent lineages. Contrary to expectations from models that interpret differentiation islands as genomic regions involved in reproductive isolation that are shielded from gene flow, patterns of sequence divergence (d_xy and relative node depth) do not support a major role of gene flow in the evolution of the differentiation landscape in these species. Instead, as predicted by models of linked selection, genome-wide variation in diversity and differentiation can be explained by variation in recombination rate and the density of targets for selection. We thus conclude that the heterogeneous landscape of differentiation in Ficedula flycatchers evolves mainly as the result of background selection and selective sweeps in genomic regions of low recombination. Our results emphasize the necessity of incorporating linked selection as a null model to identify genome regions involved in adaptation and speciation.

How a haemosporidian parasite of bats gets around: the genetic structure of a parasite, vector and host compared

Parasite population structure is often thought to be largely shaped by that of its host. In the case of a parasite with a complex life cycle, two host species, each with their own patterns of demography and migration, spread the parasite. However, the population structure of the parasite is predicted to resemble only that of the most vagile host species. In this study, we tested this prediction in the context of a vector-transmitted parasite. We sampled the haemosporidian parasite Polychromophilus melanipherus across its European range, together with its bat fly vector Nycteribia schmidlii and its host, the bent-winged bat Miniopterus schreibersii. Based on microsatellite analyses, the wingless vector, and not the bat host, was identified as the least structured population and should therefore be considered the most vagile host. Genetic distance matrices were compared for all three species based on a mitochondrial DNA fragment. Both host and vector populations followed an isolation-by-distance pattern across the Mediterranean, but not the parasite. Mantel tests found no correlation between the parasite and either the host or vector populations. We therefore found no support for our hypothesis; the parasite population structure matched neither vector nor host. Instead, we propose a model where the parasite's gene flow is represented by the added effects of host and vector dispersal patterns.