Genome-wide analysis of histone modifications can contribute to the identification of candidate cis-regulatory regions in the threespine stickleback fish.

BACKGROUND: Cis-regulatory mutations often underlie phenotypic evolution. However, because identifying the locations of promoters and enhancers in non-coding regions is challenging, we have fewer examples of identified causative cis-regulatory mutations that underlie naturally occurring phenotypic variations than of causative amino acid-altering mutations. Because cis-regulatory elements have epigenetic marks of specific histone modifications, we can detect cis-regulatory elements by mapping and analyzing them. Here, we investigated histone modifications and chromatin accessibility with cleavage under targets and tagmentation (CUT&Tag) and assay for transposase-accessible chromatin-sequencing (ATAC-seq). RESULTS: Using the threespine stickleback (Gasterosteus aculeatus) as a model, we confirmed that the genes for which nearby regions showed active marks, such as H3K4me1, H3K4me3, and high chromatin accessibility, were highly expressed. In contrast, the expression levels of genes for which nearby regions showed repressive marks, such as H3K27me3, were reduced, suggesting that our chromatin analysis protocols overall worked well. Genomic regions with peaks of histone modifications showed higher nucleotide diversity within and between populations. By comparing gene expression in the gills of the marine and stream ecotypes, we identified several insertions and deletions (indels) with transposable element fragments in the candidate cis-regulatory regions. CONCLUSIONS: Thus, mapping and analyzing histone modifications can help identify cis-regulatory elements and accelerate the identification of causative mutations in the non-coding regions underlying naturally occurring phenotypic variations.

Genomic Approaches Are Improving Taxonomic Representation in Genetic Studies of Speciation.

Until recently, our understanding of the genetics of speciation was limited to a narrow group of model species with a specific set of characteristics that made genetic analysis feasible. Rapidly advancing genomic technologies are eliminating many of the distinctions between laboratory and natural systems. In light of these genomic developments, we review the history of speciation genetics, advances that have been gleaned from model and non-model organisms, the current state of the field, and prospects for broadening the diversity of taxa included in future studies. Responses to a survey of speciation scientists across the world reveal the ongoing division between the types of questions that are addressed in model and non-model organisms. To bridge this gap, we suggest integrating genetic studies from model systems that can be reared in the laboratory or greenhouse with genomic studies in related non-models where extensive ecological knowledge exists.

Gudgeon fish with and without genetically determined countershading coexist in heterogeneous littoral environments of an ancient lake.

Countershading, characterized by a darker dorsal surface and lighter ventral surface, is common among many animals. This dorsoventral pigment polarity is often thought to be adaptive coloration for camouflage. By contrast, noncountershaded (melanistic) morphs often occur within a species due to genetic color polymorphism in terrestrial animals. However, the polymorphism with either countershaded or melanistic morphs is poorly known in wild aquatic animals. This study explored the genetic nature of diverged color morphs of a lineage of gudgeon fish (genus Sarcocheilichthys) in the ancient Lake Biwa and propose this system as a novel model for testing hypotheses of functional aspects of countershading and its loss in aquatic environments. This system harbors two color morphs that have been treated taxonomically as separate species; Sarcocheilichthys variegatus microoculus which occurs throughout the littoral zone and Sarcocheilichthys biwaensis which occurs in and around rocky areas. First, we confirmed that the divergence of dorsoventral color patterns between the two morphs is under strict genetic control at the levels of chromatophore distribution and melanin-related gene expression under common garden rearing. The former morph displayed sharp countershading coloration, whereas the latter morph exhibited a strong tendency toward its loss. The crossing results indicated that this divergence was likely controlled by a single locus in a two-allele Mendelian inheritance pattern. Furthermore, our population genomic and genome-wide association study analyses detected no genome-wide divergence between the two morphs, except for one region near a locus that may be associated with the color divergence. Thus, these morphs are either in a state of intraspecific color polymorphism or two incipient species. Evolutionary forces underlying this polymorphism appear to be associated with heterogeneous littoral environments in this lake. Future ecological genomic research will provide insight into adaptive functions of this widespread coloration, including the eco-evolutionary drivers of its loss, in the aquatic world.

Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates.

How early stages of speciation in free-spawning marine invertebrates proceed is poorly understood. The Western Pacific abalones, Haliotis discus, H. madaka, and H. gigantea, occur in sympatry with shared breeding season and are capable of producing viable F1 hybrids in spite of being ecologically differentiated. Population genomic analyses revealed that although the three species are genetically distinct, there is evidence for historical and ongoing gene flow among these species. Evidence from demographic modeling suggests that reproductive isolation among the three species started to build in allopatry and has proceeded with gene flow, possibly driven by ecological selection. We identified 27 differentiation islands between the closely related H. discus and H. madaka characterized by high FST and dA, but not high dXY values, as well as high genetic diversity in one H. madaka population. These genomic signatures suggest differentiation driven by recent ecological divergent selection in presence of gene flow outside of the genomic islands of differentiation. The differentiation islands showed low polymorphism in H. gigantea, and both high FST, dXY, and dA values between H. discus and H. gigantea, as well as between H. madaka and H. gigantea. Collectively, the Western Pacific abalones appear to occupy the early stages speciation continuum, and the differentiation islands associated with ecological divergence among the abalones do not appear to have acted as barrier loci to gene flow in the younger divergences but appear to do so in older divergences.

Multiple paths to the same destination: Influence of gene flow on convergent evolution.

