Intercontinental karyotype-environment parallelism supports a role for a chromosomal inversion in local adaptation in a seaweed fly

Genomic Islands of Divergence Between Drosophila yakuba Subspecies are Predominantly Driven by Chromosomal Inversions and the Recombination Landscape

During the early stages of local adaptation and speciation, genetic differences tend to accumulate at certain regions of the genome leading to the formation of genomic islands of divergence (GIDs). This pattern may be due to selection and/or difference in the rate of recombination. Here, we investigate the possible causes of GIDs in Drosophila yakuba mayottensis, and reconfirm using field collection its association with toxic noni (Morinda citrifolia) fruits on the Mayotte island. Population genomics revealed lack of genetic structure on the island and identified 23 GIDs distinguishing D. y. mayottensis from generalist mainland populations of D. y. yakuba. The GIDs were enriched with gene families involved in the metabolism of lipids, sugars, peptides and xenobiotics, suggesting a role in host shift. We assembled a new genome for D. y. mayottensis and identified five novel chromosomal inversions. Twenty one GIDs (~99% of outlier windows) fell in low recombining regions or subspecies-specific inversions. However, only two GIDs were in collinear, normally recombining regions suggesting a signal of hard selective sweeps. Unlike D. y. mayottensis, D. sechellia, the only other noni-specialist, is known to be homosequential with its generalist relatives. Thus, whereas structural variation may disproportionally shape GIDs in some species, striking parallel adaptations can occur between species despite distinct genomic architectures.

Genetics in the Ocean's Twilight Zone: Population Structure of the Glacier Lanternfish Across Its Distribution Range

The mesopelagic zone represents one of the few habitats that remains relatively untouched from anthropogenic activities. Among the many species inhabiting the north Atlantic mesopelagic zone, glacier lanternfish (Benthosema glaciale) is the most abundant and widely distributed. This species has been regarded as a potential target for a dedicated fishery despite the scarce knowledge of its population genetic structure. Here, we investigated its genetic structure across the North Atlantic and into the Mediterranean Sea using 121 SNPs, which revealed strong differentiation among three main groups: the Mediterranean Sea, oceanic samples, and Norwegian fjords. The Mediterranean samples displayed less than half the genetic variation of the remaining ones. Very weak or nearly absent genetic structure was detected among geographically distinct oceanic samples across the North Atlantic, which contrasts with the low motility of the species. In contrast, a longitudinal gradient of differentiation was observed in the Mediterranean Sea, where genetic connectivity is known to be strongly shaped by oceanographic processes such as current patterns and oceanographic discontinuities. In addition, 12 of the SNPs, in linkage disequilibrium, drove a three clusters' pattern detectable through Principal Component Analysis biplot matching the genetic signatures generally associated with large chromosomal rearrangements, such as inversions. The arrangement of this putative inversion showed frequency differences between open-ocean and more confined water bodies such as the fjords and the Mediterranean, as it was fixed in the latter for the second most common arrangement of the fjord's samples. However, whether genetic differentiation was driven by local adaptation, secondary contact, or a combination of both factors remains undetermined. The major finding of this study is that B. glaciale in the North Atlantic-Mediterranean is divided into three major genetic units, information that should be combined with demographic properties to outline the management of this species prior to any eventual fishery attempt.

Genomics of natural populations: gene conversion events reveal selected genes within the inversions of Drosophila pseudoobscura

When adaptive phenotypic variation or quantitative trait loci map within an inverted segment of a chromosome, researchers often despair because the suppression of crossing over will prevent the discovery of selective target genes that established the rearrangement. If an inversion polymorphism is old enough, then the accumulation of gene conversion tracts offers the promise that quantitative trait loci or selected loci within inversions can be mapped. The inversion polymorphism of Drosophila pseudoobscura is a model system to show that gene conversion analysis is a useful tool for mapping selected loci within inversions. D. pseudoobscura has over 30 different chromosomal arrangements on the third chromosome (Muller C) in natural populations and their frequencies vary with changes in environmental habitats. Statistical tests of five D. pseudoobscura gene arrangements identified outlier genes within inverted regions that had potentially heritable variation, either fixed amino acid differences or differential expression patterns. We use genome sequences of the inverted third chromosome (Muller C) to infer 98,443 gene conversion tracts for a total coverage of 142 Mb or 7.2× coverage of the 19.7 Mb chromosome. We estimated gene conversion tract coverage in the 2,668 genes on Muller C and tested whether gene conversion coverage was similar among arrangements for outlier vs non-outlier loci. Outlier genes had lower gene conversion tract coverage among arrangements than the non-outlier genes suggesting that selection removes exchanged DNA in the outlier genes. These data support the hypothesis that the third chromosome in D. pseudoobscura captured locally adapted combinations of alleles prior to inversion mutation events.

