The Genomic Landscape of Divergence Across the Speciation Continuum in Island-Colonising Silvereyes (Zosterops lateralis)

Genomic Landscape of Divergence in Ballan Wrasse (Labrus bergylta)

The architecture underpinning genomic divergence is still a largely uncharted territory and likely case-dependent. Here, we investigated genome-wide variation in Ballan wrasse, a northeastern Atlantic fish species that displays two sympatric colour morphs, spotty and plain, that have been suggested to represent subspecies. We produced a chromosome-level reference genome and thereafter investigated genomic divergence among 152 individuals including both morphs, from two localities in Spain and Norway each and one in France. Differences between morphs dominated in Spain in accordance with sympatric divergence, whereas in Norway allopatric differentiation was prominent and repeated genomic signals of local divergence were found. Chromosomes had large low-recombining areas shared across all populations. Within the Spanish morphs, these areas contained large islands of divergence, totalling ~11% of the genome, and showed high morph specificity and strong selection. The same regions showed frequent admixture in the French morphs and no differentiation in Norway. In contrast, divergent regions observed between sampling localities in Norway were shorter and found throughout the genome. High inbreeding and lower diversity were observed in the Norwegian samples, consistent with the proposed recolonisation bottleneck and subsequent drift. Several genomic regions were significantly associated with morphs and contained tens of genes of diverse functions, suggesting that colouration is unlikely to be the sole driver of divergence. Our results do not support the hypothesis of shared larger genomic features underlying intraspecific colour divergence. Instead, we observe gradual accumulation of differences into low-recombining regions, likely when additional factors like assortative mating and/or lack of gene flow favour their development.

Genomic patterns of strain-specific genetic structure, linkage, and selection across fall armyworm populations

BACKGROUND

Molecular genetic approaches have become vital to understanding the evolutionary processes that act on insect pest populations. From mapping the development of resistance to monitoring and predicting pest movement, genomic tools can inform and enhance pest management programs. Here, we used whole genome sequencing population genomics to unravel novel patterns of population structure, linkage, and selection across the genome of a notorious agricultural pest, the fall armyworm.

RESULTS

Our data strongly support the existence of two genetically distinct strains of fall armyworm in North America, which have previously been referred to as the C-strain and the R-strain. Although these strains have diverged genetically, we find that differentiation is not uniform across the genome. The Z-chromosome appears to drive divergence between strains with high levels of linkage observed across this chromosome. We also show that a region of the Z-chromosome containing a circadian clock gene implicated in allochronic reproductive isolation is under strain-specific selection. Our data indicates that strains differ in their geographic distributions and exhibit distinct patterns of geographic sub-structuring indicative of unique dispersal patterns. We provide the first evidence for nuclear genomic differentiation between the two major overwintering populations of fall armyworm in the US. Finally, our data reveal population-specific patterns of selection on genomic regions containing putative insecticide resistance alleles, which could relate to their biogeography.

CONCLUSIONS

Our results support the existence of the fall armyworm as a pest dyad in the US, with genetically-distinct strains differing in their population structure, dispersal patterns, and genomic signatures of selection on regions likely involved reproductive isolation and insecticide resistance. These differences should be considered when devising and implementing management strategies.

Allochrony in Atlantic Lumpfish: Genomic and Otolith Shape Divergence Between Spring and Autumn Spawners

Allochrony is a form of reproductive isolation characterized by differences in the timing of spawning and may play a crucial role in the genetic and phenotypic divergence within species. The Atlantic lumpfish (Cyclopterus lumpus) is known to spawn in spring and autumn. However, the role of allochrony on the genomic structure of this species has not been addressed. Here, by combining whole genome sequencing data and otolith shape of 64 specimens, we explore the evolutionary drivers of divergence in Atlantic lumpfish, focusing on spring and autumn spawners sampled at two well-separated spawning grounds along the Norwegian coast. Overall, we identified pronounced genomic and morphologic differences between the two spawning groups. Genomic differences between the two groups were concentrated in three chromosomes, with a region of chromosome 1 encompassing the same single nucleotide polymorphisms (SNPs) driving differential season spawning for both localities, suggesting parallel responses. The functional analysis of the SNPs in this region revealed genes associated with responses to environmental stressors, possibly adaptations to seasonal variations at high latitudes. The morphological analysis of otoliths supported these findings, showing differences compatible with adaptations to seasonal light availability. The presence of genomic islands of divergence, alongside a general lack of differentiation across the mitochondrial genome, suggest recent and rapid selection processes potentially modulated by ongoing gene flow. This study underscores the importance of considering temporal genetic structures, particularly for species with bimodal spawning time, in conservation and management strategies to prevent overexploitation and optimize breeding programs.

