Genome-wide patterns of genetic diversity, population structure and demographic history in mānuka (Leptospermum scoparium) growing on indigenous Māori land

Abstract

Leptospermum scoparium J. R. Forst et G. Forst, known as mānuka by Māori, the indigenous people of Aotearoa (New Zealand), is a culturally and economically significant shrub species, native to New Zealand and Australia. Chemical, morphological and phylogenetic studies have indicated geographical variation of mānuka across its range in New Zealand, and genetic differentiation between New Zealand and Australia. We used pooled whole genome re-sequencing of 76 L. scoparium and outgroup populations from New Zealand and Australia to compile a dataset totalling ∼2.5 million SNPs. We explored the genetic structure and relatedness of L. scoparium across New Zealand, and between populations in New Zealand and Australia, as well as the complex demographic history of this species. Our population genomic investigation suggests there are five geographically distinct mānuka gene pools within New Zealand, with evidence of gene flow occurring between these pools. Demographic modelling suggests three of these gene pools have undergone expansion events, whilst the evolutionary histories of the remaining two have been subjected to contractions. Furthermore, mānuka populations in New Zealand are genetically distinct from populations in Australia, with coalescent modelling suggesting these two clades diverged ∼9–12 million years ago. We discuss the evolutionary history of this species and the benefits of using pool-seq for such studies. Our research will support the management and conservation of mānuka by landowners, particularly Māori, and the development of a provenance story for the branding of mānuka based products.

References

[1] Gudgeon WE . Maori religion. J Polyn Soc. 1905; 14: 107-30.

[2] Grey G . Maori Lore: The Traditions of the Maori People, with the more Important of their Legends . By authority, J. Mackay. South Africa: Government Printer, 1904.

[3] Riley M. Viking Sevenseas NZ, 1994.

[4] Best E. Maori medical lore. Notes on sickness and disease among the Maori people of New Zealand, and their treatment of the sick; together with some account of various beliefs, superstitions and rites pertaining to sickness, and the treatment thereof, as collected from the Tuhoe tribe. Part I. J Polyn Soc. 1904; 13: 213-37.

[5] Crowe A . A field guide to the native edible plants of New Zealand, 3rd edn . Auckland: Godwit Publishing Limited, 2001.

[6] Colenso W. On the geographical and economic botany of the North Island of New Zealand. Trans N Z Inst. 1868; 1: 1-58.

[7] Beattie H. Nature-Lore of the Southern Maori . New Zealand: New Zealand Institute; 1920.

[8] Best E . Fishing Methods and Devices of the Maori . New Zealand: Ams Pr Inc; 1929.

[9] Hiroa TR . The Coming of the Maori. Maori Purposes Fund Board . Whitcombe & Tombs: Christchurch, New Zealand; 1949.

[10] Palmer J. The Maori Kotaha. J Polyn Soc. 1957; 66: 175-91.

[11] Stephens J, Molan PC, Clarkson BD . A review of Leptospermum scoparium (Myrtaceae) in New Zealand. N Z J Bot. 2005; 43: 431-49.

[12] Ministry for Primary Industries . Apiculture: Ministry for Primary Industries 2018 Apiculture Monitoring Programme . New Zealand: Wellington; 2018.

[13] Thompson J . Revision of the genus leptospermum (Myrtaceae). Telopea. 1989; 3: 301-449.

[14] Bean A . The genus leptospermum Forst. Et Forst. F.(Myrtaceae) in northern Australia and Malesia. Austrobaileya. 1992; 3: 643-59.

[15] Bean A . Three new species of leptospermum (Myrtaceae) from Queensland and northern New South Wales. Telopea. 2004; 10: 831-8.

[16] Lyne A. Leptospermum namadgiensis (Myrtaceae), a new species from the Australian Capital Territory-New South Wales border area. Telopea. 1993; 5: 319-24.

[17] Lyne A, Crisp M . Leptospermum jingera (Myrtaceae-Leptospermoideae): a new species from north-eastern Victoria. Aust Syst Bot. 1996; 9: 301-6.

[18] Cheeseman TF . VI. The flora of Rarotonga, the chief island of the cook group. In: Transactions of the Linnean Society of London. 2nd series. Botany.Vol. 6. 1925, 261-313.

[19] Allan HH . Zealand. Flora of New Zealand. Volume I . New Zealand: Manaaki Whenua Press, 1961.

[20] Ronghua Y, Mark A, Wilson J . Aspects of the ecology of the indigenous shrub Leptospermum scoparium (Myrtaceae) in New Zealand. N Z J Bot. 1984; 22: 483-507.

[21] Douglas MH, van Klink JW, Smallfield BM . et al. Essential oils from New Zealand manuka: triketone and other chemotypes of Leptospermum scoparium . Phytochemistry. 2004; 65: 1255-64.

