遺傳漂變及天擇間的交互作用一直是族群遺傳學家深感興趣的議題，兩者不僅能影響族群間的分歧及種化，還能影響物種本身既有的遺傳及表觀遺傳變異，使物種快速回應外在環境壓力產生適應性。臺灣杉（Taiwania cryptomerioides Hayata）目前僅分布於臺灣、中國雲緬邊境及越南北部的中高海拔山區霧林帶。我們利用AFLP及MSAP兩種顯性分子標誌工具分別評估8個臺灣杉現存族群的遺傳及表觀遺傳變異，一共得到1413個AFLP基因座及462個MSAP基因座，而且臺灣杉臺灣親緣支系群相較於亞洲大陸支系群擁有較高的多型性基因座百分比及遺傳、表觀遺傳多樣性。pairwise FST及STRUCTURE的結果顯示臺灣杉臺灣族群與亞洲大陸族群已有顯著的遺傳及表觀遺傳分化，且臺灣親緣支系群在親緣關係上與中國及越南族群相距較遠，可明顯辨認出兩個歸群（K=2）。在可能受潛在天擇作用的基因座檢測中，臺灣杉臺灣親緣支系群內因中性遺傳漂變導致臺灣族群產生分歧，使族群內並無發現任何可能受天擇作用的適應性基因座。在臺灣杉三個親緣支系群內共找到32個可能受天擇作用的outlier loci與環境因子有顯著相關，包含季節性溫度（BIO4）、季節性雨量（BIO15）及NDVI，這些環境變因不僅影響臺灣杉遺傳及表觀遺傳變異，也與遺傳漂變一同影響臺灣杉親緣支系群間的適應性分歧與族群間的分化。因此，藉由天擇保留下來或遺傳漂變所累積的遺傳及表觀遺傳變異，皆有助於臺灣杉適應因氣候變遷或環境快速改變所帶來的選汰壓力，維持臺灣杉族群的存續。

The interactions of neutral drift and naturel selection have received much attention from population geneticists because they may increase the level of differentiation between populations and potentially facilitate biological speciation. Taiwania cryptomerioides Hayata whose populations are now only distributed within mountainous cloud forest regions of northern Vietnam, Yunnan-Myanmar border of China, and Taiwan Island, ranging from 1600 – 2500 m a.s.l. We surveyed genetic and epigenetic variations of eight extant populations of Taiwania through two dominant markers, AFLP and MSAP, and then we finally obtained 1413 AFLP loci and 462 MSAP loci. The results suggested that large proportions of polymorphism loci and relatively high levels of genetic and epigenetic diversity in populations of Taiwanese lineages of Taiwania than those of Chinese and Vietnamese lineages. We also found there were significant levels of genetic and epigenetic differentiation between Taiwanese and mainland Asia lineages revealed by pairwise FST and STRUCTURE. Neighbor-Net phylogenetic network demonstrated the close relationship between Yunnan-Myanmar and Vietnamese lineages differentiated from Taiwanese populations, both results suggested Taiwania populations can be devided into two distinctly distinguishing clusters when K=2. In DFDIST and BAYESCAN analyses, no potential selective outliers were found within Taiwanese populations due to strongly neutral drift induced divergence. 32 loci were identified as potentially selective ouliers associated strongly with environmental variables including seasonary temperature (BIO4), seasonary precipitation (BIO15), and NDVI in either global or pairwise lineage comparisons. Not only had these important environmental variables influenced genetic and epigenetic variations, but they also had driven adaptive divergence and population differentiation of Taiwania with random drift. Therefore, the standing genetic and epigenetic variations either determined by selection or fixed by neutral drift, they both may ensure the capacity of Taiwania for adaptations to rapidly varing environment or climate change, and ability to survive and reproduce in the future.