Benchmarking the Mantel test and derived methods for testing association between distance matrices

Testing the association between objects is central in ecology, evolution, and quantitative sciences in general. Two types of variables can describe the relationships between objects: point variables (measured on individual objects), and distance variables (measured between pairs of objects). The Mantel test and derived methods have been extensively used for distance variables. Yet, these methods have been criticized due to low statistical power and inflated type I error when spatial autocorrelation is present. Here, we assessed the statistical power between different types of tested variables and the type I error rate over a wider range of autocorrelation intensities than previously assessed, both on univariate and multivariate data. We also illustrated the performance of distance matrix statistics through computational simulations of genetic diversity. We show that the Mantel test and derived methods are not affected by inflated type I error when spatial autocorrelation affects only one variable when investigating correlations, or when either the response or the explanatory variable(s) is affected by spatial autocorrelation while investigating causal relationships. As previously noted, with autocorrelation affecting more variables, inflated type I error could be reduced by modifying the significance threshold. Additionally, the Mantel test has no problem of statistical power when the hypothesis is formulated in terms of distance variables. We highlight that transformation of variable types should be avoided because of the potential information loss and modification of the tested hypothesis. We propose a set of guidelines to help choose the appropriate method according to the type of variables and defined hypothesis.

Circadian clock- and temperature-associated genes contribute to overall genomic differentiation along elevation in lichenized fungi

Circadian regulation is linked to local environmental adaptation, and many species with broad climatic niches display variation in circadian genes. Here, we hypothesize that lichenizing fungi occupying different climate zones tune their metabolism to local environmental conditions with the help of their circadian systems. We study two species of the genus Umbilicaria occupying similar climatic niches (Mediterranean and the cold temperate) in different continents. Using homology to Neurospora crassa genes, we identify gene sets associated with circadian rhythms (11 core, 39 peripheral genes) as well as temperature response (37 genes). Nucleotide diversity of these genes is significantly correlated with mean annual temperature, minimum temperature of the coldest month and mean temperature of the coldest quarter. Furthermore, we identify altitudinal clines in allele frequencies in several non-synonymous substitutions in core clock components, for example, white collar-like, frh-like and various ccg-like genes. A dN/dS approach revealed a few significant peripheral clock- and temperature-associated genes (e.g. ras-1-like, gna-1-like) that may play a role in fine-tuning the circadian clock and temperature-response machinery. An analysis of allele frequency changes demonstrated the strongest evidence for differentiation above the genomic background in the clock-associated genes in U. pustulata. These results highlight the likely relevance of the circadian clock in environmental adaptation, particularly frost tolerance, of lichens. Whether or not the fungal clock modulates the symbiotic interaction within the lichen consortium remains to be investigated. We corroborate the finding of genetic variation in clock components along altitude-not only latitude-as has been reported in other species.

The effect of divergent and parallel selection on the genomic landscape of divergence

While the role of selection in divergence along the speciation continuum is theoretically well understood, defining specific signatures of selection in the genomic landscape of divergence is empirically challenging. Modelling approaches can provide insight into the potential role of selection on the emergence of a heterogenous genomic landscape of divergence. Here, we extend and apply an individual-based approach that simulates the phenotypic and genotypic distributions of two populations under a variety of selection regimes, genotype-phenotype maps, modes of migration, and genotype-environment interactions. We show that genomic islands of high differentiation and genomic valleys of similarity may respectively form under divergent and parallel selection between populations. For both types of between-population selection, negative and positive frequency-dependent selection within populations generated genomic islands of higher magnitude and genomic valleys of similarity, respectively. Divergence rates decreased under strong dominance with divergent selection, as well as in models including genotype-environment interactions under parallel selection. For both divergent and parallel selection models, divergence rate was higher under an intermittent migration regime between populations, in contrast to a constant level of migration across generations, despite an equal number of total migrants. We highlight that interpreting a particular evolutionary history from an observed genomic pattern must be done cautiously, as similar patterns may be obtained from different combinations of evolutionary processes. Modelling approaches such as ours provide an opportunity to narrow the potential routes that generate the genomic patterns of specific evolutionary histories.

