Global investigation of the co‐evolution of MIRNA genes and microRNA targets during soybean domestication

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome‐wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA–target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA–target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA–target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA–target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co‐evolution of MIRs and miRNA targets during soybean domestication.     

Comparative population genomics reveals the domestication history of the peach,Prunus persica,and human influences on perennial fruit crops

Background

Recently, many studies utilizing next generation sequencing have investigated plant evolution and domestication in annual crops. Peach, Prunus persica, is a typical perennial fruit crop that has ornamental and edible varieties. Unlike other fruit crops, cultivated peach includes a large number of phenotypes but few polymorphisms. In this study, we explore the genetic basis of domestication in peach and the influence of humans on its evolution.

Results

We perform large-scale resequencing of 10 wild and 74 cultivated peach varieties, including 9 ornamental, 23 breeding, and 42 landrace lines. We identify 4.6 million SNPs, a large number of which could explain the phenotypic variation in cultivated peach. Population analysis shows a single domestication event, the speciation of P. persica from wild peach. Ornamental and edible peach both belong to P. persica, along with another geographically separated subgroup, Prunus ferganensis.We identify 147 and 262 genes under edible and ornamental selection, respectively. Some of these genes are associated with important biological features. We perform a population heterozygosity analysis in different plants that indicates that free recombination effects could affect domestication history. By applying artificial selection during the domestication of the peach and facilitating its asexual propagation, humans have caused a sharp decline of the heterozygote ratio of SNPs.

Conclusions

Our analyses enhance our knowledge of the domestication history of perennial fruit crops, and the dataset we generated could be useful for future research on comparative population genomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0415-1) contains supplementary material, which is available to authorized users.     

A draft genome of sweet cherry (Prunus avium L.) reveals genome‐wide and local effects of domestication

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety “Big Star*” and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome‐wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.     

Genetic basis of kernel nutritional traits during maize domestication and improvement

The nutritional traits of maize kernels are important for human and animal nutrition, and these traits have undergone selection to meet the diverse nutritional needs of humans. However, our knowledge of the genetic basis of selecting for kernel nutritional traits is limited. Here, we identified both single and epistatic quantitative trait loci (QTLs) that contributed to the differences of oil and carotenoid traits between maize and teosinte. Over half of teosinte alleles of single QTLs increased the values of the detected oil and carotenoid traits. Based on the pleiotropism or linkage information of the identified single QTLs, we constructed a trait–locus network to help clarify the genetic basis of correlations among oil and carotenoid traits. Furthermore, the selection features and evolutionary trajectories of the genes or loci underlying variations in oil and carotenoid traits revealed that these nutritional traits produced diverse selection events during maize domestication and improvement. To illustrate more, a mutator distance–relative transposable element (TE) in intron 1 of DXS2, which encoded a rate‐limiting enzyme in the methylerythritol phosphate pathway, was identified to increase carotenoid biosynthesis by enhancing DXS2 expression. This TE occurs in the grass teosinte, and has been found to have undergone selection during maize domestication and improvement, and is almost fixed in yellow maize. Our findings not only provide important insights into evolutionary changes in nutritional traits, but also highlight the feasibility of reintroducing back into commercial agricultural germplasm those nutritionally important genes hidden in wild relatives.     

Plant domestication slows pest evolution

Agricultural practices such as breeding resistant varieties and pesticide use can cause rapid evolution of pest species, but it remains unknown how plant domestication itself impacts pest contemporary evolution. Using experimental evolution on a comparative phylogenetic scale, we compared the evolutionary dynamics of a globally important economic pest – the green peach aphid (Myzus persicae) – growing on 34 plant taxa, represented by 17 crop species and their wild relatives. Domestication slowed aphid evolution by 13.5%, maintained 10.4% greater aphid genotypic diversity and 5.6% higher genotypic richness. The direction of evolution (i.e. which genotypes increased in frequency) differed among independent domestication events but was correlated with specific plant traits. Individual‐based simulation models suggested that domestication affects aphid evolution directly by reducing the strength of selection and indirectly by increasing aphid density and thus weakening genetic drift. Our results suggest that phenotypic changes during domestication can alter pest evolutionary dynamics.     

Exploring a Tomato Landraces Collection for Fruit-Related Traits by the Aid of a High-Throughput Genomic Platform

