DNA fragility in the parallel evolution of pelvic reduction in stickleback fish

Evolutionary feedback from the environment shapes mechanisms that generate genome variation

Darwin recognized that 'a grand and almost untrodden field of inquiry will be opened, on the causes and laws of variation.' However, because the Modern Synthesis assumes that the intrinsic probability of any individual mutation is unrelated to that mutation's potential adaptive value, attention has been focused on selection rather than on the intrinsic generation of variation. Yet many examples illustrate that the term 'random' mutation, as widely understood, is inaccurate. The probabilities of distinct classes of variation are neither evenly distributed across a genome nor invariant over time, nor unrelated to their potential adaptive value. Because selection acts upon variation, multiple biochemical mechanisms can and have evolved that increase the relative probability of adaptive mutations. In effect, the generation of heritable variation is in a feedback loop with selection, such that those mechanisms that tend to generate variants that survive recurring challenges in the environment would be captured by this survival and thus inherited and accumulated within lineages of genomes. Moreover, because genome variation is affected by a wide range of biochemical processes, genome variation can be regulated. Biochemical mechanisms that sense stress, from lack of nutrients to DNA damage, can increase the probability of specific classes of variation. A deeper understanding of evolution involves attention to the evolution of, and environmental influences upon, the intrinsic variation generated in gametes, in other words upon the biochemical mechanisms that generate variation across generations. These concepts have profound implications for the types of questions that can and should be asked, as omics databases become more comprehensive, detection methods more sensitive, and computation and experimental analyses even more high throughput and thus capable of revealing the intrinsic generation of variation in individual gametes. These concepts also have profound implications for evolutionary theory, which, upon reflection it will be argued, predicts that selection would increase the probability of generating adaptive mutations, in other words, predicts that the ability to evolve itself evolves.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Evolution repeats itself in replicate long-term studies in the wild

The extent to which evolution is repeatable remains debated. Here, we study changes over time in the frequency of cryptic color-pattern morphs in 10 replicate long-term field studies of a stick insect, each spanning at least a decade (across 30 years of total data). We find predictable "up-and-down" fluctuations in stripe frequency in all populations, representing repeatable evolutionary dynamics based on standing genetic variation. A field experiment demonstrates that these fluctuations involve negative frequency-dependent natural selection (NFDS). These fluctuations rely on demographic and selective variability that pushes populations away from equilibrium, such that they can reliably move back toward it via NFDS. Last, we show that the origin of new cryptic forms is associated with multiple structural genomic variants such that which mutations arise affects evolution at larger temporal scales. Thus, evolution from existing variation is predictable and repeatable, but mutation adds complexity even for traits evolving deterministically under natural selection.

Dissection of core promoter syntax through single nucleotide resolution modeling of transcription initiation

Our understanding of how the DNA sequences of cis-regulatory elements encode transcription initiation patterns remains limited. Here we introduce CLIPNET, a deep learning model trained on population-scale PRO-cap data that accurately predicts the position and quantity of transcription initiation with single nucleotide resolution from DNA sequence. Interpretation of CLIPNET revealed a complex regulatory syntax consisting of DNA-protein interactions in five major positions between -200 and +50 bp relative to the transcription start site, as well as more subtle positional preferences among different transcriptional activators. Transcriptional activator and core promoter motifs occupy different positions and play distinct roles in regulating initiation, with the former driving initiation quantity and the latter initiation position. We identified core promoter motifs that explain initiation patterns in the majority of promoters and enhancers, including DPR motifs and AT-rich TBP binding sequences in TATA-less promoters. Our results provide insights into the sequence architecture governing transcription initiation.

Loci associated with cave-derived traits concentrate in specific regions of the Mexican cavefish genome

A major goal of modern evolutionary biology is connecting phenotypic evolution with its underlying genetic basis. The Mexican cavefish (Astyanax mexicanus), a characin fish species comprised of a surface ecotype and a cave-derived ecotype, is well suited as a model to study the genetic mechanisms underlying adaptation to extreme environments. Here we map 206 previously published quantitative trait loci (QTL) for cave-derived traits in A. mexicanus to the newest version of the surface fish genome assembly, AstMex3. This analysis revealed that QTL cluster in the genome more than expected by chance, and this clustering is not explained by the distribution of genes in the genome. To investigate whether certain characteristics of the genome facilitate phenotypic evolution, we tested whether genomic characteristics, such as highly mutagenic CpG sites, are reliable predictors of the sites of trait evolution but did not find any significant trends. Finally, we combined the QTL map with previously collected expression and selection data to identify a list of 36 candidate genes that may underlie the repeated evolution of cave phenotypes, including rgrb which is predicted to be involved in phototransduction. We found this gene has disrupted exons in all non-hybrid cave populations but intact reading frames in surface fish. Overall, our results suggest specific "evolutionary hotspots" in the genome may play significant roles in driving adaptation to the cave environment in Astyanax mexicanus and demonstrate how this compiled dataset can facilitate our understanding of the genetic basis of repeated evolution in the Mexican cavefish.

Dissecting the genetic architecture of quantitative traits using genome-wide identity-by-descent sharing

Additive and dominance genetic variances underlying the expression of quantitative traits are important quantities for predicting short-term responses to selection, but they are notoriously challenging to estimate in most non-model wild populations. Specifically, large-sized or panmictic populations may be characterized by low variance in genetic relatedness among individuals which, in turn, can prevent accurate estimation of quantitative genetic parameters. We used estimates of genome-wide identity-by-descent (IBD) sharing from autosomal SNP loci to estimate quantitative genetic parameters for ecologically important traits in nine-spined sticklebacks (Pungitius pungitius) from a large, outbred population. Using empirical and simulated datasets, with varying sample sizes and pedigree complexity, we assessed the performance of different crossing schemes in estimating additive genetic variance and heritability for all traits. We found that low variance in relatedness characteristic of wild outbred populations with high migration rate can impair the estimation of quantitative genetic parameters and bias heritability estimates downwards. On the other hand, the use of a half-sib/full-sib design allowed precise estimation of genetic variance components and revealed significant additive variance and heritability for all measured traits, with negligible dominance contributions. Genome-partitioning and QTL mapping analyses revealed that most traits had a polygenic basis and were controlled by genes at multiple chromosomes. Furthermore, different QTL contributed to variation in the same traits in different populations suggesting heterogeneous underpinnings of parallel evolution at the phenotypic level. Our results provide important guidelines for future studies aimed at estimating adaptive potential in the wild, particularly for those conducted in outbred large-sized populations.