Microbiome Dynamics and Functional Composition in Coelopa frigida (Diptera, Coelopidae): Insights into Trophic Specialization of Kelp Flies

Coelopidae (Diptera), known as kelp flies, exhibit an ecological association with beached kelp and other rotting seaweeds. This unique trophic specialization necessitates significant adaptations to overcome the limitations of an algal diet. We aimed to investigate whether the flies' microbiome could be one of these adaptive mechanisms. Our analysis focused on assessing composition and diversity of adult and larval microbiota of the kelp fly Coelopa frigida. Feeding habits of the larvae of this species have been subject of numerous studies, with debates whether they directly consume kelp or primarily feed on associated bacteria. By using a 16S rRNA metabarcoding approach, we found that the larval microbiota displayed considerably less diversity than adults, heavily dominated by only four operational taxonomic units (OTUs). Phylogenetic placement recovered the most dominant OTU of the larval microbiome, which is the source of more than half of all metabarcoding sequence reads, as an undescribed genus of Orbaceae (Gammaproteobacteria). Interestingly, this OTU is barely found among the 15 most abundant taxa of the adult microbiome, where it is responsible for less than 2% of the metabarcoding sequence reads. The other three OTUs dominating the larval microbiome have been assigned as Psychrobacter (Gammaproteobacteria), Wohlfahrtiimonas (Gammaproteobacteria), and Cetobacterium (Fusobacteriota). Moreover, we also uncovered a distinct shift in the functional composition between the larval and adult stages, where our taxonomic profiling suggests a significant decrease in functional diversity in larval samples. Our study offers insights into the microbiome dynamics and functional composition of Coelopa frigida.

Understanding the chromosomal evolution in cuckoos (Aves, Cuculiformes): a journey through unusual rearrangements

The Cuculiformes are a family of over 150 species that live in a range of habitats, such as forests, savannas, and deserts. Here, bacterial artificial chromosome (BAC) probes (75 from chicken and 14 from zebra finch macrochromosomes 1-10 +ZW and for microchromosomes 11-28 (except 16)) were used to investigate chromosome homologies between chicken and the squirrel cuckoo (Piaya cayana). In addition, repetitive DNA probes were applied to characterize the chromosome organization and to explore the role of these sequences in the karyotype evolution of P. cayana. We also applied BAC probes for chicken chromosome 17 and Z to the guira cuckoo (Guira guira) to test whether this species has an unusual Robertsonian translocation between a microchromosome and the Z chromosome, recently described in the smooth-billed ani (Crotophaga ani). Our results revealed extensive chromosome reorganization with inter- and intrachromosomal rearrangements in P. cayana, including a conspicuous chromosome size and heterochromatin polymorphism on chromosome pair 20. Furthermore, we confirmed that the Z-autosome Robertsonian translocation found in C. ani is also found in G. guira, not P. cayana. These findings suggest that this translocation occurred prior to the divergence between C. ani and G. guira, but after the divergence with P. cayana.

Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii)

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the closure of all U.S. black abalone fisheries since 1993. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS remains unknown. To address these uncertainties, we sequenced and analysed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Outside the inversion, genetic differentiation between sites is minimal and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Demographic inference does indicate a severe population bottleneck beginning just 15 generations in the past, but this decline is short lived, with present-day size far exceeding the pre-bottleneck status quo. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of population genetic structure, uniform diversity and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

Investigating structural variant, indel and single nucleotide polymorphism differentiation between locally adapted Atlantic salmon populations

