Cryptic introgression into the kidney saxifrage (Saxifraga hirsuta) from its more abundant sympatric congener Saxifraga spathularis,and the potential risk of genetic assimilation

Background and Aims Although hybridization can play a positive role in plant evolution, it has been shown that excessive unidirectional hybridization can result in replacement of a species’ gene pool, and even the extinction of rare species via genetic assimilation. This study examines levels of introgression between the common Saxifraga spathularis and its rarer congener S. hirsuta, which have been observed to hybridize in the wild where they occur sympatrically.Methods Seven species-specific single nucleotide polymorphisms (SNPs) were analysed in 1025 plants representing both species and their hybrid, S. × polita, from 29 sites across their ranges in Ireland. In addition, species distribution modelling was carried out to determine whether the relative abundance of the two parental species is likely to change under future climate scenarios.Key Results Saxifraga spathularis individuals tended to be genetically pure, exhibiting little or no introgression from S. hirsuta, but significant levels of introgression of S. spathularis alleles into S. hirsuta were observed, indicating that populations exhibiting S. hirsuta morphology are more like a hybrid swarm, consisting of backcrosses and F₂s. Populations of the hybrid, S. × polita, were generally comprised of F₁s or F₂s, with some evidence of backcrossing. Species distribution modelling under projected future climate scenarios indicated an increase in suitable habitats for both parental species.Conclusions Levels of introgression observed in this study in both S. spathularis and S. hirsuta would appear to be correlated with the relative abundance of the species. Significant introgression of S. spathularis alleles was detected in the majority of the S. hirsuta populations analysed and, consequently, ongoing introgression would appear to represent a threat to the genetic integrity of S. hirsuta, particularly in areas where the species exists sympatrically with its congener and where it is greatly outnumbered.     

Evolutionary insights into the North American Necturus beyeri complex (Amphibia: Caudata) based on molecular genetic and morphological analyses

Necturus beyeri (Caudata: Proteidae), as conceived by some, contains paedomorphic salamanders distributed from the Ochlockonee drainage of Florida to the Angelina drainage of Texas. Because these salamanders differ in color pattern and karyotype across their geographic range, we performed a phylogeographic analysis that included representatives from all major drainages as well as of all congeners. The mitochondrially encoded ND2 gene was used to infer phylogenetic relationships using Bayesian inference. Morphometrics of head shape were analyzed and included as an independent data set. Our work suggests that Necturus comprises 11 lineages. A basal split within the genus separates an ancestor of two Atlantic Coastal Plain species (Necturus lewisi and Necturus punctatus) from the ancestor of nine distinct Gulf Coastal Plain lineages. One lineage is consistent with Necturus alabamensis, a species currently recognized in the Black Warrior drainage of Alabama. Two lineages comprise Necturus maculosus, as historically recognized, and six lineages comprise N. beyeri, as recognized by some, each of which occupies a unique drainage. Both of these species are demonstrated to be paraphyletic. Head morphometrics show the same patterns as the mtDNA. Overall, lineages within Necturus exhibit an east‐to‐west progression of appearance on the phylogenetic trees. This pattern corroborates biogeographic hypotheses based on previous karyological work. Within N. beyeri, this progression separates a pattern class of two eastern lineages lacking bold spotting and possessing relatively small mean body lengths from a pattern class of four western lineages possessing bold spotting and larger mean body sizes. Thus, the two eastern lineages of N. beyeri are similar in color pattern and body size to N. punctatus either through retention of the ancestral color pattern and size for the genus or through convergent selection in eastern streams of the Gulf Coastal Plain.     

Using genetic evaluation to guide conservation of remnant Juniperus communis (Cupressaceae) populations