Investigating the role of genetic variation in vgll3 and six6 in the domestication of gilthead seabream (Sparus aurata Linnaeus) and European seabass (Dicentrarchus labrax Linnaeus)

Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age-at-maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth-related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes ('gene-level genotyping'). The second dataset examined 300-600 k SNPs located in the chromosomes of the two genes, derived from a meta-analysis of a Pool-Seq experiment involving farmed and wild populations distributed widely across the Mediterranean ('chromosome-level genotyping'). The gene-level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome-level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function.

Drivers of genomic landscapes of differentiation across a Populus divergence gradient

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

Maintenance of local adaptation despite gene flow in a coastal songbird

Adaptation to local environments is common in widespread species and the basis of ecological speciation. The song sparrow (Melospiza melodia) is a widespread, polytypic passerine that occurs in shrubland habitats throughout North America. We examined the population structure of two parapatric subspecies that inhabit different environments: the Atlantic song sparrow (M. m. atlantica), a coastal specialist, and the eastern song sparrow (M. m. melodia), a shrubland generalist. These populations lacked clear mitochondrial population structure, yet coastal birds formed a distinct nuclear genetic cluster. We found weak overall genomic differentiation between these subspecies, suggesting either recent divergence, extensive gene flow, or a combination thereof. There was a steep genetic cline at the transition to coastal habitats, consistent with isolation by environment, not isolation by distance. A phenotype under divergent selection, bill size, varied with the amount of coastal ancestry in transitional areas, but larger bill size was maintained in coastal habitats regardless of ancestry, further supporting a role for selection in the maintenance of these subspecies. Demographic modeling suggested a divergence history of limited gene flow followed by secondary contact, which has emerged as a common theme in adaptive divergence across taxa.

A Novel Mutator-Like Transposable Elements With Unusual Structure and Recent Transpositions in Barley (Hordeum vulgare)

Mutator-like transposable elements (MULEs) represent a unique superfamily of DNA transposons as they can capture host genes and cause higher frequency of mutations in some eukaryotes. Despite their essential roles in plant evolution and functional genomics, MULEs are not fully understood yet in many important crops including barley (Hordeum vulgare). In this study, we analyzed the barley genome and identified a new mutator transposon Hvu_Abermu. This transposon is present at extremely high copy number in barley and shows unusual structure as it contains three open reading frames (ORFs) including one ORF (ORF1) encoding mutator transposase protein and one ORF (ORFR) showing opposite transcriptional orientation. We identified homologous sequences of Hvu_Abermu in both monocots and dicots and grouped them into a large mutator family named Abermu. Abermu transposons from different species share significant sequence identity, but they exhibit distinct sequence structures. Unlike the transposase proteins which are highly conserved between Abermu transposons from different organisms, the ORFR-encoded proteins are quite different from distant species. Phylogenetic analysis indicated that Abermu transposons shared closer evolutionary relationships with the maize MuDR transposon than other reported MULEs. We also found phylogenetic incongruence for the Abermu transposons identified in rice and its wild species implying the possibility of horizontal transfer of transposon. Further comparison indicated that over 200 barley genes contain Abermu-related sequences. We analyzed the barley pan genomes and detected polymorphic Hvu_Abermu transposons between the sequenced 23 wild and cultivated barley genomes. Our efforts identified a novel mutator transposon and revealed its recent transposition activity, which may help to develop genetic tools for barley and other crops.

Speciation dynamics and extent of parallel evolution along a lake-stream environmental contrast in African cichlid fishes

Understanding the dynamics of speciation is a central topic in evolutionary biology. Here, we investigated how morphological and genomic differentiation accumulated along the speciation continuum in the African cichlid fish Astatotilapia burtoni. While morphological differentiation was continuously distributed across different lake-stream population pairs, we found that there were two categories with respect to genomic differentiation, suggesting a “gray zone” of speciation at ~0.1% net nucleotide divergence. Genomic differentiation was increased in the presence of divergent selection and drift compared to drift alone. The quantification of phenotypic and genetic parallelism in four cichlid species occurring along a lake-stream environmental contrast revealed parallel and antiparallel components in rapid adaptive divergence, and morphological convergence in species replicates inhabiting the same environments. Furthermore, we show that the extent of parallelism was higher when ancestral populations were more similar. Our study highlights the complementary roles of divergent selection and drift on speciation and parallel evolution.