Genomic structural variants (SVs) are now recognized as an integral component of intraspecific polymorphism and are known to contribute to evolutionary processes in various organisms. However, they are inherently difficult to detect and genotype from readily available short-read sequencing data, and therefore remain poorly documented in wild populations. Salmonid species displaying strong interpopulation variability in both life history traits and habitat characteristics, such as Atlantic salmon (Salmo salar), offer a prime context for studying adaptive polymorphism, but the contribution of SVs to fine-scale local adaptation has yet to be explored. Here, we performed a comparative analysis of SVs, single nucleotide polymorphisms (SNPs) and small indels (<50 bp) segregating in the Romaine and Puyjalon salmon, two putatively locally adapted populations inhabiting neighboring rivers (Québec, Canada) and showing pronounced variation in life history traits, namely growth, fecundity, and age at maturity and smoltification. We first catalogued polymorphism using a hybrid SV characterization approach pairing both short- (16X) and long-read sequencing (20X) for variant discovery with graph-based genotyping of SVs across 60 salmon genomes, along with characterization of SNPs and small indels from short reads. We thus identified 115,907 SVs, 8,777,832 SNPs and 1,089,321 short indels, with SVs covering 4.8 times more base pairs than SNPs. All three variant types revealed a highly congruent population structure and similar patterns of F ST and density variation along the genome. Finally, we performed outlier detection and redundancy analysis (RDA) to identify variants of interest in the putative local adaptation of Romaine and Puyjalon salmon. Genes located near these variants were enriched for biological processes related to nervous system function, suggesting that observed variation in traits such as age at smoltification could arise from differences in neural development. This study therefore demonstrates the feasibility of large-scale SV characterization and highlights its relevance for salmonid population genomics.

Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii)

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the species' designation as critically endangered. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS is unknown. To address these uncertainties, we sequenced and analyzed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Genetic divergence between sites is minimal, and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Despite this, however, demographic inference confirms a severe population bottleneck beginning around the time of WS onset, highlighting the temporal offset that may occur between a population collapse and its potential impact on genetic diversity. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of genetic structure, uniform diversity, and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

A supergene in seaweed flies modulates male traits and female perception

Supergenes, tightly linked sets of alleles, offer some of the most spectacular examples of polymorphism persisting under long-term balancing selection. However, we still do not understand their evolution and persistence, especially in the face of accumulation of deleterious elements. Here, we show that an overdominant supergene in seaweed flies, Coelopa frigida, modulates male traits, potentially facilitating disassortative mating and promoting intraspecific polymorphism. Across two continents, the Cf-Inv(1) supergene strongly affected the composition of male cuticular hydrocarbons (CHCs) but only weakly affected CHC composition in females. Using gas chromatography-electroantennographic detection, we show that females can sense male CHCs and that there may be differential perception between genotypes. Combining our phenotypic results with RNA-seq data, we show that candidate genes for CHC biosynthesis primarily show differential expression for Cf-Inv(1) in males but not females. Conversely, candidate genes for odorant detection were differentially expressed in both sexes but showed high levels of divergence between supergene haplotypes. We suggest that the reduced recombination between supergene haplotypes may have led to rapid divergence in mate preferences as well as increasing linkage between male traits, and overdominant loci. Together this probably helped to maintain the polymorphism despite deleterious effects in homozygotes.

Ecology, the pace-of-life, epistatic selection and the maintenance of genetic variation in life-history genes

Evolutionary genetics has long struggled with understanding how functional genes under selection remain polymorphic in natural populations. Taking as a starting point that natural selection is ultimately a manifestation of ecological processes, we spotlight an underemphasized and potentially ubiquitous ecological effect that may have fundamental effects on the maintenance of genetic variation. Negative frequency dependency is a well-established emergent property of density dependence in ecology, because the relative profitability of different modes of exploiting or utilizing limiting resources tends to be inversely proportional to their frequency in a population. We suggest that this may often generate negative frequency-dependent selection (NFDS) on major effect loci that affect rate-dependent physiological processes, such as metabolic rate, that are phenotypically manifested as polymorphism in pace-of-life syndromes. When such a locus under NFDS shows stable intermediate frequency polymorphism, this should generate epistatic selection potentially involving large numbers of loci with more minor effects on life-history (LH) traits. When alternative alleles at such loci show sign epistasis with a major effect locus, this associative NFDS will promote the maintenance of polygenic variation in LH genes. We provide examples of the kind of major effect loci that could be involved and suggest empirical avenues that may better inform us on the importance and reach of this process.

Genomic responses to parallel temperature gradients in the eelgrass Zostera marina in adjacent bays

The extent of parallel genomic responses to similar selective pressures depends on a complex array of environmental, demographic, and evolutionary forces. Laboratory experiments with replicated selective pressures yield mixed outcomes under controlled conditions and our understanding of genomic parallelism in the wild is limited to a few well-established systems. Here, we examine genomic signals of selection in the eelgrass Zostera marina across temperature gradients in adjacent embayments. Although we find many genomic regions with signals of selection within each bay there is very little overlap in signals of selection at the SNP level, despite most polymorphisms being shared across bays. We do find overlap at the gene level, potentially suggesting multiple mutational pathways to the same phenotype. Using polygenic models we find that some sets of candidate SNPs are able to predict temperature across both bays, suggesting that small but parallel shifts in allele frequencies may be missed by independent genome scans. Together, these results highlight the continuous rather than binary nature of parallel evolution in polygenic traits and the complexity of evolutionary predictability.

