Identifying evolutionarily significant units for conservation of the endangered Malus sieversii using genome‐wide RADseq data

Malus sieversii, a wild progenitor of the domesticated apple, is an endangered species and is assigned second conservation priority by the China Plant Red Data Book. It is urgent to carry out in situ conservation of this species, but previous studies have not identified evolutionarily significant units (ESUs) for conservation management. In this study, we investigated the genetic diversity and relationships of six M. sieversii populations from China using integrated analysis of microsatellite (nSSR) data, genome‐wide SNPs and previous results in order to propose a reasonable conservation management. The results showed that levels of genetic diversity were inconsistently reflected by our nSSR and previous studies, suggesting that indices of genetic diversity are not effective to identify priority conservation areas for M. sieversii. Based on the selection criteria of ESUs for endangered species conservation, ESUs should reflect lineage divergence, geographical separation and different adaptive variation. Our phylogenetic tree based on genome‐wide SNPs yielded a clear relationship of divergent lineages among M. sieversii populations, leading to new different from those of previous studies. Three independent lineages, including the pairs of populations Huocheng‐Yining, Gongliu‐Xinyuan and Tuoli‐Emin, were identified. The geographic distances between populations among the different phylogenetic lineages were much greater than those within the same phylogenetic lineage. A cluster analysis on environmental variables showed that the three independent lineages inhabit different environmental conditions, suggesting that they may have adapted to different environments. Based on the results, we propose that three independent ESUs should be recognized as conservation units for M. sieversii in China.     

Fish conservation in the land of steppe and sky: Evolutionarily significant units of threatened salmonid species in Mongolia mirror major river basins

Mongolia's salmonids are suffering extensive population declines; thus, more comprehensive fisheries management and conservation strategies are required. To assist with their development, a better understanding of the genetic structure and diversity of these threatened species would allow a more targeted approach for preserving genetic variation and ultimately improve long‐term species recoveries. It is hypothesized that the unfragmented river basins that have persisted across Mongolia provide unobstructed connectivity for resident salmonid species. Thus, genetic structure is expected to be primarily segregated between major river basins. We tested this hypothesis by investigating the population structure for three salmonid genera (Hucho, Brachymystax and Thymallus) using different genetic markers to identify evolutionarily significant units (ESUs) and priority rivers to focus conservation efforts. Fish were assigned to separate ESUs when the combined evidence of mitochondrial and nuclear data indicated genetic isolation. Hucho taimen exhibited a dichotomous population structure forming two ESUs, with five priority rivers. Within the Brachymystax genus, there were three B. lenokESUs and one B. tumensisESU, along with six priority rivers. While B. tumensiswas confirmed to display divergent mtDNA haplotypes, haplotype sharing between these two congeneric species was also identified. For T. baicalensis,only a single ESU was assigned, with five priority rivers identified plus Lake Hovsgol. Additionally, we confirmed that T. nigrescens from Lake Hovsgol is a synonym of T. baicalensis. Across all species, the most prominent pattern was strong differentiation among major river basins with low differentiation and weak patterns of isolation by distance within river basins, which corroborated our hypothesis of high within‐basin connectivity across Mongolia. This new genetic information provides authorities the opportunity to distribute resources for management between ESUs while assigning additional protection for the more genetically valuable salmonid rivers so that the greatest adaptive potential within each species can be preserved.     

Genomewide scan for adaptive differentiation along altitudinal gradient in the Andrew's toad Bufo andrewsi

Recent studies of humans, dogs and rodents have started to discover the genetic underpinnings of high altitude adaptations, yet amphibians have received little attention in this respect. To identify possible signatures of adaptation to altitude, we performed a genome scan of 15 557 single nucleotide polymorphisms (SNPs) obtained with restriction site‐associated DNA sequencing of pooled samples from 11 populations of Andrew's toad (Bufo andrewsi) from the edge of the Tibetan Plateau, spanning an altitudinal gradient from 1690 to 2768 m.a.s.l. We discovered significant geographic differentiation among all sites, with an average F_ST = 0.023 across all SNPs. Apart from clear patterns of isolation by distance, we discovered numerous outlier SNPs showing strong associations with variation in altitude (1394 SNPs), average annual temperature (1859 SNPs) or both (1051 SNPs). Levels and patterns of genetic differentiation in these SNPs were consistent with the hypothesis that they have been subject to directional selection and reflect adaptation to altitudinal variation among the study sites. Genes with footprints of selection were significantly enriched in binding and metabolic processes. Several genes potentially related to high altitude adaptation were identified, although the identity and functional significance of most genomic targets of selection remain unknown. In general, the results provide genomic support for results of earlier common garden and low coverage genetic studies that have uncovered substantial adaptive differentiation along altitudinal and latitudinal gradients in amphibians.     

