Processed data supporting the manuscript "Cutting the sap: first molecular phylogeny of twig-girdler longhorn beetles (Coleoptera: Cerambycidae: Lamiinae: Onciderini) suggests shifts in host plant attack behaviors contributed to morphological evolution"

Processed data supporting the manuscript: Cutting the sap: first molecular phylogeny of twig-girdler longhorn beetles (Coleoptera: Cerambycidae: Lamiinae: Onciderini) suggests shifts in host plant attack behaviors contributed to morphological evolution

By: Diego de S. Souza 1, 2, Rowan L. K. French 3, José O. Silva Júnior 4, Eugenio H. Nearns 5, Luciane Marinoni 4, Ian P. Swift 6, Kelly B. Miller 7, Felix A. H. Sperling 2 & Marcela L. Monné 1

1 Department of Entomology, National Museum, Federal University of Rio de Janeiro, Rio de Janeiro, Rio de Janeiro, Brazil.

2 Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada.

3 Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada.

4 Department of Zoology, Federal University of Paraná, Curitiba, Paraná, Brazil.

5 National Museum of Natural History, Smithsonian Institution, Washington, DC, USA.

6 California State Collection of Arthropods, Sacramento, California, USA.

7 Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA.

Corresponding author: Diego de S. Souza, [e-mail]. Current affiliation: Field Museum of Natural History, Chicago, Illinois, USA.

List of Contents: 

Onciderini_concat_matrix.phy

Concatenated matrix (cox1, Wg and CPS) used for the phylogenetic analyses of Onciderini (Coleoptera: Cerambycidae: Lamiinae: Onciderini). 

PartitionFinder_AICc_best_scheme.txt

Results from PartitionFinder v2.1.1, containing the best partitioning scheme for the concatenated matrix of Onciderini, identified using the corrected Akaike Information Criterion (AICc), with model definitions for use in the phylogenetic analyses.

RAxML_Onciderini_concat_matrix (zip file)

- Onciderini_concat_matrix.phy: concatenated matrix (cox1, Wg and CPS) used in the RAxML phylogenetic analyses of Onciderini (Coleoptera: Cerambycidae: Lamiinae: Onciderini).

- Partitions_AICc_RAxML.txt: partitioning scheme used in the RAxML analysis as predefined by PartitionFinder v2.1.1 using the corrected Akaike Information Criterion (AICc).

- RAxML_bestTree.Onciderini_concat_matrix_ML: best-scoring maximum likelihood tree inferred by RAxML for the concatenated matrix of Onciderini.

- RAxML_bipartitions.Onciderini_concat_matrix_final: bipartitions (clades) of the maximum likelihood tree inferred by RAxML with support values estimated from 1,000 pseudoreplicates.

- RAxML_bipartitionsBranchLabels.Onciderini_concat_matrix_final: final maximum likelihood tree inferred by RAxML for the concatenated matrix of Onciderini, with labeled branches showing bootstrap support values.

- RAxML_bootstrap.Onciderini_concat_matrix_bootstrap: bootstrap trees generated from a non-parametric bootstrap analysis in RAxML based on 1,000 pseudoreplicates.

- RAxML_info.Onciderini_concat_matrix_bootstrap: log file containing details of the bootstrap analysis, including the settings and parameters used in the non-parametric bootstrap runs in RAxML.

- RAxML_info.Onciderini_concat_matrix_final: log file summarizing the RAxML analysis, including settings and convergence statistics for the final maximum likelihood tree.

- RAxML_info.Onciderini_concat_matrix_ML: log file containing details of the maximum likelihood tree search, including the parameters and models applied during the maximum likelihood analysis conducted by RAxML.

- RAxML_log.Onciderini_concat_matrix_ML: log file of the maximum likelihood tree search for the concatenated matrix of Onciderini.

- RAxML_parsimonyTree.Onciderini_concat_matrix_ML: parsimony starting tree used by RAxML during the maximum likelihood analysis for the concatenated matrix of Onciderini.

- RAxML_result.Onciderini_concat_matrix_ML: maximum likelihood tree inferred by RAxML from the concatenated matrix of Onciderini, summarizing the tree topology and likelihood score for the best tree obtained.

BI_AICc_Onciderini_concat_matrix (zip file)

- BI_AICc_Onciderini_concat_matrix.nex: nexus file containing the concatenated matrix of Onciderini used for Bayesian Inference (BI), including the best-fit model scheme identified by PartitionFinder and MCMC parameters for running the analysis in MrBayes.