Rapid adaptation in a fast-changing world: Emerging insights from insect genomics

Many researchers have questioned the ability of biota to adapt to rapid anthropogenic environmental shifts. Here, we synthesize emerging genomic evidence for rapid insect evolution in response to human pressure. These new data reveal diverse genomic mechanisms (single locus, polygenic, structural shifts; introgression) underpinning rapid adaptive responses to a variety of anthropogenic selective pressures. While the effects of some human impacts (e.g. pollution; pesticides) have been previously documented, here we highlight startling new evidence for rapid evolutionary responses to additional anthropogenic processes such as deforestation. These recent findings indicate that diverse insect assemblages can indeed respond dynamically to major anthropogenic evolutionary challenges. Our synthesis also emphasizes the critical roles of genomic architecture, standing variation and gene flow in maintaining future adaptive potential. Broadly, it is clear that genomic approaches are essential for predicting, monitoring and responding to ongoing anthropogenic biodiversity shifts in a fast-changing world.

Inversions and parallel evolution

Local adaptation leads to differences between populations within a species. In many systems, similar environmental contrasts occur repeatedly, sometimes driving parallel phenotypic evolution. Understanding the genomic basis of local adaptation and parallel evolution is a major goal of evolutionary genomics. It is now known that by preventing the break-up of favourable combinations of alleles across multiple loci, genetic architectures that reduce recombination, like chromosomal inversions, can make an important contribution to local adaptation. However, little is known about whether inversions also contribute disproportionately to parallel evolution. Our aim here is to highlight this knowledge gap, to showcase existing studies, and to illustrate the differences between genomic architectures with and without inversions using simple models. We predict that by generating stronger effective selection, inversions can sometimes speed up the parallel adaptive process or enable parallel adaptation where it would be impossible otherwise, but this is highly dependent on the spatial setting. We highlight that further empirical work is needed, in particular to cover a broader taxonomic range and to understand the relative importance of inversions compared to genomic regions without inversions.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

The Relative Power of Structural Genomic Variation versus SNPs in Explaining the Quantitative Trait Growth in the Marine Teleost Chrysophrys auratus

Background: Genetic diversity provides the basic substrate for evolution. Genetic variation consists of changes ranging from single base pairs (single-nucleotide polymorphisms, or SNPs) to larger-scale structural variants, such as inversions, deletions, and duplications. SNPs have long been used as the general currency for investigations into how genetic diversity fuels evolution. However, structural variants can affect more base pairs in the genome than SNPs and can be responsible for adaptive phenotypes due to their impact on linkage and recombination. In this study, we investigate the first steps needed to explore the genetic basis of an economically important growth trait in the marine teleost finfish Chrysophrys auratus using both SNP and structural variant data. Specifically, we use feature selection methods in machine learning to explore the relative predictive power of both types of genetic variants in explaining growth and discuss the feature selection results of the evaluated methods. Methods: SNP and structural variant callers were used to generate catalogues of variant data from 32 individual fish at ages 1 and 3 years. Three feature selection algorithms (ReliefF, Chi-square, and a mutual-information-based method) were used to reduce the dataset by selecting the most informative features. Following this selection process, the subset of variants was used as features to classify fish into small, medium, or large size categories using KNN, naïve Bayes, random forest, and logistic regression. The top-scoring features in each feature selection method were subsequently mapped to annotated genomic regions in the zebrafish genome, and a permutation test was conducted to see if the number of mapped regions was greater than when random sampling was applied. Results: Without feature selection, the prediction accuracies ranged from 0 to 0.5 for both structural variants and SNPs. Following feature selection, the prediction accuracy increased only slightly to between 0 and 0.65 for structural variants and between 0 and 0.75 for SNPs. The highest prediction accuracy for the logistic regression was achieved for age 3 fish using SNPs, although generally predictions for age 1 and 3 fish were very similar (ranging from 0–0.65 for both SNPs and structural variants). The Chi-square feature selection of SNP data was the only method that had a significantly higher number of matches to annotated genomic regions of zebrafish than would be explained by chance alone. Conclusions: Predicting a complex polygenic trait such as growth using data collected from a low number of individuals remains challenging. While we demonstrate that both SNPs and structural variants provide important information to help understand the genetic basis of phenotypic traits such as fish growth, the full complexities that exist within a genome cannot be easily captured by classical machine learning techniques. When using high-dimensional data, feature selection shows some increase in the prediction accuracy of classification models and provides the potential to identify unknown genomic correlates with growth. Our results show that both SNPs and structural variants significantly impact growth, and we therefore recommend that researchers interested in the genotype–phenotype map should strive to go beyond SNPs and incorporate structural variants in their studies as well. We discuss how our machine learning models can be further expanded to serve as a test bed to inform evolutionary studies and the applied management of species.

