INTRODUCTION
⌅The diversity of organisms is distributed unequally in the world. Species distributions are influenced by a wide range of factors, including biogeographic history, topography, or climate (Fischer 1960Fischer A.G. 1960. Latitudinal variations in organic diversity. Evolution 14: 64-81. ; Brown 2001Brown J. 2001. Mammals on mountainsides: elevational pattern of diversity. Global Ecology and Biogeography 10: 101-109. ; Rozzi & al. 2008Rozzi R., Armesto J.J., Goffinet B., Buck W., Massardo F., Silander J., Arroyo M.K.T., Russell S., Anderson C.B., Cavieres L.A. & Callicott J.B. 2008. Changing lenses to assess biodiversity: patterns of species richness in sub-Antarctic plants and implications for global conservation. Frontiers in Ecology and the Environment 6: 131-137.). These factors can result in species-rich regions, including hotspots of biodiversity where rare and endemic species occur. One of the 25 biodiversity hotspots of the world is the Mesoamerican region that extends from northern Costa Rica and Nicaragua to central Mexico (Myers & al. 2000Myers N., Mittermeier R.A., Mittermeier C.G., da Fonseca G.A.B. & Kent J. 2000. Biodiversity hotspots for conservation priorities. Nature 403: 853-858. ). The high levels of biodiversity in Mexico have been explained by the country’s geographic position between tropical and temperate regions and its importance as a refugial area for temperate species during Pleistocene glaciations (Mittermeier 1988Mittermeier R. 1988. Primate diversity and the tropical forest. Case studies from Brazil and Madagascar and the importance of the megadiversity countries. In Wilson E.O. (ed.), Biodiversity: 145-154. National Academy Press, Washington.; Ramamoorthy & al. 1993Ramamoorthy T.P., Bye R., Lot A. & Fa J. 1993. Biological Diversity of Mexico: Origins and Distribution. Oxford University Press, Oxford. ; Myers & al. 2000Myers N., Mittermeier R.A., Mittermeier C.G., da Fonseca G.A.B. & Kent J. 2000. Biodiversity hotspots for conservation priorities. Nature 403: 853-858. ). Mexico has a wide range of ecosystems from rainforest to dry deserts. We focus on the xerophytic scrub, which is among the most widespread ecosystems in Mexico, covering about 40% of the country and harbouring high biodiversity. This vegetation is adapted to face aridity and mainly consists of low trees, shrubs and succulent plants (Rzedowski 1978Rzedowski J. 1978. Vegetación de México. Editoral Limusa, México.). Lichen communities in xerophytic scrub frequently consist of crustose lichens, as well as, members of the lichen-forming fungi family Parmeliaceae F.Berchtold & J.Presl. In Parmeliaceae, species of the genera Usnea Dill. ex Adans., Parmotrema A.Massal., Hypotrachyna (Vain.) Hale, and Xanthoparmelia (Vain.) Hale are particularly common (Lücking & al. 2016Lücking R., Villaseñor J.L., Herrera-Campos M.A., Pérez-Pérez R.M., Egan S.R., Esslinger H.L. & Nash III T.H. 2016. Phylogenetic structure of metacommunities in Mexican Parmeliaceae (lichenized Ascomycota: Lecanorales). Bibliotheca Lichenologica 110: 27-54. ).
The genus Xanthoparmelia is the most diverse genus of lichen-forming fungi with about 800 species worldwide (Blanco & al. 2004Blanco O., Crespo A., Elix J.A., Hawksworth D.L. & Lumbsch, H.T. 2004. A molecular phylogeny and a new classification of parmelioid lichens containing Xanthoparmelia-type lichenan (Ascomycota: Lecanorales). Taxon 53: 959-975. ; Thell & al. 2012Thell A., Crespo A., Divakar P.K., Kärnefelt I., Leavitt S.D., Lumbsch H.T. & Seaward M.R.D. 2012. A review of the lichen family Parmeliaceae - history, phylogeny and current taxonomy. Nordic Journal of Botany 30: 641-664. ; Jaklitsch & al. 2016Jaklitsch W.M., Baral H.O., Lücking R. & Lumbsch H.T. 2016. Ascomycota. In Frey W. (ed.), Syllabus of Plant Families - Adolf Engler’s Syllabus der Pflanzenfamilien: 1-150. Stuttgart, Gebr. Borntraeger Verlagsbuchhandlung.). Xanthoparmelia is chemically diverse and the presence or absence of secondary metabolites has been widely used, in conjunction with morphological characters, to delimit species (Hale 1990Hale Jr. M.E. 1990. A synopsis of the lichen genus Xanthoparmelia (Vainio) Hale (Ascomycotina, Parmeliaceae). Smithsonian Contributions to Botany 74: 1-250. ; Elix 1994Elix J.A. 1994. Xanthoparmelia. Flora of Australia 55: 201-308.; Nash & al. 2016Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641.). Xanthoparmelia diversified dramatically during the Miocene (Kraichak & al. 2015Kraichak E., Divakar P.K., Crespo A., Leavitt S.D., Nelsen M.P., Lücking R. & Lumbsch H.T. 2015. A tale of two hyper-diversities: diversification dynamics of the two largest families of lichenized fungi. Scientific Reports 5: e10028. ). Recent studies show that Xanthoparmelia could have originated in the African continent during the early Miocene, eventually spreading to Australia and South Africa -both current centers of diversity for this genus. In contrast, the Holarctic has been more recently colonized and diversification happened mainly during the late Miocene and early Pliocene, resulting in a lower diversity of Xanthoparmelia in this region (Leavitt & al. 2018Leavitt S.D., Kirika P.M., Amo de Paz G., Huang J.P., Hur J.S., Elix J.A., Grewe F., Divakar P.K. & Lumbsch H.T. 2018. Assessing phylogeny and historical biogeography of the largest genus in lichen-forming fungi, Xanthoparmelia (Parmeliaceae, Ascomycota). Lichenologist 50: 299-312.). However, in Mexico there exists a high Xanthoparmelia diversity and a high chemosyndromic variation has been recognized inside of the genus (Culberson & al. 1979Culberson C., Nash T.H. III & Johnson A. 1979. 3-α-hydroxibarbatic acid, a new depside in chemosyndrome of some Xanthoparmeliae with α-orcinol depsides. Bryologist 82: 154-161.). Due to this, deep studies are needed to clarify many of the phylogenetic clades within the group (Nash & al. 2016Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641.). Nevertheless, at the date, there only exists a single molecular study of Xanthoparmelia that included five taxa from xerophytic scrub (Barcenas-Peña & al. 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ), despite the fact that this is the most extensive ecosystem in the country and the preferred habitat for Xanthoparmelia species. Consequently, the xerophytic scrub is an ecosystem that requires a deep molecular study of Xanthoparmelia species (Barcenas-Peña & al. 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ).
In Mexico, Xanthoparmelia species are very abundant due to the presence of numerous exposed rocky substrates, and c. 75 species are currently accepted based mainly on phenotypical characters (Nash & al. 2016Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641.; Barcenas-Peña & al. 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ). Forty six Xanthoparmelia species are recorded from xerophytic scrub vegetation in the country (Lücking & al. 2016Lücking R., Villaseñor J.L., Herrera-Campos M.A., Pérez-Pérez R.M., Egan S.R., Esslinger H.L. & Nash III T.H. 2016. Phylogenetic structure of metacommunities in Mexican Parmeliaceae (lichenized Ascomycota: Lecanorales). Bibliotheca Lichenologica 110: 27-54. ; Nash & al. 2016Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641.; Barcenas-Peña & al. 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ) mainly from the Sonoran Desert in the north of the country, one of the most studied taxonomically regions (morphological and chemically) (Nash & Elix 2004Nash III T.H. & Elix J.A. 2004. Xanthoparmelia. In Nash III T.H., Ryan B.D., Diederich P., Gries C. & Bungartz F. (eds.), Lichen Flora of the Greater Sonoran Desert Region vol. 2: 566-605. Lichens Unlimited, Arizona State University, Tempe.; Nash & al. 2016Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641.). However, the central area of the country still requires attention. In addition, in Mexico there are currently few studies of lichen fungi at the molecular level. For instance, only a work about Xanthoparmelia mexicana (Gyeln.) Hale group from xerophytic scrub in Mexico has been done (Barcenas-Peña & al 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ). Here we studied the diversity of Xanthoparmelia species in xerophytic scrub vegetation of the central part of Mexico using morphology, secondary chemistry, and molecular sequence data to understand the phylogenetic diversity and relationship among the species.
MATERIAL AND METHODS
⌅Study area
⌅All Xanthoparmelia specimens were collected on rocks from different localities throughout arid regions of xerophytic scrub in the central part of Mexico: Zacatecas, Aguas Calientes, San Luis Potosi, Jalisco, Guanajuato, Querétaro, Hidalgo, Estado de Mexico and Mexico City (Fig. 1). Xerophytic scrub occupies approximately 40% of the country’s surface and is distributed from the Baja California Peninsula, the Coastal Plain and the lower Sierra de Sonora. Likewise, it is characteristic of the Altiplano from Chihuahua and Coahuila to Jalisco, Guanajuato, Hidalgo and Estado de Mexico, extending to Puebla and Oaxaca. It also constitutes the vegetation of a part of the northeastern Coastal Plain of Coahuila and Tamaulipas. These areas show average temperatures between 12°C and 26°C and an average annual rainfall of 100 ml to 400 ml (Rzedowski 1978Rzedowski J. 1978. Vegetación de México. Editoral Limusa, México.).
Anatomical studies
⌅For all specimens collected (about 440, deposited at F and MEXU; Appendix 1) the morphology and chemistry were assessed. Morphological characters, such as, shape and size of the thallus and lobules, isidia shape, lower surface and medulla color were studied according to Hale (1990)Hale Jr. M.E. 1990. A synopsis of the lichen genus Xanthoparmelia (Vainio) Hale (Ascomycotina, Parmeliaceae). Smithsonian Contributions to Botany 74: 1-250. and Nash & al. (2016)Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641. using a Zeiss Stemi 2000-C stereoscope. Ascomatal anatomy, ascospore, conidia shape and size were studied using a Zeiss Axioscope. Secondary metabolites were identified using spot test with 10% KOH, KC, C, PD and high-performance thin layer chromatography (HPTLC) using solvent system C following established methods (Culberson & Johnson 1982Culberson C. & Johnson A. 1982. Substitution of methyl tert.-butyl ether for diethyl ether in standardized thin layer chromatographic method for lichen products. Journal of Chromatography B 238: 438-487.; Arup & al. 1993Arup U., Ekman S., Lindblom L. & Mattsson J-E. 1993. High performance thin layer chromatography (HPTLC), an improved technique for screening lichen substances. Lichenologist 25: 61-71.; Lumbsch 2002Lumbsch H.T. 2002. Analysis of phenolic products in lichens for identification and taxonomy. In Kranner I., Beckett R. & Varma A. (eds.), Protocols in Lichenology Culturing, biochemistry, ecophysiology and use in biomonitoring: 281-295. Springer, Berlin.; Orange & al. 2010Orange A., James P.W. & White F.J. 2010. Microchemical methods for the identification of lichens, second edition with additions and corrections. British Lichen Society, London.).
Taxon sampling
⌅Since a large number of specimens had lichenicolous fungi, only selected specimens of each species were used to obtain DNA. A total of 29 specimens of 18 species were selected for molecular analyses from xerophytic scrub vegetation representing the range of morphological and chemical variation, and supplemented with sequences of 20 specimens of 4 species from Puebla and Oaxaca (Mexico) of a previous analysis (Barcenas-Peña & al. 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ) (Appendix 2). In addition, sequences of 138 specimens of 27 species from different parts of the world were downloaded from GenBank (Appendix 2). Four species that have previously been shown to be distantly related to the Xanthoparmelia species from North America were used as outgroups, including X. crespoae Elix, Louwhoff & M.C. Molina, X. filarszkyana (Gyeln.) Hale, X. substrigosa (Hale) Hale and X. lithophiloides (Kurok.) Elix (Leavitt & al. 2018Leavitt S.D., Kirika P.M., Amo de Paz G., Huang J.P., Hur J.S., Elix J.A., Grewe F., Divakar P.K. & Lumbsch H.T. 2018. Assessing phylogeny and historical biogeography of the largest genus in lichen-forming fungi, Xanthoparmelia (Parmeliaceae, Ascomycota). Lichenologist 50: 299-312.). Altogether, a total of 187 specimens of 38 species were included in this study (Appendix 2).
Molecular methods
⌅From the 29 specimens sampled from Mexico, total genomic DNA was extracted from thallus fragments following the manufacturers’ instructions using the ZR Fungal/Bacterial DNA Miniprep Kit (Zymo Research Corp., Irvine, CA). DNA sequences were generated for three markers using the polymerase chain reaction (PCR): the nuclear ribosomal internal transcribed spacer region (ITS), a region of the mitochondrial small subunit rDNA (mtSSU), and a region of the nuclear large subunit rDNA (nuLSU). PCR reactions contained 6.25 µl of MyTaq™ Red DNA Polymerase (Bioline, Taunton, MA, USA), 5.25 µl of H2O, 0.25 µl of forward and reverse primers (10 μM), and 0.5 µl of template DNA (10X), for a total reaction volume of 12.5 µl. The ITS region was amplified using primers ITS1F (Gardes & Bruns 1993Gardes M. & Bruns T. 1993. ITS primers with enhanced specificity for basidiomycetes-Application to the identification of mycorrhizae and rusts. Molecular Ecology 2: 113-118. ) and ITS4 (White & al. 1990White T.J., Bruns T.D., Lee S. & Taylor J. 1990. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In Innis M.A., Gelfand D.H., Sninsky J.J. & White T.J. (eds.), PCR Protocols: 315-322. Academic Press, San Diego.); mtSSU using primers mrSSU1 and mrSSU3R (Zoller & al. 1999Zoller S., Scheidegger C. & Sperisen C. 1999. PCR primers for the amplication of mitochondrial small subunit ribosomal DNA of lichen-forming ascomycetes. Lichenologist 31: 511-516.), and nuLSU rDNA using primers AL2R (Mangold & al. 2008Mangold A., Martín M.P., Lücking R. & Lumbsch H.T. 2008. Molecular phylogeny suggests synonymy of Thelotremataceae within Graphidaceae (Ascomycota: Ostropales). Taxon 57: 476-486.) and LR6 (Vilgalys & Hester 1990Vilgalys R. & Hester M. 1990. Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species. Journal of Bacteriology 172: 4238-4246. ). PCR products were sequenced using the same primers used for amplification and ABI PRISM 3730 DNA Analyzer (Applied Biosystems) at the Pritzker Laboratory for Molecular Systematics and Evolution at The Field Museum, Chicago, Illinois, USA.
Sequence alignment and phylogenetic analysis
⌅ITS, mtSSU and nuLSU sequences were aligned independently using the ‘auto’ option with FFT-NS-i algorithm, in Mafft v7 (Katoh & Standley 2013Katoh K. & Standley D.M. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular Biology and Evolution 30: 772-780. ), with the remaining parameters set to default values. Ambiguous positions of each alignment were removed using options for a “less stringent” selection on Gblocks 0.91b (Castresana 2000Castresana J. 2000. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Molecular Biology and Evolution 17: 540-552. ). SequenceMatrix software (Vaidya & al. 2011Vaidya G., Lohman D.J. & Meier R. 2011. SequenceMatrix: concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics 27: 171-180. ) was used to concatenate all three alignments. Phylogenetic analyses were performed using maximum likelihood (ML) and Bayesian analyses (BA). ML trees were calculated with RAxML-HPC2 on XSEDE 8.2.10 (Stamatakis 2014Stamatakis A. 2014. RAxML Version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30: 1312-1313. ) on the Cipres Science Gateway (Miller & al. 2010Miller M.A., Pfeifer W. & Schwartz T. 2010. CIPRES Science Gateway survey. Webside: http://www.phylo.org [accessed: 24 Mar. 2020]) using GTR+G+I substitution model with 1000 bootstrap pseudoreplicates, with the data partitioned by loci. For the BA, substitution models for each locus were estimated using jModelTest-2.1.9 (Guindon & Gascuel 2003Guindon S. & Gascuel O. 2003. A simple, fast and accurate method to estimate large phylogenies by maximum-likelihood. Systematic Biology 52: 696-704. ; Darriba & al. 2012Darriba D., Taboada G.L., Doallo R. & Posada D. 2012. jModelTest 2: more models, new heuristics and parallel computing. Nature Methods 9: 772. ), which recommended for ITS locus the TIM2ef+I+G model, for mtSSU locus the F81+I model and for the nuLSU locus the TIM2ef+I+G model. Due the TIM2ef substitution models are not implemented in MrBayes were replaced by the GTR model (Ronquist & Huelsenbeck 2003Ronquist F. & Huelsenbeck J.P. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574. ). The proportion of invariable sites (I) and gamma distributed rates (G) defined in jModeltest were conserved in both cases. Two parallel Markov chain Monte Carlo (MCMC) runs were performed in MrBayes 3.2.6 (Huelsenbeck & Ronquist 2001Huelsenbeck J.P. & Ronquist F. 2001. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754-755. ; Ronquist & Huelsenbeck 2003Ronquist F. & Huelsenbeck J.P. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574. ), each using 10,000,000 generations which were sampled every 100 steps. A 50% majority rule consensus tree was generated from the combined sampled trees (149,965) of both runs after discarding the first 25% trees as burn-in. Convergence diagnostic as PSRF was reasonably close to 1.0 for all parameters and average deviation of split frequencies was below 0.01 (Gelman & Rubin 1992Gelman A. & Rubin D.B. 1992. Inference from iterative simulation using multiple sequences. Statistical Science 7: 457-511.). Tree files were visualised with FigTree 1.4.2 (Rambaut 2014Rambaut A. 2014. FigTree 1.2.2. Webside: http://tree.bio.ed.ac.uk/software/fgtree/ [accessed: 24 Mar. 2020].). The ITS, mtSSU and nuLSU sequences are deposited in GenBank (Appendix 2).
RESULTS AND DISCUSSION
⌅Phylogeny
⌅Our phylogenetic analyses recovered nine well supported clades of Xanthoparmelia (Fig. 2). Seven of them (clades 1 to 7) include Mexican specimens of the same species forming monophyletic clusters, the rest of the Mexican specimens are distributed all over the phylogenetic tree. In total our phylogenetic analysis revealed 18 species-level lineages from xerophytic scrublands in Mexico (Fig. 2, Appendix 2). Clades 8 and 9 are X. mexicana groups from USA-Spain and USA respectively. The seven major lineages of Mexican Xanthoparmelia samples from xerophytic scrublands included: 1) X. mexicana s. str., 2) the X. moctezumensis T.H.Nash group, 3) the X. ajoensis (T.H.Nash) Egan group, 4) the X. lavicola (Gyeln.) Hale group, 5) the X. pedregalensis Barcenas Peña, Lumbsch & S.D.Leav. group, 6) the X. hypomelaena (Hale) Hale group, and 7) the X. subramigera (Gyeln.) Hale group.
The 18 species recognized in this analysis have been previously reported in Mexico in morphological and chemical studies (Nash & Elix 2004Nash III T.H. & Elix J.A. 2004. Xanthoparmelia. In Nash III T.H., Ryan B.D., Diederich P., Gries C. & Bungartz F. (eds.), Lichen Flora of the Greater Sonoran Desert Region vol. 2: 566-605. Lichens Unlimited, Arizona State University, Tempe.; Nash & al. 2004Nash III T.H., Herrera-Campos M.A. & Elix J.A. 2004. Preliminary treatment of Xanthoparmelia for Mexico. Symbolae Botanicae Upsaliensis 34: 289-326.; Nash & al. 2016Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641.; Barcenas-Peña & al. 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ). However, many of these species did not have sequences available in Genbank including Xanthoparmelia hypomelaena, X. joranadia (T.H.Nash) Egan, X. neotaractica Hale, X. planilobata (Gyeln.) Hale, X. subtasmanica Elix & T.H.Nash, X. tuckeriana Elix & T.H.Nash and X. tucsonensis (T.H. Nash) Egan (Appendix 2). The ML and BA phylogeny (Fig. 2) show that the phenotype-based taxonomy in the genus requires revision, with numerous species identified using morphological and chemical characters not forming monophyletic groups. In addition, with exception of X. pedregalensis none of these species are type materials from Mexico. In the X. mexicana group from the USA and Spain (clade 8) there are not samples from Mexico. As well, within clade 8 samples of X. dierythra (Hale) Hale, X. lineola (E.C.Berry) Hale, X. mexicana, and X. plittii (Gyeln.) Hale are included and do not form monophyletic clades. These species are currently distinguished based on their reproduction and secondary chemistry (Hale 1990Hale Jr. M.E. 1990. A synopsis of the lichen genus Xanthoparmelia (Vainio) Hale (Ascomycotina, Parmeliaceae). Smithsonian Contributions to Botany 74: 1-250. ). However, the substances used as diagnostic characters (norstictic, salazinic, and stictic acids) belong to the same chemosyndrome, which has, in some cases, been interpreted as intraspecific variation (Lumbsch 1998aLumbsch H.T. 1998a. The use of metabolic data in lichenology at the species and subspecific levels. Lichenologist 30: 357-367. , 1998bLumbsch H.T. 1998b. The taxonomic use of metabolic data in lichen-forming fungi. In Frisvad J.C., Bridge P.D. & Arora D.K. (eds.), Chemical Fungal Taxonomy: 345-387. Marcel Dekker, New York.). Further, one sample agreeing with the current circumscription of the vagrant X. chlorochroa (Tuck.) Hale clustered in this group. In fact, all phenotypes (=currently accepted species) in this clade also clustered in other clades in the phylogenetic tree. This strongly indicates that the characters used for the circumscription of species in this clade need re-evaluation, considering that a recent study from Australian Xanthoparmelia species has observed that closely related samples may exhibit distinct chemical profiles. Additionally, that in the evolution of secondary metabolite composition can be rapid, which may result in convergence between distantly related samples (Autumn & al. 2020Autumn K., Barcenas-Peña A., Kish-Levine S., Huang J-P & Lumbsch T.H. 2020. Repeated colonization between arid and seasonal wet habitats, frequent transition among substrate preferences, and chemical diversity in Western Australian Xanthoparmelia lichens. Frontiers in Ecology and Evolution 8: 129. ). A similar pattern is found in the clade of Mexican X. mexicana specimens (clade 1 of Mexican specimens in Fig. 2), in which samples identified as X. chlorochroa, X. dierythra, X. lineola, and X. mexicana are found and none of them form a monophyletic group. Even though clade 1 is a little below the supported values (bootstrap values above 75%), we still considered it a good example that reflects the need of re-evaluation of species. Additionally, the presence of several well-supported clades within clade 1 suggests that even after the segregation of X. pedregalensis (Barcenas-Peña & al. 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ), X. mexicana in Mexico is not well understood. However, given that the species was described from central Mexico, we regard clade 1 as X. mexicana whereas specimens currently accepted as X. mexicana from other parts of the world could belong to other taxa.
The Xanthoparmelia mexicana group from the USA (clade 9) includes samples of the nonisidiate X. cumberlandia (Gyeln.) Hale, X. maricopensis T.H.Nash & Elix (also containing hyposalazinic acid), and the chemically different X. psoromifera (Hale) Hale (with psoromic acid). Psoromic acid has been found to be inconsistent with monophyletic groups in the distantly related genus Cladonia (Pino-Bodas & al. 2012Pino-Bodas R., Rosa Burgaz A., Martin M.P. & Lumbsch H.T. 2012. Species delimitations in the Cladonia cariosa group (Cladoniaceae, Ascomycota). Lichenologist 44: 121-135. ).
The Xanthoparmelia moctezumensis clade also includes a sample of X. ajoensis. However, the two species contain the closely-related depsides 3-a-hydroxybarbatic or diffractaic acids, respectively. This is consistent with the hypothesis that morphologically identical specimens with closely-related substances often represent variation within a single species (Feige & Lumbsch 1995Feige G.B. & Lumbsch H.T. 1995. Some types of chemical variation in lichens. Cryptogamic Botany 5: 31-35.). However, the majority of specimens phenotypically identified as X. ajoensis are found in a separate clade, the X. ajoensis group (clade 3) with unresolved relationships to the X. moctezumensis group (clade 2). Additional studies with an extended sampling of specimens and genetic markers of these two groups is required to better understand the phylogenetic relationships in this part of the phylogenetic tree.
The Xanthoparmelia lavicola group (clade 4) is represented by three samples of this species that was originally described from central Mexico. It was strongly supported as a monophyletic group but its phylogenetic relationships remain unresolved. Another strongly supported clade in the ML analysis was the X. pedregalensis group (clade 5) that is closely related to X. neotaractica, X. sublaevis (Cout.) Hale, X. lineola and X. coloradoensis (Gyeln.) Hale. X. pedregalensis has been treated in detail elsewhere (Barcenas & al. 2018Barcenas-Peña A., Leavitt S.D., Huang J.P., Grewe F. & Lumbsch H.T. 2018. Phylogenetic study and taxonomic revision of the Xanthoparmelia mexicana group, including the description of a new species (Parmeliaceae, Ascomycota). MycoKeys 40: 13-28. ). The X. hypomelaena group (clade 6) included two Mexican specimens-here again the relationships remain unresolved. All studied samples of X. subramigera from Mexico and Kenya clustered together (clade 7), indicating that this species indeed has a larger distributional range and is a well-delimited taxon.
A number of Mexican samples did not cluster in one of the seven clades discussed above including Xanthoparmelia californica Hale, X. conspersa (Ehrh. ex Ach) Hale, X. coloradoensis, X. cumberlandia, X. joranadia, X. planilobata, X. tuckeriana, and X. tucsonensis, as well as samples of X. lineola, X. neotaractica, and X. subtasmanica, none of them forming monophyletic groups. Samples of some of these species from other localities out of Mexico did not cluster with the Mexican specimens, indicating that the delimitation of these species is poorly understood. Since none of these species’ names are based on type materials from Mexico, they might represent distinct lineages that might require formal description. This, however, has to wait until we have gained a better understanding of delimitation of species in Holarctic species of Xanthoparmelia. Since the Holarctic species of the genus are hypothesized to have originated and diversified relatively recently, species delimitation with multi-gene data sets have been shown to be difficult (Leavitt & al. 2011Leavitt S.D., Johnson L.A., Goward T. & St. Clair L.L. 2011. Species delimitation in taxonomically difficult lichen-forming fungi: an example from morphologically and chemically diverse Xanthoparmelia (Parmeliaceae) in North America. Molecular Phylogenetics and Evolution 60: 317-332. , 2013Leavitt S.D., Lumbsch H.T., Stenroos S. & St Clair L.L. 2013. Pleistocene speciation in North American lichenized fungi and the impact of alternative species circumscriptions and rates of molecular evolution on divergence estimates. Plos One 8: e85240.) and hence reduced genomic data sets, such as RADseq might help elucidate species delimitations (Grewe & al. 2017Grewe F., Huang J.P., Leavitt S.D. & Lumbsch H.T. 2017. Reference-based RADseq resolves robust relationships among closely related species of lichen-forming fungi using metagenomic DNA. Scientific Reports 7: 9884. , 2018Grewe F., Lagostina E., Wu H., Printzen C. & Lumbsch H.T. 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. ). In addition, our collection sampling in Mexico will need to be extended, since the described diversity of Xanthoparmelia species based on phenotypical characters in the Mexican xerophytic scrub vegetation is higher (44 species) (Nash & al. 2016Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641.; Lücking & al. 2009Lücking R., Rivas Plata E., Chaves J.L., Umaña L. & Sipman H.J.M. 2009. How many tropical lichens are there. . . really? Bibliotheca Lichenologica 100: 399-418.) than found in this work (23 species, 18 included in the phylogenetic analyses).
Taxonomy
⌅We identified 23 species from the Mexican xerophytic scrub based on their morphology, mainly presence of isidia and lower surface colour, as well as secondary substances (Table 1, Fig. 3). Both isidiate and not isidiate species were found in almost the same proportion, while species with pale to brown lower surface were more frequent. Additionally, we see a high variety of secondary substances in species morphologically similar (Table 1). Nevertheless, nearly 44 Xanthoparmelia species were found in the Mexican xerophytic scrub mainly from the northern part of Mexico (Nash & Elix 2004Nash III T.H. & Elix J.A. 2004. Xanthoparmelia. In Nash III T.H., Ryan B.D., Diederich P., Gries C. & Bungartz F. (eds.), Lichen Flora of the Greater Sonoran Desert Region vol. 2: 566-605. Lichens Unlimited, Arizona State University, Tempe.; Nash & al. 2004Nash III T.H., Herrera-Campos M.A. & Elix J.A. 2004. Preliminary treatment of Xanthoparmelia for Mexico. Symbolae Botanicae Upsaliensis 34: 289-326.; Nash & al. 2016Nash III T.H., Herrera-Campos M.A., Esslinger T.L. & Elix J.A. 2016. Xanthoparmelia in Mexico. Bibliotheca Lichenologica 110: 621-641.). While in the present study we collected in the central part of Mexico and the diversity of species found is high (23 putative species-level lineages). Additional studies that can include all of the diversity of the country are necessary. Incorporation of molecular analyses using different genetic markers to the morphological and chemical study towards a better species delimitation is also necessary.
Species | Chemistry | Isidiate | Lower surface color |
---|---|---|---|
X. ajoensis (T.H.Nash) Egan [G3] | Diffractaic and barbatic | + | Pale tan or brown |
X. californica Hale * | Norstictic and connorstictic | - | Pale to medium brown |
X. coloradoensis (Gyeln.) Hale * | Consalazinic, norstictic, protocetraric and salazinic | - | Brown |
X. conspersa (Ehrh. ex Ach.) Hale * | Stictic, constictic, cryptostictic and norstictic | + | Black |
X. cumberlandia (Gyeln.) Hale * | Stictic and constictic and norstictic | - | Pale brown or brown |
X. hypomelaena (Hale) Hale [G6] | Fumarprotocetraric | - | Black |
X. jornadia (T.H.Nash) Hale * | Lecanoric | + | Pale brown |
X. lavicola (Gyeln.) Hale [G4] | Psoromic | + | Pale to medium brown |
X. lineola (E.C.Berry) Hale * | Salazinic and consalazinic | - | Pale to medium brown |
X. mexicana (Gyeln.) Hale s. str. [G1] | Salazinic, consalazinic and norstictic | + | Pale to medium brown |
X. moctezumensis T.H.Nash [G2] | 3-α-hydroxybarbatic, barbatic, baeomycesic and squamatic | + | Pale brown |
X. neoconspersa (Gyeln.) Hale | Stictic, constictic, cryptostictic and norstictic | - | Black |
X. neotaractica Hale * | Stictic, norstictic and constictic | - | Pale brown to brown |
X. nigropsoromifera (T.H.Nash) Egan | Psoromic, 2’-O-demethylpsoromic | - | Black |
X. novomexicana (Gyeln.) Hale | Fumarprotocetraric and protocetraric | - | Pale to medium brown |
X. pedregalensis Barcenas Peña, Lumbsch & S.D.Leav. [G5] | Salazinic and norstictic | + | Tan to brown |
X. planilobata (Gyeln.) Hale * | Stictic, constictic, cryptostictic and norstictic | - | Black |
X. plittii (Gyeln.) Hale | Stictic, constictic, cryptostictic and norstictic | + | Pale to dark brown |
X. subramigera (Gyeln.) Hale [G7] | Succinprotocetraric, fumarprotocetraric and protocetraric | + | Pale to medium brown |
X. subtasmanica Elix & T.H.Nash * | Salazinic and consalazinic | - | Black |
X. tinctina (Maheu & A.Gillet) Hale | Salazinic and consalazinic | + | Black |
X. tuckeriana Elix & T.H.Nash * | Fumarprotocetraric and protocetraric | - | Ivory to pale brown |
X. tucsonensis (T.H.Nash) Egan * | Diffractaic and barbatic | - | Pale brown |