During its evolution and domestication Solanum lycopersicum has undergone various genetic ‘bottlenecks’ and extreme inbreeding of limited genotypes. In Europe the tomato found a secondary centre for diversification, which resulted in a wide array of fruit shape variation given rise to a range of landraces that have been cultivated for centuries. Landraces represent a reservoir of genetic diversity especially for traits such as abiotic stress resistance and high fruit quality. Information about the variation present among tomato landrace populations is still limited. A collection of 123 genotypes from different geographical areas was established with the aim of capturing a wide diversity. Eighteen morphological traits were evaluated, mainly related to the fruit. About 45% of morphological variation was attributed to fruit shape, as estimated by the principal component analysis, and the dendrogram of relatedness divided the population in subgroups mainly on the basis of fruit weight and locule number. Genotyping was carried out using the tomato array platform SolCAP able to interrogate 7,720 SNPs. In the whole collection 87.1% markers were polymorphic but they decreased to 44–54% when considering groups of genotypes with different origin. The neighbour-joining tree analysis clustered the 123 genotypes into two main branches. The STRUCTURE analysis with K = 3 also divided the population on the basis of fruit size. A genomic-wide association strategy revealed 36 novel markers associated to the variation of 15 traits. The markers were mapped on the tomato chromosomes together with 98 candidate genes for the traits analyzed. Six regions were evidenced in which candidate genes co-localized with 19 associated SNPs. In addition, 17 associated SNPs were localized in genomic regions lacking candidate genes. The identification of these markers demonstrated that novel variability was captured in our germoplasm collection. They might also provide a viable indirect selection tool in future practical breeding programs.     

Animal evolution during domestication: the domesticated fox as a model

We review the evolution of domestic animals, emphasizing the effect of the earliest steps of domestication on its course. Using the first domesticated species, the dog (Canis familiaris), for illustration, we describe the evolutionary peculiarities during the historical domestication, such as the high level and wide range of diversity. We suggest that the process of earliest domestication via unconscious and later conscious selection of human‐defined behavioral traits may accelerate phenotypic variations. The review is based on the results of a long‐term experiment designed to reproduce early mammalian domestication in the silver fox (Vulpes vulpes) selected for tameability or amenability to domestication. We describe changes in behavior, morphology and physiology that appeared in the fox during its selection for tameability, which were similar to those observed in the domestic dog. Based on the data of the fox experiment and survey of relevant data, we discuss the developmental, genetic and possible molecular genetic mechanisms underlying these changes. We ascribe the causative role in evolutionary transformation of domestic animals to the selection for behavior and to the neurospecific regulatory genes it affects.     

利用多组学手段解析番茄育种过程中代谢物变化的机制

番茄(Solanum lycopersicum)在育种过程中经历了驯化、改良、分化和渐渗等不同阶段, 在这一选择过程中番茄的果重和风味等均发生了显著改变, 但是目前对于番茄育种过程中代谢物的变化及其遗传基础却不是十分清楚。近期, 中国农业科学院深圳农业基因组研究所黄三文研究组与华中农业大学罗杰研究组利用多组学(变异组、转录组及代谢组)手段系统解析了番茄育种过程中代谢物的变化。结果表明, 在番茄驯化过程中有46个甾醇糖基生物碱类物质(SGAs)含量逐渐降低, 并获得了7个与其中44个物质显著相关的遗传位点。因此, 在番茄育种过程中通过优异等位位点的组合可以显著降低SGAs的含量; 同时发现在番茄以果重为目标的选择过程中, 控制果重基因周围其它基因的“搭车效应”是引起许多代谢物变化的重要遗传因素, 及在育种过程中对某一性状的选择会对其它性状产生重要影响。该研究首次利用多组学手段系统解析了选择对作物代谢物的影响, 为番茄品质改良奠定了良好的理论基础。     

Genome‐wide analysis of polymorphisms identified domestication‐associated long low‐diversity region carrying important rice grain size/weight quantitative trait loci