Optical Mapping: Detecting Genomic Resistance Cassettes in MRSA

Methicillin-resistant Staphylococcus aureus (MRSA) is a multidrug-resistant bacterium with a global presence in healthcare facilities as well as community settings. The resistance of MRSA to beta-lactam antibiotics can be attributed to a mobile genetic element called the staphylococcal cassette chromosome mec (SCCmec), ranging from 23 to 68 kilobase pairs in length. The mec gene complex contained in SCCmec allows MRSA to survive in the presence of penicillin and other beta-lactam antibiotics. We demonstrate that optical mapping (OM) is able to identify the bacterium as S. aureus, followed by an investigation of the presence of kilobase pair range SCCmec elements by examining the associated OM-generated barcode patterns. By employing OM as an alternative to traditional DNA sequencing, we showcase its potential for the detection of complex genetic elements such as SCCmec in MRSA. This approach holds promise for enhancing our understanding of antibiotic resistance mechanisms and facilitating the development of targeted interventions against MRSA infections.

Higher evolutionary dynamics of gene copy number for Drosophila glue genes located near short repeat sequences

BACKGROUND

During evolution, genes can experience duplications, losses, inversions and gene conversions. Why certain genes are more dynamic than others is poorly understood. Here we examine how several Sgs genes encoding glue proteins, which make up a bioadhesive that sticks the animal during metamorphosis, have evolved in Drosophila species.

RESULTS

We examined high-quality genome assemblies of 24 Drosophila species to study the evolutionary dynamics of four glue genes that are present in D. melanogaster and are part of the same gene family - Sgs1, Sgs3, Sgs7 and Sgs8 - across approximately 30 millions of years. We annotated a total of 102 Sgs genes and grouped them into 4 subfamilies. We present here a new nomenclature for these Sgs genes based on protein sequence conservation, genomic location and presence/absence of internal repeats. Two types of glue genes were uncovered. The first category (Sgs1, Sgs3x, Sgs3e) showed a few gene losses but no duplication, no local inversion and no gene conversion. The second group (Sgs3b, Sgs7, Sgs8) exhibited multiple events of gene losses, gene duplications, local inversions and gene conversions. Our data suggest that the presence of short "new glue" genes near the genes of the latter group may have accelerated their dynamics.

CONCLUSIONS

Our comparative analysis suggests that the evolutionary dynamics of glue genes is influenced by genomic context. Our molecular, phylogenetic and comparative analysis of the four glue genes Sgs1, Sgs3, Sgs7 and Sgs8 provides the foundation for investigating the role of the various glue genes during Drosophila life.

Mutation and evolution: Conceptual possibilities

Although random mutation is central to models of evolutionary change, a lack of clarity remains regarding the conceptual possibilities for thinking about the nature and role of mutation in evolution. We distinguish several claims at the intersection of mutation, evolution, and directionality and then characterize a previously unrecognized category: complex conditioned mutation. Empirical evidence in support of this category suggests that the historically famous fluctuation test should be revisited, and new experiments should be undertaken with emerging experimental techniques to facilitate detecting mutation rates within specific loci at an ultra-high, individual base pair resolution.

Morphological description of spontaneous pelvic fin loss in a neotropical cichlid fish

Pelvic fins are a characteristic structure of the vertebrate Bauplan. Yet, pelvic fin loss has occurred repeatedly across a wide diversity of other lineages of tetrapods and at least 48 times in teleost fishes. This pelvic finless condition is often associated with other morphological features such as body elongation, loss of additional structures, and bilateral asymmetry. However, despite the remarkable diversity in the several thousand cichlid fish species, none of them are characterized by the complete absence of pelvic fins. Here, we examined the musculoskeletal structure and associated bilateral asymmetry in Midas cichlids (Amphilophus cf. citrinellus) that lost their pelvic fins spontaneously in the laboratory. Due to this apparent mutational loss of the pelvic girdle and fins, the external and internal anatomy are described in a series of "normal" Midas individuals and their pelvic finless sibling tankmates. First, other traits associated with teleost pelvic fin loss, the genetic basis of pelvic fin loss, and the potential for pleiotropic effects of these genes on other traits in teleosts were all reviewed. Using these traits as a guide, we investigated whether other morphological differences were associated with the pelvic girdle/fin loss. The mean values of the masses of muscle of the pectoral fin, fin ray numbers in the unpaired fins, and oral jaw tooth numbers did not differ between the two pelvic fin morphotypes. However, significant differences in meristic values of the paired traits assessed were observed for the same side of the body between morphotypes. Notably, bilateral asymmetry was found exclusively for the posterior lateral line scales. Finally, we found limited evidence of pleiotropic effects, such as lateral line scale numbers and fluctuating asymmetry between the Midas pelvic fin morphotypes. The fast and relatively isolated changes in the Midas cichlids suggest minor but interesting pleiotropic effects could accompany loss of cichlid pelvic fins.

Genome architecture underlying salinity adaptation in the invasive copepod Eurytemora affinis species complex: A review