Population genomics of wild and laboratory zebrafish (Danio rerio)

Understanding a wider range of genotype–phenotype associations can be achieved through ecological and evolutionary studies of traditional laboratory models. Here, we conducted the first large‐scale geographic analysis of genetic variation within and among wild zebrafish (Danio rerio) populations occurring in Nepal, India, and Bangladesh, and we genetically compared wild populations to several commonly used lab strains. We examined genetic variation at 1832 polymorphic EST‐based single nucleotide polymorphisms (SNPs) and the cytb mitochondrial gene in 13 wild populations and three lab strains. Natural populations were subdivided into three major mitochondrial DNA clades with an average among‐clade sequence divergence of 5.8%. SNPs revealed five major evolutionarily and genetically distinct groups with an overall F_ST of 0.170 (95% CI 0.105–0.254). These genetic groups corresponded to discrete geographic regions and appear to reflect isolation in refugia during past climate cycles. We detected 71 significantly divergent outlier loci (3.4%) and nine loci (0.5%) with significantly low F_ST values. Valleys of reduced heterozygosity, consistent with selective sweeps, surrounded six of the 71 outliers (8.5%). The lab strains formed two additional groups that were genetically distinct from all wild populations. An additional subset of outlier loci was consistent with domestication selection within lab strains. Substantial genetic variation that exists in zebrafish as a whole is missing from lab strains that we analysed. A combination of laboratory and field studies that incorporates genetic variation from divergent wild populations along with the wealth of molecular information available for this model organism provides an opportunity to advance our understanding of genetic influences on phenotypic variation for a vertebrate species.     

Population assignment reveals low migratory connectivity in a weakly structured songbird

Understanding migratory connectivity is essential for determining the drivers behind population dynamics and for implementing effective conservation strategies for migratory species. Genetic markers provide a means to describe migratory connectivity; however, they can be uninformative for species with weak population genetic structure, which has limited their application. Here, we demonstrated a genomic approach to describing migratory connectivity in the prothonotary warbler, Protonotaria citrea, a Neotropical songbird of conservation concern. Using 26,189 single nucleotide polymorphisms (SNPs), we revealed regional genetic structure between the Mississippi River Valley and the Atlantic Seaboard with overall weak genetic differentiation among populations (F_ST = 0.0055; 95% CI: 0.0051–0.0059). Genetic variation had a stronger association with geographic rather than environmental factors, with each explaining 14.5% and 8.2% of genetic variation, respectively. By varying the numbers of genomic markers used in population assignment models with individuals of known provenance, we identified a maximum assignment accuracy (89.7% to site, 94.3% to region) using a subset of 600 highly differentiated SNPs. We then assigned samples from nonbreeding sites to breeding region and found low migratory connectivity. Our results highlight the importance of filtering markers for informative loci in models of population assignment. Quantifying migratory connectivity for weakly structured species will be useful for expanding studies to a wider range of migratory species across taxonomic groups and may contribute to a deeper understanding of the evolution of migratory strategies.     

Cunningham's skinks show low genetic connectivity and signatures of divergent selection across its distribution

Establishing corridors of connecting habitat has become a mainstay conservation strategy to maintain gene flow and facilitate climate‐driven range shifts. Yet, little attention has been given to ascertaining the extent to which corridors will benefit philopatric species, which might exhibit localized adaptation. Measures of genetic connectivity and adaptive genetic variation across species’ ranges can help fill this knowledge gap. Here, we characterized the spatial genetic structure of Cunningham's skink (Egernia cunninghami), a philopatric species distributed along Australia's Great Dividing Range, and assessed evidence of localized adaptation. Analysis of 4,274 SNPs from 94 individuals sampled at four localities spanning 500 km and 4° of latitude revealed strong genetic structuring at neutral loci (mean F_ST ± SD = 0.603 ± 0.237) among the localities. Putatively neutral SNPs and those under divergent selection yielded contrasting spatial patterns, with the latter identifying two genetically distinct clusters. Given low genetic connectivity of the four localities, we suggest that the natural movement rate of this species is insufficient to keep pace with spatial shifts to its climate envelope, irrespective of habitat availability. In addition, our finding of localized adaptation highlights the risk of outbreeding depression should the translocation of individuals be adopted as a conservation management strategy.     