Rice grain size and weight are major determinants of grain quality and yield and so have been under rigorous selection since domestication. However, the genetic basis for contrasting grain size/weight trait among Indian germplasms and their association with domestication‐driven evolution is not well understood. In this study, two long (LGG) and two short grain (SGG) genotypes were resequenced. LGG (LGR and PB 1121) differentiated from SGG (Sonasal and Bindli) by 504 439 single nucleotide polymorphisms (SNPs) and 78 166 insertion‐and‐deletion polymorphisms. The LRK gene cluster was different and a truncation mutation in the LRK8 kinase domain was associated with LGG. Phylogeny with 3000 diverse rice accessions revealed that the four sequenced genotypes belonged to the japonica group and were at the edge of the clades indicating them to be the potential source of genetic diversity available in Indian rice germplasm. Six SNPs were significantly associated with grain size/weight and the top four of these could be validated in mapping a population, suggesting this study as a valuable resource for high‐throughput genotyping. A contiguous long low‐diversity region (LDR) of approximately 6 Mb carrying a major grain weight quantitative trait loci (harbouring OsTOR gene) was identified on Chromosome 5. This LDR was identified as an evolutionary important site with significant positive selection and multiple selection sweeps, and showed association with many domestication‐related traits, including grain size/weight. The aus population retained more allelic variations in the LDR than the japonica and indica populations, suggesting it to be one of the divergence loci. All the data and analyses can be accessed from the RiceSzWtBase database.     

Identification of genetic variants associated with maize flowering time using an extremely large multi‐genetic background population

Flowering time is one of the major adaptive traits in domestication of maize and an important selection criterion in breeding. To detect more maize flowering time variants we evaluated flowering time traits using an extremely large multi‐ genetic background population that contained more than 8000 lines under multiple Sino‐United States environments. The population included two nested association mapping (NAM) panels and a natural association panel. Nearly 1 million single‐nucleotide polymorphisms (SNPs) were used in the analyses. Through the parallel linkage analysis of the two NAM panels, both common and unique flowering time regions were detected. Genome wide, a total of 90 flowering time regions were identified. One‐third of these regions were connected to traits associated with the environmental sensitivity of maize flowering time. The genome‐wide association study of the three panels identified nearly 1000 flowering time‐associated SNPs, mainly distributed around 220 candidate genes (within a distance of 1 Mb). Interestingly, two types of regions were significantly enriched for these associated SNPs – one was the candidate gene regions and the other was the approximately 5 kb regions away from the candidate genes. Moreover, the associated SNPs exhibited high accuracy for predicting flowering time.     

Using Association Mapping in Teosinte to Investigate the Function of Maize Selection-Candidate Genes

Background

Large-scale screens of the maize genome identified 48 genes that show the putative signature of artificial selection during maize domestication or improvement. These selection-candidate genes may act as quantitative trait loci (QTL) that control the phenotypic differences between maize and its progenitor, teosinte. The selection-candidate genes appear to be located closer in the genome to domestication QTL than expected by chance.

Methods and Findings

As a step toward defining the traits controlled by these genes, we performed phenotype-genotype association mapping in teosinte for 32 of the 48 plus three other selection-candidate genes. Our analyses assayed 32 phenotypic traits, many of which were altered during maize domestication or improvement. We observed several significant associations between SNPs in the selection-candidate genes and trait variation in teosinte. These included two associations that surpassed the Bonferroni correction and five instances where a gene significantly associated with the same trait in both of our association mapping panels. Despite these significant associations, when compared as a group the selection-candidate genes performed no better than randomly chosen genes.

Conclusions

Our results suggest association analyses can be helpful for identifying traits under the control of selection-candidate genes. Indeed, we present evidence for new functions for several selection-candidate genes. However, with the current set of selection-candidate genes and our association mapping strategy, we found very few significant associations overall and no more than we would have found with randomly chosen genes. We discuss possible reasons that a large number of significant genotype-phenotype associations were not discovered.     

Genomic variation at the tips of the adaptive radiation of Darwin's finches