With climate change, habitat salinity is shifting rapidly throughout the globe. In addition, many destructive freshwater invaders are recent immigrants from saline habitats. Recently, populations of the copepod Eurytemora affinis species complex have invaded freshwater habitats multiple times independently from saline estuaries on three continents. This review discusses features of this species complex that could enhance their evolutionary potential during rapid environmental change. Remarkably, across independent freshwater invasions, natural selection has repeatedly favored the same alleles far more than expected. This high degree of parallelism is surprising, given the expectation of nonparallel evolution for polygenic adaptation. Factors such as population structure and the genome architecture underlying critical traits under selection might help drive rapid adaptation and parallel evolution. Given the preponderance of saline-to-freshwater invasions and climate-induced salinity change, the principles found here could provide invaluable insights into mechanisms operating in other systems and the potential for adaptation in a changing planet.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

Replication stress in mammalian embryo development, differentiation, and reprogramming

Duplicating a genome of 3 billion nucleotides is challenged by a variety of obstacles that can cause replication stress and affect the integrity of the genome. Recent studies show that replication fork slowing and stalling is prevalent in early mammalian development, resulting in genome instability and aneuploidy, and constituting a barrier to development in human reproduction. Genome instability resulting from DNA replication stress is a barrier to the cloning of animals and to the reprogramming of differentiated cells to induced pluripotent stem cells, as well as a barrier to cell transformation. Remarkably, the regions most impacted by replication stress are shared in these different cellular contexts, affecting long genes and flanking intergenic areas. In this review we integrate our knowledge of DNA replication stress in mammalian embryos, in programming, and in reprogramming, and we discuss a potential role for fragile sites in sensing replication stress and restricting cell cycle progression in health and disease.

Understanding natural selection and similarity: Convergent, parallel and repeated evolution

Parallel and convergent evolution offer some of the most compelling evidence for the significance of natural selection in evolution, as the emergence of similar adaptive solutions is unlikely to occur by random chance alone. However, these terms are often employed inconsistently, leading to misinterpretation and confusion, and recently proposed definitions have unintentionally diminished the emphasis on the evolution of similar adaptive solutions. Here, I examine various conceptual frameworks and definitions related to parallel and convergent evolution and propose a consolidated framework that enhances our comprehension of these evolutionary patterns. The primary aim of this framework is to harmonize the concepts of parallel and convergent evolution together with natural selection and the idea of similarity. Both concepts involve the evolution of similar adaptive solutions as a result of environmental challenges. The distinction lies in ancestral phenotypes. Parallel evolution takes place when the ancestral phenotypes (before selection) of the lineages are similar. Convergent evolution happens when the lineages have distinct ancestral phenotypes (before selection). Because an ancestral-based distinction will inevitably lead to cases where uncertainty in the distinction may arise, the framework includes a general term, repeated evolution, which can be used as a term applying to the evolution of similar phenotypes and genotypes as well as similar responses to environmental pressures. Based on the argument that genetic similarity may frequently arise without selection, the framework posits that the similarity of genetic sequences is not of great interest unless linked to the actions of natural selection or to the origins (mutation, standing genetic variation, gene flow) and locations of the similar sequences.

Revisiting the Role of Genetic Variation in Adaptation

AbstractEvolutionary biologists have thought about the role of genetic variation during adaptation for a very long time-before we understood the organization of the genetic code, the provenance of genetic variation, and how such variation influenced the phenotypes on which natural selection acts. Half a century after the discovery of the structure of DNA and the unraveling of the genetic code, we have a rich understanding of these problems and the means to both delve deeper and widen our perspective across organisms and natural populations. The 2022 Vice Presidential Symposium of the American Society of Naturalists highlighted examples of recent insights into the role of genetic variation in adaptive processes, which are compiled in this special section. The work was conducted in different parts of the world, included theoretical and empirical studies with diverse organisms, and addressed distinct aspects of how genetic variation influences adaptation. In our introductory article to the special section, we discuss some important recent insights about the generation and maintenance of genetic variation, its impacts on phenotype and fitness, its fate in natural populations, and its role in driving adaptation. By placing the special section articles in the broader context of recent developments, we hope that this overview will also serve as a useful introduction to the field.

Whole-genome Comparisons Identify Repeated Regulatory Changes Underlying Convergent Appendage Evolution in Diverse Fish Lineages

Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 percomorph fish species that span over 100 million years of evolution and either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA, a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, syngnathid, and synbranchid species with caudal fin reduction, and alters caudal fin development when targeted by genome editing. Our study illustrates a broadly applicable strategy for mapping phenotypes to genotypes across a tree of vertebrate species and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes across 100 million years of fish evolution.

Developmental Selection and the Perception of Mutation Bias

The notion that mutations are random relative to their fitness effects is central to the Neo-Darwinian view of evolution. However, a recent interpretation of the patterns of mutation accumulation in the genome of Arabidopsis thaliana has challenged this notion, arguing for the presence of a targeted DNA repair mechanism that causes a nonrandom association of mutation rates and fitness effects. Specifically, this mechanism was suggested to cause a reduction in the rates of mutations on essential genes, thus lowering the rates of deleterious mutations. Central to this argument were attempts to rule out selection at the population level. Here, we offer an alternative and parsimonious interpretation of the patterns of mutation accumulation previously attributed to mutation bias, showing how they can instead or additionally be caused by developmental selection, that is selection occurring at the cellular level during the development of a multicellular organism. Thus, the depletion of deleterious mutations in A. thaliana may indeed be the result of a selective process, rather than a bias in mutation. More broadly, our work highlights the importance of considering development in the interpretation of population-genetic analyses of multicellular organisms, and it emphasizes that efforts to identify mechanisms involved in mutational biases should explicitly account for developmental selection.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

The power of "controllers": Transposon-mediated duplicated genes evolve towards neofunctionalization