• Many critically endangered plant species exist in small, genetically depauperate or inbred populations, making assisted gene flow interventions necessary for long‐term population viability. However, before such interventions are implemented, conservation practitioners must consider the genetic and demographic status of extant populations, which are strongly affected by species’ life‐history traits. In northwestern Europe, Juniperus communis, a dioecious, wind‐pollinated and bird‐dispersed gymnosperm, has been declining for the past century and largely exists in small, isolated and senescent populations.
• To provide useful recommendations for a recovery plan involving translocation of plants, we investigated genetic diversity and structure of populations in Belgium using four microsatellite and five plastid single‐nucleotide polymorphism (SNP) markers.
• We detected no clonality in the populations, suggesting predominantly sexual reproduction. Populations exhibited high genetic diversity (H_e = 0.367–0.563) and low to moderate genetic differentiation (F_ST ≤ 0.133), with no clear geographic structure. Highly positive inbreeding coefficients (F_IS = 0.221–0.507) were explained by null alleles, population substructuring and biparental inbreeding. No isolation by distance was observed among distant populations, but isolation at close geographic proximity was found. Patterns were consistent with high historical gene flow through pollen and seed dispersal at both short and long distances. We also tested four pre‐germination treatments among populations to improve germination rates; however, germination rates remained low and only cold‐stratification treatments induced germination in some populations.
• To bolster population regeneration, introductions of cuttings from several source populations are recommended, in combination with in situ management practices that improve seedling survival and with ex situ propagation.

     

Habitat of the endangered salt marsh harvest mouse (Reithrodontomys raviventris) in San Francisco Bay

Understanding habitat associations is vital for conservation of at‐risk marsh‐endemic wildlife species, particularly those under threat from sea level rise. We modeled environmental and habitat associations of the marsh‐endemic, Federally endangered salt marsh harvest mouse (Reithrodontomys raviventris, RERA) and co‐occurrence with eight associated small mammal species from annual trap data, 1998–2014, in six estuarine marshes in North San Francisco Bay, California. Covariates included microhabitat metrics of elevation and vegetation species and cover; and landscape metrics of latitude–longitude, distance to anthropogenic features, and habitat patch size. The dominant cover was pickleweed (Salicornia pacifica) with 86% mean cover and 37 cm mean height, and bare ground with about 10% mean cover. We tested 38 variants of Bayesian network (BN) models to determine covariates that best account for presence of RERA and of all nine small mammal species. Best models had lowest complexity and highest classification accuracy. Among RERA presence models, three best BN models used covariates of latitude–longitude, distance to paved roads, and habitat patch size, with 0% error of false presence, 20% error of false nonpresence, and 20% overall error. The all‐species presence models suggested that within the pickleweed marsh environment, RERA are mostly habitat generalists. Accounting for presence of other species did not improve prediction of RERA. Habitat attributes compared between RERA and the next most frequently captured species, California vole (Microtus californicus), suggested substantial habitat overlap, with RERA habitat being somewhat higher in marsh elevation, greater in percent cover of the dominant plant species, closer to urban areas, further from agricultural areas, and, perhaps most significant, larger in continuous size of marsh patch. Findings will inform conservation management of the marsh environment for RERA by identifying best microhabitat elements, landscape attributes, and adverse interspecific interactions.     

A landscape genetics approach reveals ecological‐based differentiation in populations of holm oak (Quercus ilex L.) at the northern limit of its range

The holm oak plays a relevant role in the functioning of Mediterranean forests. In the area north of Garda Lake, Italian Prealps, holm oak populations are at the northernmost edge of their distribution. Being peripheral, these populations are of particular interest for ecological, evolutionary and conservation studies. Through an explicit individual‐based landscape genetics approach, we addressed the following questions: (1) are levels of genetic variation reduced in these marginal populations compared with central populations?; (2) despite the narrow geographical scale, do individual‐based analyses have some power to detect genetic differentiation?; (3) do environmental and/or climatic factors exert a role in shaping patterns of genetic variation and differentiation? Through a Bayesian method, we identified three clusters whose genetic variability can be considered to be of the same order as that recorded in central Quercus ilex populations. Although being geographically very close (< 20 km), the differentiation was statistically significant (P < 0.05) with global F _st and Φ _Pt values of 0.019 and 0.038, respectively. Geography and phylogeography could not be invoked to explain this differentiation. A redundancy discriminant analysis revealed that relevant eco‐pedological and climatic features, such as soil depth, aspect, elevation and humidity, were correlated with the observed pattern of differentiation. Toblino was ecologically separated from the other clusters, as it lies on deep soil with subhumid conditions. The differentiation of the Brione–Ranzo–Val Busa cluster appeared to be related to superficial soils and drier conditions, whereas the Nanzone–Padaro cluster was differentiated mainly according to its mid‐elevation. Coupling spatial and genetic information on a local scale proved to be effective to investigate the evolutionary and demographic history of peripheral populations. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Balancing genetic uniqueness and genetic variation in determining conservation and translocation strategies: a comprehensive case study of threatened dwarf galaxias,Galaxiella pusilla (Mack) (Pisces: Galaxiidae)