Environmentally associated chromosomal structural variation influences fine-scale population structure of Atlantic Salmon (Salmo salar)

Chromosomal rearrangements (e.g., inversions, fusions, and translocations) have long been associated with environmental variation in wild populations. New genomic tools provide the opportunity to examine the role of these structural variants in shaping adaptive differences within and among wild populations of non-model organisms. In Atlantic Salmon (Salmo salar), variations in chromosomal rearrangements exist across the species natural range, yet the role and importance of these structural variants in maintaining adaptive differences among wild populations remains poorly understood. We genotyped Atlantic Salmon (n = 1429) from 26 populations within a highly genetically structured region of southern Newfoundland, Canada with a 220K SNP array. Multivariate analysis, across two independent years, consistently identified variation in a structural variant (translocation between chromosomes Ssa01 and Ssa23), previously associated with evidence of trans-Atlantic secondary contact, as the dominant factor influencing population structure in the region. Redundancy analysis suggested that variation in the Ssa01/Ssa23 chromosomal translocation is strongly correlated with temperature. Our analyses suggest environmentally mediated selection acting on standing genetic variation in genomic architecture introduced through secondary contact may underpin fine-scale local adaptation in Placentia Bay, Newfoundland, Canada, a large and deep embayment, highlighting the importance of chromosomal structural variation as a driver of contemporary adaptive divergence.

A large chromosomal inversion shapes gene expression in seaweed flies (Coelopa frigida)

Inversions often underlie complex adaptive traits, but the genic targets inside them are largely unknown. Gene expression profiling provides a powerful way to link inversions with their phenotypic consequences. We examined the effects of the Cf-Inv(1) inversion in the seaweed fly Coelopa frigida on gene expression variation across sexes and life stages. Our analyses revealed that Cf-Inv(1) shapes global expression patterns, most likely via linked variation, but the extent of this effect is variable, with much stronger effects in adults than larvae. Furthermore, within adults, both common as well as sex-specific patterns were found. The vast majority of these differentially expressed genes mapped to Cf-Inv(1). However, genes that were differentially expressed in a single context (i.e., in males, females, or larvae) were more likely to be located outside of Cf-Inv(1). By combining our findings with genomic scans for environmentally associated SNPs, we were able to pinpoint candidate variants in the inversion that may underlie mechanistic pathways that determine phenotypes. Together the results of this study, combined with previous findings, support the notion that the polymorphic Cf-Inv(1) inversion in this species is a major factor shaping both coding and regulatory variation resulting in highly complex adaptive effects.

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis

Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.

Pre-introduction introgression contributes to parallel differentiation and contrasting hybridization outcomes between invasive and native marine mussels

Non-native species experience novel selection pressures in introduced environments and may interbreed with native lineages. Species introductions therefore provide opportunities to investigate repeated patterns of adaptation and introgression across replicated contact zones. Here, we investigate genetic parallelism between multiple introduced populations of the invasive marine mussel, Mytilus galloprovincialis, in the absence (South Africa and California) and presence of hybridization with a native congener (Mytilus planulatus in Batemans Bay and Sydney Harbour, Australia). Repeatability in post-introduction differentiation from native-range populations varied between genetically distinct Atlantic and Mediterranean lineages, with Atlantic-derived introductions displaying high differentiation (maxF_ST  > 0.4) and parallelism at outlier loci. Identification of long noncoding RNA transcripts (lncRNA) additionally allowed us to clarify that parallel responses are largely limited to protein-coding loci, with lncRNAs likely evolving under evolutionary constraints. Comparisons of independent hybrid zones revealed differential introgression most strongly in Batemans Bay, with an excess of M. galloprovincialis ancestry and resistance to introgression at loci differentiating parental lineages (M. planulatus and Atlantic M. galloprovincialis). Additionally, contigs putatively introgressed with divergent alleles from a closely related species, Mytilus edulis, showed stronger introgression asymmetries compared with genome-wide trends and also diverged in parallel in both Atlantic-derived introductions. These results suggest that divergent demographic histories experienced by introduced lineages, including pre-introduction introgression, influence contemporary admixture dynamics. Our findings build on previous investigations reporting contributions of historical introgression to intrinsic reproductive architectures shared between marine lineages and illustrate that interspecific introgression history can shape differentiation between colonizing populations and their hybridization with native congeners.