Since the discovery of the first transposon by Dr. Barbara McClintock, the prevalence and diversity of transposable elements (TEs) have been gradually recognized. As fundamental genetic components, TEs drive organismal evolution not only by contributing functional sequences (e.g., regulatory elements or "controllers" as phrased by Dr. McClintock) but also by shuffling genomic sequences. In the latter respect, TE-mediated gene duplications have contributed to the origination of new genes and attracted extensive interest. In response to the development of this field, we herein attempt to provide an overview of TE-mediated duplication by focusing on common rules emerging across duplications generated by different TE types. Specifically, despite the huge divergence of transposition machinery across TEs, we identify three common features of various TE-mediated duplication mechanisms, including end bypass, template switching, and recurrent transposition. These three features lead to one common functional outcome, namely, TE-mediated duplicates tend to be subjected to exon shuffling and neofunctionalization. Therefore, the intrinsic properties of the mutational mechanism constrain the evolutionary trajectories of these duplicates. We finally discuss the future of this field including an in-depth characterization of both the duplication mechanisms and functions of TE-mediated duplicates.

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Selection-driven trait loss in independently evolved cavefish populations

Laboratory studies have demonstrated that a single phenotype can be produced by many different genotypes; however, in natural systems, it is frequently found that phenotypic convergence is due to parallel genetic changes. This suggests a substantial role for constraint and determinism in evolution and indicates that certain mutations are more likely to contribute to phenotypic evolution. Here we use whole genome resequencing in the Mexican tetra, Astyanax mexicanus, to investigate how selection has shaped the repeated evolution of both trait loss and enhancement across independent cavefish lineages. We show that selection on standing genetic variation and de novo mutations both contribute substantially to repeated adaptation. Our findings provide empirical support for the hypothesis that genes with larger mutational targets are more likely to be the substrate of repeated evolution and indicate that features of the cave environment may impact the rate at which mutations occur.

Dynamic alternative DNA structures in biology and disease

Repetitive elements in the human genome, once considered 'junk DNA', are now known to adopt more than a dozen alternative (that is, non-B) DNA structures, such as self-annealed hairpins, left-handed Z-DNA, three-stranded triplexes (H-DNA) or four-stranded guanine quadruplex structures (G4 DNA). These dynamic conformations can act as functional genomic elements involved in DNA replication and transcription, chromatin organization and genome stability. In addition, recent studies have revealed a role for these alternative structures in triggering error-generating DNA repair processes, thereby actively enabling genome plasticity. As a driving force for genetic variation, non-B DNA structures thus contribute to both disease aetiology and evolution.

Adaptability of Bony Armor Elements of the Threespine Stickleback Gasterosteus aculeatus (Teleostei: Gasterosteidae): Ecological and Evolutionary Insights from Symmetry Analyses

Differentiation in the defensive armor of the threespine stickleback, Gasterosteus aculeatus, is caused by predator-driven divergent selection. Most studies considered armor traits related to swimming behavior, hence combining pre- and post-capture responses to gape-limited predators. Here, we focus exclusively on the defensive complex (DC), the post-capture predator defense. This complex consists of a series of bony elements surrounding the anterior part of the abdomen. Relaxation from predation pressure not only drives reduction of bony elements but is also expected to increase asymmetry in the DC. To test this hypothesis, we used four Austrian freshwater populations that differed distinctly in the formation of the DC. We found significant left–right asymmetries in the DC in the population with a distinctly reduced DC and, surprisingly, also in the population with a significantly enhanced DC. These populations occur in vastly different habitats (stream and lake) characterized by distinct regimes of gape-limited predators (none vs. many). Apparently, both a shift to very low and very high pressure by gape-limited predators can boost asymmetry. We conclude that greater asymmetries in the two populations at the opposite ends of the predatory gradient result from an ongoing process of adaptation to decreased or increased environmental stress.

Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review

Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.

Prevalent Introgression Underlies Convergent Evolution in the Diversification of Pungitius Sticklebacks

New mutations and standing genetic variations contribute significantly to repeated phenotypic evolution in sticklebacks. However, less is known about the role of introgression in this process. We analyzed taxonomically and geographically comprehensive genomic data from Pungitius sticklebacks to decipher the extent of introgression and its consequences for the diversification of this genus. Our results demonstrate that introgression is more prevalent than suggested by earlier studies. Although gene flow was generally bidirectional, it was often asymmetric and left unequal genomic signatures in hybridizing species, which might, at least partly, be due to biased hybridization and/or population size differences. In several cases, introgression of variants from one species to another was accompanied by transitions of pelvic and/or lateral plate structures-important diagnostic traits in Pungitius systematics-and frequently left signatures of adaptation in the core gene regulatory networks of armor trait development. This finding suggests that introgression has been an important source of genetic variation and enabled phenotypic convergence among Pungitius sticklebacks. The results highlight the importance of introgression of genetic variation as a source of adaptive variation underlying key ecological and taxonomic traits. Taken together, our study indicates that introgression-driven convergence likely explains the long-standing challenges in resolving the taxonomy and systematics of this small but phenotypically highly diverse group of fish.

Noncanonical DNA structures are drivers of genome evolution

In addition to the canonical right-handed double helix, other DNA structures, termed 'non-B DNA', can form in the genomes across the tree of life. Non-B DNA regulates multiple cellular processes, including replication and transcription, yet its presence is associated with elevated mutagenicity and genome instability. These discordant cellular roles fuel the enormous potential of non-B DNA to drive genomic and phenotypic evolution. Here we discuss recent studies establishing non-B DNA structures as novel functional elements subject to natural selection, affecting evolution of transposable elements (TEs), and specifying centromeres. By highlighting the contributions of non-B DNA to repeated evolution and adaptation to changing environments, we conclude that evolutionary analyses should include a perspective of not only DNA sequence, but also its structure.

