台灣油杉(Keteleeria davidiana var. formosana)族群於過去AFLP (amplified fragment length polymorphisms)研究中發現，因為分歧天擇作用而發生在地適應性演化。本研究透過ddRADseq (double digest Restriction Associated DNA sequencing)以檢測天擇作用如何造成台灣油杉族群基因組的分歧。收集台灣油杉五個族群共62株樣本及10株與其同屬的近源種鐵堅油杉(Keteleeria davidiana)。ddRADseq獲得15,094個單一核甘酸多形性(single nucleotide polymorphism, SNP)。BayeScan的結果找到61個FST outliers。LFMM及BayEnv則分別發現288跟221個與溫度、降水、地貌、植被豐富度具關連性的基因座。其中包含57個FST outliers，顯示台灣油杉族群受到與棲地環境差異有關的分歧天擇作用。然而台灣油杉兩兩族群間L形的FST頻率分布圖，顯示大部分基因座的分歧程度低。鐵堅油杉與台灣油杉的FST頻率分布則有右移的趨勢。此外TreeMix及f3 statixtic發現台灣油杉存在由南部往北部族群的基因交流。然而台灣油杉與鐵堅油杉並沒有發現基因交流，加上分歧的時間較長，遺傳漂變可能才是造成兩個近緣物種遺傳分歧的主因。

In a previous study, divergent selection and local adaptation in populations of Keteleeria davidiana var. formosana has been identified based on survey of genetic variation using AFLP. In the present study, we aimed to investigate how population genomic divergence can be shaped by natural selection in the populations of K. davidiana var. formosana using method of ddRADseq. We collected 62 samples from five populations K. davidiana var. formosana and 10 from its closely related species Keteleeria davidiana. ddRADseq generated 15,094 SNPs. Sixty-one FST outliers were identified by BayeScan. LFMM and BayEnv identified 288 and 221 loci potentially under selection, respectively. Of these 57 were found to be associated with temperature, precipitation, topography and vegetation abundance, indicating ecological divergent in K. davidiana var. formosana populations. However , many loci showed low levels of divergence between K. davidiana var. formosana populations, resulting in L shape pairwise FST distributions. Pairwise FST distributions shifted to right when compared populations of K. davidiana var. formosana to K. davidiana. Results of TreeMix and f3 statistic suggest significant gene flow was mainly from south to north populations between K. davidiana var. formosana. However, no significant gene flow was found between K. davidiana var. formosana and K. davidiana, probably because they have been diverged for a long time, and genetic drift is most probably responsible for the genetic divergence between these two closely related species.

