Genomic architecture of habitat‐related divergence and signature of directional selection in the body shapes of Gnathopogon fishes

Evolution of ecomorphologically relevant traits such as body shapes is important to colonize and persist in a novel environment. Habitat‐related adaptive divergence of these traits is therefore common among animals. We studied the genomic architecture of habitat‐related divergence in the body shape of Gnathopogon fishes, a novel example of lake–stream ecomorphological divergence, and tested for the action of directional selection on body shape differentiation. Compared to stream‐dwelling Gnathopogon elongatus, the sister species Gnathopogon caerulescens, exclusively inhabiting a large ancient lake, had an elongated body, increased proportion of the caudal region and small head, which would be advantageous in the limnetic environment. Using an F₂ interspecific cross between the two Gnathopogon species (195 individuals), quantitative trait locus (QTL) analysis with geometric morphometric quantification of body shape and restriction‐site associated DNA sequencing‐derived markers (1622 loci) identified 26 significant QTLs associated with the interspecific differences of body shape‐related traits. These QTLs had small to moderate effects, supporting polygenic inheritance of the body shape‐related traits. Each QTL was mostly located on different genomic regions, while colocalized QTLs were detected for some ecomorphologically relevant traits that are proxy of body and caudal peduncle depths, suggesting different degree of modularity among traits. The directions of the body shape QTLs were mostly consistent with the interspecific difference, and QTL sign test suggested a genetic signature of directional selection in the body shape divergence. Thus, we successfully elucidated the genomic architecture underlying the adaptive changes of the quantitative and complex morphological trait in a novel system.     

Satellite DNA restriction site variability as a molecular basis of taxonoprint method: evidence from the study of Caucasian rock lizards

The restriction site distribution in satellite DNA of 17 Caucasian rock lizard species of the genus Lacerta (Darevskia gen. nov.), (Squamata, Lacertidae) was analyzed. The distribution patterns were shown to reflect the degree of satellite DNA evolutionary divergence, which could be revealed by taxonprint method, i.e., through the analysis of genomic DNA with a set of restriction endonucleases and subsequent computer-aided treatment. Thus, the taxonprint method offers an opportunity to examine the satellite DNA divergence in closely related species and infer their phylogeny of the species studied without reserting to costly and labor-consuming procedures. This is the advantage of using this technique at the early stages of genomic DNA phylogenetic analysis for rapid and effective estimation of relationships between closely related species as well as in the cases when DNA cloning and sequencing are too expensive or not feasible.     

Detection,Identification and Phylogenetic Analysis of Indian Ralstonia solanacearum Isolates by Comparison of Partial hrpB Gene Sequences

A species‐specific Polymerase Chain Reaction (sPCR) method was developed to identify and detect isolates of Ralstonia solanacearum, the cause of bacterial wilt disease in chilli. PCR primers for R. solanacearum were identified by alignment of hrpB gene sequences and selection of sequences specific for R. solanacearum at their 3′ ends. The primers were shown to be specific for R. solanacearum, as no PCR product was obtained when genomic DNA from other bacterial species including closely related Ralstonia species, were used as test species. Lone pair of primers (RshrpBF and RshrpBR) was designed using hrpB gene sequence, unique to R. solanacearum which amplified a predicted PCR product of 810 bp from 20 different isolates. Phylogenetic analysis was also attempted to understand the evolutionary divergence of Indian R. solanacearum isolates. Based on phylogenetic analysis, Indian isolates showed homology with the standard reference isolates from other countries but, interestingly, one new isolate showed complete evolutionary divergence by forming an out‐group.     

Patterns of divergence across the geographic and genomic landscape of a butterfly hybrid zone associated with a climatic gradient

Hybrid zones are a valuable tool for studying the process of speciation and for identifying the genomic regions undergoing divergence and the ecological (extrinsic) and nonecological (intrinsic) factors involved. Here, we explored the genomic and geographic landscape of divergence in a hybrid zone between Papilio glaucus and Papilio canadensis. Using a genome scan of 28,417 ddRAD SNPs, we identified genomic regions under possible selection and examined their distribution in the context of previously identified candidate genes for ecological adaptations. We showed that differentiation was genomewide, including multiple candidate genes for ecological adaptations, particularly those involved in seasonal adaptation and host plant detoxification. The Z chromosome and four autosomes showed a disproportionate amount of differentiation, suggesting genes on these chromosomes play a potential role in reproductive isolation. Cline analyses of significantly differentiated genomic SNPs, and of species‐diagnostic genetic markers, showed a high degree of geographic coincidence (81%) and concordance (80%) and were associated with the geographic distribution of a climate‐mediated developmental threshold (length of the growing season). A relatively large proportion (1.3%) of the outliers for divergent selection were not associated with candidate genes for ecological adaptations and may reflect the presence of previously unrecognized intrinsic barriers between these species. These results suggest that exogenous (climate‐mediated) and endogenous (unknown) clines may have become coupled and act together to reinforce reproductive isolation. This approach of assessing divergence across both the genomic and geographic landscape can provide insight about the interplay between the genetic architecture of reproductive isolation and endogenous and exogenous selection.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009 ; Meyer 2011 ; Nosil 2012 ). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009 ), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009 ), or pollinator‐mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012 ), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011 ). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006 ), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011 ). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. ( 2014 ) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early‐stage divergence in Midas cichlids.     