Inside the supergene of the bird with four sexes

The white-throated sparrow (Zonotrichia albicollis) offers unique opportunities to understand the adaptive value of supergenes, particularly their role in alternative phenotypes. In this species, alternative plumage morphs segregate with a nonrecombining segment of chromosome 2, which has been called a 'supergene'. The species mates disassortatively with respect to the supergene; that is, each breeding pair consists of one individual with it and one without it. This species has therefore been called the "bird with four sexes". The supergene segregates with a behavioral phenotype; birds with it are more aggressive and less parental than birds without it. Here, we review our efforts to identify the genes inside the supergene that are responsible for the behavioral polymorphism. The gene ESR1, which encodes estrogen receptor α, differs between the morphs and predicts both territorial and parental behavior. Variation in the regulatory regions of ESR1 causes an imbalance in expression of the two alleles, and the degree to which this imbalance favors the supergene allele predicts territorial singing. In heterozygotes, knockdown of ESR1 causes a phenotypic switch, from more aggressive to less aggressive. We recently showed that another gene important for social behavior, vasoactive intestinal peptide (VIP), is differentially expressed between the morphs and predicts territorial singing. We hypothesize that ESR1 and VIP contribute to behavior in a coordinated way and could represent co-adapted alleles. Because the supergene contains more than 1000 individual genes, this species provides rich possibilities for discovering alleles that work together to mediate life-history trade-offs and maximize the fitness of alternative complex phenotypes.

Distribution modelling of an introduced species: do adaptive genetic markers affect potential range?

Biological invasions have increased in the last few decades mostly due to anthropogenic causes such as globalization of trade. Because invaders sometimes cause large economic losses and ecological disturbances, estimating their origin and potential geographical ranges is useful. Drosophila subobscura is native to the Old World but was introduced in the New World in the late 1970s and spread widely. We incorporate information on adaptive genetic markers into ecological niche modelling and then estimate the most probable geographical source of colonizers; evaluate whether the genetic bottleneck experienced by founders affects their potential distribution; and finally test whether this species has spread to all its potential suitable habitats worldwide. We find the environmental space occupied by this species in its native and introduced distributions are notably the same, although the introduced niche has shifted slightly towards higher temperature and lower precipitation. The genetic bottleneck of founding individuals was a key factor limiting the spread of this introduced species. We also find that regions in the Mediterranean and north-central Portugal show the highest probability of being the origin of the colonizers. Using genetically informed environmental niche modelling can enhance our understanding of the initial colonization and spread of invasive species, and also elucidate potential areas of future expansions worldwide.

Consequences of Single-Locus and Tightly Linked Genomic Architectures for Evolutionary Responses to Environmental Change

Genetic and genomic architectures of traits under selection are key factors influencing evolutionary responses. Yet, knowledge of their impacts has been limited by a widespread assumption that most traits are controlled by unlinked polygenic architectures. Recent advances in genome sequencing and eco-evolutionary modeling are unlocking the potential for integrating genomic information into predictions of population responses to environmental change. Using eco-evolutionary simulations, we demonstrate that hypothetical single-locus control of a life history trait produces highly variable and unpredictable harvesting-induced evolution relative to the classically applied multilocus model. Single-locus control of complex traits is thought to be uncommon, yet blocks of linked genes, such as those associated with some types of structural genomic variation, have emerged as taxonomically widespread phenomena. Inheritance of linked architectures resembles that of single loci, thus enabling single-locus-like modeling of polygenic adaptation. Yet, the number of loci, their effect sizes, and the degree of linkage among them all occur along a continuum. We review how linked architectures are often associated, directly or indirectly, with traits expected to be under selection from anthropogenic stressors and are likely to play a large role in adaptation to environmental disturbance. We suggest using single-locus models to explore evolutionary extremes and uncertainties when the trait architecture is unknown, refining parameters as genomic information becomes available, and explicitly incorporating linkage among loci when possible. By overestimating the complexity (e.g., number of independent loci) of the genomic architecture of traits under selection, we risk underestimating the complexity (e.g., nonlinearity) of their evolutionary dynamics.