Genetic markers are widely used to define and manage populations of threatened species based on the notion that populations with unique lineages of mtDNA and well‐differentiated nuclear marker frequencies should be treated separately. However, a danger of this approach is that genetic uniqueness might be emphasized at the cost of genetic diversity, which is essential for adaptation and is potentially boosted by mixing geographically separate populations. Here, we re‐explore the issue of defining management units, focussing on a detailed study of Galaxiella pusilla, a small freshwater fish of national conservation significance in Australia. Using a combination of microsatellite and mitochondrial markers, 51 populations across the species range were surveyed for genetic structure and diversity. We found an inverse relationship between genetic differentiation and genetic diversity, highlighting a long‐term risk of deliberate isolation of G. pusilla populations based on protection of unique lineages. Instead, we adopt a method for identifying genetic management units that takes into consideration both uniqueness and genetic variation. This produced a management framework to guide future translocation and re‐introduction efforts for G. pusilla, which contrasted to the framework based on a more traditional approach that may overlook important genetic variation in populations.     

The molecular phylogeny of the sea star Echinaster (Asteroidea: Echinasteridae) provides insights for genus taxonomy

In this study, we used sequences of two mitochondrial genes, cytochrome c oxidase I (COI) and 16S rRNA, and one nuclear gene, 28S rRNA, to test the monophyly of the sea star genus Echinaster, and understand the phylogenetic relationships among species and subgenera within this genus. Phylogenetic analyses based on Bayesian inference and maximum likelihood methods revealed three clades with high values of genetic divergence among them (K2P distances for COI over 23%). One of the clades grouped all Echinaster (Othilia) species, and the other two clades included Echinaster (non‐Othilia) species and Henricia species, respectively. Although the relationships among Henricia, Othilia, and Echinaster could not be completely clarified, the Othilia clade was a well‐supported group with shared diagnostic morphological characters. Moreover, the approximately unbiased test applied to the phylogenetic reconstruction rejected the hypothesis of the genus Echinaster as a monophyletic group. According to these results, we suggest the revalidation of Othilia as a genus instead of a subgenus within Echinaster. Our study clarifies important points about the phylogenetic relationships among species of Echinaster. Other important systematic questions about the taxonomic classification of Echinaster and Henricia still remain open, but this molecular study provides bases for future research on the topic.     

Molecular analysis of hybridisation between wild and domestic cats (Felis silvestris) in Portugal: implications for conservation

The endangered European wildcat (Felis silvestris silvestris) is represented, today, by fragmented and declining populations whose genetic integrity is considered to be seriously threatened by crossbreeding with widespread free-ranging domestic cats. Extensive and recent hybridisation has been described in Hungary and Scotland, in contrast with rare introgression of domestic alleles in Italy and Germany. In Portugal, the wildcat is now listed as VULNERABLE in the Red Book of Portuguese Vertebrates. Nevertheless, genetic diversity of populations and the eventual interbreeding with domestic cats remain poorly studied. We surveyed genetic variation at 12 autosomal microsatellites for 34 wild and 64 domestic cats collected across Portugal. Wild and domestic cats were significantly differentiated both at allele frequencies and sizes (F _ST=0.11, R _ST = 0.18, P < 0.001). Population structure and admixture analyses performed using Bayesian approaches also showed evidence of two discrete groups clustering wild and domestic populations. Results did not show significant genetic divergence among Northern, Central and Southern wildcats. Six morphologically identified wildcats were significantly assigned to the domestic cluster, revealing some discrepancy between phenotypic and genetic identifications. We detected four hybrids (approximately 14%) using a consensus analysis of different Bayesian model-based software. These hybrids were identified throughout all sampled areas, suggesting that hybridisation is of major concern for the appropriate implementation of wildcat conservation strategies in Portugal.     

Populations at risk: conservation genetics of kangaroo mice (Microdipodops) of the Great Basin Desert