Adaptive radiation unfolds as selection acts on the genetic variation underlying functional traits. The nature of this variation can be revealed by studying the tips of an ongoing adaptive radiation. We studied genomic variation at the tips of the Darwin's finch radiation; specifically focusing on polymorphism within, and variation among, three sympatric species of the genus Geospiza. Using restriction site‐associated DNA (RAD‐seq), we characterized 32 569 single‐nucleotide polymorphisms (SNPs), from which 11 outlier SNPs for beak and body size were uncovered by a genomewide association study (GWAS). Principal component analysis revealed that these 11 SNPs formed four statistically linked groups. Stepwise regression then revealed that the first PC score, which included 6 of the 11 top SNPs, explained over 80% of the variation in beak size, suggesting that selection on these traits influences multiple correlated loci. The two SNPs most strongly associated with beak size were near genes associated with beak morphology across deeper branches of the radiation: delta‐like 1 homologue (DLK1) and high‐mobility group AT‐hook 2 (HMGA2). Our results suggest that (i) key adaptive traits are associated with a small fraction of the genome (11 of 32 569 SNPs), (ii) SNPs linked to the candidate genes are dispersed throughout the genome (on several chromosomes), and (iii) micro‐ and macro‐evolutionary variation (roots and tips of the radiation) involve some shared and some unique genomic regions.     

Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.)

Although upland cotton (Gossypium hirsutism L.) originated in the tropics, this early maturity cotton can be planted as far north as 46°N in China due to the accumulation of numerous phenotypic and physiological adaptations during domestication. However, how the genome of early maturity cotton has been altered by strong human selection remains largely unknown. Herein, we report a cotton genome variation map generated by the resequencing of 436 cotton accessions. Whole‐genome scans for sweep regions identified 357 putative selection sweeps covering 4.94% (112 Mb) of the upland cotton genome, including 5184 genes. These genes were functionally related to flowering time control, hormone catabolism, ageing and defence response adaptations to environmental changes. A genome‐wide association study (GWAS) for seven early maturity traits identified 307 significant loci, 22.48% (69) of which overlapped with putative selection sweeps that occurred during the artificial selection of early maturity cotton. Several previously undescribed candidate genes associated with early maturity were identified by GWAS. This study provides insights into the genetic basis of early maturity in upland cotton as well as breeding resources for cotton improvement.     

Genetic and molecular regulation of fruit and plant domestication traits in tomato and pepper

Tomato and pepper are two Solanaceous fruit crops that display an enormous diversity in fruit morphology. In this review, we will present an overview of the history of tomato and pepper and discuss key plant traits that were specifically selected during domestication of the two species. The traits discussed are fruit weight, shape, colour, ripening, pungency and plant architecture. We will review these characteristics as well as the genetic loci or genes that control these features, questioning whether mutations at orthologous loci occurred independently in these two species or whether unique plant and fruit features resulted in selection at different genes.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

The predictability of genomic changes underlying a recent host shift in Melissa blue butterflies

Despite accumulating evidence that evolution can be predictable, studies quantifying the predictability of evolution remain rare. Here, we measured the predictability of genome‐wide evolutionary changes associated with a recent host shift in the Melissa blue butterfly (Lycaeides melissa). We asked whether and to what extent genome‐wide patterns of evolutionary change in nature could be predicted (i) by comparisons among instances of repeated evolution and (ii) from SNP × performance associations in a laboratory experiment. We delineated the genetic loci (SNPs) most strongly associated with host use in two L. melissa lineages that colonized alfalfa. Whereas most SNPs were strongly associated with host use in none or one of these lineages, we detected a an approximately twofold excess of SNPs associated with host use in both lineages. Similarly, we found that host‐associated SNPs in nature could also be partially predicted from SNP × performance (survival and weight) associations in a laboratory rearing experiment. But the extent of overlap, and thus degree of predictability, was somewhat reduced. Although we were able to predict (to a modest extent) the SNPs most strongly associated with host use in nature (in terms of parallelism and from the experiment), we had little to no ability to predict the direction of evolutionary change during the colonization of alfalfa. Our results show that different aspects of evolution associated with recent adaptation can be more or less predictable and highlight how stochastic and deterministic processes interact to drive patterns of genome‐wide evolutionary change.     

A genome-wide scan for signatures of directional selection in domesticated pigs

Background

Animal domestication involved drastic phenotypic changes driven by strong artificial selection and also resulted in new populations of breeds, established by humans. This study aims to identify genes that show evidence of recent artificial selection during pig domestication.

Results