Whole-genome comparisons identify repeated regulatory changes underlying convergent appendage evolution in diverse fish lineages

Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 wild fish species that either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences (CONDELs) in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA , a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, s yngnathid, and synbranchid species with caudal fin reduction, and which alters caudal fin development when targeted by genome editing. Our study illustrates a general strategy for mapping phenotypes to genotypes across a tree of vertebrate species, and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes during 100 million years of fish evolution.

Repeated genetic divergence plays a minor role in repeated phenotypic divergence of lake-stream stickleback

Recent studies have shown that the repeated evolution of similar phenotypes in response to similar ecological conditions (here "parallel evolution") often occurs through mutations in the same genes. However, many previous studies have focused on known candidate genes in a limited number of systems. Thus, the question of how often parallel phenotypic evolution is due to parallel genetic changes remains open. Here, we used quantitative trait locus (QTL) mapping in F2 intercrosses between lake and stream threespine stickleback (Gasterosteus aculeatus) from four independent watersheds on Vancouver Island, Canada to determine whether the same QTL underlie divergence in the same phenotypes across, between, and within watersheds. We find few parallel QTL, even in independent crosses from the same watershed or for phenotypes that have diverged in parallel. These findings suggest that different mutations can lead to similar phenotypes. The low genetic repeatability observed in these lake-stream systems contrasts with the higher genetic repeatability observed in other stickleback systems. We speculate that differences in evolutionary history, gene flow, and/or the strength and direction of selection might explain these differences in genetic parallelism and emphasize that more work is needed to move beyond documenting genetic parallelism to identifying the underlying causes.

A Phylogenomic Assessment of Processes Underpinning Convergent Evolution in Open-Habitat Chats

Insights into the processes underpinning convergent evolution advance our understanding of the contributions of ancestral, introgressed, and novel genetic variation to phenotypic evolution. Phylogenomic analyses characterizing genome-wide gene tree heterogeneity can provide first clues about the extent of ILS and of introgression and thereby into the potential of these processes or (in their absence) the need to invoke novel mutations to underpin convergent evolution. Here, we were interested in understanding the processes involved in convergent evolution in open-habitat chats (wheatears of the genus Oenanthe and their relatives). To this end, based on whole-genome resequencing data from 50 taxa of 44 species, we established the species tree, characterized gene tree heterogeneity, and investigated the footprints of ILS and introgression within the latter. The species tree corroborates the pattern of abundant convergent evolution, especially in wheatears. The high levels of gene tree heterogeneity in wheatears are explained by ILS alone only for 30% of internal branches. For multiple branches with high gene tree heterogeneity, D-statistics and phylogenetic networks identified footprints of introgression. Finally, long branches without extensive ILS between clades sporting similar phenotypes provide suggestive evidence for the role of novel mutations in the evolution of these phenotypes. Together, our results suggest that convergent evolution in open-habitat chats involved diverse processes and highlight that phenotypic diversification is often complex and best depicted as a network of interacting lineages.

Z-DNA and Z-RNA: Methods-Past and Future

A quote attributed to Yogi Berra makes the observation that "It's tough to make predictions, especially about the future," highlighting the difficulties posed to an author writing a manuscript like the present. The history of Z-DNA shows that earlier postulates about its biology have failed the test of time, both those from proponents who were wildly enthusiastic in enunciating roles that till this day still remain elusive to experimental validation and those from skeptics within the larger community who considered the field a folly, presumably because of the limitations in the methods available at that time. If anything, the biological roles we now know for Z-DNA and Z-RNA were not anticipated by anyone, even when those early predictions are interpreted in the most favorable way possible. The breakthroughs in the field were made using a combination of methods, especially those based on human and mouse genetic approaches informed by the biochemical and biophysical characterization of the Zα family of proteins. The first success was with the p150 Zα isoform of ADAR1 (adenosine deaminase RNA specific), with insights into the functions of ZBP1 (Z-DNA-binding protein 1) following soon after from the cell death community. Just as the replacement of mechanical clocks by more accurate designs changed expectations about navigation, the discovery of the roles assigned by nature to alternative conformations like Z-DNA has forever altered our view of how the genome operates. These recent advances have been driven by better methodology and by better analytical approaches. This article will briefly describe the methods that were key to these discoveries and highlight areas where new method development is likely to further advance our knowledge.

Genes that are Used Together are More Likely to be Fused Together in Evolution by Mutational Mechanisms: A Bioinformatic Test of the Used-Fused Hypothesis

Cases of parallel or recurrent gene fusions in evolution as well as in genetic disease and cancer are difficult to explain, because unlike point mutations, they can require the repetition of a similar configuration of multiple breakpoints rather than the repetition of a single point mutation. The used-together-fused-together hypothesis holds that genes that are used together repeatedly and persistently in a specific context are more likely to undergo fusion mutation in the course of evolution for mechanistic reasons. This hypothesis offers to explain gene fusion in both evolution and disease under one umbrella. Using bioinformatic data, we tested this hypothesis against alternatives, including that all gene pairs can fuse by random mutation, but among pairs thus fused, those that had interacted previously are more likely to be favored by selection. Results show that across multiple measures of gene interaction, human genes whose orthologs are fused in one or more species are more likely to interact with each other than random pairs of genes of the same genomic distance between pair members; that an overlap exists between genes that fused in the course of evolution in non-human species and genes that undergo fusion in human cancers; and that across six primate species studied, fusions predominate over fissions and exhibit substantial evolutionary parallelism. Together, these results support the used-together-fused-together hypothesis over its alternatives. Multiple implications are discussed, including the relevance of mutational mechanisms to the evolution of genome organization, to the distribution of fitness effects of mutation, to evolutionary parallelism and more.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-022-09579-9.

Detection of Z-DNA Structures in Supercoiled Genome