Elucidation of the genetic mechanisms of convergent evolution, the evolution of similar or the same phenotypes in phylogenetically independent lineages, helps predict how populations will respond to the same selective pressures. Convergent evolution can be caused by either the fixation of identical-by-descent alleles, independent mutations at the same gene, or mutations in different genes controlling the same trait. To what extent does the fixation of identical-by-descent alleles lead to convergent evolution in isolated populations where inflow of adaptive alleles from other populations is limited? In a From the Cover article in this issue of Molecular Ecology, Kemppainen et al. (2021) compared the genetic basis for the reduction of pelvic structures in three isolated freshwater populations of nine-spined stickleback (Pungitius pungitius) from Northern Europe. The authors used quantitative trait loci (QTL) mapping to reveal that the pelvic reduction in these three populations was caused by mutations at different genetic loci. In contrast to studies in three-spined stickleback (Gasterosteus aculeatus), where independently derived Pitx1 mutations were shown to be responsible for plate reduction across multiple freshwater populations, Kemppainen et al. (2021) found Pitx1 to be the candidate causative gene for only one population of P. pungitius. This study highlights the importance of genetic studies of convergent evolution, not only in the presence of gene flow but also in its absence for a better understanding of the genetic architecture of convergent evolution.

Genome-wide patterns of divergence and introgression after secondary contact between Pungitius sticklebacks.

Speciation is a continuous process. Although it is known that differential adaptation can initiate divergence even in the face of gene flow, we know relatively little about the mechanisms driving complete reproductive isolation and the genomic patterns of divergence and introgression at the later stages of speciation. Sticklebacks contain many pairs of sympatric species differing in levels of reproductive isolation and divergence history. Nevertheless, most previous studies have focused on young species pairs. Here, we investigated two sympatric stickleback species, Pungitius pungitius and P. sinensis, whose habitats overlap in eastern Hokkaido; these species show hybrid male sterility, suggesting that they may be at a late stage of speciation. Our demographic analysis using whole-genome sequence data showed that these species split 1.73 Ma and came into secondary contact 37 200 years ago after a period of allopatry. This long period of allopatry might have promoted the evolution of intrinsic incompatibility. Although we detected on-going gene flow and signatures of introgression, overall genomic divergence was high, with considerable heterogeneity across the genome. The heterogeneity was significantly associated with variation in recombination rate. This sympatric pair provides new avenues to investigate the late stages of the stickleback speciation continuum. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Ecosystem size predicts the probability of speciation in migratory freshwater fish.

Predicting speciation is a fundamental goal of research in evolutionary ecology. The probability of speciation is often positively correlated with ecosystem size. Although the mechanisms driving this correlation are generally difficult to identify, a shared geographical and ecological context provides a suitable condition to study the mechanisms that promote speciation in large ecosystems by reducing the number of factors to be considered. Here, we determined the correlation between speciation and ecosystem size, and discuss the underlying mechanisms of this relationship, using a probable parallel ecotype formation for freshwater fish. Our population genetic analysis revealed that speciation of the landlocked goby, Rhinogobius sp. YB, of the Ryukyu Archipelago, Japan, from its migratory ancestor, R. brunneus, occurred in parallel across five islands. Logistic regression analysis showed that speciation probability could be predicted using island size. The results suggest that ecosystem size predicts the occurrence of adaptation and reproductive isolation, probably through its association with three possible factors: divergent selection strength, population persistence, and occurrence probability of habitat separation.

A key metabolic gene for recurrent freshwater colonization and radiation in fishes.

Colonization of new ecological niches has triggered large adaptive radiations. Although some lineages have made use of such opportunities, not all do so. The factors causing this variation among lineages are largely unknown. Here, we show that deficiency in docosahexaenoic acid (DHA), an essential ω-3 fatty acid, can constrain freshwater colonization by marine fishes. Our genomic analyses revealed multiple independent duplications of the fatty acid desaturase gene Fads2 in stickleback lineages that subsequently colonized and radiated in freshwater habitats, but not in close relatives that failed to colonize. Transgenic manipulation of Fads2 in marine stickleback increased their ability to synthesize DHA and survive on DHA-deficient diets. Multiple freshwater ray-finned fishes also show a convergent increase in Fads2 copies, indicating its key role in freshwater colonization.

Phylogeny, hybridization, and life history evolution of Rhinogobius gobies in Japan, inferred from multiple nuclear gene sequences.

Rhinogobius fishes (Gobiidae) are distributed widely in East and Southeast Asia, and represent the most species-rich group of freshwater gobies with diversified life histories (i.e., amphidromous, fluvial, and lentic). To reveal their phylogenetic relationships and life history evolution patterns, we sequenced six nuclear and three mitochondrial DNA (mtDNA) loci from 18 species, mainly from the mainland of Japan and the Ryukyu Archipelago. Our phylogenetic tree based on nuclear genes resolved three major clades, including several distinct subclades. The mtDNA and nuclear DNA phylogenies showed large discordance, which strongly suggested mitochondrial introgression through large-scale interspecific hybridization in these regions. On the basis of the molecular dating using geological data as calibration points, the hybridization occurred in the early to middle Pleistocene. Reconstruction of the ancestral states of life history traits based on nuclear DNA phylogeny suggests that the evolutionary change from amphidromous to freshwater life, accompanied by egg size change, occurred independently in at least three lineages. One of these lineages showed two life history alterations, i.e., from amphidromous (small egg) to fluvial (large egg) to lentic (small egg). Although more inclusive analysis using species outside Japan should be further conducted, the present results suggest the importance of the life history evolution associated with high adaptability to freshwater environments in the remarkable species diversification in this group. Such life history divergences may have contributed to the development of reproductive isolation.