The Great Basin Desert of western North America has experienced frequent habitat alterations due to a complex biogeographic history and recent anthropogenic impacts, with the more recent alterations likely resulting in the decline of native fauna and flora. Dark (Microdipodops megacephalus) and pallid (M. pallidus) kangaroo mice are ecological specialists found within the Great Basin Desert and are potentially ideal organisms for assessing ecosystem health and inferring the biogeographic history of this vulnerable region. Herein, newly acquired nuclear‐encoded microsatellite loci were utilized to assess patterns of variation within and among spatially discrete groups of kangaroo mice and to evaluate gene flow, demographic trends, and genetic integrity. Results confirm that there are at least three genetically distinct units within M. megacephalus and two such units within M. pallidus. The three units of M. megacephalus appear to have different demographic histories, with effectively no gene flow among them since their divergence. Similarly, the two units of M. pallidus also appear to have experienced different demographic histories, with effectively no gene exchange. Contemporary effective population sizes of all groups within Microdipodops appear to be low (<500), suggesting that each genetic lineage may have difficulty coping with changing environmental pressures and hence may be at risk of extirpation. Results of this study indicate that each Microdipodops group should be recognized, and therefore managed, as a separate unit in an effort to conserve these highly specialized taxa that contribute to the diversity of the Great Basin Desert ecosystem.     

Geographic variation in the Pine Barrens Treefrog (Hyla andersonii): concordance of genetic,morphometric and acoustic signal data

Delimiting species is important to every subfield in biology. Templeton's cohesion species concept uses genetic and ecological exchangeability to identify sets of populations that ought to be considered as the same species, and the lack of exchangeability helps determine which populations can be grouped as evolutionarily significant units (ESU) in conservation science. However, previous work assessing genetic and ecological interchangeability among populations has been limited in scope. Here, we provide a method for assessing exchangeability that incorporates multiple, independent lines of multivariate evidence in genetic, behavioural and morphological data. We use this approach to assess exchangeability across three disjunct groups of populations of the Pine Barrens Treefrog (Hyla andersonii) from the eastern United States. This species is considered threatened by each state in which it occurs and conservation management of this taxon requires a clearer understanding of how populations in these three regions may differ from one another. We find a strikingly concordant pattern in which the first axis of variation for each of the three types of data distinguishes populations along a latitudinal gradient and the second axis distinguishes the set of populations occurring in the Carolinas from those occurring in the New Jersey and Florida/Alabama regions. We know of no comparable data set that displays such concordance among different types of data across so large a geographic range. The overlap in trait values (i.e. exchangeability) between neighbouring regions, however, is substantial in all three types of data, which supports continued consideration of this taxon as a single species.     

Speciation and genetic divergence of three species of charr from ancient Lake El'gygytgyn (Chukotka) and their phylogenetic relationships with other representatives of the genus Salvelinus

The diversity of phenotypically different and often reproductively isolated lacustrine forms of charrs of the genus Salvelinus represents a substantial problem for taxonomists and evolutionary biologists. Based on the analysis of variability of ten microsatellite loci and two fragments of mitochondrial DNA (control region and cyt‐b gene), the evolutionary history of three charr species from Lake El'gygytgyn was reconstructed, and phylogenetic relationships between the main representatives of the genus were revealed. Three species from Lake El'gygytgyn were strongly reproductively isolated. Long‐finned charr described previously as Salvethymus svetovidovi, an ancient endemic form in the lake, originated 3.5 Mya (95% Bayesian credible intervals: 1.7, 6.1). Placement of this species in the phylogenetic tree of Salvelinus was not determined strictly, but it should be located in the basal part of the clade Salvelinus alpinus – S. malma species complex. The origin of small‐mouth charr S. elgyticus and Boganida charr S. boganidae in Lake El'gygytgyn was related to allopatric speciation. Their ancestors were represented by two glacial lineages of Taranets charr S. alpinus taranetzi from Asia. In Lake El'gygytgyn, these lineages entered into secondary contact postglacially. A revision of the main phylogenetic groups within the Salvelinus alpinus – S. malma complex is conducted. The Boganida charrs from Lakes El'gygytgyn and Lama (Taimyr) belong to different phylogenetic groups of Arctic charr and should not be regarded as a single species S. boganidae. Using the charrs from Lakes El'gygytgyn and Lama as a case study, we show that a model of sympatric speciation, which seemed more probable based on previous empirical evidence, was rejected by other data.     

