Improving the estimation of genetic distances from Next-Generation Sequencing data

Circumpolar Population Structure, Diversity and Recent Evolutionary History of the Bearded Seal in Relation to Past and Present Icescapes

The Arctic environment plays a critical role in the global climate system and marine biodiversity. The region's ice-covered expanses provide essential breeding and feeding grounds for a diverse assemblage of marine species, who have adapted to thrive in these harsh conditions and consequently are under threat from global warming. The bearded seal (Erignathus barbatus), including two subspecies (E. barbatus nauticus-Pacific and E. barbatus barbatus-Atlantic), is an ice-obligate Arctic species using sea ice for many aspects of its life history, rendering it particularly vulnerable to sea ice loss. It is one of the least studied and hence enigmatic of the Arctic marine mammals, with little knowledge regarding genetic structure, diversity, adaptations, and demographic history, consequently hampering management and conservation efforts. Here, we sequenced 70 whole nuclear genomes from across most of the species' circumpolar range, finding significant genetic structure between the Pacific and the Atlantic subspecies, which diverged during the Penultimate Glacial Period (~200 KYA). Remarkably, we found fine-scale genetic structure within both subspecies, with at least two distinct populations in the Pacific and three in the Atlantic. We hypothesise sea-ice dynamics and bathymetry had a prominent role in shaping bearded seal genetic structure and diversity. Our analyses of highly differentiated genomic regions can be used to complement the health, physiological, and behavioural research needed to conserve this species. In addition, we provide recommendations for management units that can be used to more specifically assess climatic and anthropogenic impacts on bearded seal populations.

Circumpolar and Regional Seascape Drivers of Genomic Variation in a Southern Ocean Octopus

Understanding how ecological, environmental and geographic features influence population genetic patterns provides crucial insights into a species' evolutionary history, as well as their vulnerability or resilience under climate change. In the Southern Ocean, population genetic variation is influenced across multiple spatial scales ranging from circum-Antarctic, which encompasses the entire continent, to regional, with varying levels of geographic separation. However, comprehensive analyses testing the relative importance of different environmental and geographic variables on genomic variation across these scales are generally lacking in the Southern Ocean. Here, we examine genome-wide single nucleotide polymorphisms of the Southern Ocean octopus Pareledone turqueti across the Scotia Sea and the Antarctic continental shelf, at depths between 102 and 1342 m, throughout most of this species' range. The circumpolar distribution of P. turqueti is biogeographically structured with a clear signature of isolation-by-geographical distance, but with long-distance genetic connectivity also detected between East and West Antarctica. Genomic variation of P. turqueti was also associated with bottom water temperature at a circumpolar scale, driven by a genotype-temperature association with the warmer sub-Antarctic Shag Rocks and South Georgia. Within the Scotia Sea, geographic distance, oxygen and fine-scale isolation-by-water depth were apparent drivers of genomic variation at regional scales. Putative positive selection of haemocyanin (oxygen transport protein), calcium ion transport and genes linked to RNA modification, detected within the Scotia Sea, suggest physiological adaptation to the regional sharp temperature gradient (~0-+2°C). Overall, we identified seascape drivers of genomic variation in the Southern Ocean at circumpolar and regional scales in P. turqueti and contextualised the role of environmental adaptations in the Southern Ocean.

Whole-Genome Resequencing Reveals Polygenic Signatures of Directional and Balancing Selection on Alternative Migratory Life Histories

Migration in animals and associated adaptations to contrasting environments are underpinned by complex genetic architecture. Here, we explore the genomic basis of facultative anadromy in brown trout (Salmo trutta), wherein some individuals migrate to sea while others remain resident in natal rivers, to better understand how alternative migratory tactics (AMTs) are maintained evolutionarily. To identify genomic variants associated with AMTs, we sequenced whole genomes for 194 individual trout from five anadromous-resident population pairs, situated above and below waterfalls, in five different Irish rivers. These waterfalls act as natural barriers to upstream migration and hence we predicted that loci underpinning AMTs should be under similar divergent selection across these replicate pairs. A sliding windows based analysis revealed a highly polygenic adaptive divergence between anadromous and resident populations, encompassing 329 differentiated genomic regions. These regions were associated with 292 genes involved in various processes crucial for AMTs, including energy homeostasis, reproduction, osmoregulation, immunity, circadian rhythm and neural function. Furthermore, examining patterns of diversity we were able to link specific genes and biological processes to putative AMT trait classes: migratory-propensity, migratory-lifestyle and residency. Importantly, AMT outlier regions possessed higher genetic diversity than the background genome, particularly in the anadromous group, suggesting balancing selection may play a role in maintaining genetic variation. Overall, the results from this study provide important insights into the genetic architecture of migration and the evolutionary mechanisms shaping genomic diversity within and across populations.

Assessment of genome complementarity in three beef-on-dairy crossbreds reveals sire-specific effects on production traits with comparable rates of genomic inbreeding reduction

BACKGROUND

Crossbreeding beef bulls with dairy cows can improve the economic value and fitness of calves not entering dairy production owing to increased meat yield and heterosis. However, outcrossing might reduce the dosage of alleles that confer local adaptation or result in a higher risk of dystocia due to increased calf size. Given the clear phenotypic differences between beef breeds, the varying phylogenetic distances between beef and dairy breeds, and the genomic variations within breeds, the attainable economic and fitness gains of calves will strongly depend on the selection of sires for crossing. Thus, the aim of this study was to assess genome complementarity between Angus (AAN), Limousin (LIM), or Simmental (SIM) beef bulls and Brown Swiss (BSW) dairy cows by quantifying genomic inbreeding reduction in F1 crosses and identifying genes potentially under BSW-specific selection that might be affected by outcrossing.

RESULTS

Low-pass sequencing data from 181 cows, 34 bulls, and 301 of their F1 progeny, and body weight and carcass composition measurements of 248 F1s were obtained. The high genomic inbreeding levels detected in the BSW cows were substantially reduced in the crossbreds, with only minor differences between the sire breeds. In the BSW cows, 585 candidate genes under selection were identified, overrepresenting genes associated with milk, meat and carcass, and production traits. Only a few genes were strongly differentiated at nonsynonymous variants between the BSW and beef breeds, including four tightly clustered genes (FAM184B, NCAPG, DCAF16, and LCORL) nearly fixed for alternate alleles in the BSW cows but mostly heterozygous or homozygous for the reference alleles in the AAN and LIM bulls. The alternate allele dosage at these genes significantly correlated with reduced carcass weight and protein mass in F1s.

CONCLUSION

Some of the few genes that were highly divergent between the BSW and beef breeds at nonsynonymous variants were likely under strong selection for reduced carcass weight in the BSW breed, potentially due to trade-offs between beef and dairy productions. As alleles with opposing effects still segregate in beef cattle, marker-assisted selection of mating pairs may be used to modulate the desired phenotypes and simultaneously decrease genomic inbreeding.

Population Genomics of Adaptive Radiations: Exceptionally High Levels of Genetic Diversity and Recombination in an Endemic Spider From the Canary Islands

The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, with ~60 endemic species having originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterised by phenotypic modifications encompassing morphological, metabolic and behavioural changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low coverage-based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterised genome-wide levels of nucleotide polymorphism, divergence and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. We also identified several candidate genomic regions that are potentially under positive selection, highlighting relevant biological processes, such as vision and nitrogen extraction as potential adaptation targets. These processes may ultimately drive species diversification in this genus. This pioneering study of spiders that are endemic to an oceanic archipelago lays the groundwork for broader population genomics analyses aimed at understanding the genetic mechanisms driving adaptive radiation in island ecosystems.

Genomic and Acoustic Biogeography of the Iconic Sulphur-crested Cockatoo Clarifies Species Limits and Patterns of Intraspecific Diversity

