發生於末次冰河時期的氣候波動，被認為對於現今生物的族群變化以及遺傳多樣性的程度上，有很深遠的影響。直至目前為止，這類相關的研究大都以歐洲或北美洲的溫帶物種為主，相較之下，對於亞熱帶以及熱帶地區的生物如何受到冰河作用的影響，我們在這方面的知識是較為缺乏的。因此，本研究以亞熱帶地區的一種稜果榕授粉蜂(Ceratosolen sp. 1)為研究對象，藉由探討稜果榕授粉蜂的族群變動歷史、長期以及近期的基因交流，以瞭解最近的冰河作用如何影響一個亞熱帶物種的族群變化和當代的遺傳結構。我們使用來自於219個體的粒腺體COI基因序列片段 (1052 bp)以及398個體在15個微衛星體基因座的基因型等遺傳資料，進行分析。族群變動歷史的重建結果顯示，稜果榕授粉蜂曾經歷過一次顯著的族群成長，發生在距今大約兩萬七千年前；接著在兩萬一千年前左右，也就是最近的冰期高峰之後，此成長幅度開始劇烈地增加，直到距今約一萬年前，族群趨於穩定。此結果說明，稜果榕授粉蜂的族群變動歷史與發生在台灣的末次冰期氣候變動具有緊密的相連性，並且對於末次冰期如何影響亞熱帶物種的議題上，提供了一個明確的例子。在族群數量的變化大小方面，我們發現冰河時期前後，稜果榕授粉蜂的族群大小差異，將近一千倍左右。此結果顯示了在最近一次冰河時期發生時，稜果榕授粉蜂與其宿主的族群是高度被壓縮以及破碎化的。冰河時期後的族群數量與分布空間的擴張也可能是造成稜果榕授粉蜂族群間並不存在遺傳分化(FST < 0.02)的部分原因。根據現生個體的來源族群推測，發現稜果榕授粉蜂具有長距離擴散的可能性以及族群間低比例的第一代遷入者。綜合我們的結果推測，現今觀察到稜果榕授粉蜂的族群結構是同時受到最近一次冰河時期結束之後氣候變動以及授粉蜂擴散能力的影響。除此之外，本研究亦觀察到稜果榕擁有多種授粉蜂，而關於此現象在演化上的意涵，也於此篇報告中一併陳述。

The climate oscillations during the last glacial period have profound influences on the demography and levels of genetic diversity of species. Molecular evidences of glacial effects on temperate species in Europe and North America have been well documented, whereas little is known regarding that on subtropical and tropical species. In this study, we investigated the demographic history, long term and ongoing gene flow of a subtropical fig-pollinating wasp (Ceratosolen sp. 1) of Ficus septicas in Taiwan and provided insights into how its demography and contemporary genetic structure were shaped by climatic changes since the most recent ice age. The analysis was conducted based on mitochondrial sequences of COI gene (1052 bp) from 219 individuals and genotypes of 398 individuals at 15 microsatellite loci. Reconstructed historical demography suggested that the population of C. sp. 1 has undergone a drastic population expansion that began around 27,000 years ago and then grew rapidly around 21,000 years ago, right after the last glacial maximum (LGM) in Taiwan; finally, the population of C. sp. 1 reached stability approximately 10,000 years ago. It revealed that the demographic history of C. sp. 1 was tightly coupled with the climate change since the LGM in Taiwan and provided an unambiguous case for the effects of the last glacial cycle on the demography of a subtropical species. The magnitude of the population size change was huge, approximately 1000 times larger than the size before expansion, which indicated that the population sizes of C. sp. 1 and its host might have been highly compressed and fragmented during the LGM. The postglacial demographic and spatial expansion might partially contribute to the lack of significant genetic differentiation among populations of C. sp. 1 as revealed by the FST statistics (< 0.02) and mismatch distribution of mithochondrial sequences, which implied high level of long term gene flow. Results of population assignment test revealed the potential of long distance dispersal and low frequency of first generation migrants among local populations as well. Our findings clearly demonstrated how recent climate changes since the LGM and dispersal ability of an organism shaped the genetic composition of a subtropical fig-pollinating wasp. In addition, we found multiple pollinators associated with F. septica in Taiwan, and the evolutionary implications of such phenomenon were addressed.