- BI_AICc_Onciderini_concat_matrix.nex_r1_r2_combined_consensus.tree: consensus tree from two combined independent Bayesian Inference (BI) runs based on the concatenated matrix of Onciderini, after discarding the first 25% of initial generations as burn-in.

- BI_AICc_Onciderini_concat_matrix.nex.run1.p: log file containing parameter values and likelihood scores from the first run of the Bayesian Inference (BI) based on the concatenated matrix of Onciderini.

- BI_AICc_Onciderini_concat_matrix.nex.run2.p: log file containing parameter values and likelihood scores from the second run of the Bayesian Inference (BI) based on the concatenated matrix of Onciderini.

BEAST2_Onciderini_BD_lognormal (zip file)

- BEAUTi_Onciderini_BD_lognormal.xml: XML file generated by BEAUTi for running BEAST2, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and a lognormal distribution.

- BEAST2_Onciderini_BD_lognormal_run[1-8].log: log files from eight independent runs of BEAST2, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and a lognormal distribution.

- TreeAnnotator_Onciderini_BD_lognormal_run1-run8_consensus.out: TreeAnnotator output file combining the results of eight BEAST2 runs based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and a lognormal distribution.

- TreeAnnotator_Onciderini_BD_lognormal_run1-run8_consensus.tre: consensus tree from eight combined BEAST2 runs, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and a lognormal distribution, after discarding the first 10% of initial generations as burn-in.

BEAST2_Onciderini_BD_exponential (zip file)

- BEAUTi_Onciderini_BD_exponential.xml: XML file generated by BEAUTi for running BEAST2, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and an exponential distribution.

- BEAST2_Onciderini_BD_exponential_[1-8].log: log files from eight independent runs of BEAST2, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and an exponential distribution.

- TreeAnnotator_Onciderini_BD_exponential_run1-run8_consensus.out: TreeAnnotator output file combining the results of eight BEAST2 runs based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and an exponential distribution.

- TreeAnnotator_Onciderini_BD_exponential_run1-run8_consensus.tre: consensus tree from eight combined BEAST2 runs, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and an exponential distribution, after discarding the first 10% of initial generations as burn-in.

BEAST2_Onciderini_BD_uniform (zip file)

- BEAUTi_Onciderini_BD_uniform.xml: XML file generated by BEAUTi for running BEAST2, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and a uniform distribution.

- BEAST2_Onciderini_BD_uniform_[1-8].log: log files from eight independent runs of BEAST2, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and a uniform distribution.

- TreeAnnotator_Onciderini_BD_uniform_run1-run8_consensus.out: TreeAnnotator output file combining the results of eight BEAST2 runs based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and a uniform distribution.

- TreeAnnotator_Onciderini_BD_uniform_run1-run8_consensus.tre: consensus tree from eight combined BEAST2 runs, based on the concatenated matrix of Onciderini, using a birth-death (BD) process model and a uniform distribution, after discarding the first 10% of initial generations as burn-in.

Comparative_analyses (zip file)

RawData (folder): raw morphometric and girdling data, plus tree that was later pruned for downstream comparative analyses; these data were used as input for the OncidHeadDimorphism-DatasetPREP-FINAL.R data cleaning script. 

- Onciderini_BD_lognormal_run1-run8_consensus.nwk: newick version of TreeAnnotator_Onciderini_BD_lognormal_run1-run8_consensus.tre (outputted as a newick file by importing the .tre file into FigTree and exporting in newick format).

- Measurements_Onciderini_Raw_Final.csv: individual-level raw morphometric data for Onciderini.

- Behav_Matrix_2_states_trimmed_2022-12-15.csv: species-level data on girdling status for Onciderini species, with all Lochmaeocles species classified as girdlers (2 behavioral states across Onciderini species). 

- Matrix_3_states_trimmed_final.csv: species-level data on girdling status for Onciderini species, with all Lochmaeocles species classified as facultative girdlers (3 behavioral states across Onciderini species). 

- Matrix_2_states_trimmed_1LochGirdler_Final.csv: species-level data on girdling status for Onciderini species, with only one Lochmaeocles species (L. tessellatus) classified as a girdler (2 behavioral states across Onciderini species). 

ProcessedData (folder): filtered data and pruned trees outputted by the OncidHeadDimorphism-DatasetPREP-FINAL.R script

- oncid_f_36spp_clean.csv: dataset of species means and log-ratios for morphometric traits in females, plus girdling data; only includes species that have girdling data and are in the phylogenetic tree

- oncid_m_42spp_clean.csv: dataset of species means and log-ratios for morphometric traits in males, plus girdling data; only includes species that have girdling data and are in the phylogenetic tree