Many highly recognizable species lack genetic data important for conservation due to neglect over their hyperabundance. This likely applies to the Sulfur-crested Cockatoo (Cacatua galerita), one of the world's most iconic parrots. The species is native to Australia, New Guinea, and some surrounding Melanesian islands of the latter. Four subspecies are currently recognised based on morphology. Australian subspecies and populations are abundant, but several factors threaten those in New Guinea and Melanesia. Genetic data from natural populations are scarce-information that is vital to identifying evolutionarily significant units (ESUs) important for modern conservation planning. We used whole-genome resequencing to investigate patterns of differentiation, evolutionary affinities, and demographic history across C. galerita's distribution range to assess whether currently recognised subspecies represent ESUs. We complement this with an assessment of bioacoustic variation across the species' distribution landscape. Our results point to C. galerita sensu lato (s.l.) comprising two species. We restrict C. galerita sensu stricto (s.s.) to populations in Australia and the Trans-Fly ecodomain of southern New Guinea. The second species, recognised here as Cacatua triton, likely occurs over much of the rest of New Guinea. Restricting further discussion of intraspecific diversity in C. triton, we show that within C. galerita s.s. two ESUs exist, which align to Cacatua galerita galerita in eastern Australia and southern New Guinea and Cacatua galerita fitzroyi in northern and north-western Australia. We suggest that the evolution of these species and ESUs are linked to Middle and Late Pleistocene glacial cycles and their effects on sea level and preferential habitats. We argue that conservation assessments need updating, protection of preferential forest and woodland habitats are important and reintroductions require careful management to avoid possible negative hybridization effects of non-complementary lineages.

Speciation in savanna birds in South America: The case of the Least Nighthawk Chordeiles pusillus (Aves: Caprimulgidae) in and out of the Amazon

The Least Nighthawk Chordeiles pusillus is widespread wherever there are savannas in the South American tropics, often in isolated patches, such as white-sands savannas in the Amazon rainforest realm. Here, we investigate genetic relationships between populations of the Least Nighthawk to understand historical processes leading to its diversification and to determine dispersal routes between northern and southern savannas by way of three hypothesized dispersal corridors by comparing samples from white-sand savannas to samples from other savannas outside of the Amazon rainforest region. We use 32 mtDNA samples from the range of C. pusillus to infer a dated phylogeny. In a subset of 17 samples, we use shotgun sequences to infer a distance-based phylogeny and to estimate individual admixture proportions. We calculate gene flow and shared alleles between white-sand and non-Amazonian populations using the ABBA-BABA test (D statistics), and Principal Component Analysis (PCA) to examine genetic structure within and between lineages. Finally, we use species distribution modelling (SDM) of conditions during the Last Glacial Maximum (LGM), currently, and in the future (2050-2080) to predict potential species occurrence under a climate change scenario. Two main clades (estimated to have diverged around 1.07 million years ago) were recovered with mtDNA sequences and Single Nucleotide Polymorphism (SNPs) and were supported by NGSadmix and PCA: one in the Amazon basin white-sand savannas, the other in the non-Amazonian savannas. Possible allele sharing between these clades was indicated by the D-statistics between northern non-Amazonian populations and the white-sand savanna population, but this was not corroborated by the admixture analyses. Dispersal among northern non-Amazonian populations may have occurred in a dry corridor between the Guianan and the Brazilian Shield, which has since moved eastward. Our data suggest that the lineages separated well before the Last Glacial Maximum, consequently dispersal could have happened at any earlier time during similar climatic conditions. Subsequently, non-Amazonian lineages became more divergent among themselves, possibly connecting and dispersing across the mouth of the Amazon River across Marajó island during favourable climatic conditions in the Pleistocene.

Discordance between taxonomy and population genomic data: An avian example relevant to the United States Endangered Species Act

Population genomics can reveal cryptic biological diversity that may impact fitness while simultaneously serving to delineate relevant conservation units. Here, we leverage the power of whole-genome resequencing for conservation by studying 433 individual lesser prairie-chicken (Tympanuchus pallidicinctus; LEPC, a federally endangered species of conservation concern in the United States) and greater prairie-chicken (Tympanuchus cupido; GRPC, a legally huntable species throughout much of its range). The genomic diversity of two formally recognized distinct population segments (DPSs) of LEPCs is similar, but they are genetically distinct. Neither DPS is depleted of its genomic diversity, neither is especially inbred, and temporal diversity is relatively stable in both conservation units. Interspecific differentiation between the two species was only slightly higher than that observed between LEPC DPSs, due largely to bidirectional introgression. The high resolution provided by our dataset identified a genomic continuum between the two species such that individuals sampled from the hybrid zone were imperfectly assigned to their presumptive species when considering only their physical characteristics. The admixture between the two species is reflected in the spectrum of individual ancestry coefficients, which has legal implications for the "take" of individuals under the Endangered Species Act. Overall, our data highlight the recurring dissonance between static policies and dynamic species boundaries that are increasingly obvious in the population genomic era.

Chromosomal inversions harbour excess mutational load in the coral, Acropora kenti, on the Great Barrier Reef

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future.

The Andes as a semi-permeable geographical barrier: Genetic connectivity between structured populations in a widespread spider

Geographical barriers like mountain ranges impede genetic exchange among populations, promoting diversification. The effectiveness of these barriers in limiting gene flow varies between lineages due to each species' dispersal modes and capacities. Our understanding of how the Andes orogeny contributes to species diversification comes from well-studied vertebrates and a few arthropods and plants, neglecting organisms unable to fly or walk long distances. Some arachnids, such as Gasteracantha cancriformis, have been hypothesized to disperse long distances via ballooning (i.e. using their silk to interact with the wind). Yet, we do not know how the environment and geography shape its genetic diversity. Therefore, we tested whether the Andes contributed to the diversification of G. cancriformis acting as an absolute or semi-permeable barrier to genetic connectivity between populations of this spider at opposite sides of the mountain range. We sampled thousands of loci across the distribution of the species and implemented population genetics, phylogenetic, and landscape genetic analyses. We identified two genetically distinct groups structured by the Central Andes, and a third less structured group in the Northern Andes that shares ancestry with the previous two. This structure is largely explained by the altitude along the Andes, which decreases in some regions, possibly facilitating cross-Andean dispersal and gene flow. Our findings support that altitude in the Andes plays a major role in structuring populations in South America, but the strength of this barrier can be overcome by organisms with long-distance dispersal modes together with altitudinal depressions.

Patterns of genetic divergence in the Rio Grande cooter (Pseudemys gorzugi), a riverine turtle inhabiting an arid and anthropogenically modified system

The lower Rio Grande and Pecos River of the southwest United States have been heavily modified by human activities, profoundly impacting the integrity of their aquatic wildlife. In this context, we focused our study on the population genomics of the Rio Grande Cooter (Pseudemys gorzugi), a freshwater turtle of increasing conservation concern, residing in these two rivers and their tributaries. The genetic data revealed two distinct populations: one in the Pecos and Black Rivers of New Mexico and another in the Rio Grande and Devils River of Texas, with admixed individuals identified at the confluence of the Rio Grande and Pecos River. In addition to having a smaller geographic range, we found lower observed heterozygosity, reduced nucleotide diversity, and a smaller effective population size (Ne) in New Mexico population. Our results depict a significant isolation-by-distance pattern across their distribution, with migration being notably infrequent at river confluences. These findings are pivotal for future conservation and restoration strategies, emphasizing the need to recognize the unique needs of each population.

High genetic load without purging in caribou, a diverse species at risk

High intra-specific genetic diversity is associated with adaptive potential, which is key for resilience to global change. However, high variation may also support deleterious alleles through genetic load, thereby increasing the risk of inbreeding depression if population sizes decrease. Purging of deleterious variation has been demonstrated in some threatened species. However, less is known about the costs of declines and inbreeding in species with large population sizes and high genetic diversity even though this encompasses many species globally that are expected to undergo population declines. Caribou is a species of ecological and cultural significance in North America with a wide distribution supporting extensive phenotypic variation but with some populations undergoing significant declines resulting in their at-risk status in Canada. We assessed intra-specific genetic variation, adaptive divergence, inbreeding, and genetic load across populations with different demographic histories using an annotated chromosome-scale reference genome and 66 whole-genome sequences. We found high genetic diversity and nine phylogenomic lineages across the continent with adaptive diversification of genes, but also high genetic load among lineages. We found highly divergent levels of inbreeding across individuals, including the loss of alleles by drift but not increased purging in inbred individuals, which had more homozygous deleterious alleles. We also found comparable frequencies of homozygous deleterious alleles between lineages regardless of nucleotide diversity. Thus, further inbreeding may need to be mitigated through conservation efforts. Our results highlight the "double-edged sword" of genetic diversity that may be representative of other species atrisk affected by anthropogenic activities.