Bossart JL, Prowell DP (1998) Genetic estimates of population structure and gene flow: limitations, lessons and new directions. Trends in Ecology and Evolution, 13, 202-206.
Brown JH, Lomolino MV (1998) Biogeography. Sinauer Associates, Sunderland, MA.
Burg TM, Gaston AJ, Winker K, Friesen VL (2006) Effects of Pleistocene glaciations on population structure of North American chestnut-backed chickadees. Molecular Ecology, 15, 2409-2419.
Chen SF, Rossiter SJ, Faulkes CG (2006) Population genetic structure and demographic history of the endemic Formosan lesser horseshoe bat (Rhinolophus monoceros). Molecular Ecology, 15, 1643-1656.
Chen CH, Chou LY (1997) The blastophagini of Taiwan (Hymenoptera: Agaonidae: Agaoninae). Journal of Taiwan Museum, 50, 113-154.
Cheng YP, Hwang SY, Lin TP (2005) Potential refugia in Taiwan revealed by the phylogeographical study of Castanopsis carlesii Hayata (Fagaceae). Molecular Ecology, 14, 2075-2085.
Drummond AJ & Rambaut A (2003a) BEAST version 1.0. Available from http://evolve.zoo.ox.ac.uk/beast/.
Drummond AJ & Rambaut A (2003b) TRACER version 1.0. Available from http://evolve.zoo.ox.ac.uk/beast/.
Drummond AJ, Rambaut A, Shapiro B, Pybus OG. (2005) Bayesian coalescent inference of past population dynamics from molecular sequences. Molecular Biology and Evolution, 22, 1185-1192.
Excoffier L (2004) Patterns of DNA sequence diversity and genetic structure after a range expansion: lessons from the infinite-island model. Molecular Ecology, 13, 853-864.
Excoffier L, Laval G, Schneider S (2005) Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online, 1, 47-50.
Excoffier L, Smouse PE (1994) Using allele frequencies and geographic subdivision to reconstruct gene trees within a species: molecular variance parsimony. Genetics, 136, 343-359.
Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics, 147, 915-925.
Fu YX, Li WH (1993) Statistical tests of neutrality of mutations. Genetics, 133, 693-709.
Galil J, Eisikowtich D (1968) On the pollination of Ficus sycomorus in East Africa. Ecology, 49, 259-269.
Gemmell NJ, and Akiyama S (1996) An efficient method for the extraction of DNA from vertebrate tissues. Trends in Genetics, 12, 338-339.
Grandi G (1927) Hyménoptères sycophiles récoltés aux iles Philippines par C.F. Baker, i. Agaonini. Philippine Journal of Science, 33, 309-329.
Goudet J (2001) FSTAT. A program to estimate and test gene diversities and fixation indices (version 2.9.3). Available from http://www.unil.ch/izea/softwares/fstat.html. Updated from Goudet (1995).
Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology, 52, 696-704.
Haine ER, Martin J, Cook JM (2006) Deep mtDNA divergences indicate cryptic species in a fig-pollinating wasp. BMC Evolutionary Biology, 6, 83.
Halliburton R (2004) Introduction to population genetics. Prentice Hall, Upper Saddle River, NJ.
Harrison RD (2003) Fig wasp dispersal and the stability of a keystone plant resource in Borneo. Proceedings of the Royal Society of London. Series B, Biological Sciences, 270, S76-S79.
Harrison RD, Rasplus JY (2006) Dispersal of fig pollinators in Asian tropical rain forests. Journal of Tropical Ecolgy, 22, 631-639.
Heilveil JS, Berlocher SH (2006) Phylogeography of postglacial range expansion in Nigronia serricornis Say (Megalopeta: Corydalidae) Molecular Ecology, 15, 1627-1641.
Herre EA, Machado CA, Bermingham E, Nason JD, Windsor DM, McCafferty SS, van Houten W, Bachmann K (1996) Molecular phylogenies of figs and their pollinator wasps. Journal of Biogeography, 23, 521-530.
Hewitt GM (1996) Some genetic consequences of ice ages, and their role in divergence and speciation. Biological Journal of the Linnean Society, 58, 247-276.
Hewitt GM (2000) The genetic legacy of the Quaternary ice ages. Nature, 405, 907-913.
Hewitt GM (2004) Genetic consequences of climatic oscillation in the Quaternary. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, 359, 183-195.
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.
Kerdelhue C, Hochberg ME, Rasplus JY. (1997) Active pollination of Ficus sur by two sympatric fig wasp species in West Africa. Biotropica, 29, 69-75.
Kumar S, Tamura K, Nei M (2004) MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Briefings in Bioinformatics, 5, 150-163.
Li WH (1997) Molecular Evolution. Sinauer Associates, Sunderland, MA.
Liew PM, Chung NJ (2001) Vertical migration of forests during the last glacial period in subtropical Taiwan. Western Pacific Earth Sciences, 1, 405-414.
Machado CA, Jousselin E, Kjellberg F, Compton SG, Herre EA (2001) Phylogenetic relationships, historical biogeography and character evolution of fig-pollinating wasps. Proceedings of the Royal Society of London. Series B, Biological Sciences, 268, 685-694.
Mantel N (1967) Detection of disease clustering and a generalized regression approach. Cancer Research, 27, 209-220.
Marshall TC, Slate J, Kruuk LEB, Pemberton JM (1998) Statistical confidence for likelihood-base paternity inference in natural populations. Molecular Ecology, 7, 639-655.
Martinez-Solano I, Teixeira J, Buckley D, Garcia-Paris M (2006) Mitochondrial DNA phylogeography of Lissotriton boscai (Caudata, Salamandridae): evidence for old, multiple refugia in an Iberian endemic. Molecular Ecology, 15, 3375-3388.
Michaloud G, Carriére S, Kobbi M (1996) Exceptions to the one:one relationship between African fig trees and their fig wasp pollinators: possible evolutionary scenarios. Journal of Biogeography, 23, 513-520.
Molbo D, Machado CA, Sevenster JG, Keller L, Herre EA (2003) Cryptic species of fig-pollinating wasps: implications for the evolution of the fig-wasp mutualism, sex allocation, and precision of adaptation. Proceedings of National Academy of Sciences of the United States of America, 100, 5867-5872.
Molbo D, Machado CA, Herre EA, Keller L (2004) Inbreeding and population structure in two pairs of cryptic wasp species. Molecular Ecology, 13, 1613-1623.
Nason, JD, Herre EA, Hamrick JL (1998) The breeding structure of a tropical keystone plant resource. Nature, 391, 685-687.
Noonan BP, Gaucher P (2005) Phylogeography and demography of Guianan harlequin toads (Atelopus): diversification within a refuge. Molecular Ecology, 14, 3017-3031.
Paetkau D, Slade R, Burden M, Estoup A (2004) Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Molecular Ecology, 13, 55-65.
Pauls SU, Lubsch HT, Haase P (2006) Phylogeography of the montane caddisfly Drusus discolor evidence for multiple refugia and periglacial survival. Molecular Ecology, 15, 2153-2169.
Pearse DE, Crandall KA (2004) Beyond FST: analysis of population genetic data for conservation. Conservation Genetics, 5, 585-602.
Pinot S, Ramstein G, Harrison SP, Prentice IC, Guiot J, Stute M, Joussaume S. (1999) Tropical paleoclimates at the last glacial maximum: comparison of Paleoclimate Modeling Intercomparison Project (PMIP) simulations and paleodata. Climate Dynamics, 15, 857-874.
Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A (2004) GENECLASS2: A software for genetic assignment and first-generation migrant detection. Journal of Heredity, 95, 536-539.
Posada D, Crandall KA (1998) MODELTEST: testing the model of DNA substitution. Bioinformatics, 14, 817-818.
Rannala B, Mountain JL (1997) Detecting immigration by using multilocus genotypes. Proceedings of National Academy of Sciences of the United States of America, 94, 9197-9201.
Rasplus JY (1994) The one-to-one specificity of the Ficus-Agaoninae mutualism: how causal? In The Biodiversity of African Plants, ed. LJG van der Maesen, XM van der Burgt, JM van Medenbach de Rooy, pp639-649. Dordrecht: Kluwer Academic.
Raymond M, Rousset F (1995) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. Journal of Heredity, 86, 248-249.
Rogers AR (1995) Genetic evidence for a Pleistocene population explosion. Evolution, 49, 608-615.
Rogers AR, Harpending H (1992) Population growth makes waves in the distribution of pairwise genetic differences. Molecular Biology and Evolution, 9, 552-569.
Rohlf FJ (1973) Algorithm 76: Hierarchical clustering using the minimum spanning tree. Computer Journal, 16, 93-95.
Rousser F (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics, 145, 1219-1228.
Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R (2003) DnaSP, DNA polymorphism analysis by the coalescent and other methods. Bioinformatics, 19, 2496-2497.
Runck AM, Cook JA (2005) Postglacial expansion of the southern red-backed vole (Clethrionomys gapperi) in North America. Molecular Ecology, 14, 1445-1456.
Schmitt T, Habel JC, Zimmermann M, Müller P (2006) Genetic differentiation of the marbled white butterfly, Melanargia galathea, accounts for glacial distribution patterns and postglacial range expansion in southeastern Europe. Molecular Ecology, 15, 1889-1901.
Schneider S, Excoffier L (1999) Estimation of demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: Application to human mitochondrial DNA. Genetics, 152, 1079-1089.
Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P (1994) Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Annals of the Entomological Society of America, 87, 651-701.
Slatkin M, Hudson RR (1991) Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics, 129, 555-562.
Swofford DL (2002) PAUP*: Phylogenetic Analysis Using Parsimony (* and Other Methods). 4.0b10. Sinauer, Sunderland, MA.
Templeton AR, Routman E, Phillips CA (1995) Separating population structure from population history: a cladistic analysis of the geographical distribution of mitochondrial DNA haplotypes in the tiger salamander, Ambystoma trigrinum. Genetics, 140, 767-782.
Tsukada M (1966) Late Pleistocene vegetation and climate in Taiwan (Formosa). Proceedings of National Academy of Sciences of the United States of America, 55, 543-584.
Tsukada M (1967) Vegetation in subtropical Formosa during the Pleistocene glactations and the Holocene. Palaeogeography, Palaeoclimatology, Palaeoecology, 3, 49-64.
Tzeng HY (2004) Taxonomic study of the genus Ficus in Taiwan. Doctor’s Dissertation, National Chung Hsing University, Taichung, Taiwan. (in Chinese).
Van Noort S, Compton SG (1996) Convergent evolution of agaonine and sycoecine (Agaonidae, Chalcidoidae) head shape in response to the constraints of host fig morphology. Journal of Biogeography, 23, 415-424.
van Vuuran GJJ, Kryger P, Greeff JM (2006) Isolation of six microsatellite loci in the pollinating fig wasp, Platyscapa awekei. Molecular Ecology Notes 6: 385-386.
Ware AB, Compton SG (1992) Breakdown of pollinator specificity in an African Fig tree. Biotropica, 24, 544-549.
Ware AB, Compton SG (1994a) Dispersal of adult female fig wasps. 1. Arrivals and departures. Entomologia Experimentalis et Applicata, 73, 221-229.
Ware AB, Compton SG (1994b) Dispersal of adult female fig wasps. 2. Movements between trees. Entomologia Experimentalis et Applicata, 73, 231-238.
Weiblen GD (2001) Phylogenetic relationships of fig wasps pollinating functionally dioecious Ficus based on mitochondrial DNA sequences and morphology. Systematic Biology, 50, 243-267.
Weiblen GD (2002) How to be a fig wasp. Annual Review of Entomology, 47, 299-230.
Weiblen GD, Bush GL (2002) Speciation in fig pollinators and parasites. Molecular Ecology, 11, 1573-1578.
Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution, 38, 1358-1370
Wiebes JT (1979) Co-evolution of figs and their insect pollinators. Annual Review of Ecology and Systematics, 10, 1-12.
Yokoyama J (2003) Cospeciation of figs and fig-wasps: a case study of endemic species pairs in the Ogasawara Islands. Population Ecology, 45, 249-256.
Zane L, Marcato S, Bargelloni L, Bortolotto E, Papetti C, Simonato M, Varotto V, Patarnello T (2006) Demographic history and population structure of the Antarctic silverfish Pleuragramma antarcticum. Molecular Ecology, 15, 4499-4511.
Zavodna M, Arens P, van Dijk PJ, Partomihardjo T, Vosman B, van Damme JMM. (2005a) Pollinating fig wasps: genetic consequences of island recolonization. Journal of Evolutionary Biology, 18, 1234-1243.
Zavodna M, Compton SG, Raja S, Gilmartin PM, van Damme JMM (2005b) Do fig wasps produce mixed paternity clutches? Journal of Insect Behavior, 18, 351-362.