DOI: 10.5061/dryad.h18931znp
title: Tropical origin, global diversification and dispersal in the pond damselflies (Coenagrionoidea) revealed by a new molecular phylogeny

description: The processes responsible for the formation of Earth’s most conspicuous
 diversity pattern, the latitudinal diversity gradient (LDG), remain
 unexplored for many clades in the Tree of Life. Here, we present a
 densely-sampled and dated molecular phylogeny for the most speciose clade
 of damselflies worldwide (Odonata: Coenagrionoidea), and investigate the
 role of time, macroevolutionary processes and biome-shift dynamics in
 shaping the LDG in this ancient insect superfamily. We used process-based
 biogeographic models to jointly infer ancestral ranges and speciation
 times, and to characterise within-biome dispersal and biome-shift dynamics
 across the cosmopolitan distribution of Coenagrionoidea. We also
 investigated temporal and biome-dependent variation in diversification
 rates. Our results uncover a tropical origin of pond damselflies and
 featherlegs ~ 105 Ma, while highligthing uncertainty of ancestral ranges
 within the tropics in deep time. Even though diversification rates have
 declined since the origin of this clade, global climate change and
 biome-shifts have slowly increased diversity in warm- and cold-temperate
 areas, where lineage turnover rates have been relatively higher. This
 study underscores the importance of biogeographic origin and time to
 diversify as important drivers of the LDG in pond damselflies and their
 relatives, while diversification dynamics have instead resulted in the
 formation of ephemeral species in temperate regions. Biome-shifts,
 although limited by tropical niche conservatism, have been the main factor
 reducing the steepness of the LDG in the last 30 Myr. With ongoing climate
 change and increasing northward range expansions of many damselfly taxa,
 the LDG may become less pronounced. Our results support recent calls to
 unify biogeographic and macroevolutionary approaches to increase our
 understanding of how latitudinal diversity gradients are formed and why
 they vary across time and among taxa.
description: Data collection processing and analyses are described in full in
 the Supplementary Material (PhyloSupporting_SysBio_R1.pdf)
description: All the code used in this study is recorded and described in the
 file COENPhylo_R1.html