The genic view of hybridization in the Anthropocene

Human impact is noticeable around the globe, indicating that a new era might have begun: the Anthropocene. Continuing human activities, including land-use changes, introduction of non-native species and rapid climate change, are altering the distributions of countless species, often giving rise to human-mediated hybridization events. While the interbreeding of different populations or species can have detrimental effects, such as genetic extinction, it can be beneficial in terms of adaptive introgression or an increase in genetic diversity. In this paper, I first review the different mechanisms and outcomes of anthropogenic hybridization based on literature from the last five years (2016-2020). The most common mechanisms leading to the interbreeding of previously isolated taxa include habitat change (51% of the studies) and introduction of non-native species (34% intentional and 19% unintentional). These human-induced hybridization events most often result in introgression (80%). The high incidence of genetic exchange between the hybridizing taxa indicates that the application of a genic view of speciation (and introgression) can provide crucial insights on how to address hybridization events in the Anthropocene. This perspective considers the genome as a dynamic collection of genetic loci with distinct evolutionary histories, giving rise to a heterogenous genomic landscape in terms of genetic differentiation and introgression. First, understanding this genomic landscape can lead to a better selection of diagnostic genetic markers to characterize hybrid populations. Second, describing how introgression patterns vary across the genome can help to predict the likelihood of negative processes, such as demographic and genetic swamping, as well as positive outcomes, such as adaptive introgression. It is especially important to not only quantify how much genetic material introgressed, but also what has been exchanged. Third, comparing introgression patterns in pre-Anthropocene hybridization events with current human-induced cases might provide novel insights into the likelihood of genetic swamping or species collapse during an anthropogenic hybridization event. However, this comparative approach remains to be tested before it can be applied in practice. Finally, the genic view of introgression can be combined with conservation genomic studies to determine the legal status of hybrids and take appropriate measures to manage anthropogenic hybridization events. The interplay between evolutionary and conservation genomics will result in the constant exchange of ideas between these fields which will not only improve our knowledge on the origin of species, but also how to conserve and protect them.

A phylogeny of white-eyes based on ultraconserved elements

White-eyes are an iconic radiation of passerine birds that have been the subject of studies in evolutionary biology, biogeography, and speciation theory. Zosterops white-eyes in particular are thought to have radiated rapidly across continental and insular regions of the Afro- and Indo-Pacific tropics, yet their phylogenetic history remains equivocal. Here, we sampled 77% of the genera and 47% of known white-eye species and sequenced thousands of ultraconserved elements to infer the phylogeny of the avian family Zosteropidae. We used concatenated maximum likelihood and species tree methods and found strong support for seven clades of white-eyes and three clades within the species-rich Zosterops radiation.

An island-hopping bird reveals how founder events shape genome-wide divergence

When populations colonize new areas, both strong selection and strong drift can be experienced due to novel environments and small founding populations, respectively. Empirical studies have predominantly focused on the phenotype when assessing the role of selection, and limited neutral-loci when assessing founder-induced loss of diversity. Consequently, the extent to which processes interact to influence evolutionary trajectories is difficult to assess. Genomic-level approaches provide the opportunity to simultaneously consider these processes. Here, we examine the roles of selection and drift in shaping genomic diversity and divergence in historically documented sequential island colonizations by the silvereye (Zosterops lateralis). We provide the first empirical demonstration of the rapid appearance of highly diverged genomic regions following population founding, the position of which are highly idiosyncratic. As these regions rarely contained loci putatively under selection, it is most likely that these differences arise via the stochastic nature of the founding process. However, selection is required to explain rapid evolution of larger body size in insular silvereyes. Reconciling our genomic data with these phenotypic patterns suggests there may be many genomic routes to the island phenotype, which vary across populations. Finally, we show that accelerated divergence associated with multiple founding steps is the product of genome-wide rather than localized differences, and that diversity erodes due to loss of rare alleles. However, even multiple founder events do not result in divergence and diversity levels seen in evolutionary older subspecies, and therefore do not provide a shortcut to speciation as proposed by founder-effect speciation models.