A genomic assessment of population structure and gene flow in an aquatic salamander identifies the roles of spatial scale,barriers, and river architecture

相似文献   

Genome scans and elusive candidate genes: detecting the variation that matters for speciation

Next‐generation sequencing is providing us with vast amounts of genetic data, yet we are currently struggling in our attempts to make sense of them. In particular, it has proven difficult to link phenotypic divergence and speciation to genome level divergence. In the current issue of Molecular Ecology, Ruegg et al. ( 2014 ) present new empirical results from two closely related bird taxa. They use a promising approach combining genome scan and candidate gene analysis. Their results suggest that we may have been looking in vain for candidate speciation genes by focusing only on genes found within genomic islands of divergence. This is because genes important in divergence and speciation may not be detected by genome scans and because features of the genomic architecture per se may have a large effect on the pattern of genome divergence.     

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal,poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab‐adapted ecotypes that exhibit partial reproductive isolation from wave‐adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co‐association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes—consistent with models of adaptation with gene flow—but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co‐association of loci during ecological speciation with ongoing gene flow.     

Genomics of speciation and introgression in Princess cichlid fishes from Lake Tanganyika

How variation in the genome translates into biological diversity and new species originate has endured as the mystery of mysteries in evolutionary biology. African cichlid fishes are prime model systems to address speciation‐related questions for their remarkable taxonomic and phenotypic diversity, and the possible role of gene flow in this process. Here, we capitalize on genome sequencing and phylogenomic analyses to address the relative impacts of incomplete lineage sorting, introgression and hybrid speciation in the Neolamprologus savoryi‐complex (the ‘Princess cichlids’) from Lake Tanganyika. We present a time‐calibrated species tree based on whole‐genome sequences and provide strong evidence for incomplete lineage sorting in the early phases of diversification and multiple introgression events affecting different stages. Importantly, we find that the Neolamprologus chromosomes show centre‐to‐periphery biases in nucleotide diversity, sequence divergence, GC content, incomplete lineage sorting and rates of introgression, which are likely modulated by recombination density and linked selection. The detection of heterogeneous genomic landscapes has strong implications on the genomic mechanisms involved in speciation. Collinear chromosomal regions can be protected from gene flow and harbour incompatibility genes if they reside in lowly recombining regions, and coupling can evolve between nonphysically linked genomic regions (chromosome centres in particular). Simultaneously, higher recombination towards chromosome peripheries makes these more dynamic, evolvable regions where adaptation polymorphisms have a fertile ground. Hence, differences in genome architecture could explain the levels of taxonomic and phenotypic diversity seen in taxa with collinear genomes and might have contributed to the spectacular cichlid diversity observed today.     

Strong reproductive isolation and narrow genomic tracts of differentiation among three woodpecker species in secondary contact

Hybrid zones allow the measurement of gene flow across the genome, producing insight into the genomic architecture of speciation. Such analysis is particularly powerful when applied to multiple pairs of hybridizing species, as patterns of genomic differentiation can then be related to age of the hybridizing species, providing a view into the build‐up of differentiation over time. We examined 33 809 single nucleotide polymorphisms (SNPs) in three hybridizing woodpecker species: Red‐breasted, Red‐naped and Yellow‐bellied sapsuckers (Sphyrapicus ruber, Sphyrapicus nuchalis and Sphyrapicus varius), two of which (ruber and nuchalis) are much more closely related than each is to the third (varius). To identify positions of SNPs on chromosomes, we developed a localization method based on comparative genomics. We found narrow clines, bimodal distributions of hybrid indices and genomic regions with decreased rates of introgression. These results suggest moderately strong reproductive isolation among species and selection against specific hybrid genotypes. We found 19 small regions of strong differentiation between species, partly shared among species pairs, but no large regions of differentiation. An association analysis revealed a single strong‐effect candidate locus associated with plumage, possibly explaining mismatch among the three species in genomic relatedness and plumage similarity. Our comparative analysis of species pairs of different age and their hybrid zones showed that moderately strong reproductive isolation can occur with little genomic differentiation, but that reproductive isolation is incomplete even with much greater genomic differentiation, implying there are long periods of time when hybridization is possible if diverging populations are in geographic contact.     

Variability of Restriction Sites in Satellite DNA as a Molecular Basis of Taxonprint Method: Evidence from the Study of Caucasian Rock Lizards