Andrew RL, Rieseberg LH (2013) Divergence is focused on few genomic regions early in speciation: incipient speciation of sunflower ecotypes. Evolution 67, 2468-2482.
Blair LM, Granka JM, Feldman MW (2014) On the stability of the Bayenv method in assessing human SNP-environment associations. Human Genomics 8, 1.
Bolnick DI, Fitzpatrick BM (2007) Sympatric speciation: models and empirical evidence. Annual Review of Ecology Evolution and Systematics 38, 459-487.
Borcard D, Gillet F, Legendre P (2011) Numerical Ecology with R. Springer Verlag, New York.
Bräutigam A, Gowik U (2010) What can next generation sequencing do for you? Next generation sequencing as a valuable tool in plant research. Plant Biology 12, 831-841.
Carrera-Hernandez JJ, Gaskin SJ (2007) Spatio temporal analysis of daily precipitation and temperature in the Basin of Mexico. Journal of Hydrology 336, 231-249.
Catchen J, Hohenlohe PA, Bassham S, Amores A, Cresko WA (2013) Stacks: an analysis tool set for population genomics. Molecular Ecology 22, 3124-3140.
Catchen JM, Amores A, Hohenlohe P, Cresko W, Postlethwait JH (2011) Stacks: building and genotyping loci de novo from short-read sequences. G3-Genes Genomes Genetics 1, 171-182.
Chattopadhyay B, Garg KM, Ramakrishnan U (2014) Effect of diversity and missing data on genetic assignment with RAD-Seq markers. BMC Research Notes 7, 841.
Coop G, Witonsky D, Di Rienzo A, Pritchard JK (2010) Using environmental correlations to identify loci underlying local adaptation. Genetics 185, 1411-1423.
Coyne JA, Orr HA (2004) Speciation Sinauer & Associates, Sunderland, Massachusetts.
Davey JW, Hohenlohe PA, Etter PD, et al. (2011) Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nature Reviews Genetics 12, 499-510.
De Mita S, Thuillet AC, Gay L, et al. (2013) Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations. Molecular Ecology 22, 1383-1399.
de Villemereuil P, Frichot E, Bazin E, Francois O, Gaggiotti OE (2014) Genome scan methods against more complex models: when and how much should we trust them? Molecular Ecology 23, 2006-2019.
Dehestani A, Kazemi Tabar S (2007) A rapid efficient method for DNA isolation from plants with gigh levels of secondary metabolites. Asian Journal of Plant Sciences 6, 977-981.
Dell'Acqua M, Zuccolo A, Tuna M, Gianfranceschi L, Pe ME (2014) Targeting environmental adaptation in the monocot model Brachypodium distachyon: a multi-faceted approach. BMC Genomics 15, 801.
Fang J-Y, Chung J-D, Chiang Y-C, et al. (2013) Divergent selection and local adaptation in disjunct populations of an endangered conifer, Keteleeria davidiana var. formosana (Pinaceae). PLoS ONE 8, e70162.
Feder JL, Egan SP, Nosil P (2012a) The genomics of speciation-with-gene-flow. Trends in Genetics 28, 342-350.
Feder JL, Flaxman SM, Egan SP, Comeault AA, Nosil P (2013) Geographic mode of speciation and genomic divergence. Annual Review of Ecology, Evolution, and Systematics 44, 73-97.
Feder JL, Gejji R, Yeaman S, Nosil P (2012b) Establishment of new mutations under divergence and genome hitchhiking. Philosophical Transactions of the Royal Society B-Biological Sciences 367, 461-474.
Feldstein J (2009) PHYLIP (Phylogeny Inference Package) version 3.69. Department of Genetics, University of Washington, Seattle, WA.
Foll M, Gaggiotti O (2008) A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics 180, 977-993.
Frichot E, Schoville SD, Bouchard G, Francois O (2013) Testing for associations between loci and environmental gradients using latent factor mixed models. Molecular Biology and Evolution 30, 1687-1699.
Gunther T, Coop G (2013) Robust identification of local adaptation from allele frequencies. Genetics 195, 205-220.
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25, 1965-1978.
Hohenlohe PA, Amish SJ, Catchen JM, Allendorf FW, Luikart G (2011) Next-generation RAD sequencing identifies thousands of SNPs for assessing hybridization between rainbow and westslope cutthroat trout. Molecular Ecology Resources 11, 117-122.
Hohenlohe PA, Bassham S, Etter PD, et al. (2010) Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags. PLoS Genetics 6, e1000862.
Huang C-L, Chen J-H, Tsang M-H, et al. (2015) Influences of environmental and spatial factors on genetic and epigenetic variations in Rhododendron oldhamii (Ericaceae). Tree Genetics & Genomes 11, 1-16.
Jeffreys H (1961) Theory of probability. Oxford University Press, Oxford UK.
Jombart T, Devillard S, Balloux F (2010) Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genetics 11, 94.
Jukes TH, Cantor CR (1969) Evolution of protein molecules. Mammalian Protein Metabolism 3, 21-132.
Kang J, Park J, Choi H, et al. (2011) Plant ABC transporters. The Arabidopsis book/American Society of Plant Biologists 9, e10153.
Klein M, Geisler M, Suh SJ, et al. (2004) Disruption of AtMRP4, a guard cell plasma membrane ABCC-type ABC transporter, leads to deregulation of stomatal opening and increased drought susceptibility. Plant Journal 39, 219-236.
Li H, Hsuan K (1996) Flora of Taiwan, 2nd edition. Committee of the Flora of Taiwan, Taipei, Taiwan, R.O.C.
Manel S, Holderegger R (2013) Ten years of landscape genetics. Trends in Ecology & Evolution 28, 614-621.
Mastretta‐Yanes A, Arrigo N, Alvarez N, et al. (2015) Restriction site‐associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference. Molecular Ecology Resources 15, 28-41.
Mayr E (1963) Animal species and evolution, Harvard University Press, Cambridge, Massachusetts.
Michel AP, Sim S, Powell THQ, et al. (2010) Widespread genomic divergence during sympatric speciation. Proceedings of the National Academy of Sciences of the United States of America 107, 9724-9729.
Naimi B, Hamm NAS, Groen TA, Skidmore AK, Toxopeus AG (2014) Where is positional uncertainty a problem for species distribution modelling? Ecography 37, 191-203.
Nielsen R, Paul JS, Albrechtsen A, Song YS (2011) Genotype and SNP calling from next-generation sequencing data. Nature Reviews Genetics 12, 443-451.
Nosil P (2008) Speciation with gene flow could be common. Molecular Ecology 17, 2103-2106.
Nosil P, Feder JL (2012) Genomic divergence during speciation: causes and consequences Introduction. Philosophical Transactions of the Royal Society B-Biological Sciences 367, 332-342.
Nosil P, Funk DJ, Ortiz-Barrientos D (2009a) Divergent selection and heterogeneous genomic divergence. Molecular Ecology 18, 375-402.
Nosil P, Harmon LJ, Seehausen O (2009b) Ecological explanations for (incomplete) speciation. Trends in Ecology & Evolution 24, 145-156.
Nystedt B, Street NR, Wetterbom A, et al. (2013) The Norway spruce genome sequence and conifer genome evolution. Nature 497, 579-584.
Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstra HE (2012) Double digest RADseq: Aan inexpensive method for de novo SNP siscovery and genotyping in model and non-model species. PLoS ONE 7, e37135.
Pickrell JK, Pritchard JK (2012) Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genetics 8, e1002967.
Reich D, Thangaraj K, Patterson N, Price AL, Singh L (2009) Reconstructing Indian population history. Nature 461, 489-494.
Ren J, Kotaka M, Lockyer M, et al. (2005) GTP cyclohydrolase II structure and mechanism. Journal of Biological Chemistry 280, 36912-36919.
Roesti M, Salzburger W, Berner D (2012) Uninformative polymorphisms bias genome scans for signatures of selection. BMC Evolutionary Biology 12, 94.
Rundle HD, Nosil P (2005) Ecological speciation. Ecology Letters 8, 336-352.
Schluter D (2001) Ecology and the origin of species. Trends in Ecology & Evolution 16, 372-380.
Schoville SD, Bonin A, Francois O, et al. (2012) Adaptive genetic variation on the landscape: methods and cases. Annual Review of Ecology, Evolution, and Systematics 43, 23-43.
Shendure J, Ji HL (2008) Next-generation DNA sequencing. Nature Biotechnology 26, 1135-1145.
Strasburg JL, Sherman NA, Wright KM, et al. (2012) What can patterns of differentiation across plant genomes tell us about adaptation and speciation? Philosophical Transactions of the Royal Society B-Biological Sciences 367, 364-373.
Turner TL, Hahn MW, Nuzhdin SV (2005) Genomic islands of speciation in Anopheles gambiae. PLoS Biology 3, 1572-1578.
Via S (2009) Natural selection in action during speciation. Proceedings of the National Academy of Sciences of the United States of America 106, 9939-9946.
Via S (2012) Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow. Philosophical Transactions of the Royal Society B-Biological Sciences 367, 451-460.
Via S, West J (2008) The genetic mosaic suggests a new role for hitchhiking in ecological speciation. Molecular Ecology 17, 4334-4345.
Wagner CE, Keller I, Wittwer S, et al. (2013) Genome-wide RAD sequence data provide unprecedented resolution of species boundaries and relationships in the Lake Victoria cichlid adaptive radiation. Molecular Ecology 22, 787-798.
Williams GC (1966) Adaptation and natural selection; a critique of some current evolutionary thought Princeton University Press, Princeton NJ.