Variability in the efficacy of sex pheromone lures for monitoring oriental fruit moth (Lepidoptera: Tortricidae)

Studies were conducted in Chile and the United States to compare the attractiveness of various commercial sex pheromone lures and two experimental lures for oriental fruit moth, Grapholita molesta (Busck), in peach orchards treated with or without sex pheromone dispensers. The experimental lures contained the three‐component sex pheromone blend of G. molesta: Z‐8‐dodecenyl acetate, E‐8‐dodecenyl acetate and Z‐8‐dodecenol (Z8‐12:OH), and the sex pheromone of codling moth, Cydia pomonella (L.), (E,E)‐8,10‐dodecadien‐1‐ol, (codlemone). Commercial lures varied in their substrate, initial loading and blend ratio of components. Significant differences in male catches were found among commercial lures in orchards treated with or without sex pheromone dispensers. Experimental lures with the addition of codlemone significantly increased the catches of G. molesta using lures loaded with 0%, 1% or 5% Z8‐12:OH in the G. molesta blend compared with the same ratio of components in just the G. molesta blend. The experimental lures were significantly more attractive than all commercial lures in the untreated orchard. However, moth catch with the experimental lures in the sex pheromone‐treated orchard was only intermediate among all of the lures tested. These findings highlight the need to develop more effective and standardized lures that can be used in trap‐based monitoring programme for this important pest.     

Spatial genetic analysis reveals high connectivity of tiger (Panthera tigris) populations in the Satpura–Maikal landscape of Central India

We investigated the spatial genetic structure of the tiger meta‐population in the Satpura–Maikal landscape of central India using population‐ and individual‐based genetic clustering methods on multilocus genotypic data from 273 individuals. The Satpura–Maikal landscape is classified as a global‐priority Tiger Conservation Landscape (TCL) due to its potential for providing sufficient habitat that will allow the long‐term persistence of tigers. We found that the tiger meta‐population in the Satpura–Maikal landscape has high genetic variation and very low genetic subdivision. Individual‐based Bayesian clustering algorithms reveal two highly admixed genetic populations. We attribute this to forest connectivity and high gene flow in this landscape. However, deforestation, road widening, and mining may sever this connectivity, impede gene exchange, and further exacerbate the genetic division of tigers in central India.     

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M‐lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C‐lineage). Effective conservation actions require reliable molecular tools to identify pure‐bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q‐values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.     

Development and characterization of EST‐SSR markers for the genus Rhododendron section Brachycalyx (Ericaceae)

Simple sequence repeat (SSR) markers were developed from expressed sequence tags (ESTs) for Rhododendron section Brachycalyx in order to elucidate its evolutionary processes and reproductive ecology. Nineteen polymorphic EST‐SSR markers were developed from EST libraries of R. amagianum and R. hyugaense. Polymorphisms for these markers were assessed using four species of section Brachycalyx. The number of alleles ranged from 1 to 14, and the observed and expected heterozygosity ranged from 0.000 to 0.931 and 0.000 to 0.904, respectively. The EST‐SSR markers developed in this study will be useful for elucidating population genetic structure and breeding systems in section Brachycalyx.     

Population genetic patterns among social groups of the endangered Central American spider monkey (Ateles geoffroyi) in a human‐dominated landscape

Spider monkeys (Genus: Ateles) are a widespread Neotropical primate with a highly plastic socioecological strategy. However, the Central American species, Ateles geoffroyi, was recently re‐listed as endangered due to the accelerated loss of forest across the subcontinent. There is inconsistent evidence that spider monkey populations could persist when actively protected, but their long‐term viability in unprotected, human‐dominated landscapes is not known. We analyzed noninvasive genetic samples from 185 individuals in 14 putative social groups on the Rivas Isthmus in southwestern Nicaragua. We found evidence of weak but significant genetic structure in the mitochondrial control region and in eight nuclear microsatellite loci plus negative spatial autocorrelation in Fst and kinship. The overall pattern suggests strong localized mating and at least historical female‐biased dispersal, as is expected for this species. Heterozygosity was significantly lower than expected under random mating and lower than that found in other spider monkey populations, possibly reflecting a recent decline in genetic diversity and a threat from inbreeding. We conclude that despite a long history of human disturbance on this landscape, spider monkeys were until recently successful at maintaining gene flow. We consider the recent decline to be further indication of accelerated anthropogenic disturbance, but also of an opportunity to conserve native biodiversity. Spider monkeys are one of many wildlife species in Central America that is threatened by land cover change, and an apt example of how landscape‐scale conservation planning could be used to ensure long‐term persistence.     