[22] Perry NB, Brennan NJ, van Klink JW . et al. Essential oils from New Zealand manuka and kanuka: chemotaxonomy of leptospermum . Phytochemistry. 1997; 44: 1485-94.

[23] Porter NG, Wilkins AL . Chemical, physical and antimicrobial properties of essential oils of Leptospermum scoparium and Kunzea ericoides. Phytochemistry. 1999; 50: 407-15.

[24] Adams CJ, Boult CH, Deadman BJ . et al. Isolation by HPLC and characterisation of the bioactive fraction of New Zealand manuka (Leptospermum scoparium) honey . Carbohydr Res. 2008; 343: 651-9.

[25] Adams CJ, Manley-Harris M, Molan PC . The origin of methylglyoxal in New Zealand manuka (Leptospermum scoparium) honey. Carbohydr Res. 2009; 344: 1050-3.

[26] Mavric E, Wittmann S, Barth G . et al. Identification and quantification of methylglyoxal as the dominant antibacterial constituent of Manuka (Leptospermum scoparium) honeys from New Zealand . Mol Nutr Food Res. 2008; 52: 483-9.

[27] Atrott J, Haberlau S, Henle T . Studies on the formation of methylglyoxal from dihydroxyacetone in Manuka (Leptospermum scoparium) honey. Carbohydr Res. 2012; 361: 7-11.

[28] Williams S, King J, Revell M . et al. Regional, annual, and individual variations in the dihydroxyacetone content of the nectar of manuka (Leptospermum scoparium) in New Zealand . J Agric Food Chem. 2014; 62: 10332-40.

[29] Stephens JMC . The Factors Responsible for the Varying Levels of UMF^® in mānuka (Leptospermum scoparium) Honey . New Zealand: The University of Waikato; 2006.

[30] Buys MH, Winkworth RC, de Lange PJ . et al. The phylogenomics of diversification on an island: applying anchored hybrid enrichment to New Zealand Leptospermum scoparium (Myrtaceae) . Bot J Linn Soc. 2019; 191: 1-17.

[31] Lamb S, Mortimer N, Smith E . et al. Focusing of relative plate motion at a continental transform fault: Cenozoic dextral displacement >700 km on New Zealand’s Alpine Fault, reversing >225 km of Late cretaceous sinistral motion . Geochem Geophys Geosyst. 2016; 17: 1197-213.

[32] Thornhill AH, Ho SY, Külheim C . et al. Interpreting the modern distribution of Myrtaceae using a dated molecular phylogeny . Mol Phylogenet Evol. 2015; 93: 29-43.

[33] Kamp P, Late J . Cretaceous-Cenozoic tectonic development of the Southwest Pacific region. Tectonophysics. 1986; 121: 225-51.

[34] Trewick S, Bland K . Fire and slice: palaeogeography for biogeography at New Zealand’s North Island/south island juncture. J R Soc N Z. 2012; 42: 153-83.

[35] McCulloch GA, Wallis GP, Waters JM . Onset of glaciation drove simultaneous vicariant isolation of alpine insects in New Zealand. Evolution. 2010; 64: 2033-43.

[36] Patterson MO, McKay R, Naish T . et al. Orbital forcing of the East Antarctic ice sheet during the Pliocene and Early Pleistocene . Nat Geosci. 2014; 7: 841-7.

[37] Mimee B, Duceppe M-O, Véronneau P-Y . et al. A new method for studying population genetics of cyst nematodes based on Pool-Seq and genomewide allele frequency analysis . Mol Ecol Resour. 2015; 15: 1356-65.

[38] Fischer MC, Rellstab C, Tedder A . et al. Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps . Mol Ecol. 2013; 22: 5594-607.

[39] Schlötterer C, Tobler R, Kofler R . et al. Sequencing pools of individuals-mining genome-wide polymorphism data without big funding . Nat Rev Genet. 2014; 15: 749.

[40] Dennenmoser S, Vamosi SM, Nolte AW . et al. Adaptive genomic divergence under high gene flow between freshwater and brackish-water ecotypes of prickly sculpin (Cottus asper) revealed by Pool-Seq . Mol Ecol. 2017; 26: 25-42.

[41] Roda F, Ambrose L, Walter GM . et al. Genomic evidence for the parallel evolution of coastal forms in the Senecio lautus complex . Mol Ecol. 2013; 22: 2941-52.

[42] Futschik A, Schlötterer C . The next generation of molecular markers from massively parallel sequencing of pooled DNA samples. Genetics. 2010; 186: 207-18.

[43] Fracassetti M, Griffin PC, Willi Y . Validation of pooled whole-genome re-sequencing in Arabidopsis lyrata . PLoS One. 2015; 10: e0140462.

[44] Lynch M, Bost D, Wilson S . et al. Population-genetic inference from pooled-sequencing data . Genome Biol Evol. 2014; 6: 1210-8.

