Evolution of Portulacineae Marked by Gene Tree Conflict and Gene Family Expansion Associated with Adaptation to Harsh Environments
Total Page:16
File Type:pdf, Size:1020Kb
Supplementary Figures Evolution of Portulacineae marked by gene tree conflict and gene family expansion associated with adaptation to harsh environments Ning Wang, Email: [email protected] Stephen A. Smith, E-mail: [email protected] Dendroscope view Limeaceae_Limeum aethiopicum Montiaceae_Phemeranthus parviflorus Basellaceae_Anredera cordifolia Anacampserotaceae_Anacampseros kurtzii Portulacaceae_Portulaca amilis Cactaceae_Leuenbergeria lychnidiflora Cactaceae_Stenocereus yunckeri Cactaceae_Maihuenia poeppigii Cactaceae_Opuntia bravoana Cactaceae_Pereskia grandifolia Talinaceae_Talinum paniculatum A Didiereaceae_Portulacaria afra PhyloPlot view Limeaceae_Limeum aethiopicum Montiaceae_Phemeranthus parviflorus Basellaceae_Anredera cordifolia Anacampserotaceae_Anacampseros kurtzii 0.008 Portulacaceae_Portulaca amilis 0.992 0.118 Cactaceae_Leuenbergeria lychnidiflora Cactaceae_Stenocereus yunckeri 0.24 0.146 0.76 Cactaceae_Maihuenia poeppigii 0.854 0.882 0.364 Cactaceae_Opuntia bravoana 0.636 Cactaceae_Pereskia grandifolia B Talinaceae_Talinum paniculatum Didiereaceae_Portulacaria afra FIG. S1. The phylogenetic network inferred using MPL method in PhyloNet. Taxa were selected from each plant family based on their gene occupancy statistics. A: network visualized in Dendroscope, and B: the same network with inheritance probabilities between hybridization lineages visualized by PhyloPlot that implemented in PhyloNetworks (Solís-Lemus et al. 2017). Anacampserotaceae Basellaceae Anacampseros A. kurtzii Talinopsis frutescens Anredera cordifolia Basella alba filamentosa Bese 400 4000 4000 3000 3000 200 2000 2000 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Portulacaceae Portulaca amilis P. cryptopetala P. grandiflora P. molokiniensis P. oleracea P. pilosa 300 500 800 800 200 200 300 150 400 400 100 100 0 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Talinaceae P. suffrutescens P. umbraticola Talinum fruticosum T. paniculatum T. portulacifolium 8000 600 200 150 300 4000 1500 100 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Cactaceae Coryphantha Ariocarpus retusus Astrophytum asterias A.myriostigma Copiapoa desertorum maiz-tablasensis 8000 8000 8000 8000 3000 4000 4000 4000 4000 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Echinocereus pectinatus Echinopsis aurea Eriosyce Ferocactus latispinus Grusonia bradtiana wagenknechtii 2000 8000 8000 4000 6000 1000 4000 4000 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 FIG. S2. Distribution of within-taxon synonymous distances (Ks) among paralogs gene pairs based on BLASTP hits. Zoomed-in Ks plots of certain species from Didiereaceae, Montiaceae, Molluginaceae and Outgroups were shown next to corresponding species. Potential WGD has been marked by arrows. Yang* represents Yang et al. (2018). Node 13: the origin of Portulacineae, Node 27: the origin of Didiereoideae, Node 17: the origin of Calyptridium, and Node 19: the MRCA of Calandrinia and Claytonia in Montiaceae. Cactaceae (continue) Gymnocalycium Leuenbergeria bleo Leuenbergeria L. lychnidiflora Lophophora williamsii mihanovichii guamacho 8000 150 3000 3000 3000 4000 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Maihuenia poeppigii Maihueniopsis Matucana aurantiaca Opuntia arenaria Opuntia bravoana conoidea 12000 4000 8000 6000 3000 6000 2000 4000 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Opuntia cochenillifera Pachycereus gatesii Peniocereus Pereskia aculeata P. grandifolia cuixmalensis 4000 300 8000 6000 1500 2000 150 4000 3000 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Pterocactus tuberosus Rhipsalis baccifera Salmiopuntia salmiana Stenocereus yunckeri Stetsonia coryne 8000 8000 2000 6000 4000 4000 4000 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Tacinga lilae Tephrocactus articulatus T. bonnieae Tunilla corrugata 8000 8000 8000 WGD in Node 13 6000 4000 4000 4000 Base Potential WGD 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 FIG. S2. Continued Didiereaceae Alluaudia dumosa Alluaudia humbertii Alluaudia procera Didi Didi Didi Didi Didi Didi 3000 2000 12000 1500 15000 15000 1500 6000 1000 0 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.01 0.4 0.8 Alluaudiopsis marnieriana Decarya madagascariensis Didierea madagascariensis Didi Didi Didi Didi Didi Didi 1200 12000 1000 8000 1000 600 6000 6000 600 4000 0 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.01 0.4 0.8 Didierea trollii Portulacaria afra Portulacaria pygmaea Didi Didi 800 10000 1000 Didi WGD in Node 27 2000 400 4000 WGD in Node 13 0 0 0 0 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Montiaceae Calandrinia grandiflora Calyptridium pygmaeum Calyptridium umbellatum Cagr Cagr Caly Caly Caly 1200 Caly 800 1200 8000 4000 600 6000 400 600 4000 0 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.01 0.4 0.8 Lewisia nevadensis Claytonia nevadensis Claytonia virginica Mont Mont Lewi Lewi Mont 3500 Mont Mont Mont 3000 2000 1500 15000 6000 1000 0 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.5 1.5 2.5 0.01 1.0 2.0 3.0 0.5 1.5 2.5 FIG. S2. Continued Montiaceae (continue) Montia chamissoi Cistanthe grandiflora Phemeranthus parviflorus Mont WGD in Node 13 1200 6000 Caly WGD in Node 17 4000 Mont 3000 Mont WGD in Node 19 600 2000 3000 1500 Cagr Potential WGD 0 0 0 0 Lewi Potential WGD 0.01 1.0 2.0 3.0 0.5 1.5 2.5 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Molluginaceae Mollugo cerviana Mollugo nudicaulis Pharnaceum exiguum Suessenguthiella caespitosa 1200 2000 3000 200 600 100 1000 1500 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Mollugo pentaphylla Mollugo verticillata Glinus dahomensis Moll 500 200 1200 Moll 300 600 6000 Moll Potential WGD 100 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.5 1.5 2.5 0.01 1.0 2.0 3.0 Outgroups Sesuvium Delosperma echinatum Anisomeria littoralis Guapira obtusata Beta vulgaris portulacastrum 600 200 PHYT2 PHYT2 400 3000 300 100 1500 200 0 0 0 0 0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 0.01 1.0 2.0 3.0 Limeum aethiopicum Stegnosperma halimifolium Lime Lime Steg 2000 2000 12000 Lime Potential WGD 15000 Steg 1000 Potential WGD 6000 Steg 1000 0 0 0 0 PHYT2 WGD as in Yang* 0.01 1.0 2.0 3.0 0.01 0.4 0.8 0.01 1.0 2.0 3.0 0.5 1.5 2.5 FIG. S2. Continued FIG. S3. The probability density Ks plots inferred form the within species paralogs and reciprocal orthologs (within the vertical panel). Considering the high number of gene duplication events mapping on the phylogenetic tree (fig. 1), the WGD probably occurred before the speciation of Montia chamissoi and Claytonia, and at approximately the same time as the divergence of Lewisia, but possibly after the divergence of Calandrina grandiflora. Ks values > 0.8 and < 0.01 are not shown.