Mitochondria,sex and variation in routine metabolic rate

Variation in the metabolic costs associated with organismal maintenance may play a key role in determining fitness, and thus these differences among individuals are likely to be subject to natural selection. Although the evolvability of maintenance metabolism depends on its underlying genetic architecture, relatively little is known about the nature of genetic variation that underlies this trait. To address this, we measured variation in routine metabolic rate (?O_2routine), an index of maintenance metabolism, within and among three populations of Atlantic killifish, Fundulus heteroclitus, including a population from a region of genetic admixture between two subspecies. Polygenic association tests among individuals from the admixed population identified 54 single nucleotide polymorphisms (SNPs) that were associated with ?O_2routine, and these SNPs accounted for 43% of interindividual variation in this trait. However, genetic associations with ?O_2routine involved different SNPs if females and males were analysed separately, and there was a sex‐dependent effect of mitochondrial genotype on variation in routine metabolism. These results imply that there are sex‐specific genetic mechanisms, and potential mitonuclear interactions, that underlie variation in ?O_2routine. Additionally, there was evidence for epistatic interactions between 17% of the possible pairs of trait‐associated SNPs, suggesting that epistatic effects on ?O_2routine are common. These data demonstrate not only that phenotypic variation in this ecologically important trait has a polygenic basis with considerable epistasis among loci, but also that these underlying genetic mechanisms, and particularly the role of mitochondrial genotype, may be sex‐specific.     

Trouble at the top? Restricted distribution and extreme population isolation in an alpine crustacean assemblage with unexpected lineage diversity

相似文献   

Genetic diagnosis of a rare myrmecochorous species,Plagiorhegma dubium (Berberidaceae): Historical genetic bottlenecks and strong spatial structures among populations

Distribution of genetic variation over time and space is relevant to demographic histories and tightly linked to ecological disturbances as well as evolutionary potential of an organism. Therefore, understanding the pattern of genetic diversity is a primary step in conservation and management projects for rare and threatened plant species. We used eight microsatellite markers to examine the level of genetic diversity, spatial structure, and demographic history of Plagiorhegma dubium, a rare myrmecochorous herb, populations sampled across northeast Asia and Siberia. We found low within‐population genetic variation associated with historical bottlenecks. Although pairwise F_ST values were not much higher than the ones found in similar life form species, STRUCTURE and PCoA revealed a clear broadscale spatial pattern of genetic structure. Bayesian clustering (best K = 6) and PCoA identified three populations that are distinctive from neighboring populations in the Korean peninsula, which suggests potential units for conservation and management plans in Korea. MIGRATE‐N and BAYESASS showed that both contemporary (0.003–0.045) and historical migration rates (2 × e^?5?4.6 × e^?4) were low. Our findings provide a good example, where genetic considerations should be integrated for conservation and management plans of rare and threatened species.     

Phylogeography,Population Structure,and Conservation of the Javan Gibbon (<Emphasis Type="Italic">Hylobates moloch</Emphasis>)

An understanding of phylogeography and population genetics is needed for a comprehensive long-term conservation management strategy. The Javan gibbon (Hylobates moloch), an Endangered species endemic to the island of Java, has been protected since 1924 but is threatened by ongoing habitat loss, habitat degradation, and the wildlife trade. We studied the phylogeography and population genetic structure of the Javan gibbon, to define the number of Evolutionary Significant Units (ESUs) in the species, and the population genetic structure in each ESU. We sampled 47 individuals, analyzing 35 for variation in mitochondrial DNA control region, 41 for variation in 8 nuclear DNA microsatellites, and 13 for variation in 45 nuclear DNA single nucleotide polymorphisms (SNPs). We found support for two ESUs across the species range: a western ESU, extending from Ujung Kulon to Gunung Gede–Pangrango, and a central ESU, extending from Gunung Masigit–Simpang–Tilu to Gunung Slamet. Analysis of molecular variance and population structure analysis indicate significant structuring in the western ESU between Ujung Kulon and Gunung Halimun–Salak–Gede–Pangrango, and little to moderate structure in the central ESU, underscoring the importance of conserving as many populations as possible to preserve the full array of genetic diversity in this species. Our results will inform future more comprehensive population genetic surveys and the conservation genetic management of the Javan gibbon. This study demonstrates the importance of genetics when designing conservation management strategies for endangered primates.     