[45] Corander J, Majander KK, Cheng L . et al. High degree of cryptic population differentiation in the Baltic Sea herring Clupea harengus . Mol Ecol. 2013; 22: 2931-40.

[46] Ferretti L, Ramos-Onsins SE, Pérez-Enciso M . Population genomics from pool sequencing. Mol Ecol. 2013; 22: 5561-76.

[47] Endler L, Betancourt AJ, Nolte V . et al. Reconciling differences in pool-GWAS between populations: a case study of female abdominal pigmentation in Drosophila melanogaster . Genetics. 2016; 202: 843-55.

[48] Bertelsen HP, Gregersen VR, Poulsen N . et al. Detection of genetic variation affecting milk coagulation properties in Danish Holstein dairy cattle by analyses of pooled whole-genome sequences from phenotypically extreme samples (pool-seq) . J Anim Sci. 2016; 94: 1365-76.

[49] Turner TL, Bourne EC, Von Wettberg EJ . et al. Population resequencing reveals local adaptation of Arabidopsis lyrata to serpentine soils . Nat Genet. 2010; 42: 260.

[50] Christe C, Stölting KN, Paris M . et al. Adaptive evolution and segregating load contribute to the genomic landscape of divergence in two tree species connected by episodic gene flow . Mol Ecol. 2017; 26: 59-76.

[51] Leblois R, Gautier M, Rohfritsch A . et al. Deciphering the demographic history of allochronic differentiation in the pine processionary moth Thaumetopoea pityocampa . Mol Ecol. 2018; 27: 264-78.

[52] Rellstab C, Zoller S, Tedder. A et al. Validation of SNP allele frequencies determined by pooled next-generation sequencing in natural populations of a non-model plant species. PLoS One. 2013; 8: e80422.

[53] Anderson EC, Skaug HJ, Barshis DJ . Next-generation sequencing for molecular ecology: a caveat regarding pooled samples. Mol Ecol. 2014; 23: 502-12.

[54] Bhatia D, Wing R, Singh K . Genotyping by sequencing, its implications and benefits. Crop Improv. 2013; 40: 101-11.

[55] Thrimawithana AH, Jones D, Hilario E . et al. A whole genome assembly of Leptospermum scoparium (Myrtaceae) for mānuka research . N Z J Crop Hortic Sci. 2019; 47: 1-28.

[56] Tribunal W. Ko Aotearoa tēnei: A Report into Claims Concerning New Zealand Law and Policy Affecting Māori Culture and Identity: Te Taumata Tuatahi (Wai 262) . Wellington, New Zealand: Waitangi Tribunal Report; 2011.

[57] Morgan ER, Perry NB, Chagné D . Science at the intersection of cultures-Māori, Pākehā and mānuka. N Z J Crop Hortic Sci. 2019; 47: 225-32.

[58] Hudson M, Garrison NA, Sterling R . et al. Rights, interests and expectations: Indigenous perspectives on unrestricted access to genomic data . Nat Rev Genet. 2020; 21: 377-84.

[59] Jombart T, Devillard S, Balloux F . Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010; 11: 94.

[60] Hivert V, Leblois R, Petit EJ. et al. Measuring genetic differentiation from Pool-seq data. Genetics. 2018; 210: 315-30.

[61] Huson DH, Bryant D . Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2005; 23: 254-67.

[62] Pickrell JK, Pritchard JK . Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012; 8: 1-17, e1002967.

[63] Cutler DJ, Jensen JD . To pool, or not to pool? Genetics . 2010; 186: 41-3.

[64] Esler AE, Astridge SJ . Tea tree (leptospermum) communities of the Waitakere range, Auckland New Zealand. N Z J Bot. 1974; 12: 485-501.

[65] Heenan PB, McGlone MS . Cenozoic formation and colonisation history of the New Zealand vascular flora based on molecular clock dating of the plastid rbcL gene. N Z J Bot. 2019; 57: 204-26.

[66] Prebble JG, Reichgelt T, Mildenhall DC . et al. Terrestrial climate evolution in the Southwest Pacific over the past 30 million years . Earth Planet Sci Lett. 2017; 459: 136-44.

[67] Kamp PJ, Green PF, White SH . Fission track analysis reveals character of collisional tectonics in New Zealand. Tectonics. 1989; 8: 169-95.

[68] McGlone M. Plant biogeography and the late Cenozoic history of New Zealand. N Z J Bot. 1985; 23: 723-49.

[69] Alloway BV, Lowe DJ, Barrell DJA . et al. Towards a climate event stratigraphy for New Zealand over the past 30 000 years (NZ-INTIMATE project) . J Quat Sci. 2007; 22: 9-35.