Development of nuclear microsatellite markers for the threatened wetland plant Geranium soboliferum var. kiusianum (Geraniaceae)

Nuclear microsatellite markers were developed for the threatened plant Geranium soboliferum var. kiusianum, which has decreased its population size as a result of loss of its wetland habitat in Kyushu, Japan. Utilizing RNA‐seq data obtained by next‐generation sequencing techniques, 10 polymorphic microsatellite markers with 3–16 alleles in a nuclear genome were developed and characterized. Two to 15 alleles were observed in G. soboliferum. These markers will be used to investigate the genetic circumstance of remnant populations of G. soboliferum var. kiusianum and their phylogenetic relationship with G. soboliferum.     

The spatial genetic structure of the White‐banded Tanager (Aves,Passeriformes) in fragmented Neotropical savannas suggests two evolutionarily significant units

Understanding the genetic structure of a species is crucial for evolutionary biology research and species conservation. The objectives of this study were to investigate the genetic structure of Neothraupis fasciata in Brazilian savannas and to assess genetic differentiation of its disjunct population in the Amazonian savannas of the state of Amapá. Population genetic structure was assessed in relation to isolation by distance and landscape variables connected with habitat heterogeneity. The influences of factors, such as habitat fragmentation and core–periphery distribution, on genetic diversity were also examined. Data were derived from a set of microsatellite loci of adult individuals from nine localities: eight distributed across the Cerrado and one in the disjunct Amazonian savanna of Amapá. Analysis revealed moderate genetic diversity and moderate population genetic structure, with at least two genetic clusters, one of which is represented exclusively by the disjunct Amapá population. The genetic structure found is not the result of significant influences by geographical distance, habitat heterogeneity, or the core–periphery effect, nor by intense biome fragmentation due to anthropic action. The disjunct Amapá population exhibited a moderate level of genetic differentiation compared to the Cerrado population, suggesting that both can be considered distinct evolutionarily significant units for conservation purposes. Abstract in Portuguese is available with online material.     

Limber pine (Pinus flexilis James) genetic map constructed by exome‐seq provides insight into the evolution of disease resistance and a genomic resource for genomics‐based breeding

Limber pine ( Pinus flexilis ) is a keystone species of high‐elevation forest ecosystems of western North America, but some parts of the geographic range have high infection and mortality from the non‐native white pine blister rust caused by Cronartium ribicola . Genetic maps can provide essential knowledge for understanding genetic disease resistance as well as local adaptation to changing climates. Exome‐seq was performed to construct high‐density genetic maps in two seed families. Composite maps positioned 9612 unigenes across 12 linkage groups ( LG s). Syntenic analysis of genome structure revealed that the majority of orthologs were positional orthologous genes ( POG s) with localization on homologous LG s among conifer species. Gene ontology ( GO) enrichment analysis showed relatively fewer constraints for POG s with putative roles in adaptation to environments and relatively more conservation for POG s with roles in basic cell function and maintenance. The mapped genes included 639 nucleotide‐binding site leucine‐rich repeat genes ( NBS ‐ LRR s) , 290 receptor‐like protein kinase genes ( RLK s), and 1014 genes with potential roles in the defense response and induced systemic resistance to attack by pathogens. Orthologous loci for resistance to rust pathogens were identified and were co‐positioned with multiple members of the R gene family, revealing the evolutionary pressure acting upon them. This high‐density genetic map provides a genomic resource and practical tool for breeding and genetic conservation programs, with applications in genome‐wide association studies ( GWASs ), the characterization of functional genes underlying complex traits, and the sequencing and assembly of the full‐length genomes of limber pine and related Pinus species.