Balancing selection via life-history trade-offs maintains an inversion polymorphism in a seaweed fly

How natural diversity is maintained is an evolutionary puzzle. Genetic variation can be eroded by drift and directional selection but some polymorphisms persist for long time periods, implicating a role for balancing selection. Here, we investigate the maintenance of a chromosomal inversion polymorphism in the seaweed fly Coelopa frigida. Using experimental evolution and quantifying fitness, we show that the inversion underlies a life-history trade-off, whereby each haplotype has opposing effects on larval survival and adult reproduction. Numerical simulations confirm that such antagonistic pleiotropy can maintain polymorphism. Our results also highlight the importance of sex-specific effects, dominance and environmental heterogeneity, whose interaction enhances the maintenance of polymorphism through antagonistic pleiotropy. Overall, our findings directly demonstrate how overdominance and sexual antagonism can emerge from a life-history trade-off, inviting reconsideration of antagonistic pleiotropy as a key part of multi-headed balancing selection processes that enable the persistence of genetic variation.

Few studies empirically pinpoint how balanced polymorphisms are maintained. “Mérot et al”. identify an inversion polymorphism that is maintained in seaweed fly populations because of antagonistic pleiotropy that mediates a classic life history tradeoff between larval survival and adult reproduction.

Weak genetic structure despite strong genomic signal in lesser sandeel in the North Sea

Sandeels are an ecologically important group of fishes; they are a key part of the food chain serving as food for marine mammals, seabirds and fish. Sandeels are further targeted by a large industrial fishery, which has led to concern about ecosystem effects. In the North Sea, the lesser sandeel Ammodytes marinus is by far the most prevalent species of sandeel in the fishery. Management of sandeel in the North Sea plus the Kattegat is currently divided into seven geographical areas, based on subtle differences in demography, population dynamics and results from simulations of larval dispersal. However, little is known about the underlying genetic population structure. In this study, we used 2,522 SNPs derived from restriction site-associated DNA sequencing (RADseq) typed in 429 fish representing four main sandeel management areas. Our main results showed (a) a lack of a clear spatially defined genetic structure across the majority of genetic markers and (b) the existence of a group of at least 13 SNPs under strong linkage disequilibrium which together separate North Sea sandeel into three haplotype clusters, suggestive of one or more structural variants in the genome. Analyses of the spatial distribution of these putative structural variants suggest at least partial reproductive isolation of sandeel in the western management area along the Scottish coast, supporting a separate management. Our results highlight the importance of the application of a large number of markers to be able to detect weak patterns of differentiation. This study contributes to increasing the genetic knowledge of this important exploited species, and results can be used to improve our understanding of population dynamics and stock structure.

Cytogenetics of the small-sized fish, Copeina guttata (Characiformes, Lebiasinidae): Novel insights into the karyotype differentiation of the family

Lebiasinidae is a small fish family composed by miniature to small-sized fishes with few cytogenetic data (most of them limited to descriptions of diploid chromosome numbers), thus preventing any evolutionary comparative studies at the chromosomal level. In the present study, we are providing, the first cytogenetic data for the red spotted tetra, Copeina guttata, including the standard karyotype, C-banding, repetitive DNA mapping by fluorescence in situ hybridization (FISH) and comparative genomic hybridization (CGH), providing chromosomal patterns and novel insights into the karyotype differentiation of the family. Males and females share diploid chromosome number 2n = 42 and karyotype composed of 2 metacentric (m), 4 submetacentric (sm) and 36 subtelocentric to acrocentric (st-a) chromosomes. Blocks of constitutive heterochromatin were observed in the centromeric and interstitial regions of several chromosomes, in addition to a remarkably large distal block, heteromorphic in size, which fully corresponded with the 18S rDNA sites in the fourth chromosomal pair. This overlap was confirmed by 5S/18S rDNA dual-color FISH. On the other hand, 5S rDNA clusters were situated in the long and short arms of the 2nd and 15th pairs, respectively. No sex-linked karyotype differences were revealed by male/female CGH experiments. The genomic probes from other two lebiasinid species, Lebiasina melanoguttata and Pyrrhulina brevis, showed positive hybridization signals only in the NOR region in the genome of C. guttata. We demonstrated that karyotype diversification in lebiasinids was accompanied by a series of structural and numeric chromosome rearrangements of different types, including particularly fusions and fissions.

Chromosomal Evolution and Evolutionary Relationships of Lebiasina Species (Characiformes, Lebiasinidae)