Whole-genome resequencing of 30 individual pigs from domesticated breeds, Landrace and Yorkshire, and 10 Asian wild boars at ~16-fold coverage was performed resulting in over 4.3 million SNPs for 19,990 genes. We constructed a comprehensive genome map of directional selection by detecting selective sweeps using an F_ST-based approach that detects directional selection in lineages leading to the domesticated breeds and using a haplotype-based test that detects ongoing selective sweeps within the breeds. We show that candidate genes under selection are significantly enriched for loci implicated in quantitative traits important to pig reproduction and production. The candidate gene with the strongest signals of directional selection belongs to group III of the metabolomics glutamate receptors, known to affect brain functions associated with eating behavior, suggesting that loci under strong selection include loci involved in behaviorial traits in domesticated pigs including tameness.

Conclusions

We show that a significant proportion of selection signatures coincide with loci that were previously inferred to affect phenotypic variation in pigs. We further identify functional enrichment related to behavior, such as signal transduction and neuronal activities, for those targets of selection during domestication in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1330-x) contains supplementary material, which is available to authorized users.     

Effect prediction of identified SNPs linked to fruit quality and chilling injury in peach [<Emphasis Type="Italic">Prunus persica</Emphasis> (L.) Batsch]

Single nucleotide polymorphisms (SNPs) are a fundamental source of genomic variation. Large SNP panels have been developed for Prunus species. Fruit quality traits are essential peach breeding program objectives since they determine consumer acceptance, fruit consumption, industry trends and cultivar adoption. For many cultivars, these traits are negatively impacted by cold storage, used to extend fruit market life. The major symptoms of chilling injury are lack of flavor, off flavor, mealiness, flesh browning, and flesh bleeding. A set of 1,109 SNPs was mapped previously and 67 were linked with these complex traits. The prediction of the effects associated with these SNPs on downstream products from the ‘peach v1.0’ genome sequence was carried out. A total of 2,163 effects were detected, 282 effects (non-synonymous, synonymous or stop codon gained) were located in exonic regions (13.04 %) and 294 placed in intronic regions (13.59 %). An extended list of genes and proteins that could be related to these traits was developed. Two SNP markers that explain a high percentage of the observed phenotypic variance, UCD_SNP_1084 and UCD_SNP_46, are associated with zinc finger (C3HC4-type RING finger) family protein and AOX1A (alternative oxidase 1a) protein groups, respectively. In addition, phenotypic variation suggests that the observed polymorphism for SNP UCD_SNP_1084 [A/G] mutation could be a candidate quantitative trait nucleotide affecting quantitative trait loci for mealiness. The interaction and expression of affected proteins could explain the variation observed in each individual and facilitate understanding of gene regulatory networks for fruit quality traits in peach.     

Classic Selective Sweeps Revealed by Massive Sequencing in Cattle

Human driven selection during domestication and subsequent breed formation has likely left detectable signatures within the genome of modern cattle. The elucidation of these signatures of selection is of interest from the perspective of evolutionary biology, and for identifying domestication-related genes that ultimately may help to further genetically improve this economically important animal. To this end, we employed a panel of more than 15 million autosomal SNPs identified from re-sequencing of 43 Fleckvieh animals. We mainly applied two somewhat complementary statistics, the integrated Haplotype Homozygosity Score (iHS) reflecting primarily ongoing selection, and the Composite of Likelihood Ratio (CLR) having the most power to detect completed selection after fixation of the advantageous allele. We find 106 candidate selection regions, many of which are harboring genes related to phenotypes relevant in domestication, such as coat coloring pattern, neurobehavioral functioning and sensory perception including KIT, MITF, MC1R, NRG4, Erbb4, TMEM132D and TAS2R16, among others. To further investigate the relationship between genes with signatures of selection and genes identified in QTL mapping studies, we use a sample of 3062 animals to perform four genome-wide association analyses using appearance traits, body size and somatic cell count. We show that regions associated with coat coloring significantly (P<0.0001) overlap with the candidate selection regions, suggesting that the selection signals we identify are associated with traits known to be affected by selection during domestication. Results also provide further evidence regarding the complexity of the genetics underlying coat coloring in cattle. This study illustrates the potential of population genetic approaches for identifying genomic regions affecting domestication-related phenotypes and further helps to identify specific regions targeted by selection during speciation, domestication and breed formation of cattle. We also show that Linkage Disequilibrium (LD) decays in cattle at a much faster rate than previously thought.     

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.