description: # Tropical origin, global diversification and dispersal in the pond
 damselflies (Coenagrionoidea) revealed by a new molecular phylogeny This
 README file was generated on 2024-01-29 by Beatriz Willink. GENERAL
 INFORMATION 1. Title of Dataset: Tropical origin, global diversification
 and dispersal in the pond damselflies \(Coenagrionoidea\) revealed by a
 new molecular phylogeny 2. Author Information A. Principal Investigator
 Contact Information Name: Beatriz Willink Institution: Stockholm
 University Address: Stockholm, Sweden Email: bwillink@gmail.com B.
 Associate or Co-investigator Contact Information Name: Jessica L. Ware
 Institution: American Museum of Natural History Address: NY, NY USA Email:
 jware@amnh.org C. Associate or Co-investigator Contact Information Name:
 Erik I. Svensson Institution: Lund University Address: Lund, Sweden Email:
 erik.svensson@biologi.lu.se 3. Date of data collection (single date,
 range, approximate date): 2014\-2021 4. Geographic location of data
 collection: A. Field sites Sweden: Lund Municipality Guyana: CEIBA
 Biological Center, Karanambu Lodge, Iwokrama River Lodge and Research
 Centre Argentina: Iguazu National Park, Ibera Wetlands Cameroon: Mefou
 Park, Nyasoso Village, Mount Kupe, Mount Mbam Minkom, Afamba River, Doume
 River, Nyong River, Mbakaou Village, Dja Faunal Reserve B. Specimen
 collections Florida State Collection of Arthropods (FSCA), Smithsonian
 Institute (SI), Brigham Young University (BYU) New Zealand: New Zealand
 Arthropod Collection (NZAC) Japan: National Institute of Advanced
 Industrial Science and Technology (AIST) Singapore: National University of
 Singapore (NUS) Netherlands: National Museum of Natural History
 Naturalis(RMNH) 5. Information about funding sources that supported the
 collection of the data: A. Swedish Research Council (VR: grant no.
 2016-03356) B. Swedish Foundation for International Cooperation in
 Research and Higher Education C. Erik Philip Sörensens Stiftelse D. Royal
 Physiographic Society of Lund E. Jörgen Lindströms Fund F. Lunds
 Djurskyddsfond SHARING/ACCESS INFORMATION 1. Licenses/restrictions placed
 on the data: CC0 1\.0 Universal \(CC0 1\.0\) Public Domain 2. Links to
 publications that cite or use the data: Willink, B., Ware, J., &
 Svensson, E. I. (2023). Tropical origin, global diversification and
 dispersal in the pond damselflies (Coenagrionoidea) revealed by a new
 molecular phylogeny. bioRxiv. https://doi.org/10.1101/2022.01.21.477207 3.
 Links to other publicly accessible locations of the data:
 https://github\.com/bwillink/Damsel\_Phylo 4. Links/relationships to
 ancillary data sets: None 5. Was data derived from another source? Yes,
 part of the data was derived from literature searches and online resources
 A. If yes, list source(s): DNA sequences: National Center for
 Biotechnology Information (NCBI) GenBank Map files: Global Administrative
 areas (GADM) Species distribution data: See complete list of sources in
 data files Table\_S4.csv and Table\_S4.BibTeX 6. Recommended citation for
 this dataset: Willink, B., Ware, J., & Svensson, E. I. (2023).
 Data from: Tropical origin, global diversification and dispersal in the
 pond damselflies (Coenagrionoidea) revealed by a new molecular phylogeny.
 Dryad Digital Repository. https://doi.org/10.5061/dryad.h18931znp DATA
 & FILE OVERVIEW 1. raw A. 16S\_allspp\_unaligned.fst B.
 COI\_allspp\_unaligned.fst C. D7\_allspp\_unaligned.fst D.
 H3\_allspp\_unaligned.fst E. PMTR\_allspp\_unaligned.fst F. biome\_amb.nex
 G. biome.nex H. present\_states.csv I. Table\_S2.csv J. Table\_S4.csv K.
 Table\_S4.BibTeX 2. processed A. mol \- 16S\_allspp.fst \-
 16S\_allspp\_sorted.fs \- 16S\_allspp\_corrected.fst \-
 16S\_allspp\_trim.fst \- 16S\_spp.fst \- 16S\_spp\_sorted.fst \- 16S.nex
 \- COI\_allspp.fst \- COI\_allspp\_trim.fst \- COI\_spp.fst \-
 COI\_spp\_sorted.fst \- COI.nex \- D7\_allspp.fst \-
 D7\_allspp\_corrected.fst \- D7\_allspp\_trim.fst \- D7\_spp.fst \-
 D7\_spp\_sorted.fst \- D7.nex \- H3\_allspp.fst \- H3\_allspp\_trim.fst \-
 H3\_spp.fst \- H3\_spp\_sorted.fst \- H3.nex \- PMTR\_allspp.best.fas \-
 PMTR\_allspp\_trim.fst \- PMTR\_spp.fst \- PMTR\_spp\_sorted.fst \-
 PMTR\_exon\_pruned.fst \- PMTR.nex \- Coen.mol.namecheck.nex \-
 Coen.mol.final.nex B. nxs \- Coen.mol.no\_parts.nex \- M1.mol.nex \-
 M1.partitions.nex \- M2.mol.nex \- M2.partitions.nex \- M3.mol.nex \-
 M3.partitions.nex \- M4.mol.nex \- M4.partitions.nex \- M5.mol.nex \-
 M5.partitions.nex \- M6.mol.nex \- M6.partitions.nex C. bg\_dating \-
 atlas \* earth25\.short\.atlas\.matrix\.1\.txt \*
 earth25\.short\.atlas\.matrix\.2\.txt \*
 earth25\.short\.atlas\.matrix\.3\.txt \*
 earth25\.short\.atlas\.matrix\.4\.txt \*
 earth25\.short\.atlas\.matrix\.5\.txt \*
 earth25\.short\.atlas\.matrix\.6\\.txt \*
 earth25\.short\.atlas\.matrix\.7\.txt \*
 earth25\.short\.atlas\.matrix\.8\.txt \*
 earth25\.short\.atlas\.matrix\.9\.txt \*
 earth25\.short\.atlas\.matrix\.10\.txt \*
 earth25\.short\.atlas\.matrix\.11\.txt \*
 earth25\.short\.atlas\.matrix\.12\.txt \*
 earth25\.short\.atlas\.matrix\.13\.txt \*
 earth25\.short\.atlas\.matrix\.14\.txt \*
 earth25\.short\.atlas\.matrix\.15\.txt \*
 earth25\.short\.atlas\.matrix\.16\.txt \*
 earth25\.short\.atlas\.matrix\.17\.txt \*
 earth25\.short\.atlas\.matrix\.18\.txt \*
 earth25\.short\.atlas\.matrix\.19\.txt \*
 earth25\.short\.atlas\.matrix\.20\.txt \*
 earth25\.med\.atlas\.matrix\.1\.txt \* earth25\.med\.atlas\.matrix\.2\.txt
 \* earth25\.med\.atlas\.matrix\.3\.txt \*
 earth25\.med\.atlas\.matrix\.4\.txt \* earth25\.med\.atlas\.matrix\.5\.txt
 \* earth25\.med\.atlas\.matrix\.6\.txt \*
 earth25\.med\.atlas\.matrix\.7\.txt \* earth25\.med\.atlas\.matrix\.8\.txt
 \* earth25\.med\.atlas\.matrix\.9\.txt \*
 earth25\.med\.atlas\.matrix\.10\.txt \*
 earth25\.med\.atlas\.matrix\.11\.txt \*
 earth25\.med\.atlas\.matrix\.12\.txt \*
 earth25\.med\.atlas\.matrix\.13\.txt \*
 earth25\.med\.atlas\.matrix\.14\.txt \*
 earth25\.med\.atlas\.matrix\.15\.txt \*
 earth25\.med\.atlas\.matrix\.16\.txt \*
 earth25\.med\.atlas\.matrix\.17\.txt \*
 earth25\.med\.atlas\.matrix\.18\.txt \*
 earth25\.med\.atlas\.matrix\.19\.txt \*
 earth25\.med\.atlas\.matrix\.20\.txt \* earth25\.short\.times\.txt \*
 earth25\.med\.times\.txt \- Coen.bg.nex \- Coen.mol.nex \- Coen.start.tre
 D.adm\_map \- merged.shp \- merged.shx \- merged.bdf \- merged.prj E.
 compare\_tree \- Coenagrionidae.tre \- Core.tre \- Ridge.tre \-
 Platynemididae.tre F. biome\_shift \- atlas\_damselfly \* tropical i.
 0\_tropical.graph.txt ii. 1\_tropical.graph.txt iii. 2\_tropical.graph.txt
 iv. 3\_tropical.graph.txt v. 4\_tropical.graph.txt vi.
 5\_tropical.graph.txt vii. 6\_tropical.graph.txt viii.
 7\_tropical.graph.txt \* warm i. 0\_warm.graph.txt ii. 1\_warm.graph.txt
 iii. 2\_warm.graph.txt iv. 3\_warm.graph.txt v. 4\_warm.graph.txt vi.
 5\_warm.graph.txt vii. 6\_warm.graph.txt viii. 7\_warm.graph.txt \* cold
 i. 0\_cold.graph.txt ii. 1\_cold.graph.txt iii. 2\_cold.graph.txt iv.
 3\_cold.graph.txt v. 4\_cold.graph.txt vi. 5\_cold.graph.txt vii.
 6\_cold.graph.txt viii. 7\_cold.graph.txt \* null i. 0\_null.graph.txt ii.
 1\_null.graph.txt iii. 2\_null.graph.txt iv. 3\_null.graph.txt v.
 4\_null.graph.txt vi. 5\_null.graph.txt vii. 6\_null.graph.txt viii.
 7\_null.graph.txt \* land i. 0\_land.graph.txt ii. 1\_land.graph.txt iii.
 2\_land.graph.txt iv. 3\_land.graph.txt v. 4\_land.graph.txt vi.
 5\_land.graph.txt vii. 6\_land.graph.txt viii. 7\_land.graph.txt \*
 epoch\_names\.txt \* epoch\_times\.txt \* region\_names\.txt \-
 biome\_dat.csv \- coen\_biome.nex G. biogeo \- backbone.tre \-
 G1.strong.1000.trees \- G1.weak.1000.trees \- Node\_data.txt \-
 Node\_data\_weak.txt \- Biogeo\_data.txt \- Biogeo\_data\_weak.txt 3.
 output A. topology \- M3.Uniform.MAP.tree B. bg\_dating \-
 Coen.g0.flat.calib.433570.params.txt \-
 Coen.g0.flat.calib.433570.states.txt \- Coen.g0.flat.calib.433570.trees \-
 Coen.g1.beta.calib.strong.673464.params.txt \-
 Coen.g1.beta.calib.strong.673464.states.txt \-
 Coen.g1.beta.calib.strong.673464.trees \-
 Coen.g1.beta.calib.weak.354635.params.txt \-
 Coen.g1.beta.calib.weak.354635.states.txt \-
 Coen.g1.beta.calib.weak.354635.trees \- G1\_beta\_strong.673464\_MAP.tree
 \- G1\_beta\_weak.354635\_MAP.tree C. node\_dating \-
 Coen.g0.fossil.calib.flat.311333.params.txt \-
 Coen.g0.fossil.calib.flat.311333.trees \-
 Coen.g0.fossil.calib.flat.311333.burned.trees \-
 G0\_flat\_yule.311333\_MAP.tree \-
 Coen.g0.fossil.calib.flat.prior.495364.params.txt \-
 Coen.g0.fossil.calib.flat.prior.495364.trees D. EBED \-
 EBED.strong2020\_EBD\_extinction\_rates.log \-
 EBED.strong2020\_EBD\_extinction\_times.log \-
 EBED.strong2020\_EBD\_speciation\_rates.log \-
 EBED.strong2020\_EBD\_speciation\_times.log \-
 EBED.weak2020\_EBD\_extinction\_rates.log \-
 EBED.weak2020\_EBD\_extinction\_times.log \-
 EBED.weak2020\_EBD\_speciation\_rates.log \-
 EBED.weak2020\_EBD\_speciation\_times.log E. HiSSE \-
 BiomeHiSSE\_strong\_r2model.log \- BiomeHiSSE\_strong\_r2.ase.tre \-
 BiomeHiSSE\_strong\_amb\_r1model.log \-
 BiomeHiSSE\_strong\_amb\_r1.ase.tre \- BiomeHiSSE\_weakmodel.log \-
 BiomeHiSSE\_weak.ase.tre \-
 BiomeHiSSE\_strong\_prior\_r1model\_combined.log \- Biome\_ML \*
 BiomeHiSSE\_strong\_ML\_r1\.out \- Null\_ML \*
 BiomeHiSSE\_strong\_null\_ML\_r1\.out F. biome\_shift \-
 run2.paleo.model.log \- run2.paleo.ase.tre \- run2.paleo.history.tsv 2.
 Relationship between files, if important: files in the "raw"
 directory are required to produced files in the "processed"
 directory, files in both directories can be used to reproduce files in the
 "output" directory\. 3. Additional related data collected that
 was not included in the current data package: None 4. Are there multiple
 versions of the dataset? Yes A. If yes, name of file(s) that was updated:
 16S\_allspp\_unaligned.fst and 16S\_allspp.fst i. Why was the file
 updated? To remove one unused sequence ii. When was the file updated?
 06-09-2023 B. If yes, name of file(s) that was updated: biome.nex and
 biome\_amb.nex i. Why was the file updated? Files were created as
 suggested during peer review ii. When was the file updated? 24-04-2023 and
 12-05-2023, respectively C. If yes, name of file(s) that was updated:
 Node\_data.txt, Node\_data\_weak.txt, Biogeo\_data.txt,
 Biogeo\_data\_weak.txt i. Why was the file updated? To remove two unused
 columns ii. When was the file updated? 07-09-2023, 13-09-2023, 08-09-2023,
 and 13-09-2023, respectively D. If yes, name of file(s) that was updated:
 all files in output directory i. Why was the file updated? these
 intermediate output files were not included in the first data submission,
 but were added after peer review ii. When was the file updated? 07-09-2023
 E. If yes, name of file(s) that was updated: all files in raw directory i.
 Why was the file updated? to correct typos in species names and add
 accession numbers for all sequence data ii. When was the file updated?
 23-01-2024 F. If yes, name of file(s) that was updated:
 G1\_beta\_strong.673464\_MAP.tree, G1\_beta\_weak.354635\_MAP.tree,
 run2.paleo.ase.tre, BiomeHiSSE\_strong\_amb\_r1.ase.tre,
 BiomeHiSSE\_strong\_r2.ase.tre, BiomeHiSSE\_weak.ase.tre,
 BiomeHiSSE\_weak.ase.tre, M3.Uniform.MAP.tree i. Why was the file updated?
 to correct typos in species names ii. When was the file updated?
 12-01-2024
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 16S\_allspp\_unaligned.fst 1. File type:
 fasta 2. Number of sequences: 621 3. Number of taxa: 621 4. Sequence
 description: third domain of the 16S ribosomal DNA 5. Notes: unaligned
 sequences
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: COI\_allspp\_unaligned.fst 1. File type:
 fasta 2. Number of sequences: 573 3. Number of taxa: 573 4. Sequence
 description: 481 bp fragment of the cytochrome oxidase subunit I 5. Notes:
 unaligned sequences
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: D7\_allspp\_unaligned.fst 1. File type:
 fasta 2. Number of sequences: 518 3. Number of taxa: 518 4. Sequence
 description: seventh hypervariable \(divergent\) region of the 28S
 ribosomal DNA 5. Notes: unaligned sequences
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: H3\_allspp\_unaligned.fst 1. File type:
 fasta 2. Number of sequences: 395 3. Number of taxa: 395 4. Sequence
 description: 328 bp coding region of the histone 3 gene 5. Notes:
 unaligned sequences
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: PMTR\_allspp\_unaligned.fst 1. File type:
 fasta 2. Number of sequences: 260 3. Number of taxa: 260 4. Sequence
 description: fragment including coding \(336 bp\) and intronic regions of
 the arginine phenyl\-methyltransferase gene 5. Notes: unaligned sequences
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: biome\_amb.nex 1. File type: nexus 2.
 Number of taxa: 669 3. Number of characters: 1 4. Symbols: 0 = tropical 1
 = warm 2 = cold 5. Missing data code: \- 6. Notes: taxa occupying multiple
 biome states are coded as ambiguous \(e\.g\. "\(0 1\)" represent
 a taxon distributed between tropical and warm biomes\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: biome.nex 1. File type: nexus 2. Number of
 taxa: 669 3. Number of characters: 1 4. Symbols: 0 = tropical 1 = warm 2 =
 cold 5. Missing data code: \- 6. Notes: taxa occupying multiple biome
 states are assigned to the coldest biome they inhabit \(e\.g\. a taxon
 distributed between tropical and warm biomes is coded as "1"\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: present\_states.csv 1. File type: csv 2.
 Number of variables: 11 3. Number of cases/rows: 669 4. Variable List: \*
 Taxon: species name coded as Genus\_species \* area1: first biogeographic
 area inhabited by taxon \* area2: second biogeographic area inhabited by
 taxon \* area3: third biogeographic area inhabited by taxon \* area4:
 fourth biogeographic area inhabited by taxon \* area5: fifth biogeographic
 area inhabited by taxon \* area6: sixth biogeographic area inhabited by
 taxon \* area7: seventh biogeographic area inhabited by taxon \* area8:
 eighth biogeographic area inhabited by taxon \* area9: nineth
 biogeographic area inhabited by taxon \* area10: tenth biogeographic area
 inhabited by taxon 5. Missing data codes: None 6. Specialized formats or
 other abbreviations used: see Table S5 for area abbreviations 7. Notes:
 area2 \- area10 fields used only for taxa inhabiting more than one
 biogeographic area
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Table\_S2.csv Molecular sequence data used
 in this study. 1. File type: csv 2. Number of variables: 9 3. Number of
 cases/rows: 669 4. Variable List: \* Family: family within Coenagrionoidea
 (Platycnemididae or Coenagrionidae) \* Taxon: species name coded as
 Genus\_species \* 16S: sequence identifier for the 16S marker (third
 domain of the 16S ribosomal DNA) \* COI: sequence identifier for the COI
 marker (481 bp fragment of the cytochrome oxidase subunit I) \* H3:
 sequence identifier for the H3 marker (328 bp coding region of the histone
 3 gene) \* D7: sequence identifier for the D7 marker (seventh
 hypervariable (divergent) region of the 28S ribosomal DNA) \* PMTR:
 sequence identifier for the PMTR marker (fragment including coding (336
 bp) and intronic regions of the arginine phenyl-methyltransferase gene) \*
 Collection: museum or other academic collection hosting the specimens used
 for sequencing, or the sequence data if previously available (see 6.) \*
 Identification: initials of researcher who identified the specimen (see
 6.) 5. Missing data codes: \-\- 6. Specialized formats or other
 abbreviations used: AIST = National Institute of Advanced Industrial
 Science and Technology, Japan BYU = Brigham Young University, USA FSCA =
 Florida State Collection of Arthropods, USA LBEE = Laboratory of Biometry
 and Evolutionary Biology, S. Charlat Lab, France LU = Lund University, E.
 I. Svensson, Sweden NUS = National University of Singapore, R. Meier,
 Singapore NCBI = National Center for Biotechnology Information(Genbank),
 USA NZAC = New Zealand Arthropod Collection, New Zealand RMNH = National
 Museum of Natural History Naturalis, The Netherlands Rutgers = Rutgers
 University, J. Ware Lab, USA USDA = United States Department of
 Agriculture, N. Von Ellenrieder, USA AM = Antonia Monteiro BW = Beatriz
 Willink EIS = Erik I. Svensson IS = Issah Seidu JB = Jackie Brown JN =
 Julio Neto KDD = Klaas D. Dijkstra KK = Kawsar Khan MSH = Melissa
 Sanchez-Herrera NVE = Nathalia von Ellenrieder OK = Oleg Kosterin RF = Ryo
 Futahashi RM = Rudolf Meier SB = Seth Bybee SC = Sylvain Charlat SF =
 Sonia Ferreira SR = Samuel Renner VNO = Viktor Nilsson-Ortman
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Table\_S4.csv Geographic distribution data
 for the species included in this study. 1. File type: csv 2. Number of
 variables: 7 3. Number of cases/rows: 669 4. Variable List: \* Family:
 Family within Coenagrionoidea (Platycnemididae or Coenagrionidae) \*
 Taxon: Species name coded as Genus\_species \* Countries: 3-letter country
 codes following ISO 3166 \* Province/Region/States: names of provinces,
 states or regions for countries two span two or more biogeographic areas
 (Australia, Brazil, Canada, China, India, Indonesia, Malaysia, Russia,
 USA) \* Areas: biogeographic areas used for dating analysis \*
 Biome\_region: biome regions used for biome shift and dispersal analysis
 \* Sources: literature and online sources of distribution data (for full
 references see Table\_S4.BibTeX) 5. Missing data codes: NA 6. Specialized
 formats or other abbreviations used: See Table S5 for Areas and Table S7
 for Biome\_region in the Supporting Material
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Table\_S4.BibTeX Literature and online
 sources of distribution data for the species included in this study. 1.
 File type: BibTeX 2. Number of references: 395 3. Notes: file is
 complementary to Table\_S4\.csv
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 16S\_allspp.fst 1. File type: fasta 2.
 Number of sequences: 621 3. Number of taxa: 621 4. Sequence description:
 third domain of the 16S ribosomal DNA 5. Notes: sequences aligned using
 MUSCLE \(github\.com/rcedgar/muscle/releases/tag/v5\.1\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 16S\_allspp\_sorted.fst 1. File type: fasta
 2. Number of sequences: 621 3. Number of taxa: 621 4. Sequence
 description: third domain of the 16S ribosomal DNA 5. Notes: sequences
 aligned and sorted alphabetically for clarity
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 16S\_allspp\_corrected.fst 1. File type:
 fasta 2. Number of sequences: 621 3. Number of taxa: 621 4. Sequence
 description: third domain of the 16S ribosomal DNA 5. Notes: sequence
 alignment corrected after visual verification of secondary structure
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 16S\_allspp\_trim.fst 1. File type: fasta
 2. Number of sequences: 621 3. Number of taxa: 621 4. Sequence
 description: third domain of the 16S ribosomal DNA 5. Notes: sequence
 alignment after trimming primers and long sequences downloaded from NCBI
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 16S\_spp.fst 1. File type: fasta 2. Number
 of sequences: 621 3. Number of taxa: 621 4. Sequence description: third
 domain of the 16S ribosomal DNA 5. Notes: trimmed sequence alignment after
 removing sequence identifier \(for concatetation by taxon name\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 16S\_spp\_sorted.fst 1. File type: fasta 2.
 Number of sequences: 621 3. Number of taxa: 621 4. Sequence description:
 third domain of the 16S ribosomal DNA 5. Notes: trimmed and sorted
 sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 16S.nex 1. File type: nexus 2. Number of
 sequences: 621 3. Number of taxa: 621 4. Number of characters: 317 5.
 Number of partitions: 1 6. Sequence description: third domain of the 16S
 ribosomal DNA 7. Notes: trimmed and sorted sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: COI\_allspp.fst 1. File type: fasta 2.
 Number of sequences: 573 3. Number of taxa: 573 4. Sequence description:
 481 bp fragment of the cytochrome oxidase subunit I 5. Notes: sequences
 aligned using MUSCLE \(github\.com/rcedgar/muscle/releases/tag/v5\.1\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: COI\_allspp\_trim.fst 1. File type: fasta
 2. Number of sequences: 573 3. Number of taxa: 573 4. Sequence
 description: 481 bp fragment of the cytochrome oxidase subunit I 5. Notes:
 sequence alignment after trimming primers and long sequences downloaded
 from NCBI
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: COI\_spp.fst 1. File type: fasta 2. Number
 of sequences: 573 3. Number of taxa: 573 4. Sequence description: 481 bp
 fragment of the cytochrome oxidase subunit I 5. Notes: trimmed sequence
 alignment after removing sequence identifier \(for concatetation by taxon
 name\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: COI\_spp\_sorted.fst 1. File type: fasta 2.
 Number of sequences: 573 3. Number of taxa: 573 4. Sequence description:
 481 bp fragment of the cytochrome oxidase subunit I 5. Notes: trimmed and
 sorted sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: COI.nex 1. File type: nexus 2. Number of
 sequences: 573 3. Number of taxa: 573 4. Number of characters: 481 5.
 Number of partitions: 1 6. Sequence description: 481 bp fragment of the
 cytochrome oxidase subunit I 7. Notes: trimmed and sorted sequence
 alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: D7\_allspp.fst 1. File type: fasta 2.
 Number of sequences: 518 3. Number of taxa: 518 4. Sequence description:
 seventh hypervariable \(divergent\) region of the 28S ribosomal DNA 5.
 Notes: sequences aligned using MUSCLE
 \(github\.com/rcedgar/muscle/releases/tag/v5\.1\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: D7\_allspp\_corrected.fst 1. File type:
 fasta 2. Number of sequences: 518 3. Number of taxa: 518 4. Sequence
 description: seventh hypervariable \(divergent\) region of the 28S
 ribosomal DNA 5. Notes: sequence alignment corrected after visual
 verification of secondary structure
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: D7\_allspp\_trim.fst 1. File type: fasta 2.
 Number of sequences: 518 3. Number of taxa: 518 4. Sequence description:
 seventh hypervariable \(divergent\) region of the 28S ribosomal DNA 5.
 Notes: sequence alignment after trimming primers and long sequences
 downloaded from NCBI
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: D7\_spp.fst 1. File type: fasta 2. Number
 of sequences: 518 3. Number of taxa: 518 4. Sequence description: seventh
 hypervariable \(divergent\) region of the 28S ribosomal DNA 5. Notes:
 trimmed sequence alignment after removing sequence identifier \(for
 concatetation by taxon name\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: D7\_spp\_sorted.fst 1. File type: fasta 2.
 Number of sequences: 518 3. Number of taxa: 518 4. Sequence description:
 seventh hypervariable \(divergent\) region of the 28S ribosomal DNA 5.
 Notes: trimmed and sorted sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: D7.nex 1. File type: nexus 2. Number of
 sequences: 518 3. Number of taxa: 518 4. Number of characters: 591 5.
 Number of partitions: 1 6. Sequence description: seventh hypervariable
 \(divergent\) region of the 28S ribosomal DNA 7. Notes: trimmed and sorted
 sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: H3\_allspp.fst 1. File type: fasta 2.
 Number of sequences: 395 3. Number of taxa: 395 4. Sequence description:
 328 bp coding region of the histone 3 gene 5. Notes: sequences aligned
 using MUSCLE \(github\.com/rcedgar/muscle/releases/tag/v5\.1\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: H3\_allspp\_trim.fst 1. File type: fasta 2.
 Number of sequences: 395 3. Number of taxa: 395 4. Sequence description:
 328 bp coding region of the histone 3 gene 5. Notes: sequence alignment
 after trimming primers and long sequences downloaded from NCBI
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: H3\_spp.fst 1. File type: fasta 2. Number
 of sequences: 395 3. Number of taxa: 395 4. Sequence description: 328 bp
 coding region of the histone 3 gene 5. Notes: trimmed sequence alignment
 after removing sequence identifier \(for concatetation by taxon name\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: H3\_spp\_sorted.fst 1. File type: fasta 2.
 Number of sequences: 395 3. Number of taxa: 395 4. Sequence description:
 328 bp coding region of the histone 3 gene 5. Notes: trimmed and sorted
 sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: H3.nex 1. File type: nexus 2. Number of
 sequences: 395 3. Number of taxa: 395 4. Number of characters: 328 5.
 Number of partitions: 1 6. Sequence description: 328 bp coding region of
 the histone 3 gene 7. Notes: trimmed and sorted sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: PMTR\_allspp\_best.fas 1. File type: fasta
 2. Number of sequences: 260 3. Number of taxa: 260 4. Sequence
 description: fragment including coding \(336 bp\) and intronic regions of
 the arginine phenyl\-methyltransferase gene 5. Notes: sequences aligned
 using PRANK \(http://wasabiapp\.org/software/prank/\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: PMTR\_allspp\_trim.fst 1. File type: fasta
 2. Number of sequences: 260 3. Number of taxa: 260 4. Sequence
 description: fragment including coding \(336 bp\) and intronic regions of
 the arginine phenyl\-methyltransferase gene 5. Notes: sequence alignment
 after trimming primers and long sequences downloaded from NCBI
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: PMTR\_spp.fst 1. File type: fasta 2. Number
 of sequences: 260 3. Number of taxa: 260 4. Sequence description: fragment
 including coding \(336 bp\) and intronic regions of the arginine
 phenyl\-methyltransferase gene 5. Notes: trimmed sequence alignment after
 removing sequence identifier \(for concatetation by taxon name\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: PMTR\_spp\_sorted.fst 1. File type: fasta
 2. Number of sequences: 260 3. Number of taxa: 260 4. Sequence
 description: fragment including coding \(336 bp\) and intronic regions of
 the arginine phenyl\-methyltransferase gene 5. Notes: trimmed and sorted
 sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: PMTR\_spp\_exon\_pruned.fst 1. File type:
 fasta 2. Number of sequences: 260 3. Number of taxa: 260 4. Sequence
 description: fragment including coding \(336 bp\) regions of the arginine
 phenyl\-methyltransferase gene 5. Notes: trimmed and sorted sequence
 alignment without introns to test for codon saturation
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: PMTR.nex 1. File type: nexus 2. Number of
 sequences: 260 3. Number of taxa: 260 4. Number of characters: 750 5.
 Number of partitions: 1 6. Sequence description: fragment including coding
 \(336 bp\) and intronic regions of the arginine phenyl\-methyltransferase
 gene 7. Notes: trimmed and sorted sequence alignment
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.mol.namecheck.nex 1. File type: nexus
 2. Number of sequences: 669 3. Number of taxa: 669 4. Number of
 characters: 2467 5. Number of partitions: 5 6. Sequence description:
 concatenation of COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.mol.final.nex 1. File type: nexus 2.
 Number of sequences: 669 3. Number of taxa: 669 4. Number of characters:
 2467 5. Number of partitions: 5 6. Sequence description: concatenation of
 COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7. Notes: updated
 taxon names in Coen\.mol\.nex to match current taxonomy
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.mol.no\_parts.nex 1. File type: nexus
 2. Number of sequences: 669 3. Number of taxa: 669 4. Number of
 characters: 2467 5. Number of partitions: 1 6. Sequence description:
 concatenation of COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7\.
 Notes: removed partition block for model comparison
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M1.mol.nex 1. File type: nexus 2. Number of
 sequences: 669 3. Number of taxa: 669 4. Number of characters: 2467 5.
 Number of partitions: 6 6. Sequence description: concatenation of
 COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7\. Notes: different
 partitions for exons and introns in PMTR\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M1.partitions.nex 1. File type: nexus 2.
 Number of partitions: 6 3. Notes: partition block for M1\.mol\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M2.mol.nex 1. File type: nexus 2. Number of
 sequences: 669 3. Number of taxa: 669 4. Number of characters: 2467 5.
 Number of partitions: 9 6. Sequence description: concatenation of
 COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7. Notes: different
 partitions for exons and introns in PMTR\.nex, different partitions for
 codon positions 1 and 2 vs 3 in all coding sequences
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M2.partitions.nex 1. File type: nexus 2.
 Number of partitions: 9 3. Notes: partition block for M2\.mol\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M3.mol.nex 1. File type: nexus 2. Number of
 sequences: 669 3. Number of taxa: 669 4. Number of characters: 2467 5.
 Number of partitions: 10 6. Sequence description: concatenation of
 COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7. Notes: different
 partitions for exons and introns in PMTR\.nex, different partitions for
 codon positions 1 and 2 vs 3 in H3\.nex and PMTR\.nex, different
 partitions for all codon positions in COI\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M3.partitions.nex 1. File type: nexus 2.
 Number of partitions: 10 3. Notes: partition block for M3\.mol\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M4.mol.nex 1. File type: nexus 2. Number of
 sequences: 669 3. Number of taxa: 669 4. Number of characters: 2467 5.
 Number of partitions: 11 6. Sequence description: concatenation of
 COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7. Notes: different
 partitions for exons and introns in PMTR\.nex, different partitions for
 codon positions 1 and 2 vs 3 in PMTR\.nex, different partitions for all
 codon positions in COI\.nex and H3\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M4.partitions.nex 1. File type: nexus 2.
 Number of partitions: 11 3. Notes: partition block for M4\.mol\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M5.mol.nex 1. File type: nexus 2. Number of
 sequences: 669 3. Number of taxa: 669 4. Number of characters: 2467 5.
 Number of partitions: 11 6. Sequence description: concatenation of
 COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7. Notes: different
 partitions for exons and introns in PMTR\.nex, different partitions for
 codon positions 1 and 2 vs 3 in H3\.nex, different partitions for all
 codon positions in COI\.nex and PMTR\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M5.partitions.nex 1. File type: nexus 2.
 Number of partitions: 11 3. Notes: partition block for M5\.mol\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M6.mol.nex 1. File type: nexus 2. Number of
 sequences: 669 3. Number of taxa: 669 4. Number of characters: 2467 5.
 Number of partitions: 12 6. Sequence description: concatenation of
 COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7. Notes: different
 partitions for exons and introns in PMTR\.nex, different partitions for
 all codon positions in all coding sequences
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M6.partitions.nex 1. File type: nexus 2.
 Number of partitions: 12 3. Notes: partition block for M6\.mol\.nex
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: earth25.short.atlas.matrix.1.txt -
 earth25.short.atlas.matrix.20.txt 1. File type: plain text 2. Description:
 25 x 25 adjacency matrices for short\-distance dispersal routes\. Each
 file represents a time epoch between 240 Ma \(matrix\.1\) and 0 Ma
 \(matrix\.20\) 3. Symbols: 0 = not adjacent 1 = adjacent 4. Notes: data
 created by Landis, 2017 Syst Biol
 \(https://doi\.org/10\.1093/sysbio/syw040\) and renamed for this study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: earth25.med.atlas.matrix.1.txt -
 earth25.med.atlas.matrix.20.txt 1. File type: plain text 2. Description:
 25 x 25 adjacency matrices for medium\-distance dispersal routes\. Each
 file represents a time epoch between 240 Ma \(matrix\.1\) and 0 Ma
 \(matrix\.20\) 3. Symbols: 0 = not adjacent 1 = adjacent 4. Notes: data
 created by Landis, 2017 Syst Biol
 \(https://doi\.org/10\.1093/sysbio/syw040\) and renamed for this study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: earth25.short.times.txt 1. File type: plain
 text 2. Description: comma\-separated vector of times defining 20 epochs
 between 240 Ma and 0 Ma 3. Notes: data created by Landis, 2017 Syst Biol
 \(https://doi\.org/10\.1093/sysbio/syw040\) and modified for this study\.
 File is complementary to earth25\.short\.atlas\.matrix\.1\.txt \-
 earth25\.short\.atlas\.matrix\.20\.txt
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: earth25.med.times.txt 1. File type: plain
 text 2. Description: comma\-separated vector of times defining 20 epochs
 between 240 Ma and 0 Ma 3. Notes: data created by Landis, 2017 Syst Biol
 \(https://doi\.org/10\.1093/sysbio/syw040\) and modified for this study\.
 File is complementary to earth25\.med\.atlas\.matrix\.1.txt \-
 earth25\.med\.atlas\.matrix\.20\.txt
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.bg.nex 1. File type: nexus 2. Number
 of taxa: 669 3. Number of characters: 1 4. Symbols: 0 = SAm\_N 1 = SAm\_E
 2 = SAm\_S 3 = NAm\_NW 4 = NAm\_NE 5 = NAm\_SE 6 = NAm\_SW 7 = Grn 8 = Eur
 9 = As\_C A = As\_E B = As\_SE C = As\_NE D = Afr\_W E = Afr\_S F = Afr\_E
 G = Afr\_N H = Aus\_W I = Aus\_E J = Ind K = Mdg L = Ant\_W M = Ant\_E N =
 Mly O = NZ 5. Missing data code: \- 6. Specialized formats or other
 abbreviations used: see Table S5 in the Supporting Material 7. Notes: taxa
 occupying multiple biogeographic areas are coded as ambiguous \(e\.g\.
 "\(0 1\)" represent a taxon distributed between SAm\_N and
 SAm\_E\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.mol.nex 1. File type: nexus 2. Number
 of sequences: 669 3. Number of taxa: 669 4. Number of characters: 2467 5.
 Number of partitions: 10 6. Sequence description: concatenation of
 COI\.nex, H3\.nex, PMTR\.nex, 16S\.nex, and D7\.nex 7. Notes: different
 partitions for exons and introns in PMTR\.nex, different partitions for
 codon positions 1 and 2 vs 3 in H3\.nex and PMTR\.nex, different
 partitions for all codon positions in COI\.nex\. File used in bigeographic
 and node dating analyses
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.start.tre 1. File type: newick 2.
 Number of trees: 1 3. Number of taxa: 669 4. Notes: starting tree for
 biogeographic and node dating analyses
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: merged.shp 1. File type: shape 2. Number of
 features: 528 3. Description: map of administrative areas used to record
 distribution ranges and biogeographic areas in this study 4. Notes: data
 downloaded from https://gadm\.org/ and edited for this study\. File
 complementary to merged\.shx, merged\.bdf, and merged\.prj
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: merged.shx 1. File type: shape index 2.
 Number of features: 528 3. Description: map of administrative areas used
 to record distribution ranges and biogeographic areas in this study 4.
 Notes: data downloaded from https://gadm\.org/ and edited for this study\.
 File complementary to merged\.shp, merged\.bdf, and merged\.prj
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: merged.dbf 1. File type: attribute 2.
 Number of attributes: 12 3. Number of cases/rows: 527 4. Attribute List:
 \* ENGTYPE\_1: Type of administrative boundary associated to the second
 level of subdivision (e.g. region, departament, province, etc) \* HASC\_1:
 Index compounded by the initials of the first and second level of
 subdivision (e.g. "AF.BD" for country = Afghanistan and province
 = Badakhshan) \* ID\_0: unique country-level ID (first level of
 subdivision) \* ID\_1: unique ID at the second level of subdivision (e.g.
 region, departament, province, etc) \* ISO: ISO-Alpha 3 code for countries
 \* NAME\_0: Name of feature at first level of subdivision (i.e. country
 name) \* NAME\_1: Name of feature at second level of subdivision (i.e.
 name of region, department, province, etc) \* OBJECTID: unique identifier
 \* SHAPE\_AREA: feature area in square kilometers * 10^-1 \* SHAPE\_LENG:
 always 0, as all features are polygons \* TYPE\_1: Alternative type of
 administrative boundary associated to the second level of subdivision
 (e.g. region, departament, province, etc) \* Area: Biogeographic area
 assing in this study 5. Description: attributes of administrative areas
 used to record distribution ranges and biogeographic areas in this study
 6. Specialized formats or other abbreviations used: see Table S5 in the
 Supporting Material 7. Notes: data downloaded from https://gadm\.org/ and
 edited for this study\. File complementary to merged\.shp, merged\.shx,
 and merged\.prj
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: merged.prj 1. File type: projection
 description 2. Description: map of administrative areas used to record
 distribution ranges and biogeographic areas in this study 3. Notes: data
 downloaded from https://gadm\.org/ and edited for this study\. File
 complementary to merged\.shp, merged\.bdf, and merged\.shx
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coenagrionidae.tre 1. File type: nexus 2.
 Number of trees: 4 3. Number of taxa: 4 4. Notes: alternative topologies
 for the main clades within Coenagrionidae
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Core.tre 1. File type: nexus 2. Number of
 trees: 4 3. Number of taxa: 9 4. Notes: alternative topologies for the
 main clades within 'core' pond damselflies \(family
 Coenagrionidae\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Ridge.tre 1. File type: nexus 2. Number of
 trees: 4 3. Number of taxa: 9 4. Notes: alternative topologies for the
 main clades within 'ridge\-face' pond damselflies \(family
 Coenagrionidae\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Platynemididae.tre 1. File type: nexus 2.
 Number of trees: 4 3. Number of taxa: 5 4. Notes: alternative topologies
 for the main clades within featherleg damselflies \(family
 Platycnemididae\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 0\_tropical.graph.txt -
 7\_tropical.graph.txt 1. File type: plain text 2. Description: 9 x 9
 availability and connectivity matrices for tropical paleobiomes\. Each
 file represents a time epoch between 106 Ma \(0\_tropical\.graph\) and 0
 Ma \(7\_tropical\.graph\) 3. Symbols: 0 = biome regional coverage less
 than 1%, not connected between regions 1 = biome regional coverage between
 1-25%, weakly connected between regions 2 = biome regional coverage
 greater than 25%, strongly connected between regions 4. Notes: data
 created by Landis et al\., 2021 Syst Biol
 \(https://doi\.org/10\.1093/sysbio/syaa045\), edited and expanded for this
 study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 0\_warm.graph.txt - 7\_warm.graph.txt 1.
 File type: plain text 2. Description: 9 x 9 availability and connectivity
 matrices for warm\-temperate paleobiomes\. Each file represents a time
 epoch between 106 Ma \(0\_warm\.graph\) and 0 Ma \(7\_warm\.graph\) 3.
 Symbols: 0 = biome regional coverage less than 1%, not connected between
 regions 1 = biome regional coverage between 1-25%, weakly connected
 between regions 2 = biome regional coverage greater than 25%, strongly
 connected between regions 4. Notes: data created by Landis et al\., 2021
 Syst Biol \(https://doi\.org/10\.1093/sysbio/syaa045\), edited and
 expanded for this study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 0\_cold.graph.txt - 7\_cold.graph.txt 1.
 File type: plain text 2. Description: 9 x 9 availability and connectivity
 matrices for cold\-temperate paleobiomes\. Each file represents a time
 epoch between 106 Ma \(0\_cold\.graph\) and 0 Ma \(7\_cold\.graph\) 3.
 Symbols: 0 = biome regional coverage less than 1%, not connected between
 regions 1 = biome regional coverage between 1-25%, weakly connected
 between regions 2 = biome regional coverage greater than 25%, strongly
 connected between regions 4. Notes: data created by Landis et al\., 2021
 Syst Biol \(https://doi\.org/10\.1093/sysbio/syaa045\), edited and
 expanded for this study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 0\_null.graph.txt - 7\_null.graph.txt 1.
 File type: plain text 2. Description: 9 x 9 availability and connectivity
 matrices for null features\. Each file represents a time epoch between 106
 Ma \(0\_null\.graph\) and 0 Ma \(7\_null\.graph\) 3. Symbols: 2 = null
 features always have regional coverage greater than 25% and are strongly
 connected between regions 4. Notes: data created by Landis et al\., 2021
 Syst Biol \(https://doi\.org/10\.1093/sysbio/syaa045\), edited and
 expanded for this study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: 0\_land.graph.txt - 7\_land.graph.txt 1.
 File type: plain text 2. Description: 9 x 9 connectivity matrices for land
 regions\. Each file represents a time epoch between 106 Ma
 \(0\_land\.graph\) and 0 Ma \(7\_land\.graph\) 3. Symbols: 0 = regions not
 connected by land 1 = regions weakly connected by land 2 = regions
 connected by land 4. Notes: data created by Landis et al\., 2021 Syst Biol
 \(https://doi\.org/10\.1093/sysbio/syaa045\), edited and expanded for this
 study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: epoch\_names.txt 1. File type: csv 2.
 Number of variables: 4 3. Number of cases/rows: 8 4. Variable List: \*
 index: epoch indices as integers from 0 to 7 \* start\_age: start time of
 each epoch in Ma \* end\_age: end time of each epoch in Ma \* name: epoch
 names 5. Missing data codes: None 6. Notes: data created by Landis et
 al\., 2021 Syst Biol \(https://doi\.org/10\.1093/sysbio/syaa045\), edited
 for this study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: epoch\_times.txt 1. File type: plain text
 2. Description: line\-separated vector of end times defining 8 epochs
 between 66 Ma and 0 Ma 3. Notes: data created by Landis et al\., 2021 Syst
 Biol \(https://doi\.org/10\.1093/sysbio/syaa045\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: region\_names.txt 1. File type: csv 2.
 Number of variables: 3 3. Number of cases/rows: 9 4. Variable List: \*
 index: region indices as integers from 1 to 9 \* region\_code:
 abbreviations used for regions \* region\_name: full region names 5.
 Missing data codes: None 6. Notes: data created by Landis et al\., 2021
 Syst Biol \(https://doi\.org/10\.1093/sysbio/syaa045\), edited for this
 study
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: biome\_dat.csv 1. File type: csv 2. Number
 of variables: 8 3. Number of cases/rows: 669 4. Variable List: \* Taxon:
 Species name coded as Genus\_species \* Region\_1: first region inhabited
 by taxon \* Region\_2: second region inhabited by taxon \* Region\_3:
 third region inhabited by taxon \* Region\_4: fourth region inhabited by
 taxon \* Biome\_1: first biome inhabited by taxon \* Biome\_2: second
 biome inhabited by taxon \* Biome\_3: third biome inhabited by taxon 5.
 Missing data codes: None 6. Specialized formats or other abbreviations
 used: see Table S7 in the Supporting Material for region abbreviations T =
 Tropical W = Warm-Temperate C = Cold-Temperate 7. Notes: Region\_2 \-
 Region\_4 fields used only for taxa inhabiting more than one region,
 Biome\_2 \- Biome\_3 fields used only for taxa inhabiting more than one
 biome
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: coen\_biome.nex 1. File type: nexus 2.
 Number of taxa: 669 3. Number of characters: 1 4. Symbols: A =
 Tropical+Aus B = Tropical+SEAs C = Tropical+EAs D = Tropical+Eur E =
 Tropical+Ind F = Tropical+Afr G = Tropical+NAm H = Tropical+CAm I =
 Tropical+SAm J = Warm+Aus K = Warm+SEAs L = Warm+EAs M = Warm+Eur N =
 Warm+Ind O = Warm+Afr P = Warm+NAm Q = Warm+CAm R = Warm+SAm S = Cold+Aus
 T = Cold+SEAs U = Cold+EAs V = Cold+Eur W = Cold+Ind X = Cold+Afr Y =
 Cold+NAm 0 = Cold+CAm Z = Cold+SAm 5. Missing data code: \- 6. Specialized
 formats or other abbreviations used: see Table S7 in the Supporting
 Material 7. Notes: taxa occupying multiple biome\-region states are coded
 as ambiguous \(e\.g\. "\(A B\)" represent a taxon distributed
 between Tropical\+Aus and Tropical\+SEAs\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: backbone.tre 1. File type: newick 2. Number
 of trees: 1 3. Number of taxa: 3 4. Description: backbone topology for
 Coeangrionoidea
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: G1.strong.1000.trees 1. File type: newick
 2. Number of trees: 1000 3. Number of taxa: 669 4. Description: randomly
 sampled posterior trees from the strongly informed biogeographic dating
 analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: G1.weak.1000.trees 1. File type: newick 2.
 Number of trees: 1000 3. Number of taxa: 669 4. Description: randomly
 sampled posterior trees from the weakly informed biogeographic dating
 analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Node\_dat.txt 1. File type: space separated
 value 2. Number of variables: 4 3. Number of cases/rows: 2674000 4.
 Variable List: \* tree: posterior tree index (1-1000) \* node: node index
 in R (1-1337) \* age: node age in Ma \* index: node index in RevBayes
 (1-1337) 5. Missing data codes: None 6. Specialized formats or other
 abbreviations used: None 7. Notes: used to match node indices from R and
 Revbayes in a posterior sample of trees from the strongly informed
 biogeographic dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Node\_dat\_weak.txt 1. File type: space
 separated value 2. Number of variables: 4 3. Number of cases/rows: 2674000
 4. Variable List: \* tree: posterior tree index (1-1000) \* node: node
 index in R (1-1337) \* age: node age in Ma \* index: node index in
 RevBayes (1-1337) 5. Missing data codes: None 6. Specialized formats or
 other abbreviations used: None 7. Notes: used to match node indices from R
 and Revbayes in a posterior sample of trees from the weakly informed
 biogeographic dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Biogeo\_dat.txt 1. File type: space
 separated value 2. Number of variables: 6 3. Number of cases/rows: 1337001
 4. Variable List: \* tree: posterior tree index (1-1000) \* node: node
 index in R (1-1337) \* age: node age in Ma \* index: node index in
 RevBayes (1-1337) \* start: start state of branch before node
 (biogeographic region) \* end: end state of branch before node
 (biogeographic region) 5. Missing data codes: None 6. Specialized formats
 or other abbreviations used: see Table S5 in the Supporting Material for
 start and end state codes 7. Notes: ancestral states sampled from the
 strongly informed biogeographic dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Biogeo\_dat\_weak.txt 1. File type: space
 separated value 2. Number of variables: 6 3. Number of cases/rows: 1337001
 4. Variable List: \* tree: posterior tree index (1-1000) \* node: node
 index in R (1-1337) \* age: node age in Ma \* index: node index in
 RevBayes (1-1337) \* start: start state of branch before node
 (biogeographic region) \* end: end state of branch before node
 (biogeographic region) 5. Missing data codes: None 6. Specialized formats
 or other abbreviations used: see Table S5 in the Supporting Material for
 start and end state symbols 7. Notes: ancestral states sampled from the
 strongly informed biogeographic dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: M3.Uniform.MAP.tree 1. File type: nexus 2.
 Number of trees: 1 3. Number of taxa: 669 4. Description: maximum a
 posteriori \(MAP\) tree of the topology inference
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g0.flat.calib.433570.params.txt 1.
 File type: tsv 2. Number of variables: 10 3. Number of cases/rows: 3001 4.
 Variable List: \* Iteration: MCMC iteration index \* Posterior: posterior
 probability of sampled parameters \* Likelihood: likelihood of sampled
 parameters \* Prior: prior probability of sampled parameters \* b\_long:
 relative weight of long-distance dispersal \* b\_med: relative weight of
 medium-distance dispersal \* b\_short: relative weight of short-distance
 dispersal \* height\_bg: expected number of disperal events (biogeographic
 clock * root age) \* height\_mol: expected number of substitutions per
 site (molecular clock * root age) \* root\_age: age (in Ma) of the most
 recent common ancestor of Coenagrionoidea 5. Missing data codes: None 6.
 Specialized formats or other abbreviations used: None 7. Notes: parameters
 sampled from the biogeographic dating analysis ignoring empirical
 paleogeography
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g0.flat.calib.433570.states.txt 1.
 File type: tsv 2. Number of variables: 2675 3. Number of cases/rows: 3001
 4. Variable List: \* Iteration: MCMC iteration index \* start\_i: starting
 state (biogeographic area) of branch "i" of a total of 1337
 branches \* end\_i: ending state (biogeographic area) of branch
 "i" of a total of 1337 branches 5. Missing data codes: None 6.
 Specialized formats or other abbreviations used: see Table S5 in the
 Supporting Material for start and end state codes 7. Notes: ancestral
 states sampled at the start and end of each branch from the biogeographic
 dating analysis ignoring empirical paleogeography
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g0.flat.calib.433570.trees 1. File
 type: tsv 2. Number of variables: 5 3. Number of cases/rows: 3001 4.
 Variable List: \* Iteration: MCMC iteration index \* Posterior: posterior
 probability of sampled parameters \* Likelihood: likelihood of sampled
 parameters \* Prior: prior probability of sampled parameters \* tau:
 newick-formatted posterior tree sampled 5. Missing data codes: None 6.
 Specialized formats or other abbreviations used: none 7. Notes: parameters
 sampled from the biogeographic dating analysis ignoring empirical
 paleogeography
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g1.beta.calib.strong.673464.params.txt
 1. File type: tsv 2. Number of variables: 10 3. Number of cases/rows: 3001
 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 b\_long: relative weight of long-distance dispersal \* b\_med: relative
 weight of medium-distance dispersal \* b\_short: relative weight of
 short-distance dispersal \* height\_bg: expected number of disperal events
 (biogeographic clock * root age) \* height\_mol: expected number of
 substitutions per site (molecular clock * root age) \* root\_age: age (in
 Ma) of the most recent common ancestor of Coenagrionoidea 5. Missing data
 codes: None 6. Specialized formats or other abbreviations used: None 7.
 Notes: parameters sampled from the strongly informed biogeographic dating
 analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g1.beta.calib.strong.673464.states.txt
 1. File type: tsv 2. Number of variables: 2675 3. Number of cases/rows:
 3001 4. Variable List: \* Iteration: MCMC iteration index \* start\_i:
 starting state (biogeographic area) of branch "i" of a total of
 1337 branches \* end\_i: ending state (biogeographic area) of branch
 "i" of a total of 1337 branches 5. Missing data codes: None 6.
 Specialized formats or other abbreviations used: see Table S5 in the
 Supporting Material for start and end state codes 7. Notes: ancestral
 states sampled at the start and end of each branch from the strongly
 informed biogeographic dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g1.beta.calib.strong.673464.trees 1.
 File type: tsv 2. Number of variables: 5 3. Number of cases/rows: 3001 4.
 Variable List: \* Iteration: MCMC iteration index \* Posterior: posterior
 probability of sampled parameters \* Likelihood: likelihood of sampled
 parameters \* Prior: prior probability of sampled parameters \* tau:
 newick-formatted posterior tree sampled 5. Missing data codes: None 6.
 Specialized formats or other abbreviations used: none 7. Notes: parameters
 sampled from the strongly informed biogeographic dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g1.beta.calib.weak.354635.params.txt
 1. File type: tsv 2. Number of variables: 10 3. Number of cases/rows: 3001
 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 b\_long: relative weight of long-distance dispersal \* b\_med: relative
 weight of medium-distance dispersal \* b\_short: relative weight of
 short-distance dispersal \* height\_bg: expected number of disperal events
 (biogeographic clock * root age) \* height\_mol: expected number of
 substitutions per site (molecular clock * root age) \* root\_age: age (in
 Ma) of the most recent common ancestor of Coenagrionoidea 5. Missing data
 codes: None 6. Specialized formats or other abbreviations used: None 7.
 Notes: parameters sampled from the weakly informed biogeographic dating
 analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g1.beta.calib.weak.354635.states.txt
 1. File type: tsv 2. Number of variables: 2675 3. Number of cases/rows:
 3001 4. Variable List: \* Iteration: MCMC iteration index \* start\_i:
 starting state (biogeographic area) of branch "i" of a total of
 1337 branches \* end\_i: ending state (biogeographic area) of branch
 "i" of a total of 1337 branches 5. Missing data codes: None 6.
 Specialized formats or other abbreviations used: see Table S5 in the
 Supporting Material for start and end state codes 7. Notes: ancestral
 states sampled at the start and end of each branch from the weakly
 informed biogeographic dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g1.beta.calib.weak.354635.trees 1.
 File type: tsv 2. Number of variables: 5 3. Number of cases/rows: 3001 4.
 Variable List: \* Iteration: MCMC iteration index \* Posterior: posterior
 probability of sampled parameters \* Likelihood: likelihood of sampled
 parameters \* Prior: prior probability of sampled parameters \* tau:
 newick-formatted posterior tree sampled 5. Missing data codes: None 6.
 Specialized formats or other abbreviations used: none 7. Notes: parameters
 sampled from the weakly informed biogeographic dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: G1\_beta\_strong.673464\_MAP.tree 1. File
 type: nexus 2. Number of trees: 1 3. Number of taxa: 669 4. Description:
 maximum a posteriori \(MAP\) tree of the strongly informed biogeographic
 dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: G1\_beta\_weak.354635\_MAP.tree 1. File
 type: nexus 2. Number of trees: 1 3. Number of taxa: 669 4. Description:
 maximum a posteriori \(MAP\) tree of the weakly informed biogeographic
 dating analysis
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g0.fossil.calib.flat.311333.params.txt
 1. File type: tsv 2. Number of variables: 12 3. Number of cases/rows:
 60001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 age\_Argia: age (in Ma) of the most recent ancestor of the genus Argia \*
 age\_Ischnura: age (in Ma) of the most recent ancestor of the genus
 Ischnura \* age\_Nehalennia: age (in Ma) of the most recent ancestor of
 the genus Nehalennia \* age\_Platycnemididae: age (in Ma) of the most
 recent ancestor of featherlegs (family Platycnemididae) \*
 extinction\_rate: extinction rate of the pure-birth process (always 0) \*
 height\_mol: expected number of substitutions per site (molecular clock \*
 root age) \* root\_age: age (in Ma) of the most recent common ancestor of
 Coenagrionoidea \* speciation\_rate: speciation rate of the pure-birth
 process 5. Missing data codes: None 6. Specialized formats or other
 abbreviations used: None 7. Notes: parameters sampled from the node dating
 analysis using fossil calibrations
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: Coen.g0.fossil.calib.flat.311333.trees 1.
 File type: tsv 2. Number of variables: 5 3. Number of cases/rows: 12001 4.
 Variable List: \* Iteration: MCMC iteration index \* Posterior: posterior
 probability of sampled parameters \* Likelihood: likelihood of sampled
 parameters \* Prior: prior probability of sampled parameters \* tau:
 newick-formatted posterior tree sampled 5. Missing data codes: None 6.
 Specialized formats or other abbreviations used: none 7. Notes: parameters
 sampled from the node dating analysis using fossil calibrations
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR:
 Coen.g0.fossil.calib.flat.311333.burned.trees 1. File type: tsv 2. Number
 of variables: 5 3. Number of cases/rows: 10001 4. Variable List: \*
 Iteration: MCMC iteration index \* Posterior: posterior probability of
 sampled parameters \* Likelihood: likelihood of sampled parameters \*
 Prior: prior probability of sampled parameters \* tau: newick-formatted
 posterior tree sampled 5. Missing data codes: None 6. Specialized formats
 or other abbreviations used: none 7. Notes: parameters sampled from the
 node dating analysis using fossil calibrations after discarding the first
 2000 trees
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: G0\_flat\_yule.311333\_MAP.tree 1. File
 type: nexus 2. Number of trees: 1 3. Number of taxa: 669 4. Description:
 maximum a posteriori \(MAP\) tree of the node dating analysis using fossil
 calibrations
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR:
 Coen.g0.fossil.calib.flat.prior.495364.params.txt 1. File type: tsv 2.
 Number of variables: 11 3. Number of cases/rows: 30001 4. Variable List:
 \* Iteration: MCMC iteration index \* Posterior: posterior probability of
 sampled parameters \* Likelihood: likelihood of sampled parameters \*
 Prior: prior probability of sampled parameters \* age\_Argia: age (in Ma)
 of the most recent ancestor of the genus Argia \* age\_Ischnura: age (in
 Ma) of the most recent ancestor of the genus Ischnura \* age\_Nehalennia:
 age (in Ma) of the most recent ancestor of the genus Nehalennia \*
 extinction\_rate: extinction rate of the pure-birth process (always 0) \*
 height\_mol: expected number of substitutions per site (molecular clock *
 root age) \* root\_age: age (in Ma) of the most recent common ancestor of
 Coenagrionoidea \* speciation\_rate: speciation rate of the pure-birth
 process 5. Missing data codes: None 6. Specialized formats or other
 abbreviations used: None 7. Notes: parameters sampled from the node dating
 analysis using fossil calibrations and ran under the prior \(i\.e\.
 without data\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR:
 Coen.g0.fossil.calib.flat.prior.495364.trees 1. File type: tsv 2. Number
 of variables: 5 3. Number of cases/rows: 6001 4. Variable List: \*
 Iteration: MCMC iteration index \* Posterior: posterior probability of
 sampled parameters \* Likelihood: likelihood of sampled parameters \*
 Prior: prior probability of sampled parameters \* tau: newick-formatted
 posterior tree sampled 5. Missing data codes: None 6. Specialized formats
 or other abbreviations used: none 7. Notes: parameters sampled from the
 node dating analysis using fossil calibrations and ran under the prior
 \(i\.e\. without data\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: EBED.strong2020\_EBD\_extinction\_rates.log
 1. File type: tsv 2. Number of variables: 24 3. Number of cases/rows:
 10001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 extinction[i]: extinction rate of the episodic birth-death process during
 time interval "i" of a total of 20 time intervals 5. Missing
 data codes: None 6. Specialized formats or other abbreviations used: None
 7. Notes: parameters sampled from the episodic birth\-death model over 20
 time intervals and using the maximum a posteriori tree of the strongly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: EBED.strong2020\_EBD\_extinction\_times.log
 1. File type: tsv 2. Number of variables: 23 3. Number of cases/rows:
 10001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 interval\_times[i]: end times (in Ma) of time interval "i" of a
 total of 19 time intervals before the present time interval 5. Missing
 data codes: None 6. Specialized formats or other abbreviations used: None
 7. Notes: parameters sampled from the episodic birth\-death model over 20
 time intervals and using the maximum a posteriori tree of the strongly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: EBED.strong2020\_EBD\_speciation\_rates.log
 1. File type: tsv 2. Number of variables: 24 3. Number of cases/rows:
 10001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 speciation[i]: speciation rate of the episodic birth-death process during
 time interval "i" of a total of 20 time intervals 5. Missing
 data codes: None 6. Specialized formats or other abbreviations used: None
 7. Notes: parameters sampled from the episodic birth\-death model over 20
 time intervals and using the maximum a posteriori tree of the strongly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: EBED.strong2020\_EBD\_speciation\_times.log
 1. File type: tsv 2. Number of variables: 23 3. Number of cases/rows:
 10001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 interval\_times[i]: end times (in Ma) of time interval "i" of a
 total of 19 time intervals before the present time interval 5. Missing
 data codes: None 6. Specialized formats or other abbreviations used: None
 7. Notes: parameters sampled from the episodic birth\-death model over 20
 time intervals and using the maximum a posteriori tree of the strongly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: EBED.weak2020\_EBD\_extinction\_rates.log
 1. File type: tsv 2. Number of variables: 24 3. Number of cases/rows:
 10001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 extinction[i]: extinction rate of the episodic birth-death process during
 time interval "i" of a total of 20 time intervals 5. Missing
 data codes: None 6. Specialized formats or other abbreviations used: None
 7. Notes: parameters sampled from the episodic birth\-death model over 20
 time intervals and using the maximum a posteriori tree of the weakly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: EBED.weak2020\_EBD\_extinction\_times.log
 1. File type: tsv 2. Number of variables: 23 3. Number of cases/rows:
 10001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 interval\_times[i]: end times (in Ma) of time interval "i" of a
 total of 19 time intervals before the present time interval 5. Missing
 data codes: None 6. Specialized formats or other abbreviations used: None
 7. Notes: parameters sampled from the episodic birth\-death model over 20
 time intervals and using the maximum a posteriori tree of the weakly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: EBED.weak2020\_EBD\_speciation\_rates.log
 1. File type: tsv 2. Number of variables: 24 3. Number of cases/rows:
 10001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 speciation[i]: speciation rate of the episodic birth-death process during
 time interval "i" of a total of 20 time intervals 5. Missing
 data codes: None 6. Specialized formats or other abbreviations used: None
 7. Notes: parameters sampled from the episodic birth\-death model over 20
 time intervals and using the maximum a posteriori tree of the weakly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: EBED.weak2020\_EBD\_speciation\_times.log
 1. File type: tsv 2. Number of variables: 23 3. Number of cases/rows:
 10001 4. Variable List: \* Iteration: MCMC iteration index \* Posterior:
 posterior probability of sampled parameters \* Likelihood: likelihood of
 sampled parameters \* Prior: prior probability of sampled parameters \*
 interval\_times[i]: end times (in Ma) of time interval "i" of a
 total of 19 time intervals before the present time interval 5. Missing
 data codes: None 6. Specialized formats or other abbreviations used: None
 7. Notes: parameters sampled from the episodic birth\-death model over 20
 time intervals and using the maximum a posteriori tree of the weakly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: BiomeHiSSE\_strong\_r2model.log 1. File
 type: tsv 2. Number of variables: 46 3. Number of cases/rows: 10001 4.
 Variable List: \* Iteration: MCMC iteration index \* Replicate\_ID: index
 for independent MCMC run (0 or 1) \* Posterior: posterior probability of
 sampled parameters \* Likelihood: likelihood of sampled parameters \*
 Prior: prior probability of sampled parameters \* extinction[1]:
 extinction rate for Tropical biome in hidden state 1 \* extinction[2]:
 extinction rate for Warm-Temperate biome in hidden state 1 \*
 extinction[3]: extinction rate for Cold-Temperate biome in hidden state 1
 \* extinction[4]: extinction rate for Tropical biome in hidden state 2 \*
 extinction[5]: extinction rate for Warm-Temperate biome in hidden state 2
 \* extinction[6]: extinction rate for Cold-Temperate biome in hidden state
 2 \* extinction\_alpha[1]: turnover rate for Tropical biome in hidden
 state 1 \* extinction\_alpha[2]: turnover rate for Warm-Temperate biome in
 hidden state 1 \* extinction\_alpha[3]: turnover rate for Cold-Temperate
 biome in hidden state 1 \* extinction\_beta: additive effect of hidden
 state 2 on turnover rate \* hidden\_rate: transition rate between hidden
 states \* Q[1]: transition rate from Tropical to Warm-Temperate biome \*
 Q[2]: transition rate from Tropical to Cold-Temperate biome \* Q[3]:
 transition rate from Warm-Temperate to Tropical biome \* Q[4]: transition
 rate from Warm-Temperate to Cold-Temperate biome \* Q[5]: transition rate
 from Cold-Temperate to Tropical biome \* Q[6]: transition rate from
 Cold-Temperate to Warm-Temperate biome \* R[1]: transition rate from
 hidden state 1 to hidden state 2 \* R[2]: transition rate from hidden
 state 1 to hidden state 2 \* rate\_01: transition rate from Tropical to
 Warm-Temperate biome \* rate\_02: transition rate from Tropical to
 Cold-Temperate biome \* rate\_10: transition rate from Warm-Temperate to
 Tropical biome \* rate\_12: transition rate from Warm-Temperate to
 Cold-Temperate biome \* rate\_02: transition rate from Cold-Temperate to
 Tropical biome \* rate\_21: transition rate from Cold-Temperate to
 Warm-Temperate biome \* rate\_category\_prior[1]: root frequency of
 Tropical-hidden state 1 \* rate\_category\_prior[2]: root frequency of
 Warm-Temperate-hidden state 1 \* rate\_category\_prior[3]: root frequency
 of Cold-Temperate-hidden state 1 \* rate\_category\_prior[4]: root
 frequency of Tropical-hidden state 2 \* rate\_category\_prior[5]: root
 frequency of Warm-Temperate-hidden state 2 \* rate\_category\_prior[6]:
 root frequency of Cold-Temperate-hidden state 2 \* speciation[1]:
 speciation rate for Tropical biome in hidden state 1 \* speciation[2]:
 speciation rate for Warm-Temperate biome in hidden state 1 \*
 speciation[3]: speciation rate for Cold-Temperate biome in hidden state 1
 \* speciation[4]: speciation rate for Tropical biome in hidden state 2 \*
 speciation[5]: speciation rate for Warm-Temperate biome in hidden state 2
 \* speciation[6]: speciation rate for Cold-Temperate biome in hidden state
 2 \* speciation\_alpha[1]: diversification rate for Tropical biome in
 hidden state 1 \* speciation\_alpha[2]: diversification rate for
 Warm-Temperate biome in hidden state 1 \* speciation\_alpha[3]:
 diversification rate for Cold-Temperate biome in hidden state 1 \*
 speciation\_beta: multiplicative effect of hidden state 2 on
 diversification rate 5. Missing data codes: None 6. Specialized formats or
 other abbreviations used: None 7. Notes: parameters sampled from a Hidden
 State Speciation and Extinction model, using the maximum a posteriori tree
 of the strongly informed biogeographic dating analysis as input, and
 assigning taxa occupying multiple biome states to the coldest biome they
 inhabit
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: BiomeHiSSE\_strong\_r2.ase.tre 1. File
 type: nexus 2. Number of trees: 1 3. Number of taxa: 669 4. Description:
 tree annotated with ancestral state reconstructions from Hidden State
 Speciation and Extinction model, using the maximum a posteriori tree of
 the strongly informed biogeographic dating analysis as input, and
 assigning taxa occupying multiple biome states to the coldest biome they
 inhabit
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: BiomeHiSSE\_strong\_amb\_r1model.log 1.
 File type: tsv 2. Number of variables: 46 3. Number of cases/rows: 10001
 4. Variable List: \* Iteration: MCMC iteration index \* Replicate\_ID:
 index for independent MCMC run (0 or 1) \* Posterior: posterior
 probability of sampled parameters \* Likelihood: likelihood of sampled
 parameters \* Prior: prior probability of sampled parameters \*
 extinction[1]: extinction rate for Tropical biome in hidden state 1 \*
 extinction[2]: extinction rate for Warm-Temperate biome in hidden state 1
 \* extinction[3]: extinction rate for Cold-Temperate biome in hidden state
 1 \* extinction[4]: extinction rate for Tropical biome in hidden state 2
 \* extinction[5]: extinction rate for Warm-Temperate biome in hidden state
 2 \* extinction[6]: extinction rate for Cold-Temperate biome in hidden
 state 2 \* extinction\_alpha[1]: turnover rate for Tropical biome in
 hidden state 1 \* extinction\_alpha[2]: turnover rate for Warm-Temperate
 biome in hidden state 1 \* extinction\_alpha[3]: turnover rate for
 Cold-Temperate biome in hidden state 1 \* extinction\_beta: additive
 effect of hidden state 2 on turnover rate \* hidden\_rate: transition rate
 between hidden states \* Q[1]: transition rate from Tropical to
 Warm-Temperate biome \* Q[2]: transition rate from Tropical to
 Cold-Temperate biome \* Q[3]: transition rate from Warm-Temperate to
 Tropical biome \* Q[4]: transition rate from Warm-Temperate to
 Cold-Temperate biome \* Q[5]: transition rate from Cold-Temperate to
 Tropical biome \* Q[6]: transition rate from Cold-Temperate to
 Warm-Temperate biome \* R[1]: transition rate from hidden state 1 to
 hidden state 2 \* R[2]: transition rate from hidden state 1 to hidden
 state 2 \* rate\_01: transition rate from Tropical to Warm-Temperate biome
 \* rate\_02: transition rate from Tropical to Cold-Temperate biome \*
 rate\_10: transition rate from Warm-Temperate to Tropical biome \*
 rate\_12: transition rate from Warm-Temperate to Cold-Temperate biome \*
 rate\_02: transition rate from Cold-Temperate to Tropical biome \*
 rate\_21: transition rate from Cold-Temperate to Warm-Temperate biome \*
 rate\_category\_prior[1]: root frequency of Tropical-hidden state 1 \*
 rate\_category\_prior[2]: root frequency of Warm-Temperate-hidden state 1
 \* rate\_category\_prior[3]: root frequency of Cold-Temperate-hidden state
 1 \* rate\_category\_prior[4]: root frequency of Tropical-hidden state 2
 \* rate\_category\_prior[5]: root frequency of Warm-Temperate-hidden state
 2 \* rate\_category\_prior[6]: root frequency of Cold-Temperate-hidden
 state 2 \* speciation[1]: speciation rate for Tropical biome in hidden
 state 1 \* speciation[2]: speciation rate for Warm-Temperate biome in
 hidden state 1 \* speciation[3]: speciation rate for Cold-Temperate biome
 in hidden state 1 \* speciation[4]: speciation rate for Tropical biome in
 hidden state 2 \* speciation[5]: speciation rate for Warm-Temperate biome
 in hidden state 2 \* speciation[6]: speciation rate for Cold-Temperate
 biome in hidden state 2 \* speciation\_alpha[1]: diversification rate for
 Tropical biome in hidden state 1 \* speciation\_alpha[2]: diversification
 rate for Warm-Temperate biome in hidden state 1 \* speciation\_alpha[3]:
 diversification rate for Cold-Temperate biome in hidden state 1 \*
 speciation\_beta: multiplicative effect of hidden state 2 on
 diversification rate 5. Missing data codes: None 6. Specialized formats or
 other abbreviations used: None 7. Notes: parameters sampled from a Hidden
 State Speciation and Extinction model, using the maximum a posteriori tree
 of the strongly informed biogeographic dating analysis as input, and
 coding taxa occupying multiple biome states as ambiguous \(e\.g\.
 "\(0 1\)" represent a taxon distributed between tropical and
 warm biomes\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: BiomeHiSSE\_strong\_amb\_r1.ase.tre 1. File
 type: nexus 2. Number of trees: 1 3. Number of taxa: 669 4. Description:
 tree annotated with ancestral state reconstructions from Hidden State
 Speciation and Extinction model, using the maximum a posteriori tree of
 the strongly informed biogeographic dating analysis as input, and coding
 taxa occupying multiple biome states as ambiguous \(e\.g\. "\(0
 1\)" represent a taxon distributed between tropical and warm biomes\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: BiomeHiSSE\_weakmodel.log 1. File type: tsv
 2. Number of variables: 46 3. Number of cases/rows: 10001 4. Variable
 List: \* Iteration: MCMC iteration index \* Replicate\_ID: index for
 independent MCMC run (0 or 1) \* Posterior: posterior probability of
 sampled parameters \* Likelihood: likelihood of sampled parameters \*
 Prior: prior probability of sampled parameters \* extinction[1]:
 extinction rate for Tropical biome in hidden state 1 \* extinction[2]:
 extinction rate for Warm-Temperate biome in hidden state 1 \*
 extinction[3]: extinction rate for Cold-Temperate biome in hidden state 1
 \* extinction[4]: extinction rate for Tropical biome in hidden state 2 \*
 extinction[5]: extinction rate for Warm-Temperate biome in hidden state 2
 \* extinction[6]: extinction rate for Cold-Temperate biome in hidden state
 2 \* extinction\_alpha[1]: turnover rate for Tropical biome in hidden
 state 1 \* extinction\_alpha[2]: turnover rate for Warm-Temperate biome in
 hidden state 1 \* extinction\_alpha[3]: turnover rate for Cold-Temperate
 biome in hidden state 1 \* extinction\_beta: additive effect of hidden
 state 2 on turnover rate \* hidden\_rate: transition rate between hidden
 states \* Q[1]: transition rate from Tropical to Warm-Temperate biome \*
 Q[2]: transition rate from Tropical to Cold-Temperate biome \* Q[3]:
 transition rate from Warm-Temperate to Tropical biome \* Q[4]: transition
 rate from Warm-Temperate to Cold-Temperate biome \* Q[5]: transition rate
 from Cold-Temperate to Tropical biome \* Q[6]: transition rate from
 Cold-Temperate to Warm-Temperate biome \* R[1]: transition rate from
 hidden state 1 to hidden state 2 \* R[2]: transition rate from hidden
 state 1 to hidden state 2 \* rate\_01: transition rate from Tropical to
 Warm-Temperate biome \* rate\_02: transition rate from Tropical to
 Cold-Temperate biome \* rate\_10: transition rate from Warm-Temperate to
 Tropical biome \* rate\_12: transition rate from Warm-Temperate to
 Cold-Temperate biome \* rate\_02: transition rate from Cold-Temperate to
 Tropical biome \* rate\_21: transition rate from Cold-Temperate to
 Warm-Temperate biome \* rate\_category\_prior[1]: root frequency of
 Tropical-hidden state 1 \* rate\_category\_prior[2]: root frequency of
 Warm-Temperate-hidden state 1 \* rate\_category\_prior[3]: root frequency
 of Cold-Temperate-hidden state 1 \* rate\_category\_prior[4]: root
 frequency of Tropical-hidden state 2 \* rate\_category\_prior[5]: root
 frequency of Warm-Temperate-hidden state 2 \* rate\_category\_prior[6]:
 root frequency of Cold-Temperate-hidden state 2 \* speciation[1]:
 speciation rate for Tropical biome in hidden state 1 \* speciation[2]:
 speciation rate for Warm-Temperate biome in hidden state 1 \*
 speciation[3]: speciation rate for Cold-Temperate biome in hidden state 1
 \* speciation[4]: speciation rate for Tropical biome in hidden state 2 \*
 speciation[5]: speciation rate for Warm-Temperate biome in hidden state 2
 \* speciation[6]: speciation rate for Cold-Temperate biome in hidden state
 2 \* speciation\_alpha[1]: diversification rate for Tropical biome in
 hidden state 1 \* speciation\_alpha[2]: diversification rate for
 Warm-Temperate biome in hidden state 1 \* speciation\_alpha[3]:
 diversification rate for Cold-Temperate biome in hidden state 1 \*
 speciation\_beta: multiplicative effect of hidden state 2 on
 diversification rate 5. Missing data codes: None 6. Specialized formats or
 other abbreviations used: None 7. Notes: parameters sampled from a Hidden
 State Speciation and Extinction model, using the maximum a posteriori tree
 of the weakly informed biogeographic dating analysis as input, and
 assigning taxa occupying multiple biome states to the coldest biome they
 inhabit
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: BiomeHiSSE\_weak.ase.tre 1. File type:
 nexus 2. Number of trees: 1 3. Number of taxa: 669 4. Description: tree
 annotated with ancestral state reconstructions from Hidden State
 Speciation and Extinction model, using the maximum a posteriori tree of
 the weakly informed biogeographic dating analysis as input, and assigning
 taxa occupying multiple biome states to the coldest biome they inhabit
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR:
 BiomeHiSSE\_strong\_prior\_r1model\_combined.log 1. File type: tsv 2.
 Number of variables: 46 3. Number of cases/rows: 10001 4. Variable List:
 \* Iteration: MCMC iteration index \* Replicate\_ID: index for independent
 MCMC run (0 or 1) \* Posterior: posterior probability of sampled
 parameters \* Likelihood: likelihood of sampled parameters \* Prior: prior
 probability of sampled parameters \* extinction[1]: extinction rate for
 Tropical biome in hidden state 1 \* extinction[2]: extinction rate for
 Warm-Temperate biome in hidden state 1 \* extinction[3]: extinction rate
 for Cold-Temperate biome in hidden state 1 \* extinction[4]: extinction
 rate for Tropical biome in hidden state 2 \* extinction[5]: extinction
 rate for Warm-Temperate biome in hidden state 2 \* extinction[6]:
 extinction rate for Cold-Temperate biome in hidden state 2 \*
 extinction\_alpha[1]: turnover rate for Tropical biome in hidden state 1
 \* extinction\_alpha[2]: turnover rate for Warm-Temperate biome in hidden
 state 1 \* extinction\_alpha[3]: turnover rate for Cold-Temperate biome in
 hidden state 1 \* extinction\_beta: additive effect of hidden state 2 on
 turnover rate \* hidden\_rate: transition rate between hidden states \*
 Q[1]: transition rate from Tropical to Warm-Temperate biome \* Q[2]:
 transition rate from Tropical to Cold-Temperate biome \* Q[3]: transition
 rate from Warm-Temperate to Tropical biome \* Q[4]: transition rate from
 Warm-Temperate to Cold-Temperate biome \* Q[5]: transition rate from
 Cold-Temperate to Tropical biome \* Q[6]: transition rate from
 Cold-Temperate to Warm-Temperate biome \* R[1]: transition rate from
 hidden state 1 to hidden state 2 \* R[2]: transition rate from hidden
 state 1 to hidden state 2 \* rate\_01: transition rate from Tropical to
 Warm-Temperate biome \* rate\_02: transition rate from Tropical to
 Cold-Temperate biome \* rate\_10: transition rate from Warm-Temperate to
 Tropical biome \* rate\_12: transition rate from Warm-Temperate to
 Cold-Temperate biome \* rate\_02: transition rate from Cold-Temperate to
 Tropical biome \* rate\_21: transition rate from Cold-Temperate to
 Warm-Temperate biome \* rate\_category\_prior[1]: root frequency of
 Tropical-hidden state 1 \* rate\_category\_prior[2]: root frequency of
 Warm-Temperate-hidden state 1 \* rate\_category\_prior[3]: root frequency
 of Cold-Temperate-hidden state 1 \* rate\_category\_prior[4]: root
 frequency of Tropical-hidden state 2 \* rate\_category\_prior[5]: root
 frequency of Warm-Temperate-hidden state 2 \* rate\_category\_prior[6]:
 root frequency of Cold-Temperate-hidden state 2 \* speciation[1]:
 speciation rate for Tropical biome in hidden state 1 \* speciation[2]:
 speciation rate for Warm-Temperate biome in hidden state 1 \*
 speciation[3]: speciation rate for Cold-Temperate biome in hidden state 1
 \* speciation[4]: speciation rate for Tropical biome in hidden state 2 \*
 speciation[5]: speciation rate for Warm-Temperate biome in hidden state 2
 \* speciation[6]: speciation rate for Cold-Temperate biome in hidden state
 2 \* speciation\_alpha[1]: diversification rate for Tropical biome in
 hidden state 1 \* speciation\_alpha[2]: diversification rate for
 Warm-Temperate biome in hidden state 1 \* speciation\_alpha[3]:
 diversification rate for Cold-Temperate biome in hidden state 1 \*
 speciation\_beta: multiplicative effect of hidden state 2 on
 diversification rate 5. Missing data codes: None 6. Specialized formats or
 other abbreviations used: None 7. Notes: parameters sampled from a Hidden
 State Speciation and Extinction model, using the maximum a posteriori tree
 of the weakly informed biogeographic dating analysis as input, and ran
 under the prior \(i\.e\. without data\)
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: BiomeHiSSE\_strong\_ML\_r1.out 1. File
 type: tsv 2. Number of variables: 3 3. Number of cases/rows: 408 4.
 Variable List: \* state: sampling generation within each stepping stone
 (300 - \* power: power to which likelihood is raised for sampling model
 parameters (1 = sampling from the posterior and 0 = sampling from the
 prior) \* likelihood: likelihood of the data given the sampled model
 parameters 5. Missing data codes: None 6. Specialized formats or other
 abbreviations used: None 7. Notes: power posterior analysis for Hidden
 State Speciation and Extinction \(HiSSE\) model wherein diversification
 rates vary between biomes \(Tropical, Warm\-Temperate and
 Cold\-Temperate\) and between two states of a hidden trait\. This file is
 used to estimate the marginal likelihood of the HiSSE model
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: BiomeHiSSE\_strong\_null\_ML\_r1.out 1.
 File type: tsv 2. Number of variables: 3 3. Number of cases/rows: 408 4.
 Variable List: \* state: sampling generation within each stepping stone
 (300 - 1000) \* power: power to which likelihood is raised for sampling
 model parameters (1 = sampling from the posterior and 0 = sampling from
 the prior) \* likelihood: likelihood of the data given the sampled model
 parameters 5. Missing data codes: None 6. Specialized formats or other
 abbreviations used: None 7. Notes: power posterior analysis for Hidden
 State Speciation and Extinction \(HiSSE\) model wherein diversification
 rates vary between two states of a hidden trait, but are equal across
 biomes\. This file is used to estimate the marginal likelihood of a
 suitable null model
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: run2.paleo.model.log 1. File type: tsv 2.
 Number of variables: 53 3. Number of cases/rows: 604 4. Variable List: \*
 Iteration: MCMC iteration index \* Posterior: posterior probability of
 sampled parameters \* Likelihood: likelihood of sampled parameters \*
 Prior: prior probability of sampled parameters \* biome\_rates[1]: base
 biome shift rate from Tropical to Warm-Temperate \* biome\_rates[2]: base
 biome shift rate from Warm-Temperate to Tropical \* biome\_rates[3]: base
 biome shift rate from Warm-Temperate to Cold-Temperate \* biome\_rates[4]:
 base biome shift rate from Cold-Temperate to Warm-Temperate \* br\_base:
 expected proportion of biome shift events relative to dispersal events \*
 br\_CW: base biome shift rate from Cold-Temperate to Warm-Temperate \*
 br\_TW: base biome shift rate from Tropical to Warm-Temperate \* br\_WC:
 base biome shift rate from Warm-Temperate to Cold-Temperate \* br\_WT:
 base biome shift rate from Warm-Temperate to Tropical \* dr\_base:
 expected proportion of dispersal events relative to biome shift events (1
 - br\_base) \* rate\_biome: biome shift rate \* rf\_simplex[1]: root
 frequency of biome-region Tropical+Aus \* rf\_simplex[2]: root frequency
 of biome-region Tropical+SE\_As \* rf\_simplex[3]: root frequency of
 biome-region Tropical+E\_As \* rf\_simplex[4]: root frequency of
 biome-region Tropical+Eur \* rf\_simplex[5]: root frequency of
 biome-region Tropical+Ind \* rf\_simplex[6]: root frequency of
 biome-region Tropical+Afr \* rf\_simplex[7]: root frequency of
 biome-region Tropical+N\_Am \* rf\_simplex[8]: root frequency of
 biome-region Tropical+C\_Am \* rf\_simplex[9]: root frequency of
 biome-region Tropical+S\_Am \* rf\_simplex[10]: root frequency of
 biome-region Warm-Temperate+Aus \* rf\_simplex[11]: root frequency of
 biome-region Warm-Temperate+SE\_As \* rf\_simplex[12]: root frequency of
 biome-region Warm-Temperate+E\_As \* rf\_simplex[13]: root frequency of
 biome-region Warm-Temperate+Eur \* rf\_simplex[14]: root frequency of
 biome-region Warm-Temperate+Ind \* rf\_simplex[15]: root frequency of
 biome-region Warm-Temperate+Afr \* rf\_simplex[16]: root frequency of
 biome-region Warm-Temperate+N\_Am \* rf\_simplex[17]: root frequency of
 biome-region Warm-Temperate+C\_Am \* rf\_simplex[18]: root frequency of
 biome-region Warm-Temperate+S\_Am \* rf\_simplex[19]: root frequency of
 biome-region Cold-Temperate+Aus \* rf\_simplex[20]: root frequency of
 biome-region Cold-Temperate+SE\_As \* rf\_simplex[21]: root frequency of
 biome-region Cold-Temperate+E\_As \* rf\_simplex[22]: root frequency of
 biome-region Cold-Temperate+Eur \* rf\_simplex[23]: root frequency of
 biome-region Cold-Temperate+Ind \* rf\_simplex[24]: root frequency of
 biome-region Cold-Temperate+Afr \* rf\_simplex[25]: root frequency of
 biome-region Cold-Temperate+N\_Am \* rf\_simplex[26]: root frequency of
 biome-region Cold-Temperate+C\_Am \* rf\_simplex[27]: root frequency of
 biome-region Cold-Temperate+S\_Am \* w\_adj[1]: marginal features for
 availability/connectivity among biomes/regions (always 0) \* w\_adj[2]:
 weak features for availability/connectivity among biomes/regions (equal to
 y and between 0 and 1) \* w\_adj[3]: strong features for
 availability/connectivity among biomes/regions (always 1) \* w\_B: weight
 of biome features \* w\_G: weight of biome-independent land (geographical)
 features \* w\_graph[1]: weight of null (biome- and land-independent)
 features \* w\_graph[2]: weight of biome-independent land (geographical)
 features \* w\_graph[3]: weight of biome features \* w\_not\_B: weight of
 null and land features combined \* w\_U: weight of null (biome- and
 land-independent) features \* y: effect of weak features relative to
 strong and marginal features 5. Missing data codes: None 6. Specialized
 formats or other abbreviations used: see Table S7 for biome\-region
 abbreviations 7. Notes: parameters sampled from a joint dispersal and
 biome shift model, using the maximum a posteriori tree of the strongly
 informed biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: run2.paleo.ase.tre 1. File type: nexus 2.
 Number of trees: 1 3. Number of taxa: 669 4. Description: tree annotated
 with ancestral state reconstructions from a joint dispersal and biome
 shift model, using the maximum a posteriori tree of the strongly informed
 biogeographic dating analysis as input
 #########################################################################
 DATA-SPECIFIC INFORMATION FOR: run2.paleo.history.tsv 1. File type: tsv 2.
 Number of variables: 11 3. Number of cases/rows: 1470376 4. Variable List:
 \* iteration: MCMC iteration index \* node\_index: RevBayes index for tree
 node (1-1337) \* branch\_start\_time: start time of branch (in Ma) \*
 branch\_end\_time: end time of branch (in Ma) \* start\_state:
 biome-region state at the start of the branch \* end\_state: biome-region
 state at the end of the branch \* transition\_time: time (in Ma) when the
 transition between biome-region states occurred \* transition\_type:
 always anagenetic (i.e. occuring along branches) \* parent\_index: node
 index of parent node \* child1\_index: node index of first child node \*
 child2\_index: node index of second child node 5. Missing data codes: None
 6. Specialized formats or other abbreviations used: see Table S7 for
 biome\-region state symbols 7. Notes: stochastic character histories
 simulated in a joint dispersal and biome shift model, using the maximum a
 posteriori tree of the strongly informed biogeographic dating analysis as
 input
 #########################################################################