- oncid_sd_35spp_clean.csv: dataset of species means for sexual dimorphism in morphometric traits, plus girdling data; only includes species that have girdling data and are in the phylogenetic tree

- oncid_girdlingbehav_allingroupspp_clean.csv: full dataset of girdling behavior for 56 Onciderini species that are in the phylogenetic tree; includes separate columns for the three alternative girdling classification schemes

- oncid_tree_behavfull_56spp.nwk: pruned phylogenetic tree for the full girdling dataset (56 species)

- oncid_tree_f_36spp.nwk: pruned phylogenetic tree for the female morphometric dataset (36 species)

- oncid_tree_m_42spp.nwk: pruned phylogenetic tree for the male morphometric dataset (42 species)

- oncid_tree_mf_43spp.nwk: pruned phylogenetic tree for all species with morphometric data for males or females; used for the stochastic character map next to the heatmap plot (Fig 4)

- oncid_tree_sd_35spp.nwk: pruned phylogenetic tree for the sexual dimorphism dataset (35 species)

FittedModels (folder): fitted models (mvgls, model comparison analyses, OUM models, simmaps) outputted by the OncidHeadDimorphism-Analysis-FINAL.R script 

- MacroModelFits_logRtraits_f_36spp-2024-10-12.Rdata: summary of model comparison results for Brownian Motion (BM), single-peak Ornstein-Uhlenbeck (OU), multipeak OU (OUM), and multi-rate Brownian motion (BMM) models (univariate and multivariate) fitted to female morphometric data across 100 stochastic character maps of girdling behavior

- MacroModelFits_logRtraits_m_42spp-2024-10-12.Rdata: summary of model comparison results for Brownian Motion (BM), single-peak Ornstein-Uhlenbeck (OU), multipeak OU (OUM), and multi-rate Brownian motion (BMM) models (univariate and multivariate) fitted to male morphometric data across 100 stochastic character maps of girdling behavior

- MacroModelFits-SDDI-35spp_2024-10-11.Rdata: summary of model comparison results for Brownian Motion (BM), single-peak Ornstein-Uhlenbeck (OU), multipeak OU (OUM), and multi-rate Brownian motion (BMM) models (univariate and multivariate) fitted to sexual dimorphism data across 100 stochastic character maps of girdling behavior

- mvgls-results-headsize-mf-2024-10-12.Rdata: fitted mvgls regression models for male and female traits (analyzed separately)

- mvgls-results-sddi-2024-10-12.Rdata: fitted mvgls regression models for sexual dimorphism

- OUM_headtraits_f_36spp-2024-10-12.Rdata: fitted OUM models and summary statistics for female head traits

- OUM_headtraits_m_42spp-2024-10-12.Rdata: fitted OUM models and summary statistics for male head traits

- OUM-SDDI-35spp-2024-10-12.Rdata: fitted OUM models and summary statistics for sexual dimorphism in head traits

- simmaps_ard_full_2state.RDS: stochastic character maps of girdling behaviour for all 56 species with girdling data, with two behavioral states (girdling or non-girdling) - all Lochmaeocles species are classified as girdlers

- simmaps_ard_full_3state.RDS: stochastic character maps of girdling behaviour for all 56 species with girdling data, with three behavioral states (girdling, non-girdling, or facultative girdling)

- simmaps_ard_full_1Loch.RDS: stochastic character maps of girdling behaviour for all 56 species with girdling data, with two behavioral states (girdling or non-girdling) - only one Lochmaeocles species (L. tessellatus) is classified as a girdler

- simmaps_ard_m.RDS: stochastic character map for the 42 species used in the analyses of male morphometric traits

- simmaps_ard_f.RDS: stochastic character map for the 36 species used in the analyses of female morphometric traits

- simmaps_ard.RDS: stochastic character map for the 35 species used in the sexual dimorphism analyses; 2 behavioral states.

Rscripts (folder): R scripts used to process data and run phylogenetic comparative analyses of head size and girdling behavior.

- OncidHeadDimorphism-DatasetPREP-FINAL.R: R script used to filter data and prune trees from the RawData folder for downstream phylogenetic comparative analyses; outputs of this script are in the ProcessedData folder.

- OncidHeadDimorphism-Analysis-FINAL.R: R script used to analyze data in the ProcessedData folder to answer questions about the origin and evolution of girdling behavior and the relationship between girdling and head size or head size sexual dimorphism; fitted models outputted by this script are in the FittedModels folder

- OncidHeadDimorphism-Plots-FINAL.R: R script used to generate plots for the manuscript