“…The models used were as follows: TVM+F+G4 for both ITS-1 and ITS2, and K2P+I for 5.8S, all included in the ribosomal gene cluster ITS1-5.8S-ITS2 alignment (3 partitions); TVM+F+I+G4 for the 28S dataset alone; TIM3+F+I+G4, TIM+F+G4, GTR+F+G4, for the 1st, 2nd and 3rd codon positions of the cox 1 dataset, respectively (3 partitions); GTR+F+G4, K2P+I, TVM+F+G4, TVM+F+I+G4, TN+F+G4, K3Pu+F+G4, TPM2+F+G4, for the ITS1, 5.8S, ITS2, 28S, cox 1 1st + 2nd + 3rd codon positions of the concatenated dataset (7 partitions). Outgroups were chosen based on previous phylogenetic studies ( Briosio-Aguilar et al, 2019 ; Rosser et al, 2020 ; Montes et al, 2021 ).…”