Z-DNAs are nucleic acid secondary structures that form within a special pattern of nucleotides and are promoted by DNA supercoiling. Through Z-DNA formation, DNA encodes information by dynamic changes in its secondary structure. A growing body of evidence indicates that Z-DNA formation can play a role in gene regulation; it can affect chromatin architecture and demonstrates its association with genomic instability, genetic diseases, and genome evolution. Many functional roles of Z-DNA are yet to be discovered highlighting the need for techniques to detect genome-wide folding of DNA into this structure. Here, we describe an approach to convert linear genome into supercoiled genome sponsoring Z-DNA formation. Applying permanganate-based methodology and high-throughput sequencing to supercoiled genome allows genome-wide detection of single-stranded DNA. Single-stranded DNA is characteristic of the junctions between the classical B-form of DNA and Z-DNA. Consequently, analysis of single-stranded DNA map provides snapshots of the Z-DNA conformation over the whole genome.

Methods to Study Z-DNA-Induced Genetic Instability

Alternative DNA structures that differ from the canonical B-DNA double helix, including Z-DNA, have received much attention recently due to their impact on DNA metabolic processes, including replication, transcription, and genome maintenance. Non-B-DNA-forming sequences can also stimulate genetic instability associated with disease development and evolution. Z-DNA can stimulate different types of genetic instability events in different species, and several different assays have been established to detect Z-DNA-induced DNA strand breaks and mutagenesis in prokaryotic and eukaryotic systems. In this chapter, we will introduce some of these methods including Z-DNA-induced mutation screening and detection of Z-DNA-induced strand breaks in mammalian cells, yeast, and mammalian cell extracts. Results from these assays should provide better insight into the mechanisms of Z-DNA-related genetic instability in different eukaryotic model systems.

The remodeling of Z-DNA in the mammalian germ line

We recently discovered a novel biological process, the scheduled remodeling of Z-DNA structures in the developing fetal mouse male germ cells [Nat. Cell Biol. 24, 1141-1153]. This process affects purine/pyrimidine dinucleotide repeat (PPR) rich sequences, which can form stable left-handed Z-DNA structures. The protein that carries out this function is identified as ZBTB43, member of a large family of ZBTB proteins. Z-DNA remodeling by ZBTB43 not only coincides with global remodeling of DNA methylation and chromatin events in the male germ line, but it also is a prerequisite for de novo DNA methylation. When ZBTB43 changes DNA structure from the left-handed zigzag shaped Z-DNA to the regular smooth right-handed B-DNA, it also generates a suitable substrate for the de novo DNA methyltransferase, DNMT3A. By instructing de novo DNA methylation at PPRs in prospermatogonia, ZBTB43 safeguards epigenomic integrity of the male gamete. PPRs are fragile sequences, sites of large deletions and rearrangements in mammalian cells, and this fragility is thought to be due to Z-DNA structure formation rather than the sequence itself. This idea is now supported by the in vivo finding that DNA double strand breaks accumulate in mutant prospermatogonia which lack ZBTB43-dependent Z-DNA remodeling. If unrepaired, double stranded DNA breaks can lead to germ line mutations. Therefore, by preventing such breaks ZBTB43 is critical for guarding genome stability between generations. Here, we discuss the significance and implications of these findings in more detail.

Adaptive sequence divergence forged new neurodevelopmental enhancers in humans

Searches for the genetic underpinnings of uniquely human traits have focused on human-specific divergence in conserved genomic regions, which reflects adaptive modifications of existing functional elements. However, the study of conserved regions excludes functional elements that descended from previously neutral regions. Here, we demonstrate that the fastest-evolved regions of the human genome, which we term "human ancestor quickly evolved regions" (HAQERs), rapidly diverged in an episodic burst of directional positive selection prior to the human-Neanderthal split, before transitioning to constraint within hominins. HAQERs are enriched for bivalent chromatin states, particularly in gastrointestinal and neurodevelopmental tissues, and genetic variants linked to neurodevelopmental disease. We developed a multiplex, single-cell in vivo enhancer assay to discover that rapid sequence divergence in HAQERs generated hominin-unique enhancers in the developing cerebral cortex. We propose that a lack of pleiotropic constraints and elevated mutation rates poised HAQERs for rapid adaptation and subsequent susceptibility to disease.

Relaxed risk of predation drives parallel evolution of stickleback behavior

The occurrence of similar phenotypes in multiple independent populations derived from common ancestral conditions (viz. parallel evolution) is a testimony of evolution by natural selection. Parallel evolution implies that populations share a common phenotypic response to a common selection pressure associated with habitat similarity. Examples of parallel evolution at genetic and phenotypic levels are fairly common, but the driving selective agents often remain elusive. Similarly, the role of phenotypic plasticity in facilitating early stages of parallel evolution is unclear. We investigated whether the relaxation of predation pressure associated with the colonization of freshwater ponds by nine-spined sticklebacks (Pungitius pungitius) likely explains the divergence in complex behaviors between marine and pond populations, and whether this divergence is parallel. Using laboratory-raised individuals exposed to different levels of perceived predation risk, we calculated vectors of phenotypic divergence for four behavioral traits between habitats and predation risk treatments. We found a significant correlation between the directions of evolutionary divergence and phenotypic plasticity, suggesting that divergence in behavior between habitats is aligned with the response to relaxation of predation pressure. Finally, we show alignment across multiple pairs of populations, and that relaxation of predation pressure has likely driven parallel evolution of behavior in this species.

A DNA Replication Mechanism Can Explain Structural Variation at the Pigeon Recessive Red Locus

For species to adapt to their environment, evolution must act upon genetic variation that is present in the population. Elucidating the molecular mechanisms that give rise to this variation is thus of crucial importance for understanding how organisms evolve. In addition to variation caused by point mutations, structural variation (deletions, duplications, inversions, translocations) is also an important source of variety. Mechanisms involving recombination, transposition and retrotransposition, and replication have been proposed for generating structural variation, and each are capable of explaining certain rearrangements. In this study, we conduct a detailed analysis of two partially overlapping rearrangements (e1 and e2 allele) in domestic rock pigeon (Columba livia) which are both associated with the recessive red phenotype. We find that a replicative mechanism is best able to explain the complex architecture of the e1 allele, and is also compatible with the simpler architecture of the e2 allele as well.