Population structure and genetic diversity of trout (<Emphasis Type="Italic">Oncorhynchus mykiss</Emphasis>) above and below natural and man-made barriers in the Russian River,California

The effects of landscape features on gene flow in threatened and endangered species play an important role in influencing the genetic structure of populations. We examined genetic variation of trout in the species Oncorhynchus mykiss at 22 microsatellite loci from 20 sites in the Russian River basin in central California. We assessed relative patterns of genetic structure and variation in fish from above and below both natural (waterfalls) and man-made (dams) barriers. Additionally, we compared sites sampled in the Russian River with sites from 16 other coastal watersheds in California. Genetic variation among the 20 sites sampled within the Russian River was significantly partitioned into six groups above natural barriers and one group consisting of all below barrier and above dam sites. Although the below-barrier sites showed moderate gene flow, we found some support for sub-population differentiation of individual tributaries in the watershed. Genetic variation at all below-barrier sites was high compared to above-barrier sites. Fish above dams were similar to those from below-barrier sites and had similar levels of genetic diversity, indicating they have not been isolated very long from below-barrier populations. Population samples from above natural barriers were highly divergent, with large F _st values, and had significantly lower genetic diversity, indicating relatively small population sizes. The origins of populations above natural barriers could not be ascertained by comparing microsatellite diversity to other California rivers. Finally, below-barrier sites farther inland were more genetically differentiated from other watersheds than below-barrier sites nearer the river’s mouth.     

Natural selection contributes to geographic patterns of thermal plasticity in Plantago lanceolata