The distribution of restriction sites in satellite DNA of 17 Caucasian rock lizard species of the genus Lacerta (Darevskia gen. nov., (Squamata, Lacertidae) was analyzed. The distribution patterns were shown to reflect the degree of satellite DNA evolutionary divergence, which could be revealed by taxonprint method, i.e., through the analysis of genomic DNA with a set of restriction endonucleases and subsequent computer-aided analysis. Thus, the taxonprint method offers an opportunity to examine the satellite DNA divergence in closely related species and infer the phylogeny of the species studied without reserting to costly and labor-consuming procedures. This is the advantage of using this technique at the early stages phylogenetic analysis of genomic DNA for rapid and effective estimation of relationships between closely related species as well as in the cases when DNA cloning and sequencing are too expensive or not feasible.     

Genomic collinearity and the genetic architecture of floral differences between the homoploid hybrid species Iris nelsonii and one of its progenitors,Iris hexagona

Hybrid speciation represents a relatively rapid form of diversification. Early models of homoploid hybrid speciation suggested that reproductive isolation between the hybrid species and progenitors primarily resulted from karyotypic differences between the species. However, genic incompatibilities and ecological divergence may also be responsible for isolation. Iris nelsonii is an example of a homoploid hybrid species that is likely isolated from its progenitors primarily by strong prezygotic isolation, including habitat divergence, floral isolation and post-pollination prezygotic barriers. Here, we used linkage mapping and quantitative trait locus (QTL) mapping approaches to investigate genomic collinearity and the genetic architecture of floral differences between I. nelsonii and one of its progenitor species I. hexagona. The linkage map produced from this cross is highly collinear with another linkage map produced between I. fulva and I. brevicaulis (the two other species shown to have contributed to the genomic makeup of I. nelsonii), suggesting that karyotypic differences do not contribute substantially to isolation in this homoploid hybrid species. Similar to other studies of the genetic architecture of floral characteristics, at least one QTL was found that explained >20% variance in each color trait, while minor QTLs were detected for each morphological trait. These QTLs will serve as hypotheses for regions under selection by pollinators.     

Genome scans for the contemporary response to selection in quantitative traits

Genome scans have been an important approach for discovering historical signatures of selection in both model and nonmodel species. An exciting new experimental design for genome scans is to measure the change in allele frequency before and after contemporary selection within a generation, from a single population. The most widely‐used methods, however, have two major limitations: they are based on testing one locus at a time, and they only have power to uncover loci that have evolved under relatively strong selection. On the other hand, complex quantitative traits are common in nature and are caused by several loci of small effect. Selection on a quantitative trait at the phenotypic level is predicted to be accompanied by subtle allele frequency changes in many loci that covary (a polygenic soft sweep), rather than a large, single‐effect allele (a selective sweep). In this issue of Molecular Ecology, Bourret et al. (2014) measure the contemporary response to natural selection across the genome in multiple cohorts of Atlantic salmon during their first year at sea. They introduce a multilocus framework based on groups of markers that covary in their genotypic distribution. While the traditional, single‐locus approach did not find evidence for repeated patterns of selection, the multivariate approach found that a group of covarying SNPs was selected for in different cohorts at one site. Their multilocus framework has potential to be a more fruitful approach for uncovering the genomic basis of adaptation in quantitative traits, although caution should be applied as the framework has yet to be validated with simulated data.     

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits — including body shape, colour, lips and jaws — is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high‐bodied benthic species (Amphilophus astorquii). A total of 453 genome‐wide informative SNPs were identified in 240 F₂ hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy‐two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL‐linked genomic regions, genes that might contribute to divergence in body shape along the benthic–limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.     

Speciation as a sieve for ancestral polymorphism

Because they are considered rare, balanced polymorphisms are often discounted as crucial constituents of genome‐wide variation in sequence diversity. Despite its perceived rarity, however, long‐term balancing selection can elevate genetic diversity and significantly affect observed divergence between species. Here, we discuss how ancestral balanced polymorphisms can be “sieved” by the speciation process, which sorts them unequally across descendant lineages. After speciation, ancestral balancing selection is revealed by genomic regions of high divergence between species. This signature, which resembles that of other evolutionary processes, can potentially confound genomic studies of population divergence and inferences of “islands of speciation.”     

Cryptic lineages of a common alpine mayfly show strong life‐history divergence

Understanding ecological divergence of morphologically similar but genetically distinct species – previously considered as a single morphospecies – is of key importance in evolutionary ecology and conservation biology. Despite their morphological similarity, cryptic species may have evolved distinct adaptations. If such ecological divergence is unaccounted for, any predictions about their responses to environmental change and biodiversity loss may be biased. We used spatio‐temporally replicated field surveys of larval cohort structure and population genetic analyses (using nuclear microsatellite markers) to test for life‐history divergence between two cryptic lineages of the alpine mayfly Baetis alpinus in the Swiss Alps. We found that the more widespread and abundant cryptic lineage represents a ‘generalist’ with at least two cohorts per year, whereas the less abundant lineage is restricted to higher elevations and represents a ‘specialist’ with a single cohort per year. Importantly, our results indicate partial temporal segregation in reproductive periods between these lineages, potentially facilitating local coexistence and reproductive isolation. Taken together, our findings emphasize the need for a taxonomic revision: widespread and apparently generalist morphospecies can hide cryptic lineages with much narrower ecological niches and distribution ranges.