Population structure based on mitochondrial data
We obtained the DNA sequences of of 633 bp from the COI gene and 428 bp from the cyt b gene and characterized five individuals randomly selected from each of the 14 population. The analysis of theCOI gene revealed two different mitochondrial haplotypes in the 66 individuals sampled from 14 populations. The frequency of the haplotype A (GenBank accession no. KP670306.1) was 28.8 % (19 individuals). While haplotype A was common in commercial populations, haplotype B (GenBank accession no. KP670307.1) was common in native and nearby greenhouses populations. The frequency of the haplotype B was 71.2% (47 individuals).
Six different haplotypes (haplotype 1-6) were detected in 68 individuals sampled from 14 populations for the cyt b gene. The frequencies of the haplotype 1, 2, 3, 4, 5 and 6 were 7.3 % (5 individuals), 5.9 % (4 individuals), 11.8% (8 individuals), 63.2 % (43 individuals), 5.9 % (4 individuals), 5.9 % (4 individuals) respectively. When the sequences of the haplotypes were checked with the BLAST program in the Gene Bank, it was found that the haplotype 4 and haplotype 5 genotypes clustered within the B. t. dalmatinus subspecies. However, haplotype 1, haplotype 2, haplotype 3 and haplotype 6 genotypes clustered within both B. t. terrestris and B. t. dalmatinus subspecies. Maximum likelihood (ML) phylogeny inferred from COI and cyt b gene regions sequences of B. terrestris individuals were given in Figure S1 and Figure S2, respectively.