[70] McGlone MS, Newnham RM, Moar NT . The vegetation cover of New Zealand during the last glacial maximum: do pollen records under-represent woody vegetation. Terra Australis. 2010; 32: 49-68.

[71] Kitzberger T, Perry GLW, Paritsis J . et al. Fire-vegetation feedbacks and alternative states: common mechanisms of temperate forest vulnerability to fire in southern South America and New Zealand . N Z J Bot. 2016; 54: 247-72.

[72] Wyse SV, Wilmshurst JM, Burns BR . et al. New Zealand forest dynamics: a review of past and present vegetation responses to disturbance, and development of conceptual forest models . N Z J Ecol. 2018; 42: 87-106.

[73] Battersby PF, Wilmshurst JM, Curran TJ . et al. Exploring fire adaptation in a land with little fire: serotiny in Leptospermum scoparium (Myrtaceae) . J Biogeogr. 2017; 44: 1306-18.

[74] Clark PU, Dyke AS, Shakun JD . et al. The last glacial maximum . Science. 2009; 325: 710-4.

[75] McGlone M, Duncan R, Heenan P . Endemism, species selection and the origin and distribution of the vascular plant flora of New Zealand. J Biogeogr. 2001; 28: 199-216.

[76] Lockhart PJ, McLenachan PA, Havell D . et al. Phylogeny, radiation, and transoceanic dispersal of New Zealand alpine buttercups: molecular evidence under split decomposition . Ann Missouri Bot Gard. 2001; 88: 458-77.

[77] Winkworth RC, Wagstaff SJ, Glenny D et al. Plant dispersal news from New Zealand. Trends Ecol Evol. 2002; 17: 514-20.

[78] Knapp M, Stöckler K, Havell D . et al. Relaxed molecular clock provides evidence for long-distance dispersal of Nothofagus (southern beech) . PLoS Biol. 2005; 3: e14.

[79] Breitwieser I, Glenny DS, Thorne A . et al. Phylogenetic relationships in Australasian Gnaphalieae (Compositae) inferred from ITS sequences . N Z J Bot. 1999; 37: 399-412.

[80] Goldberg J, Morgan-Richards M, Trewick SA . Intercontinental island hopping: colonization and speciation of the grasshopper genus Phaulacridium (Orthoptera: Acrididae) in Australasia. Zool Anz. 2015; 255: 71-9.

[81] Goldberg J, Trewick SA, Paterson AM . Evolution of New Zealand’s terrestrial fauna: a review of molecular evidence. Philos Trans R Soc B: Biol Sci. 2008; 363: 3319-34.

[82] Koot EM, Morgan-Richards M, Trewick SA . An alpine grasshopper radiation older than the mountains, on Kā Tiritiri o te Moana (Southern Alps) of Aotearoa (New Zealand). Mol Phylogenet Evol. 2020; 147: 1-19, 106783.

[83] Sanmartín I, Wanntorp L, Winkworth RC . West wind drift revisited: testing for directional dispersal in the southern hemisphere using event-based tree fitting. J Biogeogr. 2007; 34: 398-416.

[84] Sanmartíin I, Ronquist F . Southern hemisphere biogeography inferred by event-based models: plant versus animal patterns. Sys Biol. 2004; 53: 216-43.

[85] Gardiner S, Bassett HC, Noiton DA . et al. A detailed linkage map around an apple scab resistance gene demonstrates that two disease resistance classes both carry the Vf gene . Theor Appl Genet. 1996; 93: 485-93.

[86] Tung P-Y, Blischak JD, Hsiao CJ . et al. Batch effects and the effective design of single-cell gene expression studies . Sci Rep. 2017; 7: 39921.

[87] Andrews S. FastQC: A Quality Control Tool for High Throughput Sequence Data [Online]. 2016. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.

[88] Langmead B, Salzberg SL . Fast gapped-read alignment with bowtie 2. Nat Methods. 2012; 9: 357.

[89] Li H, Handsaker B, Wysoker A . et al. The sequence alignment/map format and SAMtools . Bioinformatics. 2009; 25: 2078-9.

[90] Li H . A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011; 27: 2987-93.

[91] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[92] R Core Team . R: A language and environment for statistical computing. 2013.

[93] Knaus BJ, Grünwald NJ . vcfr: a package to manipulate and visualize variant call format data in R. Mol Ecol Resour. 2017; 17: 44-53.

[94] Jombart T . Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008; 24: 1403-5.

[95] Dray S, Dufour A-B . The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007; 22: 1-20.

[96] Excoffier L, Dupanloup I, Huerta-Sánchez E et al. Robust demographic inference from genomic and SNP data. PLoS Genet. 2013; 9: 1-17, e1003905.

[97] DeGiorgio M, Huber CD, Hubisz MJ et al. SweepFinder2: increased sensitivity, robustness and flexibility . Bioinformatics. 2016; 32: 1895-7.