Research Article |
Corresponding author: Lucie Vaníčková ( luci.vanickova@gmail.com ) Academic editor: Anthony Clarke
© 2015 Lucie Vaníčková, Radka Břízová, Antonio Pompeiano, Sunday Ekesi, Marc De Meyer.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Vaníčková L, Břízová R, Pompeiano A, Ekesi S, De Meyer M (2015) Cuticular hydrocarbons corroborate the distinction between lowland and highland Natal fruit fly (Tephritidae, Ceratitis rosa) populations. In: De Meyer M, Clarke AR, Vera MT, Hendrichs J (Eds) Resolution of Cryptic Species Complexes of Tephritid Pests to Enhance SIT Application and Facilitate International Trade. ZooKeys 540: 507-524. https://doi.org/10.3897/zookeys.540.9619
|
The cuticular hydrocarbons (CHs) and morphology of two Ceratitis rosa Karsch (Diptera: Tephritidae) populations, putatively belonging to two cryptic taxa, were analysed. The chemical profiles were characterised by two-dimensional gas chromatography with mass spectrometric detection. CHs of C. rosa that originated from the lowlands and highlands of Kenya comprised of n-alkanes, monomethylalkanes, dimethylalkanes and unsaturated hydrocarbons in the range of the carbon backbone from C14 to C37. Hydrocarbons containing C29, C31, C33 and C35 carbon atoms predominated in these two populations. 2-Methyltriacontane was the predominant compound in both populations. Quantitative differences in the distribution of hydrocarbons of different chain lengths, mainly the C22, C32, C33 and C34 compounds of these two populations, were observed despite indistinct qualitative differences in these hydrocarbons. Morphological analyses of male legs confirmed that the flies belong to different morphotypes of C. rosa previously labelled as R1 and R2 for lowland and highland populations, respectively. A statistical analysis of the CH compositions of the putative R1 and R2 species showed distinct interspecific identities, with several CHs specific for each of the lowland and highland populations. This study supports a hypothesis that the taxon C. rosa consists of at least two biological species.
Ceratitis rosa , cryptic species, chemotaxonomy, GC×GC/MS, integrative taxonomy
Sexual selection within populations can play an important role in speciation when divergence in mating signals and corresponding mate preference occur along different evolutionary trajectories in different populations (
In species of the fruit fly genus Ceratitis, courtship generally includes visual, auditory, tactile and olfactory cues (
The Natal fruit fly, Ceratitis rosa Karsch (Diptera, Tephritidae), is a polyphagous species attacking a wide range of fruits on the African mainland. It has invaded some Indian Ocean islands, where it displaced the similarly introduced C. capitata (
The chemical analyses of the cuticular hydrocarbon profiles of these putative species found significant differences between the A, F2 and R2 genotypes and characterised chemotaxonomic markers to distinguish these groups (
The literature provides conflicting information regarding developmental physiology and climatic niche for C. rosa. Some studies indicate that C. rosa might be more tolerant of colder and wetter conditions than C. capitata (
The purpose of the present study was, therefore, to identify the chemical constituents of the CHs and to analyse their variation between two populations of C. rosa (one highland and one lowland - based on morphological differentiation) originating from Kenya. These two populations were chosen for this study because they had previously been shown to be sexually incompatible (Ekesi et al. unpublished data), as well as having distinct male-borne volatile profiles (Kalinová et al. unpublished data). Additional to inter-population differences, we also evaluated sexual dimorphism in CHs composition within each population.
Pupae of two laboratory populations of C. rosa were obtained from the International Centre of Insect Physiology and Ecology (ICIPE, Nairobi, Kenya). The source colonies were established in 2012 and came from one lowland locality [Mwajamba, Msambweni, Coast Province, 04°18.21'S; 39°29.88'E, host fruit Psidium guajava (Myrtaceae), altitude 106 m, average temperature 28.1 °C] and one highland locality [Kithoka, Meru, Central Province, 00°05.59'N; 37°40.40'E, host fruit Mangifera indica (Anacardiaceae), altitude 1425 m, average temperature 21.5 °C] in Kenya (see
The extraction of the cuticular hydrocarbons of 20-day-old virgin males (N = 10) and females (N = 10) of the R1 and R2 morphotypes (resulting in N = 20 for R1 and N = 20 for R2) followed the methodology described in
Two-dimensional gas chromatography with time-of-flight mass spectrometric detection (GC×GC/MS) was used for the quantification and identification of CH profiles. The analyses were performed on a LECO Pegasus 4D instrument (LECO Corp., St. Joseph, MI, USA) equipped with a non-moving quad-jet cryomodulator. A DB-5 column (J&W Scientific, Folsom, CA, USA; 30 m × 250 µm i.d. × 0.25 µm film) was used for GC in the first dimension. The second-dimension analysis was performed on a polar BPX-50 column (SGE Inc., Austin, TX, USA; 2 m × 100 µm i.d. × 0.1 µm film). Helium was used as a carrier gas at a constant flow of 1 mL min-1. The temperature program for the primary GC oven was as follows: 150 °C for 2 min, then 150–300 °C at 5 °C min-1, and finally a 10 min hold at 320 °C. The program in the secondary oven was 10 °C higher than in the primary one and was operated in an iso-ramping mode. The modulation period, the hot-pulse duration and the cool time between the stages were set to 3.0, 0.4 and 1.1 sec, respectively. The transfer line to the TOFMS was operated at 260 °C. The source temperature was 250 °C with a filament bias voltage of −70 eV. The data-acquisition rate was 100 Hz (scans/sec) for the mass range of 29–400 amu. The detector voltage was 1750V. For each sample, 1µL was injected in splitless mode. The inlet temperature was 200 °C. The purge time was 60 sec at a flow of 60 mL min-1. The data were processed and consecutively visualized on 2D and 3D chromatograms using LECO ChromaTOFTM software. The n-alkane standard (C8–C38; Sigma-Aldrich) was co-injected with authentic samples to determine the retention indices (RI) of the analytes. The hydrocarbons were identified by a comparison of their mass spectra fragmentation patterns and RI (
Male specimens were shipped to the Royal Museum for Central Africa (RMCA), Tervuren, Belgium, where identifications were confirmed by M. De M. based on the pilosity and coloration of mid tibia (
The relative peak areas of 46 CH compounds (as identified by the GC×GC/MS in the deconvoluted total-ion chromatogram mode) were calculated in 10 replicate specimens for each sex of the two species (N = 40). Following
A heat map was used to visualise the complex data sets organised as matrices. Heat maps make it possible to identify differences in the relative amounts of CHs between populations, with different compounds tending to form small clusters according to their quantities. To achieve this, the heat map performed two actions on a matrix of chromatographic peak areas. First, it reordered the rows and columns so that rows and columns with similar profiles were closer to one another, causing these profiles to be more visible to the eye. Second, each entry in the data matrix was displayed in a different colour, making it possible to view the patterns graphically. The dendrograms were created using correlation-based distances and the Ward method of agglomeration was applied in the present analysis (
To examine the differences between the two populations and sexes further, the percentage contribution of each compound to the average dissimilarity between the aforementioned factors was calculated with similarity percentage analysis (SIMPER) (
The GC×GC/MS analyses identified 46 peaks. The chain-length of the carbon backbones ranged from C14 to C37. The hydrocarbon profiles of the males and females included 5 n-alkanes, 19 methylbranched alkanes, 19 unsaturated alkanes, squalene, 1 aldehyde and 1 unidentified compound. The heat map characterised differences in the relative amounts of CHs between the C. rosa flies originating from highland and coastal regions (Figure
A heat map of the 46 cuticular hydrocarbons (columns, CH1-46) and the two Ceratitis rosa populations (rows, f-female, m-male) from the GC×GC/MS data set. The dendrograms are created using correlation-based distances and the Ward method of hierarchical clustering (P < 0.05). Putative morphotypes (R1 for the coastal population and R2 for the highland population) are depicted in the row dendrogram.
A comparison of the average abundance of important cuticular hydrocarbons between two morphotypes of Ceratitis rosa [coastal R1, highland R2]. The compounds are listed in the order of their contribution (δi) to the average dissimilarity 5(δi) between the two groups, with a cut-off when the cumulative percent contribution (∑δi%) to δi reaches 70%. The numbering of the compounds corresponds to Figure
No. | Compound | RI | Abundance | δi | δi /SD(δi) | % contr. diss. | ∑δi% | |
---|---|---|---|---|---|---|---|---|
R1 male | R2 male | |||||||
15 | 2-MeC28 | 2865 | 1.037 | 1.731 | 0.016 | 3.384 | 0.104 | 15 |
11 | C22:1 | 2182 | 0.425 | 0.995 | 0.015 | 1.592 | 0.096 | 11 |
35 | diMeC31 | 3297 | 0.883 | 0.447 | 0.010 | 2.884 | 0.065 | 35 |
26 | diMeC28 | 3105 | 1.575 | 1.147 | 0.010 | 4.198 | 0.064 | 26 |
29 | 3-MeC31 | 3178 | 1.093 | 0.797 | 0.007 | 3.241 | 0.044 | 29 |
12 | C27:1 | 2622 | 0.198 | 0.419 | 0.006 | 1.531 | 0.039 | 12 |
16 | diMeC26 | 2902 | 0.613 | 0.869 | 0.006 | 2.443 | 0.039 | 16 |
30 | diMeC29 | 3205 | 0.395 | 0.196 | 0.005 | 2.623 | 0.030 | 30 |
37 | MeC33 | 3331 | 0.721 | 0.906 | 0.004 | 1.826 | 0.029 | 37 |
36 | C34:1 | 3308 | 0.815 | 0.650 | 0.004 | 1.389 | 0.028 | 36 |
38 | C34:1 | 3342 | 0.206 | 0.370 | 0.004 | 1.431 | 0.026 | 38 |
23 | 2-MeC30 | 3064 | 2.045 | 1.882 | 0.004 | 1.633 | 0.026 | 23 |
2 | unknown | 1402 | 0.816 | 0.649 | 0.004 | 2.035 | 0.025 | 2 |
27 | 7-/9-MeC31 | 3142 | 1.044 | 0.882 | 0.004 | 1.768 | 0.024 | 27 |
1 | C14 | 1400 | 0.934 | 0.771 | 0.004 | 1.826 | 0.024 | 1 |
34 | C33:1 | 3291 | 0.378 | 0.261 | 0.004 | 1.134 | 0.024 | 34 |
28 | MeC31 | 3152 | 0.317 | 0.227 | 0.004 | 1.695 | 0.023 | 28 |
No. | Compound | RI | Abundance | δi | δi /SD(δi) | % contr.diss. | ∑δi% | |
R1 female | R2 female | |||||||
11 | C22:1 | 2182 | 0.435 | 1.415 | 0.022 | 2.485 | 0.133 | 11 |
15 | 2-MeC28 | 2865 | 1.186 | 1.768 | 0.013 | 2.817 | 0.079 | 15 |
29 | 3-MeC31 | 3178 | 1.078 | 0.599 | 0.011 | 4.145 | 0.065 | 29 |
26 | diMeC28 | 3105 | 1.534 | 1.191 | 0.008 | 2.284 | 0.047 | 26 |
34 | C33:1 | 3291 | 0.352 | 0.112 | 0.007 | 1.789 | 0.044 | 34 |
28 | MeC31 | 3152 | 0.384 | 0.066 | 0.007 | 3.160 | 0.043 | 28 |
33 | C33:1 | 3280 | 1.353 | 1.444 | 0.005 | 1.337 | 0.029 | 33 |
30 | diMeC29 | 3205 | 0.362 | 0.148 | 0.005 | 2.727 | 0.029 | 30 |
27 | 7-/9-MeC31 | 3142 | 1.043 | 0.835 | 0.005 | 1.875 | 0.029 | 27 |
36 | C34:1 | 3308 | 0.705 | 0.639 | 0.005 | 1.373 | 0.029 | 36 |
1 | C14 | 1400 | 0.732 | 0.925 | 0.005 | 1.791 | 0.028 | 1 |
42 | C35:2 | 3460 | 1.263 | 1.223 | 0.004 | 1.381 | 0.027 | 42 |
2 | unknown | 1402 | 0.623 | 0.811 | 0.004 | 1.831 | 0.027 | 2 |
35 | diMeC31 | 3297 | 0.803 | 0.629 | 0.004 | 1.384 | 0.026 | 35 |
16 | diMeC26 | 2902 | 0.689 | 0.858 | 0.004 | 1.512 | 0.026 | 16 |
38 | C34:1 | 3342 | 0.199 | 0.318 | 0.004 | 2.264 | 0.024 | 38 |
24 | C31:1 | 3082 | 0.414 | 0.354 | 0.004 | 1.426 | 0.024 | 24 |
No. | Compound | RI | Abundance | δi | δi /SD(δi) | % contr. diss. | ∑δi% | |
R1 male | R1 female | |||||||
33 | C33:1 | 3280 | 1.026 | 1.353 | 0.008 | 1.680 | 0.071 | 33 |
11 | C22:1 | 2182 | 0.425 | 0.435 | 0.005 | 1.118 | 0.049 | 11 |
24 | C31:1 | 3082 | 0.196 | 0.414 | 0.005 | 1.633 | 0.045 | 24 |
36 | C34:1 | 3308 | 0.815 | 0.705 | 0.005 | 1.387 | 0.043 | 36 |
1 | C14 | 1400 | 0.934 | 0.732 | 0.005 | 2.207 | 0.043 | 1 |
2 | unknown | 1402 | 0.816 | 0.623 | 0.004 | 2.217 | 0.041 | 2 |
15 | 2-MeC28 | 2865 | 1.037 | 1.186 | 0.004 | 1.220 | 0.040 | 15 |
40 | C34:2 | 3371 | 0.311 | 0.242 | 0.004 | 1.376 | 0.034 | 40 |
42 | C35:2 | 3460 | 1.352 | 1.263 | 0.004 | 1.500 | 0.034 | 42 |
34 | C33:1 | 3291 | 0.378 | 0.352 | 0.003 | 1.160 | 0.031 | 34 |
35 | diMeC31 | 3297 | 0.883 | 0.803 | 0.003 | 1.759 | 0.031 | 35 |
28 | MeC31 | 3152 | 0.317 | 0.384 | 0.003 | 1.215 | 0.030 | 28 |
22 | C31:1 | 3047 | 0.155 | 0.258 | 0.003 | 1.958 | 0.027 | 22 |
32 | 3-MeC32 | 3262 | 1.031 | 1.125 | 0.003 | 1.557 | 0.027 | 32 |
31 | C33:1 | 3240 | 1.406 | 1.516 | 0.003 | 1.614 | 0.027 | 31 |
19 | MeC29 | 2960 | 0.477 | 0.586 | 0.003 | 1.350 | 0.026 | 19 |
12 | C27:1 | 2622 | 0.198 | 0.101 | 0.003 | 0.915 | 0.026 | 12 |
13 | MeC26 | 2649 | 0.189 | 0.112 | 0.003 | 1.194 | 0.025 | 13 |
38 | C34:1 | 3342 | 0.206 | 0.199 | 0.003 | 1.058 | 0.024 | 38 |
26 | diMeC28 | 3105 | 1.575 | 1.534 | 0.003 | 1.439 | 0.024 | 26 |
16 | diMeC26 | 2902 | 0.613 | 0.689 | 0.002 | 1.200 | 0.023 | 16 |
No. | Compound | RI | Abundance | δi | δi /SD(δi) | % contr. diss. | ∑δi% | |
R2 male | R2 female | |||||||
11 | C22:1 | 2182 | 0.995 | 1.415 | 0.013 | 1.324 | 0.093 | 11 |
33 | C33:1 | 3280 | 1.095 | 1.444 | 0.008 | 1.904 | 0.058 | 33 |
31 | C33:1 | 3240 | 1.357 | 1.673 | 0.007 | 2.177 | 0.052 | 31 |
12 | C27:1 | 2622 | 0.419 | 0.115 | 0.007 | 1.604 | 0.050 | 12 |
34 | C33:1 | 3291 | 0.261 | 0.112 | 0.006 | 1.486 | 0.042 | 34 |
24 | C31:1 | 3082 | 0.141 | 0.354 | 0.005 | 1.491 | 0.038 | 24 |
35 | diMeC31 | 3297 | 0.447 | 0.629 | 0.005 | 1.451 | 0.033 | 35 |
29 | 3-MeC31 | 3178 | 0.797 | 0.599 | 0.005 | 2.054 | 0.033 | 29 |
42 | C35:2 | 3460 | 1.275 | 1.223 | 0.004 | 1.211 | 0.030 | 42 |
2 | unknown | 1402 | 0.649 | 0.811 | 0.004 | 1.605 | 0.028 | 2 |
32 | 3-MeC32 | 3262 | 0.942 | 1.070 | 0.004 | 1.422 | 0.028 | 32 |
1 | C14 | 1400 | 0.771 | 0.925 | 0.004 | 1.521 | 0.027 | 1 |
36 | C34:1 | 3308 | 0.650 | 0.639 | 0.004 | 1.574 | 0.027 | 36 |
23 | 2-MeC30 | 3064 | 1.882 | 2.025 | 0.004 | 1.581 | 0.027 | 23 |
21 | C31:1 | 3029 | 0.227 | 0.066 | 0.004 | 2.473 | 0.027 | 21 |
37 | MeC33 | 3331 | 0.906 | 0.804 | 0.004 | 1.348 | 0.026 | 37 |
15 | 2-MeC28 | 2865 | 1.731 | 1.768 | 0.003 | 1.085 | 0.025 | 15 |
41 | C34:2 | 3377 | 0.355 | 0.277 | 0.003 | 1.587 | 0.024 | 41 |
40 | C34:2 | 3371 | 0.348 | 0.482 | 0.003 | 1.816 | 0.023 | 40 |
19 | MeC29 | 2960 | 0.369 | 0.489 | 0.003 | 1.580 | 0.023 | 19 |
The CH profiles of the virgin males and females differed qualitatively. SIMPER analyses, comparing conspecific males and females, revealed sex-specific compounds. In females the most abundant compounds were docosene (C22:1, RI 2182, CH11), hentriacontene (C31:1, RI 3082, CH24), 3-methyldotriacontane (3-MeC32, RI 3272, CH32) and tritriacontene (C33:1, RI 3280, CH33) (Table
Different patterns of CHs were detected between the two populations when constructing the heat map (Figure
Significant quantitative differences in the chemical CH profiles of the two populations of C. rosa have been demonstrated and complementary morphological analyses have confirmed that these two populations belong to two different morphotypes/genotypes, previously labelled by
The characteristic compounds of the lowland R1 type, diMeC28 and 3-MeC31, were present in higher relative amounts, whereas the highland R2 flies were characterised by high amounts of C22:1 and 2-MeC28. The compounds found in the present study correspond to the estimated chain lengths of the CH clusters identified in our earlier work for C. rosa, C. anonae, C. fasciventris and C. capitata, where the C. rosa R2 type could be determined based on the presence of even methylbranched hydrocarbons and the absence of odd methylbranched CHs when compared with the other three Ceratitis species (
The intraspecific variation in the CH profiles between the two types reported here might be a result of several different factors, such as the effects of temperature, the social context and diet (
In C. rosa, we found that the differences in cuticular hydrocarbon profiles between the two populations were greater than those between the sexes, although there was still a significant quantitative sexual dimorphism. Our findings are in agreement with studies conducted on Drosophila sp., where differences between D. montana populations were found to be considerably greater than those between the sexes (
It is important to note that the two populations of C. rosa studied here originate from different host plants, nevertheless they were reared during two generations on identical laboratory diet. The identified differences in the abundance of the CH between the populations and between the sexes may be, in addition to temperature and reproductive isolation factors, a result of the effects of host plants from which they originated (
Our data on cuticular hydrocarbon profiles, along with the previously published studies on morphology, genetics and sexual compatibility suggest that there exist two different entities, almost certainly unique biological species, within the taxa C. rosa from Kenya. In order to determine whether the different entities observed are consistent, the study needs to be extended to other populations of the two entities throughout their geographic and host ranges.
The funding was provided through the Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Prague (RVO: 61388963) and through research contracts 16106 and 16965 as a part of the FAO/IAEA Coordinated Research Project Resolution of Cryptic Species Complexes of Tephritid Pests to Overcome Constraints to SIT Application and International Trade.