Standard Paper Resolving the phylogenetic relationship between Parmotrema crinitum and Parmotrema perlatum populations Ayoub Stelate1,2,4 , Ruth Del-Prado2, David Alors3 , Hikmat Tahiri4, Pradeep K. Divakar2 and Ana Crespo2 1Department of Experimental Plant Biology, Charles University Faculty of Science, Viničná 5, 128 44 Prague, Czech Republic; 2Departamento de Farmacología, Farmacognosia y Botánica, Facultad de Farmacia, Universidad Complutense de Madrid, Madrid 28040, Spain; 3Microalgal Biotechnology Laboratory, Jacob Blaustein Institutes for Desert Research, Ben Gurion University of the Negev, Sede-Boker Campus 8499000, Israel and 4Department of Biology, Faculty of Sciences, Mohammed V University, Av Ibn Battouta 4, BP 1014 RP. Rabat, Morocco Abstract The widespread species Parmotrema crinitum (Ach.) M. Choisy and Parmotrema perlatum (Huds.) M. Choisy are mainly distinguished by their reproductive strategies. While P. crinitum propagates by isidia, P. perlatum produces soredia. In this study, we aim to evaluate the phylogenetic relationship between both species and to critically examine their species boundaries. To this purpose, 46 samples belonging to P. crinitum and P. perlatum were used in our analysis, including 22 for which we studied the morphology and chemistry, before extracting their DNA. We used 35 sequences of the internal transcribed spacer region of nuclear ribosomal DNA (ITS) of Parmotrema perlatum from Europe and Africa (20 of which were newly generated), and 11 of Parmotrema crinitum from Europe, North America and North Africa (two newly generated). Additionally, 28 sequences of several species from Parmotrema were included in the ITS dataset. The ITS data matrix was analyzed using different approaches, such as traditional phylogeny (maximum likelihood and Bayesian analyses), genetic distances, auto- matic barcode gap discovery (ABGD) and the coalescent-based method poisson tree processes (PTP), in order to test congruence among results. Our results indicate that all samples referred to P. crinitum and P. perlatum nested in a well-supported monophyletic clade, but phylogenetic relationships among them remain unresolved. Delimitations inferred from PTP, ABGD and genetic distance analyses were comparable and suggested that P. crinitum and P. perlatum belong to the same lineage. Interestingly, two samples of P. perlatum separate in a different monophyletic clade, which is supported as a different lineage by all the analyses. Key words: genetic diversity, ITS, lichen, phylogenetic analyses, species delimitation (Accepted 25 February 2022) Introduction Lichenized fungi form mutualistic relationships with photo- autotrophic organisms (photobionts), mainly green algae (Trebouxiophyceae and Ulvophyceae) and/or cyanobacteria. The lichen symbiosis has been highly successful within fungi, espe- cially Ascomycota, with c. 19 400 currently accepted species (Lücking et al. 2017) and an estimated diversity of more than 28 000 species (Lücking et al. 2009; Leavitt et al. 2013). Additionally, lichens are commonly used to assess environmental disturbance, serving as bioindicators of air pollution, forest age and health, and climate change (Nimis et al. 2002; Crespo et al. 2004; Giordani & Brunialti 2015; Sujetovienė 2015; Sancho et al. 2019; Abas 2021). Recognizing phylogenetic relationships and delimiting species in lichens are crucial for ecological and conservation studies, asses- sing biotic diversity, and identifying factors driving diversification. They are also important for future investigations because phylogen- etic differences may not be fully reflected in the phenotype. Traditionally, to infer taxonomic boundaries in lichen-forming fungi, thin-layer chromatography (Culberson 1972), morphology and the expression of signature secondary metabolites, and isola- tion and identification of lichen substances have been used (Huneck & Yoshimura 1996; Huneck 1999). However, these char- acters are highly variable and their homology has proved difficult to assess. For example, it has been shown that apothecial form and spore wall thickness have changed in parallel within Pertusaria. The Pertusaria-type ascus is plesiomorphic within the Pertusariaceae and thus cannot be used to circumscribe Pertusaria (Lumbsch & Schmitt 2001). However, with the use of molecular tools, it has been shown that phenotype-based tax- onomy may not reflect natural groups, including cases in which morphologically distinct forms formerly recognized as distinct species are shown to represent a polymorphic species (see e.g. Boluda et al. 2019), and other cases where several morphologically cryptic species are masked within a single nominal taxon (Del-Prado et al. 2016). Within lichen-forming Ascomycetes, Parmeliaceae (Lecanorales) constitutes one of the largest and best-studied families (Crespo Author for correspondence: Ayoub Stelate. E-mail: Ayoub.stelate@gmail.com Cite this article: Stelate A, Del-Prado R, Alors D, Tahiri H, Divakar PK and Crespo A (2022) Resolving the phylogenetic relationship between Parmotrema crinitum and Parmotrema perlatum populations. Lichenologist 54, 183–194. https://doi.org/10.1017/ S0024282922000147 © The Author(s), 2022. Published by Cambridge University Press on behalf of the British Lichen Society. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited. The Lichenologist (2022), 54, 183–194 doi:10.1017/S0024282922000147 https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://orcid.org/0000-0001-8929-3046 https://orcid.org/0000-0001-7288-9521 https://orcid.org/0000-0002-0300-0124 https://orcid.org/0000-0002-5271-0157 mailto:Ayoub.stelate@gmail.com https://doi.org/10.1017/S0024282922000147 https://doi.org/10.1017/S0024282922000147 https://creativecommons.org/licenses/by/4.0/ https://doi.org/10.1017/S0024282922000147 et al. 2007, 2011; Thell et al. 2012; Divakar et al. 2015). This family is usually characterized morphologically by a certain type of ascoma ontogeny and the presence of an ascomatal structure known as a cupulate exciple (Henssen & Jahns 1974). Parmeliaceae also com- prises species which are frequently used in biomonitoring studies of atmospheric pollution, such as Parmelia sulcata, Flavoparmelia caperata, Parmotrema perlatum and Punctelia subrudecta (Hawksworth & Rose 1970; Crespo et al. 2004; De La Cruz et al. 2018). The application of phylogenetic analysis based on molecular (DNA) characters to delimit species allows us to determine a poster- iori which types of phenotypic characters are good predictors of phylogenetic species and demonstrate how these characters evolve in this family and in lichenized fungi in general. These molecular data have led to the recognition of morphologically cryptic species, such as Parmelia serrana (Molina et al. 2004), P. barrenoae (Divakar et al. 2005a), P. encryptata (Molina et al. 2011a) and Melanelixia californica (Divakar et al. 2010), and conversely also to the union of species traditionally regarded as morphologically distinct (see e.g. Boluda et al. 2019). Recently, coalescent-based spe- cies delimitation approaches have shown to be well suited to critic- ally evaluate species boundaries in Parmeliaceae, as well as lichens in general (Leavitt et al. 2015). Furthermore, these methods can accurately display relationships when incomplete lineage sorting and gene tree heterogeneity hide phylogenetic relationships among species (Knowles & Carstens 2007; Camargo et al. 2012). Commonly used methods to critically evaluate species boundaries include the poisson tree processes (PTP) model (Zhang et al. 2013), the automatic barcode gap discovery (ABGD) (Puillandre et al. 2012), the general mixed Yule coalescent model (GMYC) (Pons et al. 2006; Monaghan et al. 2009), and SpedeSTEM (Ence & Carstens 2011). Parmotrema (Massalongo 1860) is one of the largest genera in the parmelioid group of the family Parmeliaceae. It includes c. 300 described species with an apparent centre of speciation in the Pacific Islands, tropical and subtropical regions of South America (Spielmann & Marcelli 2020). The species of the genus are characterized by a pored epicortex, large thalli with broad lobes, a broad, naked marginal zone on the lower surface, and large, thick-walled, ellipsoid ascospores, sublageniform or filiform conidia (Elix 1993), and (commonly) marginal cilia (Hale 1974). Reproductive strategies vary among Parmotrema taxa. Sexual repro- duction is restricted to characteristic fungal fruiting bodies (ascomata) producing ascospores. Ascospores are dispersed independently of the photosynthesizing partner (photobiont) and require reacquisition of the appropriate photobiont partner to re-establish the lichenized condition. Species within Parmotrema also commonly propagate asexually by means of vegetative diaspores, either isidia or soredia. These specialized vegetative reproductive propagules contain both fungal and algal symbionts, eliminating the need for independent acquisition of the appropriate photobiont. In a molecular phylogeny of parmotremoid lichens (Blanco et al. 2005), a single sample of both Parmotrema perlatum (Huds.) M. Choisy and Parmotrema crinitum (Ach.) M. Choisy, originally from Portugal, were included and these formed a well- supported monophyletic group, indicating that these specimens could be conspecific. Therefore, a critical evaluation of species boundaries is necessary. Parmotrema crinitum is characterized by the presence of cor- alloid branched, apically ciliate isidia or often eciliate isidia (Divakar & Upreti 2005), and the stictic acid complex in the medulla. According to Elix (1994), P. crinitum is a cosmopolitan species that is widespread in temperate, tropical regions and even in high humidity, sub-boreal forests (Elix 1994; Kurokawa & Lai 2001). Many European countries have reported the presence of P. crinitum (Jablońska et al. 2009) as have some Asian countries such as Japan (Yoshimura 1974), China (Wei 1991) and Taiwan (Wang-Yang & Lai 1976). Parmotrema perlatum is a greenish grey foliose lichen, saxico- lous or corticolous, loosely adnate, with rounded lobes. It can be recognized by its broad lobes, irregularly branched (7.0 cm), lat- erally overlapping, with frequent black cilia (0.20–)0.50–1.00 × 0.02(–0.10) mm, evenly distributed but less common in the lobe apices. Its soralia are marginal and linear, sometimes widely dis- tributed or subcontinuous, concolorous with the thallus, and the medulla is white. The thallus contains atranorin, stictic acid, hypostictic acid, menegazziaic acid and norstictic acid (Spielmann & Marcelli 2009). Parmotrema perlatum is a wide- spread species in temperate regions of the Northern and Southern Hemispheres: Asia (Hale 1965; Kurokawa 1991; Kurokawa & Lai 2001), Europe (Hale 1965), Africa (Hale 1965; Swinscow & Krog 1988), North America (Hale 1965; Brodo et al. 2001; Nash & Elix 2002), Central America (Hale 1965), South America (Hale 1965), Australia (Kantvilas 2019) and New Zealand (Blanchon 2013). Morphologically, P. crinitum and P. perlatum can easily be separated based on the characteristics associated with their differ- ent dispersal strategies. Parmotrema crinitum has apically ciliate isidia or often eciliate isidia, and P. perlatum has marginal soralia. Parmotrema perlatum is one of the most common and widely uti- lized lichens in the Ayurvedic system of medicine and has been overexploited for uses in traditional medicine in India (Kumar & Upreti 2001). In India, this species is currently considered threatened (Divakar & Upreti 2005). The aim of this study was to evaluate the phylogenetic relationships between P. crinitum and P. perlatum, and elucidate the possible monophyly between both species. Additionally, PTP, ABGD and genetic distance analyses were included to crit- ically assess species boundaries. The ITS data matrix consisted of a total of 74 sequences of Parmotrema species, including 11 of P. crinitum and 35 of P. perlatum, plus two outgroup sequences of Crespoa carneopruinata. Morphological and chem- ical features of all samples were critically examined. Material and Methods Taxon sampling Sequence data of the ITS locus were analyzed from 76 specimens, of which 11 sequences were of Parmotrema crinitum (two indivi- duals were newly sequenced and nine sequences downloaded from GenBank) and 35 of Parmotrema perlatum (20 individuals were sequenced in this study and 15 sequences downloaded from GenBank). The new specimens were collected from distant geographical regions throughout the species distributions (Table 1). Additionally, 30 ITS sequences from 30 specimens, belonging to 14 different species (6 were newly sequenced and 24 sequences downloaded from GenBank), are included in this study to test the monophyly of both species within the genus. This total included Crespoa carneopruinata which was selected as outgroup since it has previously been shown to be the sister group of Parmotrema (Divakar et al. 2015). Details of the material studied, including GenBank Accession numbers, are shown in Table 1. Species recognition was based mainly on morphological and chemical characters. 184 Ayoub Stelate et al. https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://doi.org/10.1017/S0024282922000147 Table 1. Parmotrema species used in this study with GenBank Accession numbers of the ITS sequences and voucher information for the specimens. Crespoa carneopruinata was included in the phylogenetic analyses as outgroup. * = newly generated sequences from specimens for which morphological and chemical characters are given in Table S1 (available online). Species Locality Voucher specimen Collector GenBank Accession number Crespoa carneopruinata 1 C. carneopruinata 2 Parmotrema clavuliferum P. crinitum 1 P. crinitum 2 P. crinitum 3 P. crinitum 4 P. crinitum 5 P. crinitum 6* P. crinitum 7 P. crinitum 8* P. crinitum 9 P. crinitum 10 P. crinitum 11 P. dilatatum 1* P. dilatatum 2* P. flavotinctum* P. grayanum P. haitiense P. hypoleucinum 1 P. hypoleucinum 2 P. hypoleucinum 3 P. hypoleucinum 4 P. mellissii* P. paulense* P. perforatum P. perlatum 1 P. perlatum 2* P. perlatum 3* P. perlatum 4* P. perlatum 5* P. perlatum 6 P. perlatum 7* P. perlatum 8 P. perlatum 9 P. perlatum 10* P. perlatum 11* P. perlatum 12 P. perlatum 13 P. perlatum 14 P. perlatum 15 P. perlatum 16* P. perlatum 17* P. perlatum 18* P. perlatum 19 P. perlatum 20 P. perlatum 21* P. perlatum 22* P. perlatum 23 P. perlatum 24 P. perlatum 25* P. perlatum 26* Costa Rica Costa Rica Gomera, Islas Canarias, Spain North Carolina, USA La Palma, Islas Canarias, Spain North Carolina, USA North Carolina, USA La Palma, Islas Canarias, Spain Madeira, Portugal Lisboa, Igreja Nova, Portugal Madeira, Portugal Tenerife, Islas Canarias, Spain Tenerife, Islas Canarias, Spain Tenerife, Islas Canarias, Spain Gomera, Islas Canarias, Spain China Pichincha, Ecuador La Palma, Islas Canarias, Spain Australia Spain Morocco Morocco Morocco Pichincha, Ecuador Medellin, Antioqia, Colombia USA Medium Atlas, Fes, Morocco Akchour, Chaouen, Morocco Akchour, Chaouen, Morocco Akchour, Chaouen, Morocco Akchour, Chaouen, Morocco Chaouen, Morocco Gomera, Islas Canarias, Spain Eregli Vakif, Zonguldak, Turkey Tenerife, Islas Canarias, Spain Chaouen, Morocco Chaouen, Morocco Lisboa, Igreja Nova, Portugal La Palma, Islas Canarias, Spain La Palma, Islas Canarias, Spain Chaouen, Morocco Gomera, Islas Canarias, Spain Gomera, Islas Canarias, Spain Gomera, Islas Canarias, Spain La Palma, Islas Canarias, Spain Tenerife, Islas Canarias, Spain Gomera, Islas Canarias, Spain Gomera, Islas Canarias, Spain Tenerife, Islas Canarias, Spain Tenerife, Islas Canarias, Spain Gomera, Islas Canarias, Spain Gomera, Islas Canarias, Spain 15171a 15171a MAF-Lich 20689 44262 MAF-Lich 16189 43711 43863 MAF-Lich 16188 MAF-Lich 20704 MAF-Lich 6061 MAF-Lich 20703 MAF-Lich 16174 MAF-Lich 16170 MAF-Lich 16169 MAF-Lich 20698 MAF-Lich 10164 54174 MAF-Lich 16190 MAF-Lich 7657 MAF-Lich 7637 MAF-Lich 16142 MAF-Lich 16141 MAF-Lich 16147 54312 MAF-Lich 20060 7983 MAF-Lich 16145 MAF-Lich 20712 MAF-Lich 20707 MAF-Lich 20711 MAF-Lich 20708 MAF-Lich 16146 MAF-Lich 20700 MAF-Lich 16192 MAF-Lich 16173 MAF-Lich 20709 MAF-Lich 20710 MAF-Lich 6965 MAF-Lich 16179 MAF-Lich 16187 MAF-Lich 16203 MAF-Lich 20699 MAF-Lich 20693 MAF-Lich 20691 MAF-Lich 16186 MAF-Lich 16162 MAF-Lich 20702 MAF-Lich 20696 MAF-Lich 16172 MAF-Lich 16166 MAF-Lich 20690 MAF-Lich 20695 R. Lücking R. Lücking A. Crespo, A. Santos J. C. Lendemer A. Crespo, A. Santos J. C. Lendemer J. C. Lendemer A. Crespo, A. Santos P. K. Divakar, M. Talavera P. K. Divakar, M. Talavera A. Crespo, A. Santos A. Crespo, A. Santos A. Crespo, A. Santos A. Crespo, R. Del-Prado, A. Santos A. R. Burgaz, M. A. Carrasco D. L. Hawksworth P. v. d. Boom A. Crespo, R. Del-Prado, P. K. Divakar S. H. Louwhoff A. Crespo, P. K. Divakar A. Crespo, P. K. Divakar, H. Tahiri, A. Crespo, P. K. Divakar, H. Tahiri A. Crespo, P. K. Divakar, H. Tahiri P. v. d. Boom P. K. Divakar M. Cole A. Crespo, P. K. Divakar, H. Tahiri A. Stelate, H. Tahiri A. Stelate, H. Tahiri A. Stelate, H. Tahiri A. Stelate, H. Tahiri A. Crespo, P. K. Divakar, H. Tahiri A. Crespo, A. Santos P. K. Divakar, R. Del-Prado A. Crespo, R. Del-Prado A. Stelate, H. Tahiri A. Stelate, H. Tahiri P. K. Divakar, M. Talavera A. Crespo, R. Del-Prado A. Crespo, R. Del-Prado A. Crespo, R. Del-Prado A. Crespo, A. Santos A. Crespo, A. Santos A. Crespo, A. Santos A. Crespo, R. Del-Prado A. Crespo, R. Del-Prado A. Crespo, A. Santos A. Crespo, A. Santos A. Crespo, R. Del-Prado A. Crespo, R. Del-Prado A. Crespo, A. Santos A. Crespo, A. Santos KP888204 EF042904 AY586577 KP943761 HM017032 KP943759 KP943760 HM017033 ON312512 AY586565 ON312511 HM017030 HM017028 HM017029 ON312500 ON312501 ON312516 HM017026 AY581055 AY586567 HM017035 HM017037 HM017036 ON312515 ON312513 AY586568 HM017044 ON312519 ON312518 ON312521 ON312523 HM017046 ON312505 HM017041 HM017042 ON312520 ON312522 AY586566 HM017040 HM017043 HM017044 ON312506 ON312510 ON312508 HM017048 HM017052 ON312503 ON312496 HM017049 HM017051 ON312507 ON312497 (Continued ) The Lichenologist 185 https://doi.org/10.1017/S0024282922000147 Published online by Cam bridge U niversity Press https://doi.org/10.1017/S0024282922000147 Chemistry and morphology Thallus morphology of all new specimens of Parmotrema crinitum and P. perlatum included in the molecular analyses was studied using a Nikon SMZ-1000 stereomicroscope to iden- tify morphological characteristics of the thallus: lobe shape, isi- dia, soralia, cilia and rhizines. This is because P. crinitum and P. perlatum are traditionally differentiated based on these fea- tures. Photographs were taken with a Nikon 105 mm f/2.8D AF Micro-Nikkor lens coupled to a Nikon D90 camera. Spot tests were carried out on the medulla with usual chem- ical reagents such as aqueous potassium hydroxide (K), Steiner’s stable paraphenylenediamine (PD) and aqueous cal- cium hypochlorite (C). Thin-layer chromatography (TLC) was carried out following Orange et al. (2010). We used TLC solvent system C (200 ml toluene/30 ml acetic acid) (Elix & Ernst-Russell 1993), with concentrated acetone extracts at 50 °C spotted onto silica gel 60 F254 aluminum sheets (Merck, Darmstadt). The aluminum sheets were dried for 10 min in an acetic acid atmosphere to maximize resolution. DNA extraction, PCR and sequencing Total DNA was extracted from a single, clean (under a dissect- ing microscope) lichen lobe using the DNeasy Plant Mini Kit (Qiagen, Barcelona) with a slight modification to the manufac- turer’s instruction (Crespo et al. 2001). Genomic DNA (5–25 ng, quantified using a quantitative PCR machine) was used for PCR amplifications of the ITS region. Standard PCR amplifications were conducted in 25 μl reaction volumes containing 12.5 μl of Master Mix (50 units/ml of Taq DNA polymerase supplied in a proprietary reaction buffer (pH 8.5), 400 μM dATP, 400 μM dGTP, 400 μM dCTP, 400 μM dTTP, 3 mM MgCl2), and 1.5 μl of each primer at 10 μM, 4.5 μl of water (H2O) and 5 μl of DNA template. Fungal nuclear internal transcribed spacer (ITS) rDNA was amplified using the primer pair ITS1F (5ʹ [CTT GGT CAT TTA GAG GAA GTA A] 3ʹ) (Gardes & Bruns 1993)/ LR1 (5ʹ [GGT TGG TTT CTT TTC CT] 3ʹ) (Vilgalys & Hester 1990). We also tried primer pair ITS1-LM (5ʹ [GAA CCT GCG GAA GGA TCA TT] 3ʹ) (Myllys et al. 2001) and ITS2-KL (5ʹ [ATG CTT AAG TTC AGC GGG TA] 3ʹ) (Lohtander et al. 1998), but we had better results with primer pair ITS1F/LR1. For any failed samples, we tried a nested PCR using the primer pair ITS1F/LR1 for the first PCR, and 5 μl of this PCR product as DNA template for the second PCR using ITS1-LM/ITS2-KL. With this nested PCR, we aimed to amplify a smaller region than with the previous primers. However, it did not yield any additional products. Amplification was run in an automatic thermocycler (Techne Progene, Jepson Bolton & Co., Watford, Herts, UK) using the following parameters: initial denaturation at 94 °C for 5 min followed by 35 cycles at 94 °C for 1 min, 54 °C (ITS1F/LR1) for 1 min, and 72 °C for 1 min 30 s, with a final extension at 72 °C for 10 min. Amplification products were visualized on 1% agarose gel stained with SYBR green DNA (Life Technologies Corporation, Grand Island, New York, USA), and subsequently purified using ExoSAP-IT (GE Healthcare, Chalfont St Giles, UK) according to the manufac- turer’s instructions. Sequencing was performed using BigDye Terminator reaction kit (ABI PRISM, Applied Biosystems, Waltham, Massachusetts, USA). Cycle sequencing reactionsTa b le 1. (C on tin ue d) Sp ec ie s Lo ca lit y Vo uc he r sp ec im en Co lle ct or G en B an k Ac ce ss io n nu m be r P. pe rl at um 27 P. pe rl at um 28 * P. pe rl at um 29 P. pe rl at um 30 * P. pe rl at um 31 * P. pe rl at um 32 * P. pe rl at um 33 * P. pe rl at um 34 * P. pe rl at um 35 P. ps eu do re tic ul at um P. re tic ul at um P. ro bu st um P. su bt in ct or iu m 1 P. su bt in ct or iu m 2* P. tin ct or um 1 P. tin ct or um 2 P. tin ct or um 3 P. tin ct or um 4 P. tin ct or um 5 P. tin ct or um 6 P. tin ct or um 7 P. tin ct or um 8 P. tin ct or um 9 P. tin ct or um 10 Te ne ri fe , Is la s Ca na ri as , Sp ai n G om er a, Is la s Ca na ri as , Sp ai n La Pa lm a, Is la s Ca na ri as , Sp ai n G om er a, Is la s Ca na ri as , Sp ai n G om er a, Is la s Ca na ri as , Sp ai n G om er a, Is la s Ca na ri as , Sp ai n G om er a, Is la s Ca na ri as , Sp ai n Li sb oa , Ig re ja N ov a, Po rt ug al Te ne ri fe , Is la s Ca na ri as , Sp ai n Ea st er n Ca pe , So ut h Af ri ca B lu f Ar ea , K en ya Po rt ug al In di a W es te rn , B lu f Ar ea , K en ya Sh iz uo ka , Ja pa n Sh iz uo ka , Ja pa n In ha ca Is la nd , M oz am bi qu e M t Ts uk ub a, Ja pa n Ch am ol i, In di a K an ag aw a, Ja pa n K an ag aw a, Ja pa n N ar a, Ja pa n O sa ka , Ja pa n In ha ca Is la nd , M oz am bi qu e M AF -L ic h 16 16 8 M AF -L ic h 20 70 1 M AF -L ic h 16 19 1 M AF -L ic h 20 69 7 M AF -L ic h 20 69 4 M AF -L ic h 20 68 8 M AF -L ic h 20 69 2 M AF -L ic h 20 70 6 M AF -L ic h 16 17 1 M AF -L ic h 10 28 7 M AF -L ic h 10 12 1 M AF -L ic h 10 16 6 G P G C0 2- 00 06 96 M AF -L ic h 20 70 5 TN S- L- Y. O . 53 99 TN S- L- Y. O . 53 75 A M AF -L ic h 18 19 3 M AF -L ic h 18 16 3 M AF -L ic h 18 17 0 TN S- L- Y. O . 59 47 TN S- L- Y. O . 54 86 TN S- L- Y. O . 58 13 TN S- L- Y. O . 54 80 M AF -L ic h 18 19 6 A. Cr es po , R. D el -P ra do A. Cr es po , A. Sa nt os A. Cr es po , R. D el -P ra do A. Cr es po , A. Sa nt os A. Cr es po , A. Sa nt os A. Cr es po , A. Sa nt os A. Cr es po , A. Sa nt os P. K. D iv ak ar , A. Ag ud o A. Cr es po , A. Sa nt os A. Cr es po , A. Sa nt os P. K. D iv ak ar P. K. D iv ak ar P. K. D iv ak ar A. Cr es po , R. D el -P ra do Y. O hm ur a Y. O hm ur a B . Ro ca -V al ie nt e A. Cr es po , P. K. D iv ak ar P. K. D iv ak ar M . Ta ke da Y. O hm ur a Y. O hm ur a N . H am ad a B . Ro ca -V al ie nt e, A. Lu m br er as H M 01 70 50 O N 31 25 04 H M 01 70 39 O N 31 25 02 O N 31 25 09 O N 31 24 98 O N 31 24 99 O N 31 25 17 H M 01 70 27 AY 64 28 28 H M 01 69 56 AY 58 65 69 AY 58 65 58 O N 31 25 14 AB 17 74 04 AB 17 74 01 K F1 29 44 7 K F1 29 41 4 K F1 29 42 2 K F1 29 46 4 K F1 29 46 3 K F1 29 46 2 K F1 29 46 0 K F1 29 45 0 186 Ayoub Stelate et al. https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://doi.org/10.1017/S0024282922000147 were performed with the same sets of primers used in the ampli- fication step. Sequencing reactions were electrophoresed on a 3730 DNA Analyzer (Applied Biosystems) at the Unidad de Genómica (Parque Científico de Madrid). Sequence alignments and phylogenetic analyses Sequence fragments generated for this study were assembled and edited using the program SeqMan v. 7 (Lasergene R, DNASTAR, Madison, Wisconsin, USA). Sequence identity was confirmed using the megaBLAST search function in GenBank. For the align- ment we prepared a matrix of 480 base pairs (bp). Then we used the program MAFFT v. 7 (Katoh & Standley 2013) with the para- meters set to default. If sequences had different lengths, only the part shared by all the sequences was used, and after manually removing ambiguously aligned nucleotide positions, we kept 451 unambiguously aligned base pairs for the final matrix that we used as input for the phylogenetic reconstruction. The align- ment was analyzed using maximum likelihood (ML) and a Bayesian Markov chain Monte Carlo (B/MCMC) approach. The ML analysis was performed using an online version of the program RAxML-HPC BlackBox v. 8.2.8 (Stamatakis 2006; Stamatakis et al. 2008), implemented on the Cipres science gate- way (https://www.phylo.org/portal2/home.action; Miller et al. 2010). We used a GTRGAMMA model, which includes a param- eter (Γ) for rate heterogeneity among sites but chose not to include a parameter to estimate the proportion of invariable sites (Stamatakis 2006; Stamatakis et al. 2008). Support values were assessed using the ‘rapid bootstrapping’ option with 1000 replicates. For the Bayesian reconstruction, we used MrBayes v. 3.2.1 (Ronquist & Huelsenbeck 2003). The analysis was performed assuming a discrete gamma distribution with six rate categories (GTR + G). The nucleotide-substitution model parameters were selected using the Akaike Information Criterion as implemented in jModelTest (Posada 2008), molecular clock not assumed. A run with 10 million generations, starting with a random tree and employing 12 simultaneous chains, was executed. Trees were saved to a file every 200th generation. The first 2 million generations (i.e. 20 000 trees) were deleted as the ‘burn-in’ of the chains. We plotted the log-likelihood scores of sample points against generations using Tracer v. 1.5 (Rambaut & Drummond 2007) and determined that stationarity had been achieved when the log-likelihood values of the sample points reached an equilib- rium value (Huelsenbeck & Ronquist 2001). The trees obtained before stationarity were discarded. Posterior probabilities (PPs) were obtained from the 50% majority-rule consensus of sampled trees after excluding the initial 20% as burn-in. Only clades that received bootstrap support ≥ 70% in the ML analyses and PP ≥ 0.95 were considered strongly supported. The phylogenetic tree was drawn using FigTree v. 1.2.3 (Rambaut 2009) and modified with CorelDRAW v. 8. Candidate species identification In order to establish candidate species limits in the phylogenetic tree, three computational approaches were used: 1) Poisson tree processes (PTP): this method does not require an ultra-metric tree, as the transition point between intra- and inter-specific branching rates is identified directly using the number of nucleotide substitutions (Zhang et al. 2013). PTP incorporates the number of substitutions in the model of spe- ciation and assumes that the probability of a substitution giv- ing rise to a speciation event follows a Poisson distribution. The branch lengths of the input tree are supposed to be gen- erated by two independent classes of Poisson events, one cor- responding to speciation and the other to coalescence. The ML phylogeny of the ITS dataset obtained with RAxML was used as the input trees to run PTP species delimitation analysis on the PTP web server (http://species.h-its.org/ptp/). We ran the PTP analysis for 100 000 MCMC generations, with a thinning value of 100 and a burn-in of 25%. Outgroup taxa were removed for species delimitation. 2) Automatic barcode gap discovery (ABGD): this is an auto- matic procedure that sorts the sequences into hypothetical species based on the barcode gap. This method automatically finds the distance where the barcode gap is located (Puillandre et al. 2012). The ABGD method was carried out on the ITS dataset using the web interface at http://wwwabi.snv.jussieu. fr/public/abgd/abgdweb.html. Default parameters were chosen using Kimura 2-parameter (K2P) distances that correct for transition rate bias (relative to transversions) in the substitu- tion process (Kimura 1980). The default for the minimum relative gap width was set to different values between 0.1 and 0.15. 3) Genetic distances: pairwise ML distances (given as the number of nucleotide substitutions per site) among the ITS rDNA sequences of Parmotrema crinitum and P. perlatum were calculated. Genetic distances were calculated with TREE- PUZZLE v. 5.2 (Schmidt et al. 2002) using the GTR model of nucleotide substitution, assuming a discrete gamma distri- bution with six rate categories. The program generates an output file which consists of a triangular matrix with all pair- wise distances between all the samples included. This matrix was visualized with Microsoft Excel 2010 and genetic distances between different specimens of P. crinitum and P. perlatum were manually identified. Candidate species were proposed based on the threshold of 0.016 substitutions per site (s/s) which separates intra- and interspecific distances in parme- lioid lichens (Del-Prado et al. 2010). Distance values in the matrix ≤ 0.016 s/s have been considered as the values between samples of a single species. By using the Excel filter, we sepa- rated values ≤ 0.016, providing for every specimen included in the analysis a group of specimens that share the values charac- terizing its species range. Results Morphological and chemical analyses Morphological and chemical characters of the specimens (P. cri- nitum and P. perlatum) newly generated and included in the phylogenetic analysis are presented in Supplementary Material Table S1, available online. Parmotrema crinitum and P. perlatum could not be distinguished based on the colour of the thalli, pres- ence of cilia (Fig. 1A & B) and colour of the lower surface (Fig. 1E & F). However, the reproductive structures allowed unequivocal separation into the two species; Parmotrema perlatum has soralia (Fig. 1C) and P. crinitum has isidia (Fig. 1D). The TLC test showed that both species contained the stictic acid complex (stictic acid, menegazzic acid and hypostictic acid) and atranorin. The spot tests of P. crinitum were similar to those of P. perlatum (Fig. 2, Supplementary Material Table S1). The Lichenologist 187 https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://www.phylo.org/portal2/ 20home.action https://www.phylo.org/portal2/ 20home.action http://species.h-its.org/ptp/ http://species.h-its.org/ptp/ http://wwwabi.snv.jussieu.fr/public/abgd/abgdweb.html http://wwwabi.snv.jussieu.fr/public/abgd/abgdweb.html http://wwwabi.snv.jussieu.fr/public/abgd/abgdweb.html https://doi.org/10.1017/S0024282922000147 However, the samples Parmotrema perlatum 34 and P. perlatum 35 lacked menegazzic acid (nos 10 and 12 in Fig. 2). The spot test was congruent with other samples of P. perlatum, and the morphology of both samples showed the same characteristics of P. perlatum, although the cilia were not as abundant as is usual in this species. Phylogenetic analysis The ITS data matrix consisted of a total of 74 sequences of Parmotrema species, including 46 of P. crinitum and P. perlatum (11 sequences of P. crinitum and 35 sequences of P. perlatum), with two sequences of Crespoa carneopruinata were included as outgroup (Table 1). The final matrix used as input for the phylo- genetic reconstruction contained 451 unambiguously aligned base pairs (bp). The tree reconstruction (Fig. 3) comprised the follow- ing species of Parmotrema, included to test the monophyly of P. crinitum and P. perlatum: Parmotrema clavuliferum, P. dilata- tum, P. flavotinctum, P. grayanum, P. haitiense, P. hypoleucinum, P. mellissii, P. paulense, P. perforatum, P. pseudoreticulatum, P. reticulatum, P. robustum, P. subtinctorium and P. tinctorum. The phylogenetic trees estimated from ML and Bayesian methods did not show any well-supported conflict; ML topologies are presented with bootstrap and B/MCMC analysis with poster- ior probability (ML bootstrap ≥ 70%; PP ≥ 0.95 in B/MCMC ana- lysis) (Fig. 3). The LnL value was −1938.068210 for ML, and Fig. 1. Main morphological characters of P. perlatum (A, C, E) and P. crinitum. (B, D, F), as indicated with arrows. A & B, black, simple cilia. C, marginal soralia concolorous with the thallus. D, isidia with apical cilia. E & F, black lower surface with brown naked zone peripherally. In colour online. 188 Ayoub Stelate et al. https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://doi.org/10.1017/S0024282922000147 −2184.428 for the Bayesian analysis. The phylogenetic analysis showed that all samples of Parmotrema included in this analysis formed a monophyletic group. Within this genus the samples of Parmotrema crinitum and P. perlatum were grouped in one monophyletic group (clade A), with the exception of two samples of P. perlatum that were grouped in a separate clade B. While clade A comprised specimens of P. crinitum and P. perlatum col- lected from various geographical regions, clade B included only two specimens of P. perlatum collected from Tenerife (Canary Islands) and Lisbon (Portugal). Although four well-supported monophyletic clusters can be recognized in clade A, the phylogenetic relationships among them remained unresolved (Fig. 3). Cluster A1 (Fig. 3) contained 15 specimens of Parmotrema perlatum, most of them collected from Morocco, except four from the Canary Islands and two from Portugal and Turkey. Cluster A2 included eight specimens of Parmotrema crinitum from different geographical regions (Table 1). Cluster A3 grouped 18 samples of Parmotrema perlatum from the Canary Islands (Tenerife, Palma, Gomera), and cluster A4 included three samples of P. crinitum from the Canary Islands. Identifying candidate species We used the same RaxML tree obtained from the phylogenetic analysis to illustrate the delimitation of putative species recog- nized by the different approaches conducted with the ITS dataset. ABGD, PTP and genetic distance analyses applied to the ITS data- set detected two candidate species that corresponded to the well- supported clades A (including P. crinitum + P. perlatum) and B (including P. perlatum) obtained in the phylogenetic analysis (Fig. 3). Fig. 2. Thin-layer chromatography profile of Parmotrema species (included in the DNA analysis) in solvent system C. P is the control, Pleurosticta acetabulum. Lanes 1 & 7, P. dilatatum. Lane 2, P. flavotinctum. Lane 3, P. mellissii. Lane 4, P. paulense. Lanes 5, 6, 8, 11 & 19, P. crinitum. Lane 9, empty. Lanes 10 & 12, Parmotrema perlatum. Lanes 13, 14, 15, 16, 17, 18 & 20, P. perlatum; these have the same TLC result as P. crinitum but were distinguished from each other based on the morph- ology. a = norstictic acid; b = atranorin; c = protocetraric acid; d = unknown; e = hypostictic acid; f = stictic acid; g = menegazzic acid. In colour online. The Lichenologist 189 https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://doi.org/10.1017/S0024282922000147 Fig. 3. Phylogenetic relationships between Parmotrema crinitum and P. perlatum. The green (thickened) line indicates the clade supported by maximum likelihood (ML) and Bayesian analyses; * = clades supported only by ML. ** = clades supported only by Bayesian analysis. Alphanumeric labels indicate clades and clusters. Species delimitation scenarios obtained from PTP, ABGD and genetic distances are indicated in columns to the right. In colour online. 190 Ayoub Stelate et al. https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://doi.org/10.1017/S0024282922000147 Discussion Previous molecular phylogenetic studies have proposed the monophyly of Parmotrema crinitum and Parmotrema perlatum (Blanco et al. 2006; Crespo et al. 2010). The primary goal of the present investigation was to evaluate the phylogenetic rela- tionship between the two species using ML and Bayesian analyses. Additionally, we used different approaches to assess the species boundary, such as genetic distances based on the threshold of 0.015–0.017 s/s established by Del-Prado et al. (2010) to measure intra- and interspecific genetic distances in parmelioid lichens, and separate ranges of intra- and interspecific divergence. This threshold was established using both phylogenetically and mor- phologically well-delimited species (for details see Del-Prado et al. (2010)). We also used the poisson tree processes (PTP) model (Zhang et al. 2013) and the automatic barcode gap discov- ery (ABGD; Puillandre et al. 2012). Furthermore, we combined ecological, biogeographical, morphological and chemical data. Most of the samples belonging to P. crinitum and P. perlatum were recovered in a well-supported monophyletic clade (clade A, Fig. 3), with the exception of two samples of P. perlatum, one from Tenerife (Canary Islands) and one from Lisbon (Portugal), that are separated in the monophyletic clade B, supported by ML and B/MCMC analysis. Morphologically, P. perlatum in clade A has more cilia compared to P. perlatum in clade B. Thin-layer chromatography (TLC) revealed that the samples in clade B (nos 10 and 12 in Fig. 2) differ chemically from the sam- ples of Parmotrema perlatum grouped in clade A. Several studies have shown that environmental factors, such as light, temperature, pH and culture media, can influence the secondary metabolism in lichens (BeGora & Fahselt 2001). However, we have samples from the same localities and same conditions as the samples in clade A. Furthermore, ABGD, PTP and genetic distance analyses supported clade B as a new sister group to clade A (Fig. 3), sug- gesting polyphyly, a common phenomenon in Parmeliaceae in general (Lumbsch & Leavitt 2011; Leavitt et al. 2016) and Parmotrema in particular (Divakar et al. 2005b; Del-Prado et al. 2016, 2019; Widhelm et al. 2016). Clade A was split into 4 groups; however, the phylogenetic relationships among them remained unresolved (Fig. 3). Previous studies have shown that phylogenetic analyses alone are insufficient to explain phylogenetic relationships within Parmeliaceae. For example, based only on maximum parsimony and Bayesian analyses the sorediate P. sulcata was shown to belong to the same clade as the isidiate P. squarrosa (Molina et al. 2004), leaving the authors unable to reach any conclusion regarding species boundaries. Similarly, phylogenetic analyses were insufficient to resolve genetic variability among Parmelia saxatilis specimens; whereas samples from distant geographical regions formed a monophyletic group, samples from neighbour- ing localities were separated (Crespo et al. 2002; Molina et al. 2011b, 2017). A study on Usnea perpusilla demonstrated that it was necessary to use a combined approach with molecular and morphological data to assess species boundaries in closely related and morphologically variable species (Wirtz et al. 2008). For this reason, coalescent-based species delimitation analyses have been applied with the goal of explaining relationships among clades and delimiting species boundaries (Parnmen et al. 2012; Leavitt et al. 2013). PTP, ABGD and genetic distance analyses supported clades A and B as putative distinct species. The values obtained were within the interspecific ranges (genetic distances as defined in Del-Prado et al. (2010)), which would support the two samples of P. perlatum in clade B as a separate species, different from clade A. Focusing on clade A, our species delimitation approaches (PTP, ABGD and genetic distance) supported P. crinitum and P. perlatum as one species. However, based only on one molecular marker (ITS) and the various phylogenetic analyses and species delimitation approaches used, can we consider Parmotrema perla- tum conspecific with P. crinitum? Similar cases have been reported in the Parmotrema perfora- tum species complex (Widhelm et al. 2016), and authors have suggested that the phylogenetic relationships between sexual and asexual populations of this species group could be more com- plex than previously assumed. They also suggest that traditional tools based on reproductive mode and secondary metabolites (Culberson & Culberson 1973) are no longer the key to identify species such as P. perforatum (Widhelm et al. 2016). A previous study on Umbilicaria (Ott et al. 2004) focused on Umbilicaria kappenii and U. antarctica, which are distinguished only by their reproductive strategies. Umbilicaria antarctica pro- pagates by thalloconidia and U. kappenii exhibits a variety of asex- ual propagules: soredia, adventive lobes and sorediate thallyles. To infer phylogenetic relationships between both species, the authors used molecular data from three loci. Results indicated that all samples morphologically referred to U. antarctica and U. kappenii form a monophyletic group and they proposed placing U. kappe- nii into synonymy with U. antarctica (Ott et al. 2004). In Parmelia, Molina et al. (2011b) rejected the previous hypothesis that P. sulcata and P. squarrosa form a monophyletic group (Molina et al. 2004) and based on phylogenetic analyses and species delimitation approaches confirmed that P. sulcata is not conspecific with P. squarrosa. In addition, P. squarrosa is a repro- ductively isolated lineage and genetic distances clearly separate this from other Parmelia species (Del-Prado et al. 2010). Within Usnea, studies have suggested that Usnea subfloridana was a secondary species derived from U. florida (Clerc 1984, 1987, 1997; Purvis et al. 1992). Usnea florida and U. subfloridana show many morphological similarities but they have different reproductive strategies. Usnea florida always displays many apothecia and produces no specialized asexual propagules. Usnea subfloridana has soralia, isidiomorphs and occasionally apothecia. Multilocus phylogenetic analyses based on sequences of the ITS, IGS and LSU regions of the nuclear ribosomal DNA and the gene coding for β-tubulin, Mcm7, RPB1 and RPB2 showed that specimens of the two morphospecies formed one monophyletic intermixed group, and not two groups correspond- ing to morphology (Articus et al. 2002; Mark et al. 2016). This topology was further confirmed with a coalescent-based species delimitation approach (Mark et al. 2016). Authors have suggested that they could be conspecific but taxonomic conclusions must await further study. Moreover, a recent study using RADseq data suggests that closely related lichen species may need genome- wide data to test their species boundaries (Grewe et al. 2018). Traditionally it was thought that asexually reproducing species in lichens, and in filamentous fungi in general, was an evolution- ary dead end (Normark et al. 2003). However, recent molecular studies have demonstrated that lineages with vegetative propaga- tion can also present high genetic diversity (e.g. Parmelia sulcata (Molina et al. 2011b), Parmotrema reticulatum (Del-Prado et al. 2016)). However, even in the absence of sexual reproduction, lichens can exchange photobionts and this process could provide opportunities for gene transfer (Piercey-Normore 2006; Nelsen & Gargas 2008). The Lichenologist 191 https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://doi.org/10.1017/S0024282922000147 While our study confirmed the monophyly of an intermixed clade of Parmotrema crinitum and P. perlatum, the taxonomic conclusion must await additional studies including more markers. The phylogenetic tree lacked ML or B/MCMC support for other widely accepted Parmotrema species, such as P. pseudoreticulatum and P. reticulatum. A more comprehensive taxon sampling and additional molecular markers will therefore be needed before making a formal taxonomic decision on the status of clade B. Acknowledgements. This study was supported by the Spanish Ministerio de Ciencia e Innovacion (PID2019-105312GB-I00) and the Santander-Universidad Complutense de Madrid (PR87/19-22637 and G/6400100/3000). Author ORCIDs. Ayoub Stelate, 0000-0001-8929-3046; David Alors, 0000-0001-7288-9521; Pradeep Kumar Divakar, 0000-0002-0300-0124; Ana Crespo, 0000-0002-5271-0157. Supplementary Material. To view Supplementary Material for this article, please visit https://doi.org/10.1017/S0024282922000147. References Abas A (2021) A systematic review on biomonitoring using lichen as the bio- logical indicator: a decade of practices, progress and challenges. Ecological Indicators 121, 107197. Articus K, Mattsson JE, Tibell L, Grube M and Wedin M (2002) Ribosomal DNA and B-tubulin data do not support the separation of the lichens Usnea florida and U. subfloridana as distinct species. Mycological Research 106, 412–418. BeGora MD and Fahselt D (2001) Usnic acid and atranorin concentrations in lichens in relation to bands of UV irradiance. Bryologist 104, 134–140. Blanchon DJ (2013) Auckland lichens. Auckland Botanical Society Journal 68, 21–27. Blanco O, Crespo A, Divakar PK, Elix JA and Lumbsch HT (2005) Molecular phylogeny of parmotremoid lichens (Ascomycota, Parmeliaceae). Mycologia 97, 150–159. Blanco O, Crespo A, Ree RH and Lumbsch HT (2006) Major clades of par- melioid lichens (Parmeliaceae, Ascomycota) and the evolution of their mor- phological and chemical diversity. Molecular Phylogenetics and Evolution 39, 52–69. Boluda CG, Rico VJ, Divakar PK, Nadyeina O, Myllys L, McMullin RT, Zamora JC, Scheidegger C and Hawksworth DL (2019) Evaluating meth- odologies for species delimitation: the mismatch between phenotypes and genotypes in lichenized fungi (Bryoria sect. Implexae, Parmeliaceae). Persoonia 42, 75–100. Brodo IM, Sharnoff SD and Sharnoff S (2001) Lichens of North America. New Haven, Connecticut: Yale University Press. Camargo A, Morando M, Avila LJ and Sites JW (2012) Species delimitation with ABC and other coalescent-based methods: a test of accuracy with simulations and an empirical example with lizards of the Liolaemus darwi- nii complex (Squamata: Liolaemidae). Evolution 66, 2834–2849. Clerc P (1984) Contribution à la révision de la systématique des Usnées (Ascomycotina, Usnea) d’Europe. 1. Usnea florida (L.) Wigg. emend. Clerc. Cryptogamie, Bryologie et Lichénologie 5, 333–360. Clerc P (1987) Systematics of the Usnea fragilescens aggregate and its distribu- tion in Scandinavia. Nordic Journal of Botany 7, 479–495. Clerc P (1997) Notes on the genus Usnea Dill. ex Adanson. Lichenologist 29, 209–215. Crespo A, Blanco O and Hawksworth DL (2001) The potential of mitochon- drial DNA for establishing phylogeny and stabilizing generic concepts in the parmelioid lichens. Taxon 50, 807–819. Crespo A, Molina MC, Blanco O, Schroeter B, Sancho LG and Hawksworth DL (2002) rDNA ITS and β-tubulin gene sequence analyses reveal two monophyletic groups within the cosmopolitan lichen Parmelia saxatilis. Mycological Research 106, 788–795. Crespo A, Divakar PK, Argüello A, Gasca C and Hawksworth DL (2004) Molecular studies on Punctelia species of the Iberian Peninsula, with an emphasis on specimens newly colonizing Madrid. Lichenologist 36, 299–308. Crespo A, Lumbsch HT, Mattsson JE, Blanco O, Divakar PK, Articus K, Wiklund E, Bawingan PA and Wedin M (2007) Testing morphology- based hypotheses of phylogenetic relationships in Parmeliaceae (Ascomycota) using three ribosomal markers and the nuclear RPB1 gene. Molecular Phylogenetics and Evolution 44, 812–824. Crespo A, Kauff F, Divakar PK, del Prado R, Pérez-Ortega S, Amo de Paz G, Ferencova Z, Blanco O, Roca-Valiente B, Núñez-Zapata J, et al. (2010) Phylogenetic generic classification of parmelioid lichens (Parmeliaceae, Ascomycota) based on molecular, morphological and chemical evidence. Taxon 59, 1735–1753. Crespo A, Divakar PK and Hawksworth DL (2011) Generic concepts in par- melioid lichens, and the phylogenetic value of characters used in their cir- cumscription. Lichenologist 43, 511–535. Culberson CF (1972) Improved conditions and new data for identification of lichen products by a standardized thin-layer chromatographic method. Journal of Chromatography A 72, 113–125. Culberson WL and Culberson CF (1973) Parallel evolution in the lichen- forming fungi. Science 180, 196–198. De La Cruz ARH, De La Cruz JKH, Tolentino DA and Gioda A (2018) Trace element biomonitoring in the Peruvian Andes metropolitan region using Flavoparmelia caperata lichen. Chemosphere 210, 849–858. Del-Prado R, Cubas P, Lumbsch HT, Divakar PK, Blanco O, Amo de Paz G, Molina MC and Crespo A (2010) Genetic distances within and among spe- cies in monophyletic lineages of Parmeliaceae (Ascomycota) as a tool for taxon delimitation. Molecular Phylogenetics and Evolution 56, 125–33. Del-Prado R, Divakar PK, Lumbsch HT and Crespo A (2016) Hidden gen- etic diversity in an asexually reproducing lichen forming fungal group. PLoS ONE 11, e0161031. Del-Prado R, Buaruang K, Lumbsch HT, Crespo A and Divakar PK (2019) DNA sequence-based identification and barcoding of a morphologically highly plastic lichen forming fungal genus (Parmotrema, Parmeliaceae) from the tropics. Bryologist 122, 281–291. Divakar PK and Upreti DK (2005) Parmelioid Lichens in India (A Revisionary Study). Dehradun: Bishen Singh Mahendra Pal Singh. Divakar PK, Molina MC, Lumbsch HT and Crespo A (2005a) Parmelia bar- renoae, a new lichen species related to Parmelia sulcata (Parmeliaceae) based on molecular and morphological data. Lichenologist 37, 37–46. Divakar PK, Blanco O, Hawksworth DL and Crespo A (2005b) Molecular phylogenetic studies on the Parmotrema reticulatum (syn. Rimelia reticulata) complex, including the confirmation of P. pseudoreticulatum. Lichenologist 37, 55–65. Divakar PK, Figueras G, Hladun N and Crespo A (2010) Molecular phylo- genetic studies reveal an undescribed species within the North American concept of Melanelixia glabra (Parmeliaceae). Fungal Diversity 42, 47–55. Divakar PK, Crespo A, Wedin M, Leavitt SD, Hawksworth DL, Myllys L, McCune B, Randlane T, Bjerke JW, Ohmura Y, et al. (2015) Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi. New Phytologist 208, 1217–1226. Elix JA (1993) Progress in the generic delimitation of Parmelia sensu lato lichens (Ascomycotina: Parmeliaceae) and a synoptic key to the Parmeliaceae. Bryologist 96, 359–383. Elix JA (1994) Flora of Australia. Volume 55. Lichens: Lecanorales. 2. Parmeliaceae. Canberra: Australian Biological Resources Study. Elix JA and Ernst-Russell KD (1993) A Catalogue of Standardized Thin Layer Chromatographic Data and Biosynthetic Relationships for Lichen Substances. 2nd Edn. Canberra: Australian National University. Ence DD and Carstens BC (2011) SpedeSTEM: a rapid and accurate method for species delimitation. Molecular Ecology Resources 11, 473–480. Gardes M and Bruns TD (1993) ITS primers with enhanced specificity for basidiomycetes – application to the identification of mycorrhizae and rusts. Molecular Ecology 2, 113–118. Giordani P and Brunialti G (2015) Sampling and interpreting lichen diversity data for biomonitoring purposes. In Upreti DK, Divakar PK, Shukla V and Bajpai R (eds), Recent Advances in Lichenology, Vol. 1. New Delhi: Springer, pp. 19–46. Grewe F, Lagostina E, Wu H, Printzen C and Lumbsch HT (2018) Population genomic analyses of RAD sequences resolves the phylogenetic relationship of the lichen-forming fungal species Usnea antarctica and Usnea aurantiacoatra. MycoKeys 43, 91–113. 192 Ayoub Stelate et al. https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://orcid.org/ https://orcid.org/0000-0001-8929-3046 https://orcid.org/0000-0001-7288-9521 https://orcid.org/0000-0002-0300-0124 https://orcid.org/0000-0002-5271-0157 https://doi.org/10.1017/S0024282922000147 https://doi.org/10.1017/S0024282922000147 https://doi.org/10.1017/S0024282922000147 Hale ME (1965) A monograph of Parmelia subgenus Amphigymnia. Contributions from the United States National Herbarium 36, 193–358. Hale ME (1974) Bulbothrix, Parmelina, Relicina, and Xanthoparmelia, four new genera in the Parmeliaceae. Phytologia 28, 479–490. Hawksworth DL and Rose L (1970) Qualitative scale of estimating sulfur dioxide air pollution in England and Wales using epiphytic lichens. Nature 227, 145–148. Henssen A and Jahns HM (1974) Lichenes. Eine Einfuhrung in die Flechtenkunde. Stuttgart: Georg Thieme Verlag. Huelsenbeck JP and Ronquist F (2001) MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17, 754–755. Huneck S (1999) The significance of lichens and their metabolites. Naturwissenschaften 86, 559–570. Huneck S and Yoshimura I (1996) Identification of Lichen Substances. New York: Springer. Jablońska A, Oset M and Kukwa M (2009) The lichen family Parmeliaceae in Poland. I. The genus Parmotrema. Acta Mycologica 44, 211–222. Kantvilas G (2019) An annotated catalogue of the lichens of Kangaroo Island, South Australia. Swainsona 32, 1–97. Katoh K and Standley DM (2013) MAFFT multiple sequence alignment soft- ware version 7: improvements in performance and usability. Molecular Biology and Evolution 30, 772–780. Kimura M (1980) A simple method for estimating evolutionary rate of base substitution through comparative studies of nucleotide sequences. Journal of Molecular Evolution 16, 111–120. Knowles LL and Carstens BC (2007) Delimiting species without monophy- letic gene trees. Systematic Biology 56, 887–895. Kumar K and Upreti DK (2001) Parmelia spp. (Lichens) in ancient medicinal plant lore of India. Economic Botany 55, 458–459. Kurokawa S (1991) Japanese species and genera of the Parmeliaceae. Journal of Japanese Botany 66, 152–159. Kurokawa S and Lai MJ (2001) Parmelioid lichen genera and species in Taiwan. Mycotaxon 77, 225–284. Leavitt SD, Esslinger TL, Spribille T, Divakar PK and Lumbsch HT (2013) Multilocus phylogeny of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota): insights on diversity, distributions, and a com- parison of species tree and concatenated topologies.Molecular Phylogenetics and Evolution 66, 138–152. Leavitt SD, Moreau CS and Lumbsch HT (2015) The dynamic discipline of species delimitation: progress toward effectively recognizing species boundar- ies in natural populations. In Upreti DK, Divakar PK, Shukla V and Bajpai R (eds), Recent Advances in Lichenology, Vol. 2. New Delhi: Springer, pp. 11–44. Leavitt SD, Divakar PK, Crespo A and Lumbsch HT (2016) A matter of time – understanding the limits of the power of molecular data for delimit- ing species boundaries. Herzogia 29, 479–492. Lohtander K, Källersjö M and Tehler A (1998) Dispersal strategies in Roccellina capensis (Arthoniales). Lichenologist 30, 341–350. Lücking R. Rivas Plata E, Chaves JL, Umaña L and Sipman HJM (2009) How many tropical lichens are there… really? Bibliotheca Lichenologica 100, 399–418. Lücking R, Hodkinson BP and Leavitt SD (2017) The 2016 classification of lichenized fungi in the Ascomycota and Basidiomycota – approaching one thousand genera. Bryologist 119, 361–416. Lumbsch HT and Leavitt SD (2011) Goodbye morphology? A paradigm shift in the delimitation of species in lichenized fungi. Fungal Diversity 50, 59–72. Lumbsch HT and Schmitt I (2001) Molecular data suggest that the lichen genus Pertusaria is not monophyletic. Lichenologist 33, 161–170. Mark K, Cornejo C, Keller C, Flück D and Scheidegger C (2016) Barcoding lichen-forming fungi using 454 pyrosequencing is challenged by artifactual and biological sequence variation. Genome 59, 685–704. Massalongo AB (1860) Esame comparativo di alcuni generi di licheni. Atti dell’ I. R. Istituto Veneto di Scienze, Lettere ed Arti, series 3 5, 247–276. Miller MA, Pfeiffer W and Schwartz T (2010) Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In Proceedings of the Gateway Computing Environments Workshop (GCE), 14 November 2010, New Orleans, Louisiana, pp. 1–8. Molina MC, Crespo A, Blanco O, Lumbsch HT and Hawksworth DL (2004) Phylogenetic relationships and species concepts in Parmelia s. str. (Parmeliaceae) inferred from nuclear ITS rDNA and β-tubulin sequences. Lichenologist 36, 37–54. Molina MC, Divakar PK, Millanes AM, Sanchez E, Del-Prado R, Hawksworth DL and Crespo A (2011a) Parmelia sulcata (Ascomycota: Parmeliaceae), a sympatric monophyletic species complex. Lichenologist 43, 585–601. Molina MC, Del-Prado R, Divakar PK, Sánchez-Mata D and Crespo A (2011b) Another example of cryptic diversity in lichen-forming fungi: the new species Parmelia mayi (Ascomycota: Parmeliaceae). Organisms, Diversity and Evolution 11, 331–342. Molina MC, Divakar PK, Goward T, Millanes AM, Lumbsch HT and Crespo A (2017) Neogene diversification in the temperate lichen-forming fungal genus Parmelia (Parmeliaceae, Ascomycota). Systematics and Biodiversity 15, 166–181. Monaghan MT, Wild R, Elliot M, Fujisawa T, Balke M, Inward DJG, Lees DC, Ranaivosolo R, Eggleton P, Barraclough TG, et al. (2009) Accelerated species inventory on Madagascar using coalescent-based models of species delineation. Systematic Biology 58, 298–311. Myllys L, Lohtander K and Tehler A (2001) β-tubulin, ITS and group I intron sequences challenge the species pair concept in Physcia aipolia and P. caesia. Mycologia 93, 335–343. Nash TH, III and Elix JA (2002) Parmotrema. In Nash TH, III, Ryan BD, Gries C and Bungartz F (eds), Lichen Flora of the Greater Sonoran Desert Region, Vol. I. Tempe, Arizona: Lichens Unlimited, Arizona State University, pp. 318–329. Nelsen MP and Gargas A (2008) Dissociation and horizontal transmission of codispersing lichen symbionts in the genus Lepraria (Lecanorales: Stereocaulaceae). New Phytologist 177, 264–275. Nimis PL, Scheidegger C and Wolseley P (2002) Monitoring with Lichens – Monitoring Lichens. Dordrecht: Kluwer Academics. Normark BB, Judson OP and Moran NA (2003) Genomic signatures of ancient asexual lineages. Biological Journal of the Linnean Society 79, 69–84. Orange A, James PW and White FJ (2010) Microchemical Methods for the Identification of Lichens. 2nd Edn. London: British Lichen Society. Ott S, Brinkmann M, Wirtz N and Lumbsch HT (2004) Mitochondrial and nuclear ribosomal DNA data do not support the separation of the Antarctic lichens Umbilicaria kappenii and Umbilicaria antarctica as distinct species. Lichenologist 36, 227–234. Parnmen S, Rangsiruji A, Mongkolsuk P, Boonpragob K, Nutakki A and Lumbsch HT (2012) Using phylogenetic and coalescent methods to under- stand the species diversity in the Cladia aggregata complex (Ascomycota, Lecanorales). PLoS ONE 7, e52245. Piercey-Normore MD (2006) The lichen-forming ascomycete Evernia meso- morpha associated with multiple genotypes of Trebouxia jamesii. New Phytologist 169, 331–344. Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD and Vogler AP (2006) Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Systematic Biology 55, 595–609. Posada D (2008) jModelTest: phylogenetic model averaging.Molecular Biology and Evolution 25, 1253–1256. Puillandre N, Lambert A, Brouillet S and Achaz G (2012) ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Molecular Ecology 21, 1864–1877. Purvis OW, Coppins BJ, Hawksworth DL, James PW and Moore DM (1992) The Lichen Flora of Great Britain and Ireland. London: British Lichen Society. Rambaut A (2009) FigTree v. 1.2.3. [WWW resource] URL http://tree.bio.ed. ac.uk/software/figtree/ [Accessed 5 November 2010]. Rambaut A and Drummond J (2007) Tracer v. 1.5. [WWW resource] URL http://beast.bio.ed.ac.uk/Tracer. Ronquist F and Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19, 1572–1574. Sancho LG, Pintado A and Green TGA (2019) Antarctic studies show lichens to be excellent biomonitors of climate change. Diversity 11, 42. Schmidt HA, Strimmer K, Vingron M and von Haeseler A (2002) TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics 18, 502–504. The Lichenologist 193 https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press http://tree.bio.ed.ac.uk/software/figtree/ http://tree.bio.ed.ac.uk/software/figtree/ http://tree.bio.ed.ac.uk/software/figtree/ http://beast.bio.ed.ac.uk/Tracer http://beast.bio.ed.ac.uk/Tracer https://doi.org/10.1017/S0024282922000147 Spielmann AA and Marcelli MP (2009) Parmotrema s.l. (Parmeliaceae, liche- nized Ascomycota) from Serra Geral slopes in central Rio Grande do Sul State, Brazil. Hoehnea 36, 551–595. Spielmann AA and Marcelli MP (2020) Type studies on Parmotrema (Parmeliaceae, Ascomycota) with salazinic acid. Plant and Fungal Systematics 65, 403–508. Stamatakis A (2006) RAxML-VI-HPC: maximum likelihood-based phylogen- etic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688–2690. Stamatakis A, Hoover P and Rougemont J (2008) A rapid bootstrap algo- rithm for the RAxML webservers. Systematic Biology 75, 758–771. Sujetovienė G (2015) Monitoring lichen as indicators of atmospheric quality. In Upreti DK, Divakar PK, Shukla V and Bajpai R (eds), Recent Advances in Lichenology, Vol. 1. New Delhi: Springer, pp. 87–118. Swinscow TDV and Krog H (1988) Macrolichens of East Africa. London: British Museum of Natural History. Thell A, Crespo A, Divakar PK, Kearnefelt I, Leavitt SD, Lumbsch HT and Seaward MRD (2012) A review of the lichen family Parmeliaceae – history, phylogeny and current taxonomy. Nordic Journal of Botany 30, 641–664. Vilgalys R and Hester M (1990) Rapid genetic identification and mapping of enzymatically amplified DNA from several Cryptococcus species. Journal of Bacteriology 172, 4238–4246. Wang-Yang JR and Lai MJ (1976) Additions and corrections to the lichen flora of Taiwan. Taiwania 21, 226–228. Wei JC (1991) An Enumeration of Lichens in China. Beijing: International Academic Publishers. Widhelm TJ, Egan RS, Bertoletti FR, Asztalos MJ, Kraichak E, Leavitt SD and LumbschHT (2016) Pickingholes in traditional species delimitations: an integra- tive taxonomic reassessment of the Parmotrema perforatum group (Parmeliaceae, Ascomycota). Botanical Journal of the Linnean Society 182, 868–884. Wirtz N, Printzen C and Lumbsch HT (2008) The delimitation of Antarctic and bipolar species of neuropogonoid Usnea (Ascomycota, Lecanorales): a cohesion approach of species recognition for the Usnea perpusilla complex. Mycological Research 112, 472–484. Yoshimura I (1974) Lichen Flora of Japan in Colour. Osaka: Hoikusha. Zhang J, Kapli P, Pavlidis P and Stamatakis A (2013) A general species delimitation method with applications to phylogenetic placements. Bioinformatics 29, 2869–2876. 194 Ayoub Stelate et al. https://doi.org/10.1017/S0024282922000147 Published online by Cambridge University Press https://doi.org/10.1017/S0024282922000147 Resolving the phylogenetic relationship between Parmotrema crinitum and Parmotrema perlatum populations Introduction Material and Methods Taxon sampling Chemistry and morphology DNA extraction, PCR and sequencing Sequence alignments and phylogenetic analyses Candidate species identification Results Morphological and chemical analyses Phylogenetic analysis Identifying candidate species Discussion Acknowledgements References