Single-molecule fluorescence imaging techniques reveal molecular mechanisms underlying deoxyribonucleic acid damage repair

Advances in single-molecule techniques have uncovered numerous biological secrets that cannot be disclosed by traditional methods. Among a variety of single-molecule methods, single-molecule fluorescence imaging techniques enable real-time visualization of biomolecular interactions and have allowed the accumulation of convincing evidence. These techniques have been broadly utilized for studying DNA metabolic events such as replication, transcription, and DNA repair, which are fundamental biological reactions. In particular, DNA repair has received much attention because it maintains genomic integrity and is associated with diverse human diseases. In this review, we introduce representative single-molecule fluorescence imaging techniques and survey how each technique has been employed for investigating the detailed mechanisms underlying DNA repair pathways. In addition, we briefly show how live-cell imaging at the single-molecule level contributes to understanding DNA repair processes inside cells.

Extreme intraspecific divergence in mitochondrial haplotypes makes the threespine stickleback fish an emerging evolutionary mutant model for mito-nuclear interactions

Mitochondrial DNA is primarily maternally inherited in most animals and evolves about 10 times faster than biparentally inherited nuclear DNA. Mitochondrial dysfunction (mt-dys) arises when interactions between the co-evolving mitochondrial and nuclear genomes are perturbed in essential processes like oxidative phosphorylation (OXPHOS). Over time mt-dys can lead to mitochondrial diseases (mt-diseases), which are surprisingly prevalent and include common diseases such as Alzheimer’s, Parkinson’s, and diabetes. Unfortunately, the strong impact that intraspecific mitochondrial and nuclear genetic variation has on mt-disease complicates its study and the development of effective treatments. Animal models have advanced our understanding of mt-disease but their relevance to human conditions is often limited by their relatively low nuclear genetic diversity. Many traditional laboratory models also typically have a single mitochondrial haplotype (mitotype), in stark contrast to over 5,000 mitotypes in humans worldwide. The threespine stickleback fish has an evolutionary history that has made it a favorable evolutionary mutant model (EMM) for studying mito-nuclear interactions and possibly mt-diseases. EMMs are species with naturally evolved states that mimic maladaptive human diseases. In threespine stickleback, a period of isolation followed by introgression of the mitochondrial genome from a sister species resulted in the maintenance of two distinct mitochondrial haplotypes which continue to segregate within many populations of wild stickleback. The existence of two mitogenomes segregating in numerous genetically diverse populations provides a unique system for exploring complex mito-nuclear dynamics. Here we provide the first complete coding region analysis of the two threespine stickleback mitotypes, whose mitogenomic divergence exceeds that of other mammalian models for mitochondrial disease and even that between ancient and modern humans. We find that divergence is not uniform across the mitogenome, but primarily impacts protein coding genes, and significantly impacts proteins in Complex I of OXPHOS. The full characterization of these highly divergent intraspecific mitotypes provides a foundation for the development of threespine stickleback as an EMM for mito-nuclear interactions.

Joint analysis of microsatellites and flanking sequences enlightens complex demographic history of interspecific gene flow and vicariance in rear-edge oak populations

Inference of recent population divergence requires fast evolving markers and necessitates to differentiate shared genetic variation caused by ancestral polymorphism and gene flow. Theoretical research shows that the use of compound marker systems integrating linked polymorphisms with different mutational dynamics, such as a microsatellite and its flanking sequences, can improve estimation of population structure and inference of demographic history, especially in the case of complex population dynamics. However, empirical application in natural populations has so far been limited by lack of suitable methods for data collection. A solution comes from the development of sequence-based microsatellite genotyping which we used to study molecular variation at 36 sequenced nuclear microsatellites in seven Quercus canariensis and four Q. faginea rear-edge populations across Algeria. We aim to decipher their taxonomic relationship, past evolutionary history and recent demographic trajectory. First, we compare the estimation of population genetics parameters and simulation-based inference of demographic history from microsatellite sequence alone, flanking sequence alone or the combination of linked microsatellite and flanking sequence variation. Second, we apply random forest approximate Bayesian computation to identify which of these sequence types is most informative. Whereas analysing microsatellite variation alone indicates recent interspecific gene flow, additional information gained by integrating nucleotide variation in flanking sequences, by reducing homoplasy, suggests ancient interspecific gene flow followed by drift in isolation instead. The weight of each polymorphism in the inference also demonstrates the value of linked variations with contrasted mutation dynamic to improve estimation of both demographic and mutational parameters.

Multivariate selection and the making and breaking of mutational pleiotropy

The role of mutations have been subject to many controversies since the formation of the Modern Synthesis of evolution in the early 1940ties. Geneticists in the early half of the twentieth century tended to view mutations as a limiting factor in evolutionary change. In contrast, natural selection was largely viewed as a “sieve” whose main role was to sort out the unfit but which could not create anything novel alone. This view gradually changed with the development of mathematical population genetics theory, increased appreciation of standing genetic variation and the discovery of more complex forms of selection, including balancing selection. Short-term evolutionary responses to selection are mainly influenced by standing genetic variation, and are predictable to some degree using information about the genetic variance–covariance matrix (G) and the strength and form of selection (e. g. the vector of selection gradients, β). However, predicting long-term evolution is more challenging, and requires information about the nature and supply of novel mutations, summarized by the mutational variance–covariance matrix (M). Recently, there has been increased attention to the role of mutations in general and M in particular. Some evolutionary biologists argue that evolution is largely mutation-driven and claim that mutation bias frequently results in mutation-biased adaptation. Strong similarities between G and M have also raised questions about the non-randomness of mutations. Moreover, novel mutations are typically not isotropic in their phenotypic effects and mutational pleiotropy is common. Here I discuss the evolutionary origin and consequences of mutational pleiotropy and how multivariate selection directly shapes G and indirectly M through changed epistatic relationships. I illustrate these ideas by reviewing recent literature and models about correlational selection, evolution of G and M, sexual selection and the fitness consequences of sexual antagonism.