Metagenomic analysis of the honey bee queen microbiome reveals low bacterial diversity and Caudoviricetes phages

In eusocial insects, the health of the queens-the colony founders and sole reproductive females-is a primary determinant for colony success. Queen failure in the honey bee Apis mellifera, for example, is a major concern of beekeepers who annually suffer colony losses, necessitating a greater knowledge of queen health. Several studies on the microbiome of honey bees have characterized its diversity and shown its importance for the health of worker bees, the female non-reproductive caste. However, the microbiome of workers differs from that of queens, which, in comparison, is still poorly studied. Thus, direct investigations of the queen microbiome are required to understand colony-level microbiome assembly, functional roles, and evolution. Here, we used metagenomics to comprehensively characterize the honey bee queen microbiome. Comparing samples from different geographic locations and breeder sources, we show that the microbiome of queens is mostly shaped by the environment experienced since early life and is predicted to play roles in the breakdown of the diet and protection from pathogens and xenobiotics. We also reveal that the microbiome of queens comprises only four candidate core bacterial species, Apilactobacillus kunkeei, Lactobacillus apis, Bombella apis, and Commensalibacter sp. Interestingly, in addition to bacteria, we show that bacteriophages infect the queen microbiome, for which Lactobacillaceae are predicted to be the main reservoirs. Together, our results provide the basis to understand the honey bee colony microbiome assemblage, can guide improvements in queen-rearing processes, and highlight the importance of considering bacteriophages for queen microbiome health and microbiome homeostasis in eusocial insects.IMPORTANCEThe queen caste plays a central role in colony success in eusocial insects, as queens lay eggs and regulate colony behavior and development. Queen failure can cause colonies to collapse, which is one of the major concerns of beekeepers. Thus, understanding the biology behind the queen's health is a pressing issue. Previous studies have shown that the bee microbiome plays an important role in worker bee health, but little is known about the queen microbiome and its function in vivo. Here, we characterized the queen microbiome, identifying for the first time the present species and their putative functions. We show that the queen microbiome has predicted nutritional and protective roles in queen association and comprises only four consistently present bacterial species. Additionally, we bring to attention the spread of phages in the queen microbiome, which increased in abundance in failing queens and may impact the fate of the colony.

Genomic evidence for West Antarctic Ice Sheet collapse during the Last Interglacial

The marine-based West Antarctic Ice Sheet (WAIS) is considered vulnerable to irreversible collapse under future climate trajectories, and its tipping point may lie within the mitigated warming scenarios of 1.5° to 2°C of the United Nations Paris Agreement. Knowledge of ice loss during similarly warm past climates could resolve this uncertainty, including the Last Interglacial when global sea levels were 5 to 10 meters higher than today and global average temperatures were 0.5° to 1.5°C warmer than preindustrial levels. Using a panel of genome-wide, single-nucleotide polymorphisms of a circum-Antarctic octopus, we show persistent, historic signals of gene flow only possible with complete WAIS collapse. Our results provide the first empirical evidence that the tipping point of WAIS loss could be reached even under stringent climate mitigation scenarios.

Microbiome composition is shaped by geography and population structure in the parasitic wasp Asobara japonica, but not in the presence of the endosymbiont Wolbachia

The microbial community composition is crucial for diverse life-history traits in many organisms. However, we still lack a sufficient understanding of how the host microbiome is acquired and maintained, a pressing issue in times of global environmental change. Here we investigated to what extent host genotype, environmental conditions, and the endosymbiont Wolbachia influence the bacterial communities in the parasitic wasp Asobara japonica. We sampled multiple wasp populations across 10 locations in their natural distribution range in Japan and sequenced the host genome (whole genome sequencing) and microbiome (16S rRNA gene). We compared the host population structure and bacterial community composition of wasps that reproduce sexually and are uninfected with Wolbachia with wasps that reproduce asexually and carry Wolbachia. The bacterial communities in asexual wasps were highly similar due to a strong effect of Wolbachia rather than host genomic structure. In contrast, in sexual wasps, bacterial communities appear primarily shaped by a combination of population structure and environmental conditions. Our research highlights that multiple factors shape the bacterial communities of an organism and that the presence of a single endosymbiont can strongly alter their compositions. This information is crucial to understanding how organisms and their associated microbiome will react in the face of environmental change.