朱麗萍、潘富俊、楊政川、王亞男（2006）依RAPD標誌推論臺灣杉屬之遺傳督樣性。林業研究季刊，28，15-28。臺中市：國立中興大學實驗林管理處。
范瑾瑜（2012）台灣油杉之族群遺傳及適應性演化。國立臺灣師範大學生命科學研究所碩士論文。
徐新武、郭幸榮、陳冠瑋、王智弘、田玉娟、林淑婷、梁亞忠、鍾年鈞、關秉宗（2007）臺灣杉芽生長之物候學。臺灣杉國際學術研討會論文集，157-169。國立臺灣大學實驗林管處。
張仲德、林登秋、林能暉（2009）臺灣地區酸沉降臨界負荷量及其對環境與經濟影響之估算。地理學報，56，39-58。臺北市：國立臺灣大學地理環境資源學系。
陳淑華、鍾年鈞、王亞男、李金玲、李宜玲、蔡摒芳（2007）台灣杉花粉之發育。臺灣杉國際學術研討會論文集，139-155。國立臺灣大學實驗林管處。
黃增泉、楊國禎、王俊能（2007）Taiwania之稱呼、原始林、遺傳變異與地理親緣關係。臺灣杉國際學術研討會論文集，81-96。國立臺灣大學實驗林管處。
蘇鴻傑（2007）台灣杉的前世今生：植群變遷與生活史。臺灣杉國際學術研討會論文集，97-117。國立臺灣大學實驗林管處。
Angers B, Castonguay E, Massicotte R (2010) Environmentally induced phenotypes and DNA methylation: how to deal with unpredictable conditions until the next generation and after. Molecular Ecology, 19, 1283-1295.
Antao T, Beaumont MA (2011) Mcheza: a workbench to detect selection using dominant markers. Bioinformatics, 15, 1717-1718.
Arthofer W (2010) tinyFLP and tinyCAT: software for automatic peak selection and scoring of AFLP data tables. Molecular Ecology Resource, 10, 385-388.
Barton K (2017) MuMIn: Multi-model inference. R package version 1.40.0. https://CRAN.R-project.org/package=MuMIn
Beaumont MA, Balding D (2004) Identifying adaptive genetic divergence among populations from genome scans. Molecular Ecology, 13, 969-980.
Bensch S, Åkesson M (2005) Ten years of AFLP in ecology and evolution: why so few animals? Molecular Ecology, 14, 2899-2914.
Blanchet FG, Legendre P, Borcard D (2008) Forward selection of explanatory variables. Ecology, 89, 2623-2632.
Bonin A, Ehrich D, Manel S (2007) Statistical analysis of amplified fragment length polymorphism data: a toolbox for molecular ecologists and evolutionists. Molecular Ecology, 16, 3737-3758.
Bossdorf O, Richards CL, Pigliucci M (2008) Epigenetics for ecologists. Ecolgoy Letters, 11, 106-115.
Borcard D, Legendre P, Drapeau P (1992) Partialling out the spatial component of ecological variation. Ecology, 73, 1045-1055.
Bryant D, Moulton V (2004) Neighbor-Net: an agglomerative method for the construction of phylogenetic network. Molecular Biology and Evolution, 21, 255-265.
Burnham KP, Anderson DR (2002) Model selection and multimodel inference: a practical information-theoretic approach. Springer, New York.
Chang CT, Wang SF, Vadeboncoeur MA, Lin TC (2014) Relating vegetation dynamics to temperature and precipitation at monthly and annual timescales in Taiwan using MODIS vegetation indices. International Journal of Remote Sensing, 35, 598-620.
Chen IC, Hill JK, Ohlemüller R, Roy DB, Thomas CD (2011) Rapid range shifts of species associated with high levels of climate warming. Science, 333, 1024.
Chou YW, Thomas PI, Ge XJ, LePage BA, Wang CN (2011) Refugia and phylogeography of Taiwania in East Asia. Journal of Biogeography, 38, 1992-2005.
Chung JD, Lin TP, Tan YC, Lin MY, Hwang SY (2004) Genetic diversity and biogeography of Cunninghamia konishii (Cupressaceae), an island species in Taiwan: a comparison with Cunninghamia lanceolata, a mainland species in China. Molecular Phylogenetics and Evolution, 33, 791-801.
Chung JD, Lin TP, Cheng YP, Hwang SY (2007) Phylogeographic study reveals the origin and evolutionary history of a Rhododendron species complex in Taiwan. Molecular Phylogenetics and Evolution, 42, 14-24.
Cokus SJ, Feng S, Zhang X et al. (2008) Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature, 452, 215.
Cubas P, Vincent C, Coen E (1999) An epigenetic mutation responsible for natural variation in floral symmetry. Nature, 401, 157-161.
Doyle JJ, Doyle JL (1987) A rapid DNA isolation procedure for small quantities of fresh lead tissue. Phytochemical Bulletin, 19, 11-15.
Dray S, Legendre P, Blanchet G (2016) packfor: forward selection with permutation (Canoco p.46). R package version 0.0-8/r136. Available online at: http:// R-Forge.R-project.org/projects/sedar/
Eckert CG, Samis KE, Lougheed SC (2008) Genetic variation across species’ geographic ranges: the central-marginal hypothesis and beyond. Molecular Ecology, 17, 1170-1188.
Ehrich D (2006) AFLPDAT: a collection of R functions for convenient handling of AFLP data. Molecular Ecology Resource, 6, 603-604.
Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analysis under linux and windows. Molecular Ecology Resource, 10, 564-567.
Evanno G, Regnaut S, Goudet J (2005) Detecting the number of cluster of individuals using the software STRUCTURE: a simulation study. Molecular Ecology, 14, 2611-2620.
Falush D, Stephens M, Pritchard JK et al. (2007) Inference of population structure: extensions to linked loci and correlated allele frequencies. Genetics, 164, 1567-1587.
Fang JY, Chung JD, Chiang YC, Chang CT, Chen CY, Hwang SY (2013) Divergent selection and local adaptation in disjunct populations of an endangered conifer, Keteleeria davidiana var. formosana (Pinaceae). PLoS ONE, 8, e70162.
Farjon A, Wang CN, Yang KC (2002) Extensions to the natural range of Taiwania cryptomerioides. Fitzroya, 5, 4-5.
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.
Francis RM (2017) pophelper: an R package and web app to analyse and visualise population structure. Molecular Ecology Resource, 17, 27-32.
Frankham R (1998) Inbreeding and extinction: Island populations. Conservation Biology, 12, 665-675.
Franks SJ (2010) Genetics, evolution, and conservation of island plants. Journal of Plant Biology, 53, 1-9.
Funk WC, Lovich RE, Hohenlohe PA et al. (2016) Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis). Molecular Ecology, 25, 2176-2194.
Hagmann J, Becker C, Müller J et al. (2015) Century-scale methylation stability in a recenty diverged Arabidopsis thaliana Lineage. PLoS Genetics, 11, e1004920.
Hayata B (1906) On Taiwania, a new genus of coniferae from the island of Formosa. Journal of the Linnean Society of London, Botany, 37, 330-331.
Hayata B (1907) On Taiwania and its affinity to other genera. Shokubutsugaku Zasshi, 21, 21-28.
Heber S, Varsani A, Kuhn S, Girg A, Kempenaers B, Briskie J (2013) The genetic rescue of two bottlenecked South Island robin populations using translocations of inbred donors. Proceedings of the Royal Society B-Biological Sciences, 280, 20122228.
Hedrick PW (2001) Conservation genetics: where are we now? Trends in Ecology & Evolution, 16, 629-636.
Herrera CM, Bazaga P (2010) Epigenetic differentiation and relationship to adaptive genetic divergence in discrete populations of the violet Viola cazorlensis. New Phytologist, 187, 867-876.
Herrera CM, Bazaga P (2016) Genetic and epigenetic divergence between disturbed and undisturbed suppopulations of a Mediterranean shrub: a 20-year field experiment. Ecology and Evolution, 6, 3832-3847.
Herrera CM, Medrano M, Bazaga P (2017) Comparative epigenetic and genetic spatial structure of the perennial herb Helleborus foetidus: Isolation by environment, isolation by distance, and functional trait divergence. American Jounal of Botany, 104, 1195-1204.
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surface for global land areas. International Journal of Climatology, 25, 1965-1978.
Holsinger KE, Lewis PO (2003) Hickory: a package for analysis of population genetic data v1.1 distributed by the authors, department of ecology and evolution biology, University of Connecticut, Storrs, USA.
Holsinger KE, Wallace LE (2004) Bayesian approaches for the analysis of population genetic structure: an example from Platanthera leucophaea (Orchidaceae). Molecular Ecology, 13, 887-894.
Hsieh YC, Chung JD, Wang CN, Chang CT, Chen CY, Hwang SY (2013) Historical connectivity, contemporary isolation and local adaptation in a widespread but discontinuously distributed species endemic to Taiwan, Rhododendron oldhamii (Ericaceae). Heredity, 111, 147-156.
Huang CL, Chang CT, Huang BH et al. (2015a) Genetic relationships and ecological divergence in Salix species and populations in Taiwan. Tree Genetics & Genomes, 11, 39.
Huang CL, Chen JH, Tsang MH, Chung JD, Chang CT, Hwang SY (2015b) Influences of environmental and spatial factors on genetic and epigenetic variations in Rhododendron oldhamii (Ericaceae). Tree Geneitcs & Genomics, 11, 823.
Huang SSF, Hwang SY, Lin TP (2002) Spatial pattern of chloroplast DNA variation of Cyclobalanopsis glauca in Taiwan and East Asia. Molecular Ecology, 11, 2349-2358.
Huang X, Li S, Ni P et al. (2017) Rapid reponse to changing environments during biological invasions: DNA methylation perspectives. Molecular Ecology, 26, 6621-6633.
Huson DH, Bryant D (2006) Application of phylogenetic network in evolutionary studies. Molecular Biology and Evolution, 23, 254-267.
Hutchison DW, Templeton AR (1999) Correlation of pairwise genetic and geographic distance measures: inferring the relative influences of gene flow and drift on the distribution of genetic variability. Evolution, 53, 1898-1914.
Hwang SY, Lin TP, Ma CS, Lin CL, Chung JD, Yang JC (2003) Postglacial population growth of Cunninghamia konishii (Cupressaceae) inferred from phylogeographical and mismatch analysis of chloroplast DNA variation. Molecular Ecology, 12, 2689-2695.
Jablonka E, Oborny B, Molnár I, Kisdi É, Hofbauer J, Czárán T (1995) The adaptive advantage of phenotypic memory in changing environments. The Philosophical Transactions of the Royal Society B, 350, 133-141.
Jakobsson M, Rosenberg NA (2007) CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics, 23, 1801-1806.
Janes JK, Miller JM, Dupuis JR et al. (2017) The K = 2 conundrum. Molecular Biology, 26, 3594-3602.
Jeffreys H (1961) Theory of probability (third edition). Oxford: Oxford University Press.
Jombart T, Devillard S, Balloux F (2010) Discriminant analysis of principle component: a new method for the analysis of genetically structured populations. BMC genetics, 11, 94.
Jombart T, Ahmed I (2011) adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics, 27, 3070-3071.
Jump AS, Huang TJ, Chou CH (2012) Rapid altitudinal migration of mountain plants in Taiwan and its implications for high altitude biodiversity. Ecography, 35, 204-210.
Kamvar ZN, Tabima JF, Grunwald NJ (2014) Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality. Frontiers in Genetics, 6, 208.
Keller I, Seehausen O (2012) Thermal adaptation and ecological speciation. Molecular Ecology, 21, 782-799.
Kilvitis HJ, Alvarez M, Foust CM, Schrey AW, Robertson M, Richards CL (2014) Ecological epigenetics. In: Ecological Genomics (eds Landry CR, Aubin-Horth N), pp. 191-210. Springer, New York.
Kimura M (1983) The Neutral Theory of Molecular Evolution. Cambridge University Press, Cambridge.
Kolbe JJ, Leal M, Schoener TW, Spiller DA, Losos JB (2012) Founder effects persist despite adaptive differentiation: a field experiment with lizards. Science, 335, 1086-1089.
Laporte V, Cuguen J, Couvet D (2000) Effective population sizes for cytoplasmic and nuclear genes in a gynodioecious species: the role of the sex determination system. Genetics, 154, 447-458.
Lenormand T (2002) Gene flow and limits to natural selection. Trends Ecology and Evolution, 17, 183-189.
Ley AC, Hardy OJ (2013) Improving AFLP analysis of large-scale patterns of genetic variation – a case study with Central African lianas Haumania spp (Marantaceae) showing interspecific gene flow. Molecular Ecology, 22, 1984-1997.
Li C-F, Chytrý M, Zelený D et al. (2013) Classification of Taiwan forest vegetation. Applied Vegetation Science, 16, 698-719.
Lin TP, Lu CS, Chung YL Yang JC (1993) Allozyme variation in four populations of Taiwania cryptomeroides in Taiwan. Silvae Genetica, 42, 278-284.
Lin TP, Chuang WJ, Huang SSF, Hwang SY (2003) Evidence for the existence of some dissociation in an otherwise strong linkage disequilibrium between mitochondrial and chloroplastic genomes in Cyclobalanopsis glauca. Molecular Ecology, 12, 2661-2668.
Liu TS, Su HJ (1983) Biosystematic studies on Taiwania and numerical evaluations of the systematics of Taxodiaceae. Taiwan Museum Special Publication No. 2.
Manel S, Poncet BN, Legendre P, Gugerli F, Holderegger R (2010) Common factors drive adaptive genetic variation at different spatial scales in Arabis alpna. Molecular Ecology, 19, 3824-3835.
Manel S, Gugerli F, Thuiller W et al. (2012) Broad-scale adative genetic variation in alpine plants is driven by temperature and precipitation. Molecular Ecology, 21, 3729-3738.
Médail F, Diadema K (2009) Glacial refugia influence plant diversity patterns in the Mediterranean Basin. Journal of Biogeography, 36, 1333-1345.
Nicotra AB, Segal DL, Hoyle GL, Schrey AW, Verhoeven KJ, Richards CL (2015) Adaptation plasticity and epigenetic variation in response to warming in an Alpine plant. Ecology and Evolution, 5, 634-647.
Nosil P, Funk DJ, Ortiz-Barrientos D (2009) Divergent selection and heterogeneous genomic divergence. Molecular Ecology, 18, 375-402.
Oksanen J, Blanchet G, Friendly M et al. (2018) vegan: Community Ecology Package. R package version 2.5-2. Available online at: https://CRAN.R-project.org/package=vegan
Palmer JD, Stein DB (1986) Conservation of chloroplast genome structure among vascular plants. Current Genetics, 10, 823-833.
Pérez-Figueroa A, Garcia-Pereira MJ, Saura M, Rolan-Alvarez E, Caballero A (2010) Comparing three different methods to detect selective loci using dominant markers. Journal of Evolutionary Biology, 23, 2267-2276.
Petit RJ, Aguinagalde I, de Beaulieu JL (2003) Glacial refugia: hotspots but not melting pots of genetic diversity. Science, 300, 1563-1565.
Pettorelli N, Vik JO, Mysterud A, Gaiilard J-M, Tucker CJ, Stenseth NC (2005) Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends in Ecology & Evolution, 20, 503-510.
Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics, 155, 945-959.
Rapp RA, Wendel JF (2005) Epigenetics and plant evolution. New phytologist. 168, 81-91.
Räsänen K, Hendry AP (2008) Disentangling interaction between adaptive divergence and gene flow when ecology drives diversification. Ecology Letters, 11, 624-636.
Richards CL, Alonso C, Becker C et al. (2017) Ecological plant epigenetics: Evidence from model and non-model species, and the way forward. Ecology Letters, 20, 1576-1590.
Schnitzer SA (2005) A mechanistic explanation for global patterns of Liana abundance and distribution. The American Naturalist, 166, 262-276.
Schrey AW, Alvarez M, Foust CM et al. (2013) Ecological epigenetics: Beyond MS-AFLP. Intergrative and Comparative Biology, 53, 340-350.
Schulz B, Eckstein RL, Durka W (2013) Scoring and analysis of methylation-sensitive amplification polymorphism for epigenetic population studies. Molecular Ecology Resources, 13, 642-653.
Schulz B, Eckstein RL, Durka W (2014) Epigenetic variation reflects dynamic habitat conditions in a rare floodplain herb. Molecular Ecology, 23, 3523-3537.
Sexton JP, Hangartner SB, Hoffmann AA (2014) Genetic isolation by environment or distance: which pattern of gene flow is most common? Evolution, 68, 1-15.
Slatkin M (1987) Gene flow and the geographic structure of natural populations. Science, 236, 787-792.
Slatkin M, Barton NH (1989) A comparison of three indirect methods for estimating average levels of gene flow. Evolution, 43, 1349-1368.
Stewart JR, Lister AM, Barnes I, Dalén L (2010) Refugia revisited: individualistic responses of species in space and time. Proceedings of the Royal Society B-Biological Sciences, 277, 661-671.
Storz JF (2005) Using genome scans of DNA polymorphism to infer adaptive population divergence. Molecular Ecology, 14, 671-688.
Stucki S, Orozco-terWengel P, Forester BR et al. (2017) High performance computation of landscape genomic models including local indicators of spatial association. Molecular Ecology Resource, 17, 1072-1089.
Thomas P, Farjon A (2011) Taiwania cryptomerioides. The IUCN Red List of Threatened Species 2011, e.T31255A9620141. available at: http://dx.doi.org/10.2305/IUCN.UK.2011-2.RLTS.T31255A9620141.en. Downloaded on 30 June 2018.
Thornthwaite CW (1948) An approach toward a rational classification of climate. Geographical Review, 38, 55-94.
Turelli M, Barton NH, Coyne JA (2001) Theory and speciation. Trends in Ecology & Evolution, 16, 330-343.
Uyeda JC, Arnold SJ, Hohenlohe PA, Mead L (2009) Drift promotes speciation by sexual selection. Evolution, 63, 583-594.
Vekemans X, Beauwens T, Lemaire M, Roldan-Ruiz I (2002) Data from amplified fragment length polymorphism (AFLP) markers show indication of size homoplasy and of a relationship between degree of homoplasy and fragment size. Molecular Ecology, 11, 139-151.
Verhoeven KJF, Jansen JJ, van Dijk PJ, Biere A (2010) Stress-induced DNA methylation changes and their heritability in asexual dandelions. New Phytologist, 185, 1108-1118.
Via S (2012) Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow. Philosophical Transactions of the Royal Society of London B: Biological Science, 367, 451-460.
Vos P, Hogers R, Bleeker M et al. (1995) AFLP : a new technique for DNA fingerprinting. Nucleic Acids Research, 23, 4407-4414.
Vuylsteke M, Peleman JD, van Eijk MJT (2007) AFLP technology for DNA fingerprinting. Nature Protocols, 2, 1387-1398.
Wang IJ (2013) Examining the full effects of landscape heterogeneity on spatial genetic variation: a multiple matrix regression approach for quantifying geographic and ecological isolation. Evolution, 67, 3403-3411.
Weigelt P, Jetz W, Kreft H (2013) Bioclimatic and physical characterization of the world’s islands. Proceedings of the National Academy of Sciences of the United States of America, 110, 15307-15312.
Whitlock MC, Mccauley DE (1999) Indirect measures of gene flow and migration: FST ≠ 1/(4Nm+1). Heredity, 82, 117-125.
Xiong LZ, Xu CG, Saghai Maroof MA, Zhang Q (1999) Patterns of cytosine methylation in an elite rice hybrid and parental lines, detected by a methylation-sensitive amplification polymorphism technique. Molecular and General Genetics, 261, 439-446.
Yan M, Xiong Y, Liu R, Deng M, Song J (2018) The application and limitation of universal chloroplast markers in discriminating East Asian evergreen oaks. Frontiers in Plant Science, 9, 569.
Zhivotovsky (1999) Estimating population structurre in diploids with multilocus dominant DNA markers. Molecular Ecology, 8, 907-913.