We present the first cytogenetic data for Lebiasina bimaculata and L. melanoguttata with the aim of (1) investigating evolutionary events within Lebiasina and their relationships with other Lebiasinidae genera and (2) checking the evolutionary relationships between Lebiasinidae and Ctenoluciidae. Both species have a diploid number 2n = 36 with similar karyotypes and microsatellite distribution patterns but present contrasting C-positive heterochromatin and CMA3+ banding patterns. The remarkable interstitial series of C-positive heterochromatin occurring in L. melanoguttata is absent in L. bimaculata. Accordingly, L. bimaculata shows the ribosomal DNA sites as the only GC-rich (CMA3+) regions, while L. melanoguttata shows evidence of a clear intercalated CMA3+ banding pattern. In addition, the multiple 5S and 18S rDNA sites in L. melanogutatta contrast with single sites present in L. bimaculata. Comparative genomic hybridization (CGH) experiments also revealed a high level of genomic differentiation between both species. A polymorphic state of a conspicuous C-positive, CMA3+, and (CGG)n band was found only to occur in L. bimaculata females, and its possible relationship with a nascent sex chromosome system is discussed. Whole chromosome painting (WCP) and CGH experiments indicate that the Lebiasina species examined and Boulengerella maculata share similar chromosomal sequences, thus supporting the relatedness between them and the evolutionary relationships between the Lebiasinidae and Ctenoluciidae families.

Intercontinental karyotype-environment parallelism supports a role for a chromosomal inversion in local adaptation in a seaweed fly

Large chromosomal rearrangements are thought to facilitate adaptation to heterogeneous environments by limiting genomic recombination. Indeed, inversions have been implicated in adaptation along environmental clines and in ecotype specialization. Here, we combine classical ecological studies and population genetics to investigate an inversion polymorphism previously documented in Europe among natural populations of the seaweed fly Coelopa frigida along a latitudinal cline in North America. We test if the inversion is present in North America and polymorphic, assess which environmental conditions modulate the inversion karyotype frequencies, and document the relationship between inversion karyotype and adult size. We sampled nearly 2000 flies from 20 populations along several environmental gradients to quantify associations of inversion frequencies to heterogeneous environmental variables. Genotyping and phenotyping showed a widespread and conserved inversion polymorphism between Europe and America. Variation in inversion frequency was significantly associated with environmental factors, with parallel patterns between continents, indicating that the inversion may play a role in local adaptation. The three karyotypes of the inversion are differently favoured across micro-habitats and represent life-history strategies likely to be maintained by the collective action of several mechanisms of balancing selection. Our study adds to the mounting evidence that inversions are facilitators of adaptation and enhance within-species diversity.

Karyotype diversity and evolutionary trends in the Asian swamp eel Monopterus albus (Synbranchiformes, Synbranchidae): a case of chromosomal speciation?

BACKGROUND

Synbranchidae or swamp eels are fishes belonging to the order Synbranchiformes that occur in both freshwater and occasionally in brackish. They are worldwide distributed in tropical and subtropical rivers of four different continents. A large degree of chromosomal variation has been found in this family, mainly through the use of conventional cytogenetic investigations. Inside this group, a still almost unexplored species under the cytogenetic point of view is the Asian swamp eel Monopterus albus, a widely distributed species throughout Asia. Here, we tested the hypothesis of chromosomal speciation, where a case of sympatric speciation may occur as the primary consequence of chromosomal rearrangements. We performed a comparative chromosomal analysis of M. albus from 22 different localities in Thailand, using distinct staining methods (C-banding, Ag-NO3, and Chromomycin A3), and FISH with repetitive DNA probes (5S rDNA, 18S rDNA, Rex1 element and microsatellite repeats).

RESULTS

This approach evidenced two contrasting karyotypes (named karyomorphs A and B) that varied concerning their 2n and repetitive DNAs distribution, where chromosomal fusions and pericentric inversions were involved in such differentiation. While the karyomorph A has 2n = 24 chromosomes, the karyomorph B has only 2n = 18, both with NF = 24. In addition, karyomorph A contains only acrocentric chromosomes, while karyomorph B contains three unique metacentric pairs. These features highlight that M. albus has already gone through a significant genomic divergence, and may include at least two cryptic species.

CONCLUSIONS

This marked chromosomal differentiation, likely linked to the lifestyle of these fishes, point to the occurrence of a chromosomal speciation scenario, in which fusions and inversions had a prominent role. This highlights the biodiversity of M. albus and justifies its taxonomic revision, since this nominal species may constitute a species complex.