An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication

Teleost fishes are ancient tetraploids descended from an ancestral whole-genome duplication that may have contributed to the impressive diversification of this clade. Whole-genome duplications can occur via self-doubling (autopolyploidy) or via hybridization between different species (allopolyploidy). The mode of tetraploidization conditions evolutionary processes by which duplicated genomes return to diploid meiotic pairing, and subsequent genetic divergence of duplicated genes (cytological and genetic rediploidization). How teleosts became tetraploid remains unresolved, leaving a fundamental gap in the interpretation of their functional evolution. As a result of the whole-genome duplication, identifying orthologous and paralogous genomic regions across teleosts is challenging, hindering genome-wide investigations into their polyploid history. Here, we combine tailored gene phylogeny methodology together with a state-of-the-art ancestral karyotype reconstruction to establish the first high-resolution comparative atlas of paleopolyploid regions across 74 teleost genomes. We then leverage this atlas to investigate how rediploidization occurred in teleosts at the genome-wide level. We uncover that some duplicated regions maintained tetraploidy for more than 60 million years, with three chromosome pairs diverging genetically only after the separation of major teleost families. This evidence suggests that the teleost ancestor was an autopolyploid. Further, we find evidence for biased gene retention along several duplicated chromosomes, contradicting current paradigms that asymmetrical evolution is specific to allopolyploids. Altogether, our results offer novel insights into genome evolutionary dynamics following ancient polyploidizations in vertebrates.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

Z-DNA is remodelled by ZBTB43 in prospermatogonia to safeguard the germline genome and epigenome

Mutagenic purine–pyrimidine repeats can adopt the left-handed Z-DNA conformation. DNA breaks at potential Z-DNA sites can lead to somatic mutations in cancer or to germline mutations that are transmitted to the next generation. It is not known whether any mechanism exists in the germ line to control Z-DNA structure and DNA breaks at purine–pyrimidine repeats. Here we provide genetic, epigenomic and biochemical evidence for the existence of a biological process that erases Z-DNA specifically in germ cells of the mouse male foetus. We show that a previously uncharacterized zinc finger protein, ZBTB43, binds to and removes Z-DNA, preventing the formation of DNA double-strand breaks. By removing Z-DNA, ZBTB43 also promotes de novo DNA methylation at CG-containing purine–pyrimidine repeats in prospermatogonia. Therefore, the genomic and epigenomic integrity of the species is safeguarded by remodelling DNA structure in the mammalian germ line during a critical window of germline epigenome reprogramming.

Meng et al. show that ZBTB43 alters Z-DNA structures to prevent deleterious double-strand breaks and promote DNA methylation at purine–pyrimidine repeats in the mouse germ line.

Attempting genetic inference from directional asymmetry during convergent hindlimb reduction in squamates

Loss and reduction in paired appendages are common in vertebrate evolution. How often does such convergent evolution depend on similar developmental and genetic pathways? For example, many populations of the threespine stickleback and ninespine stickleback (Gasterosteidae) have independently evolved pelvic reduction, usually based on independent mutations that caused reduced Pitx1 expression. Reduced Pitx1 expression has also been implicated in pelvic reduction in manatees. Thus, hindlimb reduction stemming from reduced Pitx1 expression has arisen independently in groups that diverged tens to hundreds of millions of years ago, suggesting a potential for repeated use of Pitx1 across vertebrates. Notably, hindlimb reduction based on the reduction in Pitx1 expression produces left-larger directional asymmetry in the vestiges. We used this phenotypic signature as a genetic proxy, testing for hindlimb directional asymmetry in six genera of squamate reptiles that independently evolved hindlimb reduction and for which genetic and developmental tools are not yet developed: Agamodon anguliceps, Bachia intermedia, Chalcides sepsoides, Indotyphlops braminus, Ophisaurus attenuatuas and O. ventralis, and Teius teyou. Significant asymmetry occurred in one taxon, Chalcides sepsoides, whose left-side pelvis and femur vestiges were 18% and 64% larger than right-side vestiges, respectively, suggesting modification in Pitx1 expression in that species. However, there was either right-larger asymmetry or no directional asymmetry in the other five taxa, suggesting multiple developmental genetic pathways to hindlimb reduction in squamates and the vertebrates more generally.

Genomic changes underlying repeated niche shifts in an adaptive radiation

In adaptive radiations, single lineages rapidly diversify by adapting to many new niches. Little is known yet about the genomic mechanisms involved, that is, the source of genetic variation or genomic architecture facilitating or constraining adaptive radiation. Here, we investigate genomic changes associated with repeated invasion of many different freshwater niches by threespine stickleback in the Haida Gwaii archipelago, Canada, by resequencing single genomes from one marine and 28 freshwater populations. We find 89 likely targets of parallel selection in the genome that are enriched for old standing genetic variation. In contrast to theoretical expectations, their genomic architecture is highly dispersed with little clustering. Candidate genes and genotype-environment correlations match the three major environmental axes predation regime, light environment, and ecosystem size. In a niche space with these three dimensions, we find that the more divergent a new niche from the ancestral marine habitat, the more loci show signatures of parallel selection. Our findings suggest that the genomic architecture of parallel adaptation in adaptive radiation depends on the steepness of ecological gradients and the dimensionality of the niche space.

Genetic basis of speciation and adaptation: from loci to causative mutations

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.