Multiple paths toward repeated phenotypic evolution in the spiny-leg adaptive radiation (Tetragnatha; Hawai'i)

The repeated evolution of phenotypes provides clear evidence for the role of natural selection in driving evolutionary change. However, the evolutionary origin of repeated phenotypes can be difficult to disentangle as it can arise from a combination of factors such as gene flow, shared ancestral polymorphisms or mutation. Here, we investigate the presence of these evolutionary processes in the Hawaiian spiny-leg Tetragnatha adaptive radiation, which includes four microhabitat-specialists or ecomorphs, with different body pigmentation and size (Green, Large Brown, Maroon, and Small Brown). We investigated the evolutionary history of this radiation using 76 newly generated low-coverage, whole-genome resequenced samples, along with phylogenetic and population genomic tools. Considering the Green ecomorph as the ancestral state, our results suggest that the Green ecomorph likely re-evolved once, the Large Brown and Maroon ecomorphs evolved twice and the Small Brown evolved three times. We found that the evolution of the Maroon and Small Brown ecomorphs likely involved ancestral hybridization events, while the Green and Large Brown ecomorphs likely evolved through novel mutations, despite a high rate of incomplete lineage sorting in the dataset. Our findings demonstrate that the repeated evolution of ecomorphs in the Hawaiian spiny-leg Tetragnatha is influenced by multiple evolutionary processes.

North African fox genomes show signatures of repeated introgression and adaptation to life in deserts

Elucidating the evolutionary process of animal adaptation to deserts is key to understanding adaptive responses to climate change. Here we generated 82 individual whole genomes of four fox species (genus Vulpes) inhabiting the Sahara Desert at different evolutionary times. We show that adaptation of new colonizing species to a hot arid environment has probably been facilitated by introgression and trans-species polymorphisms shared with older desert resident species, including a putatively adaptive 25 Mb genomic region. Scans for signatures of selection implicated genes affecting temperature perception, non-renal water loss and heat production in the recent adaptation of North African red foxes (Vulpes vulpes), after divergence from Eurasian populations approximately 78 thousand years ago. In the extreme desert specialists, Rueppell's fox (V. rueppellii) and fennec (V. zerda), we identified repeated signatures of selection in genes affecting renal water homeostasis supported by gene expression and physiological differences. Our study provides insights into the mechanisms and genetic underpinnings of a natural experiment of repeated adaptation to extreme conditions.

Climate-induced range shifts drive adaptive response via spatio-temporal sieving of alleles

Quaternary climate fluctuations drove many species to shift their geographic ranges, in turn shaping their genetic structures. Recently, it has been argued that adaptation may have accompanied species range shifts via the "sieving" of genotypes during colonisation and establishment. However, this has not been directly demonstrated, and knowledge remains limited on how different evolutionary forces, which are typically investigated separately, interacted to jointly mediate species responses to past climatic change. Here, through whole-genome re-sequencing of over 1200 individuals of the carnation Dianthus sylvestris coupled with integrated population genomic and gene-environment models, we reconstruct the past neutral and adaptive landscape of this species as it was shaped by the Quaternary glacial cycles. We show that adaptive responses emerged concomitantly with the post-glacial range shifts and expansions of this species in the last 20 thousand years. This was due to the heterogenous sieving of adaptive alleles across space and time, as populations expanded out of restrictive glacial refugia into the broader and more heterogeneous range of habitats available in the present-day inter-glacial. Our findings reveal a tightly-linked interplay of migration and adaptation under past climate-induced range shifts, which we show is key to understanding the spatial patterns of adaptive variation we see in species today.

Epigenetic and Genetic Population Structure is Coupled in a Marine Invertebrate

Delineating the relative influence of genotype and the environment on DNA methylation is critical for characterizing the spectrum of organism fitness as driven by adaptation and phenotypic plasticity. In this study, we integrated genomic and DNA methylation data for two distinct Olympia oyster (Ostrea lurida) populations while controlling for within-generation environmental influences. In addition to providing the first characterization of genome-wide DNA methylation patterns in the oyster genus Ostrea, we identified 3,963 differentially methylated loci between populations. Our results show a clear coupling between genetic and epigenetic patterns of variation, with 27% of variation in interindividual methylation differences explained by genotype. Underlying this association are both direct genetic changes in CpGs (CpG-SNPs) and genetic variation with indirect influence on methylation (mQTLs). When comparing measures of genetic and epigenetic population divergence at specific genomic regions this relationship surprisingly breaks down, which has implications for the methods commonly used to study epigenetic and genetic coupling in marine invertebrates.

Contrasting Phylogeographic Patterns of Mitochondrial and Genome-Wide Variation in the Groundwater Amphipod Crangonyx islandicus That Survived the Ice Age in Iceland

The analysis of phylogeographic patterns has often been based on mitochondrial DNA variation, but recent analyses dealing with nuclear DNA have in some instances revealed mito-nuclear discordances and complex evolutionary histories. These enigmatic scenarios, which may involve stochastic lineage sorting, ancestral hybridization, past dispersal and secondary contacts, are increasingly scrutinized with a new generation of genomic tools such as RADseq, which also poses additional analytical challenges. Here, we revisited the previously inconclusive phylogeographic history, showing the mito-nuclear discordance of an endemic groundwater amphipod from Iceland, Crangonyx islandicus, which is the only metazoan known to have survived the Pleistocene beneath the glaciers. Previous studies based on three DNA markers documented a mitochondrial scenario with the main divergence occurring between populations in northern Iceland and an ITS scenario with the main divergence between the south and north. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to clarify this mito-nuclear discordance by applying several statistical methods while estimating the sensitivity to different analytical approaches (data-type, differentiation indices and base call uncertainty). A majority of nuclear markers and methods support the ITS divergence. Nevertheless, a more complex scenario emerges, possibly involving introgression led by male-biased dispersal among northern locations or mitochondrial capture, which may have been further strengthened by natural selection.

Sweepstakes reproductive success via pervasive and recurrent selective sweeps

Highly fecund natural populations characterized by high early mortality abound, yet our knowledge about their recruitment dynamics is somewhat rudimentary. This knowledge gap has implications for our understanding of genetic variation, population connectivity, local adaptation, and the resilience of highly fecund populations. The concept of sweepstakes reproductive success, which posits a considerable variance and skew in individual reproductive output, is key to understanding the distribution of individual reproductive success. However, it still needs to be determined whether highly fecund organisms reproduce through sweepstakes and, if they do, the relative roles of neutral and selective sweepstakes. Here, we use coalescent-based statistical analysis of population genomic data to show that selective sweepstakes likely explain recruitment dynamics in the highly fecund Atlantic cod. We show that the Kingman coalescent (modelling no sweepstakes) and the Xi-Beta coalescent (modelling random sweepstakes), including complex demography and background selection, do not provide an adequate fit for the data. The Durrett-Schweinsberg coalescent, in which selective sweepstakes result from recurrent and pervasive selective sweeps of new mutations, offers greater explanatory power. Our results show that models of sweepstakes reproduction and multiple-merger coalescents are relevant and necessary for understanding genetic diversity in highly fecund natural populations. These findings have fundamental implications for understanding the recruitment variation of fish stocks and general evolutionary genomics of high-fecundity organisms.

Whole-genome Analysis Reveals Contrasting Relationships Among Nuclear and Mitochondrial Genomes Between Three Sympatric Bat Species

Understanding mechanisms involved in speciation can be challenging, especially when hybridization or introgression blurs species boundaries. In bats, resolving relationships of some closely related groups has proved difficult due subtle interspecific variation both in morphometrics and molecular data sets. The endemic South American Histiotus bats, currently considered a subgenus of Eptesicus, harbor unresolved phylogenetic relationships and of those is a trio consisting of two closely related species: Eptesicus (Histiotus) macrotus and Eptesicus (Histiotus) montanus, and their relationship with a third, Eptesicus (Histiotus) magellanicus. The three sympatric species bear marked resemblance to each other, but can be differentiated morphologically. Furthermore, previous studies have been unable to differentiate the species from each other at a molecular level. In order to disentangle the phylogenetic relationships of these species, we examined the differentiation patterns and evolutionary history of the three Eptesicus (H.) species at the whole-genome level. The nuclear DNA statistics between the species suggest strong gene flow and recent hybridization between E. (H.) montanus and E. (H.) macrotus, whereas E. (H.) magellanicus shows a higher degree of isolation. In contrast, mitochondrial DNA shows a closer relationship between E. (H.) magellanicus and E. (H.) montanus. Opposing patterns in mtDNA and nuclear markers are often due to differences in dispersal, and here it could be both as a result of isolation in refugia during the last glacial maximum and female philopatry and male-biased dispersal. In conclusion, this study shows the importance of both the nuclear and mitochondrial DNA in resolving phylogenetic relationships and species histories.

Fidelity varies in the symbiosis between a gutless marine worm and its microbial consortium

BACKGROUND

Many animals live in intimate associations with a species-rich microbiome. A key factor in maintaining these beneficial associations is fidelity, defined as the stability of associations between hosts and their microbiota over multiple host generations. Fidelity has been well studied in terrestrial hosts, particularly insects, over longer macroevolutionary time. In contrast, little is known about fidelity in marine animals with species-rich microbiomes at short microevolutionary time scales, that is at the level of a single host population. Given that natural selection acts most directly on local populations, studies of microevolutionary partner fidelity are important for revealing the ecological and evolutionary processes that drive intimate beneficial associations within animal species.

RESULTS

In this study on the obligate symbiosis between the gutless marine annelid Olavius algarvensis and its consortium of seven co-occurring bacterial symbionts, we show that partner fidelity varies across symbiont species from strict to absent over short microevolutionary time. Using a low-coverage sequencing approach that has not yet been applied to microbial community analyses, we analysed the metagenomes of 80 O. algarvensis individuals from the Mediterranean and compared host mitochondrial and symbiont phylogenies based on single-nucleotide polymorphisms across genomes. Fidelity was highest for the two chemoautotrophic, sulphur-oxidizing symbionts that dominated the microbial consortium of all O. algarvensis individuals. In contrast, fidelity was only intermediate to absent in the sulphate-reducing and spirochaetal symbionts with lower abundance. These differences in fidelity are likely driven by both selective and stochastic forces acting on the consistency with which symbionts are vertically transmitted.

CONCLUSIONS

We hypothesize that variable degrees of fidelity are advantageous for O. algarvensis by allowing the faithful transmission of their nutritionally most important symbionts and flexibility in the acquisition of other symbionts that promote ecological plasticity in the acquisition of environmental resources. Video Abstract.

Speciation of pelagic zooplankton: Invisible boundaries can drive isolation of oceanic ctenophores

The study of evolution and speciation in non-model systems provides us with an opportunity to expand our understanding of biodiversity in nature. Connectivity studies generally focus on species with obvious boundaries to gene flow, but in open-ocean environments, such boundaries are difficult to identify. Due to the lack of obvious boundaries, speciation and population subdivision in the pelagic environment remain largely unexplained. Comb jellies (Phylum Ctenophora) are mostly planktonic gelatinous invertebrates, many of which are considered to have freely interbreeding distributions worldwide. It is thought that the lobate ctenophore Bolinopsis infundibulum is distributed throughout cooler northern latitudes and B. vitrea warmer. Here, we examined the global population structure for species of Bolinopsis with genetic and morphological data. We found distinct evolutionary patterns within the genus, where B. infundibulum had a broad distribution from northern Pacific to Atlantic waters despite many physical barriers, while other species were geographically segregated despite few barriers. Divergent patterns of speciation within the genus suggest that oceanic currents, sea-level, and geological changes over time can act as either barriers or aids to dispersal in the pelagic environment. Further, we used population genomic data to examine evolution in the open ocean of a distinct lineage of Bolinopsis ctenophores from the North Eastern Pacific. Genetic information and morphological observations validated this as a separate species, Bolinopsis microptera, which was previously described but has recently been called B. infundibulum. We found that populations of B. microptera from California were in cytonuclear discordance, which indicates a secondary contact zone for previously isolated populations. Discordance at this scale is rare, especially in a continuous setting.

The Genomic Landscapes of Desert Birds Form over Multiple Time Scales

Spatial models show that genetic differentiation between populations can be explained by factors ranging from geographic distance to environmental resistance across the landscape. However, genomes exhibit a landscape of differentiation, indicating that multiple processes may mediate divergence in different portions of the genome. We tested this idea by comparing alternative geographic predctors of differentiation in ten bird species that co-occur in Sonoran and Chihuahuan Deserts of North America. Using population-level genomic data, we described the genomic landscapes across species and modeled conditions that represented historical and contemporary mechanisms. The characteristics of genomic landscapes differed across species, influenced by varying levels of population structuring and admixture between deserts, and the best-fit models contrasted between the whole genome and partitions along the genome. Both historical and contemporary mechanisms were important in explaining genetic distance, but particularly past and current environments, suggesting that genomic evolution was modulated by climate and habitat There were also different best-ftit models across genomic partitions of the data, indicating that these regions capture different evolutionary histories. These results show that the genomic landscape of differentiation can be associated with alternative geographic factors operating on different portions of the genome, which reflect how heterogeneous patterns of genetic differentiation can evolve across species and genomes.

Phylogenomics reveals extensive introgression and a case of mito-nuclear discordance in the killifish genus Kryptolebias

Introgression is a widespread evolutionary process leading to phylogenetic inconsistencies among distinct parts of the genomes, particularly between mitochondrial and nuclear-based phylogenetic reconstructions (e.g., mito-nuclear discordances). Here, we used mtDNA and genome-wide nuclear sites to provide the first phylogenomic-based hypothesis on the evolutionary relationships within the killifish genus Kryptolebias. In addition, we tested for evidence of past introgression in the genus given the multiple reports of undergoing hybridization between its members. Our mtDNA phylogeny generally agreed with the relationships previously proposed for the genus. However, our reconstruction based on nuclear DNA revealed an unknown lineage - Kryptolebias sp. 'ESP' - as the sister group of the self-fertilizing mangrove killifishes, K. marmoratus and K. hermaphroditus. All individuals sequenced of Kryptolebias sp. 'ESP' had the same mtDNA haplotype commonly observed in K. hermaphroditus, demonstrating a clear case of mito-nuclear discordance. Our analysis further confirmed extensive history of introgression between Kryptolebias sp. 'ESP' and K. hermaphroditus. Population genomics analyses indicate no current gene flow between the two lineages, despite their current sympatry and history of introgression. We also confirmed introgression between other species pairs in the genus that have been recently reported to form hybrid zones. Overall, our study provides a phylogenomic reconstruction covering most of the Kryptolebias species, reveals a new lineage hidden in a case of mito-nuclear discordance, and provides evidence of multiple events of ancestral introgression in the genus. These findings underscore the importance of investigating different genomic information in a phylogenetic framework, particularly in taxa where introgression is common as in the sexually diverse mangrove killifishes.

Population Genomics of Variegated Toad-Headed Lizard Phrynocephalus versicolor and Its Adaptation to the Colorful Sand of the Gobi Desert

The variegated toad-headed agama, Phrynocephalus versicolor, lives in the arid landscape of the Chinese Gobi Desert. We analyzed populations from three different locations which vary in substrate color and altitude: Heishankou (HSK), Guazhou County (GZ), and Ejin Banner (EJN). The substrate color is either light-yellow (GZ-y), yellow (EJN-y), or black (HSK-b); the corresponding lizard population colors largely match their substrate in the degree of melanism. We assembled the P. versicolor genome and sequenced over 90 individuals from the three different populations. Genetic divergence between populations corresponds to their geographic distribution. We inferred the genetic relationships among these populations and used selection scans and differential expression to identify genes that show signatures of selection. Slc2a11 and akap12, among other genes, are highly differentiated and may be responsible for pigment adaptation to substrate color in P. versicolor.

distAngsd: Fast and accurate inference of genetic distances for Next Generation Sequencing data

Commonly used methods for inferring phylogenies were designed before the emergence of high throughput sequencing and can generally not accommodate the challenges associated with noisy, diploid sequencing data. In many applications, diploid genomes are still treated as haploid through the use of ambiguity characters; while the uncertainty in genotype calling - arising as a consequence of the sequencing technology - is ignored. In order to address this problem we describe two new probabilistic approaches for estimating genetic distances: distAngsd-geno and distAngsd-nuc, both implemented in a software suite named distAngsd. These methods are specifically designed for next generation sequencing data, utilize the full information from the data, and take uncertainty in genotype calling into account. Through extensive simulations, we show that these new methods are markedly more accurate and have more stable statistical behaviors than other currently available methods for estimating genetic distances - even for very low depth data with high error rates.

Population dynamics and genetic connectivity in recent chimpanzee history

Knowledge on the population history of endangered species is critical for conservation, but whole-genome data on chimpanzees (Pan troglodytes) is geographically sparse. Here, we produced the first non-invasive geolocalized catalog of genomic diversity by capturing chromosome 21 from 828 non-invasive samples collected at 48 sampling sites across Africa. The four recognized subspecies show clear genetic differentiation correlating with known barriers, while previously undescribed genetic exchange suggests that these have been permeable on a local scale. We obtained a detailed reconstruction of population stratification and fine-scale patterns of isolation, migration, and connectivity, including a comprehensive picture of admixture with bonobos (Pan paniscus). Unlike humans, chimpanzees did not experience extended episodes of long-distance migrations, which might have limited cultural transmission. Finally, based on local rare variation, we implement a fine-grained geolocalization approach demonstrating improved precision in determining the origin of confiscated chimpanzees.

Fast and accurate estimation of multidimensional site frequency spectra from low-coverage high-throughput sequencing data

BACKGROUND

The site frequency spectrum summarizes the distribution of allele frequencies throughout the genome, and it is widely used as a summary statistic to infer demographic parameters and to detect signals of natural selection. The use of high-throughput low-coverage DNA sequencing data can lead to biased estimates of the site frequency spectrum due to high levels of uncertainty in genotyping.

RESULTS

Here we design and implement a method to efficiently and accurately estimate the multidimensional joint site frequency spectrum for large numbers of haploid or diploid individuals across an arbitrary number of populations, using low-coverage sequencing data. The method maximizes a likelihood function that represents the probability of the sequencing data observed given a multidimensional site frequency spectrum using genotype likelihoods. Notably, it uses an advanced binning heuristic paired with an accelerated expectation-maximization algorithm for a fast and memory-efficient computation, and can generate both unfolded and folded spectra and bootstrapped replicates for haploid and diploid genomes. On the basis of extensive simulations, we show that the new method requires remarkably less storage and is faster than previous implementations whilst retaining the same accuracy. When applied to low-coverage sequencing data from the fungal pathogen Neonectria neomacrospora, results recapitulate the patterns of population differentiation generated using the original high-coverage data.

CONCLUSION

The new implementation allows for accurate estimation of population genetic parameters from arbitrarily large, low-coverage datasets, thus facilitating cost-effective sequencing experiments in model and non-model organisms.

Population structure in Neotropical plants: Integrating pollination biology, topography and climatic niches

Animal pollinators mediate gene flow among plant populations, but in contrast to well-studied topographic and (Pleistocene) environmental isolating barriers, their impact on population genetic differentiation remains largely unexplored. Comparing how these multifarious factors drive microevolutionary histories is, however, crucial for better resolving macroevolutionary patterns of plant diversification. Here we combined genomic analyses with landscape genetics and niche modelling across six related Neotropical plant species (424 individuals across 33 localities) differing in pollination strategy to test the hypothesis that highly mobile (vertebrate) pollinators more effectively link isolated localities than less mobile (bee) pollinators. We found consistently higher genetic differentiation (FST ) among localities of bee- than vertebrate-pollinated species with increasing geographical distance, topographic barriers and historical climatic instability. High admixture among montane populations further suggested relative climatic stability of Neotropical montane forests during the Pleistocene. Overall, our results indicate that pollinators may differentially impact the potential for allopatric speciation, thereby critically influencing diversification histories at macroevolutionary scales.

Population Structure and Diversity in European Honey Bees (Apismellifera L.)-An Empirical Comparison of Pool and Individual Whole-Genome Sequencing

BACKGROUND

Whole-genome sequencing has become routine for population genetic studies. Sequencing of individuals provides maximal data but is rather expensive and fewer samples can be studied. In contrast, sequencing a pool of samples (pool-seq) can provide sufficient data, while presenting less of an economic challenge. Few studies have compared the two approaches to infer population genetic structure and diversity in real datasets. Here, we apply individual sequencing (ind-seq) and pool-seq to the study of Western honey bees (Apis mellifera).

METHODS

We collected honey bee workers that belonged to 14 populations, including 13 subspecies, totaling 1347 colonies, who were individually (139 individuals) and pool-sequenced (14 pools). We compared allele frequencies, genetic diversity estimates, and population structure as inferred by the two approaches.

RESULTS

Pool-seq and ind-seq revealed near identical population structure and genetic diversities, albeit at different costs. While pool-seq provides genome-wide polymorphism data at considerably lower costs, ind-seq can provide additional information, including the identification of population substructures, hybridization, or individual outliers.

CONCLUSIONS

If costs are not the limiting factor, we recommend using ind-seq, as population genetic structure can be inferred similarly well, with the advantage gained from individual genetic information. Not least, it also significantly reduces the effort required for the collection of numerous samples and their further processing in the laboratory.

The evolutionary pathways for local adaptation in mountain hares

Understanding the evolution of local adaptations is a central aim of evolutionary biology and key for the identification of unique populations and lineages of conservation relevance. By combining RAD sequencing and whole‐genome sequencing, we identify genetic signatures of local adaptation in mountain hares (Lepus timidus) from isolated and distinctive habitats of its wide distribution: Ireland, the Alps and Fennoscandia. Demographic modelling suggested that the split of these mountain hares occurred around 20 thousand years ago, providing the opportunity to study adaptive evolution over a short timescale. Using genome‐wide scans, we identified signatures of extreme differentiation among hares from distinct geographic areas that overlap with area‐specific selective sweeps, suggesting targets for local adaptation. Several identified candidate genes are associated with traits related to the uniqueness of the different environments inhabited by the three groups of mountain hares, including coat colour, ability to live at high altitudes and variation in body size. In Irish mountain hares, a variant of ASIP, a gene previously implicated in introgression‐driven winter coat colour variation in mountain and snowshoe hares (L. americanus), may underlie brown winter coats, reinforcing the repeated nature of evolution at ASIP moulding adaptive seasonal colouration. Comparative genomic analyses across several hare species suggested that mountain hares’ adaptive variants appear predominantly species‐specific. However, using coalescent simulations, we also show instances where the candidate adaptive variants have been introduced via introgressive hybridization. Our study shows that standing adaptive variation, including that introgressed from other species, was a crucial component of the post‐glacial dynamics of species.

A holobiont view of island biogeography: Unravelling patterns driving the nascent diversification of a Hawaiian spider and its microbial associates

The diversification of a host lineage can be influenced by both the external environment and its assemblage of microbes. Here, we use a young lineage of spiders, distributed along a chronologically arranged series of volcanic mountains, to investigate how their associated microbial communities have changed as the spiders colonized new locations. Using the stick spider Ariamnes waikula (Araneae, Theridiidae) on the island of Hawai'i, and outgroup taxa on older islands, we tested whether each component of the "holobiont" (spider hosts, intracellular endosymbionts and gut microbial communities) showed correlated signatures of diversity due to sequential colonization from older to younger volcanoes. To investigate this, we generated ddRAD data for the host spiders and 16S rRNA gene amplicon data from their microbiota. We expected sequential colonizations to result in a (phylo)genetic structuring of the host spiders and in a diversity gradient in microbial communities. The results showed that the host A. waikula is indeed structured by geographical isolation, suggesting sequential colonization from older to younger volcanoes. Similarly, the endosymbiont communities were markedly different between Ariamnes species on different islands, but more homogeneous among A. waikula populations on the island of Hawai'i. Conversely, the gut microbiota, which we suspect is generally environmentally derived, was largely conserved across all populations and species. Our results show that different components of the holobiont respond in distinct ways to the dynamic environment of the volcanic archipelago. This highlights the necessity of understanding the interplay between different components of the holobiont, to properly characterize its evolution.

Pervasive hybridization with local wild relatives in Western European grapevine varieties

Grapevine (Vitis vinifera L.) diversity richness results from a complex domestication history over multiple historical periods. Here, we used whole-genome resequencing to elucidate different aspects of its recent evolutionary history. Our results support a model in which a central domestication event in grapevine was followed by postdomestication hybridization with local wild genotypes, leading to the presence of an introgression signature in modern wine varieties across Western Europe. The strongest signal was associated with a subset of Iberian grapevine varieties showing large introgression tracts. We targeted this study group for further analysis, demonstrating how regions under selection in wild populations from the Iberian Peninsula were preferentially passed on to the cultivated varieties by gene flow. Examination of underlying genes suggests that environmental adaptation played a fundamental role in both the evolution of wild genotypes and the outcome of hybridization with cultivated varieties, supporting a case of adaptive introgression in grapevine.

The genetic impact of an Ebola outbreak on a wild gorilla population

Background

Numerous Ebola virus outbreaks have occurred in Equatorial Africa over the past decades. Besides human fatalities, gorillas and chimpanzees have also succumbed to the fatal virus. The 2004 outbreak at the Odzala-Kokoua National Park (Republic of Congo) alone caused a severe decline in the resident western lowland gorilla (Gorilla gorilla gorilla) population, with a 95% mortality rate. Here, we explore the immediate genetic impact of the Ebola outbreak in the western lowland gorilla population.

Results

Associations with survivorship were evaluated by utilizing DNA obtained from fecal samples from 16 gorilla individuals declared missing after the outbreak (non-survivors) and 15 individuals observed before and after the epidemic (survivors). We used a target enrichment approach to capture the sequences of 123 genes previously associated with immunology and Ebola virus resistance and additionally analyzed the gut microbiome which could influence the survival after an infection. Our results indicate no changes in the population genetic diversity before and after the Ebola outbreak, and no significant differences in microbial community composition between survivors and non-survivors. However, and despite the low power for an association analysis, we do detect six nominally significant missense mutations in four genes that might be candidate variants associated with an increased chance of survival.

Conclusion

This study offers the first insight to the genetics of a wild great ape population before and after an Ebola outbreak using target capture experiments from fecal samples, and presents a list of candidate loci that may have facilitated their survival.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08025-y.

Against the Odds: Hybrid Zones between Mangrove Killifish Species with Different Mating Systems

Different mating systems are expected to affect the extent and direction of hybridization. Due to the different levels of sexual conflict, the weak inbreeder/strong outbreeder (WISO) hypothesis predicts that gametes from self-incompatible (SI) species should outcompete gametes from self-compatible (SC) ones. However, other factors such as timing of selfing and unilateral incompatibilities may also play a role on the direction of hybridization. In addition, differential mating opportunities provided by different mating systems are also expected to affect the direction of introgression in hybrid zones involving outcrossers and selfers. Here, we explored these hypotheses with a unique case of recent hybridization between two mangrove killifish species with different mating systems, Kryptolebias ocellatus (obligately outcrossing) and K. hermaphroditus (predominantly self-fertilizing) in two hybrid zones in southeast Brazil. Hybridization rates were relatively high (~20%), representing the first example of natural hybridization between species with different mating systems in vertebrates. All F1 individuals were sired by the selfing species. Backcrossing was small, but mostly asymmetrical with the SI parental species, suggesting pattern commonly observed in plant hybrid zones with different mating systems. Our findings shed light on how contrasting mating systems may affect the direction and extent of gene flow between sympatric species, ultimately affecting the evolution and maintenance of hybrid zones.

Museomics dissects the genetic basis for adaptive seasonal colouration in the least weasel

Dissecting the link between genetic variation and adaptive phenotypes provides outstanding opportunities to understand fundamental evolutionary processes. Here, we use a museomics approach to investigate the genetic basis and evolution of winter coat colouration morphs in least weasels (Mustela nivalis), a repeated adaptation for camouflage in mammals with seasonal pelage colour moults across regions with varying winter snow. Whole-genome sequence data was obtained from biological collections and mapped onto a newly assembled reference genome for the species. Sampling represented two replicate transition zones between nivalis and vulgaris colouration morphs in Europe, which typically develop white or brown winter coats, respectively. Population analyses showed that the morph distribution across transition zones is not a by-product of historical structure. Association scans linked a 200 kb genomic region to colouration morph, which was validated by genotyping museum specimens from inter-morph experimental crosses. Genotyping the wild populations narrowed down the association to pigmentation gene MC1R and pinpointed a candidate amino acid change co-segregating with colouration morph. This polymorphism replaces an ancestral leucine residue by lysine at the start of the first extracellular loop of the protein in the vulgaris morph. A selective sweep signature overlapped the association region in vulgaris, suggesting that past adaptation favoured winter-brown morphs and can anchor future adaptive responses to decreasing winter snow. Using biological collections as valuable resources to study natural adaptations, our study showed a new evolutionary route generating winter colour variation in mammals and that seasonal camouflage can be modulated by changes at single key genes.

Whole-genome analysis of giraffe supports four distinct species

Species is the fundamental taxonomic unit in biology and its delimitation has implications for conservation. In giraffe (Giraffa spp.), multiple taxonomic classifications have been proposed since the early 1900s.¹ However, one species with nine subspecies has been generally accepted,² likely due to limited in-depth assessments, subspecies hybridizing in captivity,^3,4 and anecdotal reports of hybrids in the wild.⁵ Giraffe taxonomy received new attention after population genetic studies using traditional genetic markers suggested at least four species.^6,7 This view has been met with controversy,⁸ setting the stage for debate.^9,10 Genomics is significantly enhancing our understanding of biodiversity and speciation relative to traditional genetic approaches and thus has important implications for species delineation and conservation.¹¹ We present a high-quality de novo genome assembly of the critically endangered Kordofan giraffe (G. camelopardalis antiquorum)¹² and a comprehensive whole-genome analysis of 50 giraffe representing all traditionally recognized subspecies. Population structure and phylogenomic analyses support four separately evolving giraffe lineages, which diverged 230-370 ka ago. These lineages underwent distinct demographic histories and show different levels of heterozygosity and inbreeding. Our results strengthen previous findings of limited gene flow and admixture among putative giraffe species^6,7,9 and establish a genomic foundation for recognizing four species and seven subspecies, the latter of which should be considered as evolutionary significant units. Achieving a consensus over the number of species and subspecies in giraffe is essential for adequately assessing their threat level and will improve conservation efforts for these iconic taxa.

Pervasive Genomic Signatures of Local Adaptation to Altitude Across Highland Specialist Andean Hummingbird Populations

Populations along steep environmental gradients are subject to differentiating selection that can result in local adaptation, despite countervailing gene flow, and genetic drift. In montane systems, where species are often restricted to narrow ranges of elevation, it is unclear whether the selection is strong enough to influence functional differentiation of subpopulations differing by a few hundred meters in elevation. We used targeted capture of 12 501 exons from across the genome, including 271 genes previously implicated in altitude adaptation, to test for adaptation to local elevations for 2 highland hummingbird species, Coeligena violifer (n = 62) and Colibri coruscans (n = 101). For each species, we described population genetic structure across the complex geography of the Peruvian Andes and, while accounting for this structure, we tested whether elevational allele frequency clines in single nucleotide polymorphisms (SNPs) showed evidence for local adaptation to elevation. Although the 2 species exhibited contrasting population genetic structures, we found signatures of clinal genetic variation with shifts in elevation in both. The genes with SNP-elevation associations included candidate genes previously discovered for high-elevation adaptation as well as others not previously identified, with cellular functions related to hypoxia response, energy metabolism, and immune function, among others. Despite the homogenizing effects of gene flow and genetic drift, natural selection on parts of the genome evidently optimizes elevation-specific cellular function even within elevation range-restricted montane populations. Consequently, our results suggest local adaptation occurring in narrow elevation bands in tropical mountains, such as the Andes, may effectively make them "taller" biogeographic barriers.

Darwin's Fancy Revised: An Updated Understanding of the Genomic Constitution of Pigeon Breeds

Through its long history of artificial selection, the rock pigeon (Columba livia Gmelin 1789) was forged into a large number of domestic breeds. The incredible amount of phenotypic diversity exhibited in these breeds has long held the fascination of scholars, particularly those interested in biological inheritance and evolution. However, exploiting them as a model system is challenging, as unlike with many other domestic species, few reliable records exist about the origins of, and relationships between, each of the breeds. Therefore, in order to broaden our understanding of the complex evolutionary relationships among pigeon breeds, we generated genome-wide data by performing the genotyping-by-sequencing (GBS) method on close to 200 domestic individuals representing over 60 breeds. We analyzed these GBS data alongside previously published whole-genome sequencing data, and this combined analysis allowed us to conduct the most extensive phylogenetic analysis of the group, including two feral pigeons and one outgroup. We improve previous phylogenies, find considerable population structure across the different breeds, and identify unreported interbreed admixture events. Despite the reduced number of loci relative to whole-genome sequencing, we demonstrate that GBS data provide sufficient analytical power to investigate intertwined evolutionary relationships, such as those that are characteristic of animal domestic breeds. Thus, we argue that future studies should consider sequencing methods akin to the GBS approach as an optimal cost-effective approach for addressing complex phylogenies.

A reference-free approach to analyse RADseq data using standard next generation sequencing toolkits

Genotyping-by-sequencing methods such as RADseq are popular for generating genomic and population-scale data sets from a diverse range of organisms. These often lack a usable reference genome, restricting users to RADseq specific software for processing. However, these come with limitations compared to generic next generation sequencing (NGS) toolkits. Here, we describe and test a simple pipeline for reference-free RADseq data processing that blends de novo elements from STACKS with the full suite of state-of-the art NGS tools. Specifically, we use the de novo RADseq assembly employed by STACKS to create a catalogue of RAD loci that serves as a reference for read mapping, variant calling and site filters. Using RADseq data from 28 zebra sequenced to ~8x depth-of-coverage we evaluate our approach by comparing the site frequency spectra (SFS) to those from alternative pipelines. Most pipelines yielded similar SFS at 8x depth, but only a genotype likelihood based pipeline performed similarly at low sequencing depth (2-4x). We compared the RADseq SFS with medium-depth (~13x) shotgun sequencing of eight overlapping samples, revealing that the RADseq SFS was persistently slightly skewed towards rare and invariant alleles. Using simulations and human data we confirm that this is expected when there is allelic dropout (AD) in the RADseq data. AD in the RADseq data caused a heterozygosity deficit of ~16%, which dropped to ~5% after filtering AD. Hence, AD was the most important source of bias in our RADseq data.

Whole exome sequencing identifies the potential for genetic rescue in iconic and critically endangered Panamanian harlequin frogs

Avoiding extinction in a rapidly changing environment often relies on a species' ability to quickly adapt in the face of extreme selective pressures. In Panamá, two closely related harlequin frog species (Atelopus varius and Atelopus zeteki) are threatened with extinction due to the fungal pathogen Batrachochytrium dendrobatidis (Bd). Once thought to be nearly extirpated from Panamá, A. varius have recently been rediscovered in multiple localities across their historical range; however, A. zeteki are possibly extinct in the wild. By leveraging a unique collection of 186 Atelopus tissue samples collected before and after the Bd outbreak in Panama, we describe the genetics of persistence for these species on the brink of extinction. We sequenced the transcriptome and developed an exome-capture assay to sequence the coding regions of the Atelopus genome. Using these genetic data, we evaluate the population genetic structure of historical A. varius and A. zeteki populations, describe changes in genetic diversity over time, assess the relationship between contemporary and historical individuals, and test the hypothesis that some A. varius populations have rapidly evolved to resist or tolerate Bd infection. We found a significant decrease in genetic diversity in contemporary (compared to historical) A. varius populations. We did not find strong evidence of directional allele frequency change or selection for Bd resistance genes, but we uncovered a set of candidate genes that warrant further study. Additionally, we found preliminary evidence of recent migration and gene flow in one of the largest persisting A. varius populations in Panamá, suggesting the potential for genetic rescue in this system. Finally, we propose that previous conservation units should be modified, as clear genetic breaks do not exist beyond the local population level. Our data lay the groundwork for genetically informed conservation and advance our understanding of how imperiled species might be rescued from extinction.

Evaluating insect-host interactions as a driver of species divergence in palm flower weevils

Plants and their specialized flower visitors provide valuable insights into the evolutionary consequences of species interactions. In particular, antagonistic interactions between insects and plants have often been invoked as a major driver of diversification. Here we use a tropical community of palms and their specialized insect flower visitors to test whether antagonisms lead to higher population divergence. Interactions between palms and the insects visiting their flowers range from brood pollination to florivory and commensalism, with the latter being species that feed on decaying-and presumably undefended-plant tissues. We test the role of insect-host interactions in the early stages of diversification of nine species of beetles sharing host plants and geographical ranges by first delimiting cryptic species and then using models of genetic isolation by environment. The degree to which insect populations are structured by the genetic divergence of plant populations varies. A hierarchical model reveals that this variation is largely uncorrelated with the kind of interaction, showing that antagonistic interactions are not associated with higher genetic differentiation. Other aspects of host use that affect plant-associated insects regardless of the outcomes of their interactions, such as sensory biases, are likely more general drivers of insect population divergence.

Divergent Selection and Primary Gene Flow Shape Incipient Speciation of a Riparian Tree on Hawaii Island

A long-standing goal of evolutionary biology is to understand the mechanisms underlying the formation of species. Of particular interest is whether or not speciation can occur in the presence of gene flow and without a period of physical isolation. Here, we investigated this process within Hawaiian Metrosideros, a hypervariable and highly dispersible woody species complex that dominates the Hawaiian Islands in continuous stands. Specifically, we investigated the origin of Metrosideros polymorpha var. newellii (newellii), a riparian ecotype endemic to Hawaii Island that is purportedly derived from the archipelago-wide M. polymorpha var. glaberrima (glaberrima). Disruptive selection across a sharp forest-riparian ecotone contributes to the isolation of these varieties and is a likely driver of newellii's origin. We examined genome-wide variation of 42 trees from Hawaii Island and older islands. Results revealed a split between glaberrima and newellii within the past 0.3-1.2 My. Admixture was extensive between lineages within Hawaii Island and between islands, but introgression from populations on older islands (i.e., secondary gene flow) did not appear to contribute to the emergence of newellii. In contrast, recurrent gene flow (i.e., primary gene flow) between glaberrima and newellii contributed to the formation of genomic islands of elevated absolute and relative divergence. These regions were enriched for genes with regulatory functions as well as for signals of positive selection, especially in newellii, consistent with divergent selection underlying their formation. In sum, our results support riparian newellii as a rare case of incipient ecological speciation with primary gene flow in trees.

Divergence, gene flow, and the origin of leapfrog geographic distributions: The history of colour pattern variation in Phyllobates poison-dart frogs

The geographic distribution of phenotypic variation among closely related populations is a valuable source of information about the evolutionary processes that generate and maintain biodiversity. Leapfrog distributions, in which phenotypically similar populations are disjunctly distributed and separated by one or more phenotypically distinct populations, represent geographic replicates for the existence of a phenotype, and are therefore especially informative. These geographic patterns have mostly been studied from phylogenetic perspectives to understand how common ancestry and divergent evolution drive their formation. Other processes, such as gene flow between populations, have not received as much attention. Here, we investigate the roles of divergence and gene flow between populations in the origin and maintenance of a leapfrog distribution in Phyllobates poison frogs. We found evidence for high levels of gene flow between neighbouring populations but not over long distances, indicating that gene flow between populations exhibiting the central phenotype may have a homogenizing effect that maintains their similarity, and that introgression between 'leapfroging' taxa has not played a prominent role as a driver of phenotypic diversity in Phyllobates. Although phylogenetic analyses suggest that the leapfrog distribution was formed through independent evolution of the peripheral (i.e. leapfrogging) populations, the elevated levels of gene flow between geographically close populations poise alternative scenarios, such as the history of phenotypic change becoming decoupled from genome-averaged patterns of divergence, which we cannot rule out. These results highlight the importance of incorporating gene flow between populations into the study of geographic variation in phenotypes, both as a driver of phenotypic diversity and as a confounding factor of phylogeographic inferences.

Combining whole-genome shotgun sequencing and rRNA gene amplicon analyses to improve detection of microbe-microbe interaction networks in plant leaves

Microorganisms from all domains of life establish associations with plants. Although some harm the plant, others antagonize pathogens or prime the plant immune system, support the acquisition of nutrients, tune plant hormone levels, or perform additional services. Most culture-independent plant microbiome research has focused on amplicon sequencing of the 16S rRNA gene and/or the internal transcribed spacer (ITS) of rRNA genomic loci, which show the relative abundance of the microbes to each other. Here, we describe shotgun sequencing of 275 wild Arabidopsis thaliana leaf microbiomes from southwest Germany, with additional bacterial 16S and eukaryotic ITS1 rRNA amplicon data from 176 of these samples. Shotgun data, which unlike the amplicon data capture the ratio of microbe to plant DNA, enable scaling of microbial read abundances to reflect the microbial load on the host. In a more cost-effective hybrid strategy, we show they also allow a similar scaling of amplicon data to overcome compositionality problems. Our wild plants were dominated by bacterial sequences, with eukaryotes contributing only a minority of reads. Microbial membership showed weak associations with both site of origin and plant genotype, both of which were highly confounded in this dataset. There was large variation among microbiomes, with one extreme comprising samples of low complexity and a high load of microorganisms typical of infected plants, and the other extreme being samples of high complexity and a low microbial load. Critically, considering absolute microbial load led to fundamentally different conclusions about microbiome assembly and the interaction networks among major taxa.

Natural variations at the Stay-Green gene promoter control lifespan and yield in rice cultivars

Increased grain yield will be critical to meet the growing demand for food, and could be achieved by delaying crop senescence. Here, via quantitative trait locus (QTL) mapping, we uncover the genetic basis underlying distinct life cycles and senescence patterns of two rice subspecies, indica and japonica. Promoter variations in the Stay-Green (OsSGR) gene encoding the chlorophyll-degrading Mg⁺⁺-dechelatase were found to trigger higher and earlier induction of OsSGR in indica, which accelerated senescence of indica rice cultivars. The indica-type promoter is present in a progenitor subspecies O. nivara and thus was acquired early during the evolution of rapid cycling trait in rice subspecies. Japonica OsSGR alleles introgressed into indica-type cultivars in Korean rice fields lead to delayed senescence, with increased grain yield and enhanced photosynthetic competence. Taken together, these data establish that naturally occurring OsSGR promoter and related lifespan variations can be exploited in breeding programs to augment rice yield.

The origin of domestication genes in goats

Goat domestication was critical for agriculture and civilization, but its underlying genetic changes and selection regimes remain unclear. Here, we analyze the genomes of worldwide domestic goats, wild caprid species, and historical remains, providing evidence of an ancient introgression event from a West Caucasian tur-like species to the ancestor of domestic goats. One introgressed locus with a strong signature of selection harbors the MUC6 gene, which encodes a gastrointestinally secreted mucin. Experiments revealed that the nearly fixed introgressed haplotype confers enhanced immune resistance to gastrointestinal pathogens. Another locus with a strong signal of selection may be related to behavior. The selected alleles at these two loci emerged in domestic goats at least 7200 and 8100 years ago, respectively, and increased to high frequencies concurrent with the expansion of the ubiquitous modern mitochondrial haplogroup A. Tracking these archaeologically cryptic evolutionary transformations provides new insights into the mechanisms of animal domestication.