Molecular phylogenetics and evolutionary history of the endemic land snail genus Everettia in northern Borneo

Supplemental
Information

Additional File 1. Species records
of Genus Everettia in Sabah for Maxent Analysis. The table consists of
the 718 records of 17 Everettia species in Sabah obtained from
BORNEENSIS collection, Universiti Malaysia Sabah. The collection lot reference
number, number of specimens in each collection lot, and the geographic
coordinates of the specimens were included in the table.

Additional
File 2. Results of ModelTest for each partition of DNA sequencing alignment and
parameter used in BEAST analysis. Each of the six partitions, namely, codons of COI, namely, 1st, 2nd and 3rd
codon positions of COI, 16S rDNA, ITS-1, and 28S rDNA, was tested for molecular evolution via ModelFinder (Kalyaanamoorthy et al.,
2017) and partition models (Chernomor et al., 2016) based on the both AIC and BIC that built into
IQ-Tree v.1.6.7 (Nguyen et al., 2015; Trifinopoulos et al. 2016). We limited
the candidate models to the six models that are available in MrBayes analysis,
namely, JC, F81, K80, HKY, SYM, and GTR.

Additional
File 3. A zipped folder consists of input files and outputs of Bayesian (BA)
and Maximum Likelihood (ML) analysis. The input files and outputs of Bayesian
(BA) and Maximum Likelihood (ML) analysis for each of two different best-fit
substitution models selected by BIC and AIC criteria, respectively. The phylogenies
for each of the analyses were summarised in the word document file.

Additional File 4. A zipped folder consists of input files and
outputs of BEAST analysis. A total of four BEAST input XML files for the combinations two different best-fit substitution models (selected
by BIC and AIC criteria) and two calibrated models (Yule model vs. Birth-Death
model). The calibrated phylogenies for each of the four analyses were
summarised in the word document file.

Additional File 5. Correlations
between bioclim variables. Bioclimatic variables were sampled 500 at random
locations in Borneo in each of the 19 climatic layers. Collinearity among the
19 climatic variables was evaluated by using pairwise Pearson’s r correlation.
After the analysis, 12 highly correlated variables (r > 0.8) were excluded
from MAXENT analysis. The seven climatic variables were used for species
distribution modelling, namely, BIO1 Annual Mean Temperature, BIO3
Isothermality, BIO4 Temperature Seasonality, BIO7 Temperature Annual Range,
BIO12 Annual Precipitation, BIO15 Precipitation Seasonality, and BIO19
Precipitation of Coldest Quarter.

Additional File 6. Concatenated DNA
Data Matrix for 16S, COI, ITS, and 28S Sequences for 73 taxa. DNA sequences
alignment in FASTA format. Position 1-501: 16S; Position 502-1059: COI;
Position 1060-1869: ITS; and Position 1870-2795: 28S.