A long‐standing debate in evolutionary biology concerns the relative importance of different evolutionary forces in explaining phenotypic diversification at large geographic scales. For example, natural selection is typically assumed to underlie divergence along environmental gradients. However, neutral evolutionary processes can produce similar patterns. We collected molecular genetic data from 14 European populations of Plantago lanceolata to test the contributions of natural selection versus neutral evolution to population divergence in temperature‐sensitive phenotypic plasticity of floral reflectance. In P. lanceolata, reflectance plasticity is positively correlated with latitude/altitude. We used population pairwise comparisons between neutral genetic differentiation (F_ST and Jost's D) and phenotypic differentiation (P_ST) to assess the contributions of geographic distance and environmental parameters of the reproductive season in driving population divergence. Data are consistent with selection having shaped large‐scale geographic patterns in thermal plasticity. The aggregate pattern of P_ST versus F_ST was consistent with divergent selection. F_ST explained thermal plasticity differences only when geographic distance was not included in the model. Differences in the extent of cool reproductive season temperatures, and not overall temperature variation, explained plasticity differences independent of distance. Results are consistent with the hypothesis that thermal plasticity is adaptive where growing seasons are shorter and cooler, that is, at high latitude/altitude.     

Microsatellite variation reveals weak genetic structure and retention of genetic variability in threatened Chinook salmon (<Emphasis Type="Italic">Oncorhynchus tshawytscha</Emphasis>) within a Snake River watershed

Pacific salmon (Oncorhynchus spp.) have been central to the development of management concepts associated with evolutionarily significant units (ESUs), yet there are still relatively few studies of genetic diversity within threatened and endangered ESUs for salmon or other species. We analyzed genetic variation at 10 microsatellite loci to evaluate spatial population structure and genetic variability in indigenous Chinook salmon (Oncorhynchus tshawytscha) across a large wilderness basin within a Snake River ESU. Despite dramatic 20th century declines in abundance, these populations retained robust levels of genetic variability. No significant genetic bottlenecks were found, although the bottleneck metric (M ratio) was significantly correlated with average population size and variability. Weak but significant genetic structure existed among tributaries despite evidence of high levels of gene flow, with the strongest genetic differentiation mirroring the physical segregation of fish from two sub-basins. Despite the more recent colonization of one sub-basin and differences between sub-basins in the natural level of fragmentation, gene diversity and genetic differentiation were similar between sub-basins. Various factors, such as the (unknown) genetic contribution of precocial males, genetic compensation, lack of hatchery influence, and high levels of current gene flow may have contributed to the persistence of genetic variability in this system in spite of historical declines. This unique study of indigenous Chinook salmon underscores the importance of maintaining natural populations in interconnected and complex habitats to minimize losses of genetic diversity within ESUs.     

Effective population size is strongly correlated with breeding pond size in the endangered California tiger salamander, <Emphasis Type="Italic">Ambystoma californiense</Emphasis>

Maintaining genetic diversity and population viability in endangered and threatened species is a primary concern of conservation biology. Genetic diversity depends on population connectivity and effective population size (N_e), both of which are often compromised in endangered taxa. While the importance of population connectivity and gene flow has been well studied, investigating effective population sizes in natural systems has received far less attention. However, N_e plays a prominent role in the maintenance of genetic diversity, the prevention of inbreeding depression, and in determining the probability of population persistence. In this study, we examined the relationship between breeding pond characteristics and N_e in the endangered California tiger salamander, Ambystoma californiense. We sampled 203 individuals from 10 breeding ponds on a local landscape, and used 11 polymorphic microsatellite loci to quantify genetic structure, gene flow, and effective population sizes. We also measured the areas of each pond using satellite imagery and classified ponds as either hydrologically-modified perennial ponds or naturally occurring vernal pools, the latter of which constitute the natural breeding habitat for A. californiense. We found no correlation between pond area and heterozygosity or allelic diversity, but we identified a strong positive relationship between breeding pond area and N_e, particularly for vernal pools. Our results provide some of the first empirical evidence that variation in breeding habitat can be associated with differences in N_e and suggest that a more complete understanding of the environmental features that influence N_e is an important component of conservation genetics and management.     

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra‐ and interspecific genetic variation at 25 climate‐related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two‐thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. F_ST outlier analysis identified six SNPs in five climate‐related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation‐with‐migration model indicated that gene flow between the species has been more restricted for climate‐related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species‐specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.     

Characterization of microsatellite markers in the threatened sand skink (Neoseps reynoldsi)

We characterized microsatellite loci for the sand skink (Neoseps reynoldsi) for future studies of genetic structure in this threatened taxon. We screened a partial genomic library enriched for microsatellites, designed primers for eight loci and assessed these markers for polymorphism across 11 populations in central Florida. Preliminary analyses indicate deviations from Hardy–Weinberg expectations for most loci, suggesting population genetic structure across the sampled populations; therefore, understanding genetic connectivity is critical for maintaining genetic variation in this species.     

Evolutionarily significant units of the critically endangered leaf frog Pithecopus ayeaye (Anura,Phyllomedusidae) are not effectively preserved by the Brazilian protected areas network

Protected areas (PAs) are essential for biodiversity conservation, but their coverage is considered inefficient for the preservation of all species. Many species are subdivided into evolutionarily significant units (ESUs) and the effectiveness of PAs in protecting them needs to be investigated. We evaluated the usefulness of the Brazilian PAs network in protecting ESUs of the critically endangered Pithecopus ayeaye through ongoing climate change. This species occurs in a threatened mountaintop ecosystem known as campos rupestres. We used multilocus DNA sequences to delimit geographic clusters, which were further validated as ESUs with a coalescent approach. Ecological niche modeling was used to estimate spatial changes in ESUs’ potential distributions, and a gap analysis was carried out to evaluate the effectiveness of the Brazilian PAs network to protect P. ayeaye in the face of climate changes. We tested the niche overlap between ESUs to gain insights for potential management alternatives for the species. Pithecopus ayeaye contains at least three ESUs isolated in distinct mountain regions, and one of them is not protected by any PA. There are no climatic niche differences between the units, and only 4% of the suitable potential area of the species is protected in present and future projections. The current PAs are not effective in preserving the intraspecific diversity of P. ayeaye in its present and future range distributions. The genetic structure of P. ayeaye could represent a typical pattern in campos rupestres endemics, which should be considered for evaluating its conservation status.     

The impact of Yangtze River discharge on the genetic structure of a population of the rock bream,Oplegnathus fasciatus

As a rocky reef species, Oplegnathus fasciatus displays significant geographic differences that may be congruent with patterns of molecular variation. We used 10 polymorphism microsatellite loci to characterize the genetic structure and test the impact of ocean currents and the Yangtze River on the population genetic connectivity of the rock bream, O. fasciatus, sampled from Jiaonan (JN), Zhoushan (ZS) and Taiwan (TW). The average number of alleles was 17.78, with a total of 160 alleles, and the average allele value in the East China Sea (ZS?=?13.70) was higher than that in the Yellow Sea (JN?=?6.70). The individuals at the ZS location showed the highest expected heterozygosity (He?=?0.86), Nei's genetic diversity (H?=?0.84) and Shannon genetic diversity (I?=?2.21). Genetic bottleneck tests showed that individuals had experienced a recent bottleneck effect in the JN location. AMOVA analysis showed that 12.19% of genetic variation existed among populations, indicating significant population genetic differentiation in the species. The values of the population pairwise differentiation index (F_st and Φ_pt) ranged from 0.01 to 0.29. The UPGMA tree showed two diverged genealogical branches corresponding to sampling localities (Clade A?=?JN, Clade B?=?ZS and TW) in the species. The AMOVA and STRUCTURE analyses suggested that O. fasciatus may be composed of two stocks in the examined range, and the freshwater discharge from the Yangtze River might act as a physical barrier to genetic exchange between north and south China stocks of the species.     

Cryptic speciation and phylogeographic relationships in the elephant ear sponge Ianthella basta (Porifera,Ianthellidae) from northern Australia

Morphological delineation of sponge species is hindered by the narrow range of fixed diagnostic characters and our limited knowledge of how much phenotypic plasticity the sponge body plan assumes in response to environmental conditions. Here, we make use of the partial mitochondrial cytochrome c oxidase subunit I (COI) gene and the second internal transcribed spacer (ITS2) region to assess the taxonomic validity of colour morphotypes observed in the elephant ear sponge Ianthella basta (Pallas, 1776) across its distribution range in northern Australia, and explore levels and patterns of genetic diversity among populations of the species collected from both sides of the Torres Strait. Molecular phylogenetic analyses revealed congruent topologies consistent with three evolutionarily significant units (ESUs) that were independent of the morphology of the sponge. ESU I includes previously morphologically and genetically delineated western Pacific specimens of I. basta (Guam), and probably corresponds to the type specimen originally described from Indonesia. ESU I occurs in almost all sampling sites across northern Australia, suggesting considerable levels of connectivity among reefs throughout the Torres Strait. ESUs II and III are each exclusively associated with a geographic region of high sponge species richness separated by Torres Strait, and probably represent the result of historical population fragmentation. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 225–235.     

The population genomic signature of environmental association and gene flow in an ecologically divergent tree species Metrosideros polymorpha (Myrtaceae)

Genomewide markers enable us to study genetic differentiation within a species and the factors underlying it at a much higher resolution than before, which advances our understanding of adaptation in organisms. We investigated genomic divergence in Metrosideros polymorpha, a woody species that occupies a wide range of ecological habitats across the Hawaiian Islands and shows remarkable phenotypic variation. Using 1659 single nucleotide polymorphism (SNP) markers annotated with the genome assembly, we examined the population genetic structure and demographic history of nine populations across five elevations and two ages of substrates on Mauna Loa, the island of Hawaii. The nine populations were differentiated into two genetic clusters distributed on the lower and higher elevations and were largely admixed on the middle elevation. Demographic modelling revealed that the two genetic clusters have been maintained in the face of gene flow, and the effective population size of the high‐altitude cluster was much smaller. A F_ST‐based outlier search among the 1659 SNPs revealed that 34 SNPs (2.05%) were likely to be under divergent selection and the allele frequencies of 21 of them were associated with environmental changes along elevations, such as temperature and precipitation. This study shows a genomic mosaic of M. polymorpha, in which contrasting divergence patterns were found. While most genomic polymorphisms were shared among populations, a small fraction of the genome was significantly differentiated between populations in diverse environments and could be responsible for the dramatic adaptation to a wide range of environments.