Genetic stability, genetic variation, and fitness performance of the genetic sexing Salaya1 strain for Bactrocera dorsalis, under long-term mass rearing conditions
BMC Genetics volume 21, Article number: 131 (2020)
A genetic sexing strain (GSS) is an essential component for pest control using the sterile insect technique (SIT). A GSS is developed using a combination of Y-autosome translocation and a selectable marker such as pupal color, resulting in heterozygous males and homozygous females that possess wild-type brown pupae (wp+) and mutant white pupae (wp) alleles, respectively. The genetic sexing Salaya1 strain developed for Bactrocera dorsalis was evaluated using a clean stream and scaled-up for subsequent production lines (e.g., initiation, injection, and release). Colony management under small- and large-scale conditions for long-term rearing may affect the sexing system, genetic background, and fitness performance of the strain. Routine monitoring was applied to study genetic stability, genetic variation, and male mating competitiveness.
The percentage of recombinants was significantly different between males (wp) and females (wp+), ranging between 0.21–0.43% and 0.01–0.04%, respectively. Using 106 bands from six ISSR markers, the genetic backgrounds of two generations (F40 and F108) of the clean stream were found to be almost identical (0.960), and between those two generations and the wild population, the similarities were 0.840 and 0.800, respectively. In addition, the sterile males performed well in competitive mating with fertile females (Relative Sterility Index = 0.67 ± 0.13). The rates of fliers calculated from both clean and release streams were higher than 0.95. Regarding the fitness of the Salaya1 strain, the fertility and pupal recovery were similar in all production lines. The sex ratio (Male/Female) distortion was also recorded.
The Salaya1 strain reared at the mass-rearing facility retained its genetic stability, genetic variation, behavior (e.g., competitive mating and flight ability), and traits related to fitness for at least 10 consecutive generations. The filter rearing system is effective at minimising the selection pressure while maintaining the genetic background and fitness performances of the clean stream. These characteristics were stable throughout the production lines. In addition, the production efficiency is comparable among the different production lines and other similar types of GSSs.
The sterile insect technique (SIT) is a widely accepted pest population suppression method due to its species specificity and environmentally-friendly nature. Pest-free areas can also be developed, maintained and protected by the integration of SIT . The SIT program involves mass rearing and sterilisation by radiation of a target insect colony in a facility. Subsequently, the sterile insects are systematically released to compete against the wild male population. Hence, the reproductive potential of these wild males is reduced. One way to improve the effectiveness and cut operational costs of the SIT is to release high quality male-only sterile insects . The release of sterile females is not desirable because ovipositing by pest damages fruits. The presence of sterile females also distracts the mating between the released sterile males and wild females during SIT application. The elimination of females lowers the cost of mass rearing and field release.
For tephritid fruit flies, the development of genetic sexing strains (GSSs) is made possible by the separation of the male and female individuals in the mass rearing. GSSs have been developed based on a classical genetic approach, which requires a combination of Y-autosome translocation (T(Y;A)) and a selectable marker (e.g. white pupae (wp), black pupae (bp), temperature-sensitive lethal (tsl) mutations) . The Y-autosome translocation results in the linkage of Y chromosome and wild-type phenotype because they are heterozygous, one other allele carrying a wild-type allele on their Y-translocated chromosome and the other mutant allele on their free autosome. However, the females show a mutant phenotype because they are homozygous for the selectable marker. Vienna 7 , Vienna 8 , and Tapachula-7  GSSs are examples of this type.
The genetic stability of the sexing phenotype is one of the major concerns in the quality management system during mass rearing. This is because recombination can occur between the chromosomal break point of the T(Y;A) and the selectable marker. Its frequency is lower if the autosomal break point is closer to the location of the selectable marker . Rearing under stress conditions also increases the recombination frequency. The recombinant individuals, if they are not filtered out, they will accumulate. This leads to genetic breakdown and gradual loss of the sexing capacity in long-term mass rearing. The accumulation of recombinant individuals can be avoided by using a filter rearing system (FRS) [4,5,6].
The principle of FRS is to maintain a small standby GSS colony under minimal stress conditions (clean stream) where recombinants are physically selected out of the production. This clean stream must be cleaned off recombinants and be regularly used to refresh and/or continuously supply the mainstream of production lines with renewed genetic material. This can be accomplished using an amplification bridge that involves two steps, initiation and injection streams, which lead to the release stream. The initiation stream is used to increase the size of the clean stream, and the injection stream is used to inject new genetic material in to the release stream. Concurrently, the release stream is used to produce a large number of eggs and sorted male pupae for sterilisation and further field release . The routine evaluation of genetic stability and biological quality of a GSS for different populations along the production streams incorporated in the FRS is essential for quality control in a mass-rearing facility. The insect quality can be used to predict the sterile male performance in the field. In addition, it can be used as a baseline to measure the effects of changes and scaling up in the production system. Furthermore, the data can be used to design, test, and improve rearing conditions [4,5,6,7].
Adaptation of the mass-reared GSS colonies is another important factors that affects the sustainability of SIT. The genetic response of the GSSs to environmental conditions may result in a deviation of genetic background and biological properties from their own ancestors, related strains, and wild populations . Routine measurement of genetic diversity, mating competitiveness, traits related to fitness and rearing efficiency (fertility, egg-to-adult recovery rate, flight ability, and sex ratio) can be implemented for use in a quality control system. Although these parameters are specific to strain or mass-rearing protocol, they can be adopted as a standard protocol for GSS development and evaluation [7, 9,10,11].
The oriental fruit fly, Bactrocera dorsalis, is a key agricultural pest, causing reduction of fruit yield and limiting trade . This fruit fly is native to far eastern Asia and was later distributed throughout Southeast Asia and the rest of Asia [13,14,15]. It has also spread outside its native range to some areas of the USA (i.e., California, Florida, and Hawaii), continental Africa (http://www.cabi.org/ISC/datasheet/17685; Accessed 26 May 2020), and recently, southern Italy . The B. dorsalis Salaya1 GSS was developed based on the T(Y;A) translocation and wp gene. Their sex-linked traits produce brown-pupae males and white-pupae females . The Salaya1 samples from F5 and F17 of the clean stream were evaluated for genetic variation. This strain was able to maintain genetic polymorphisms and mating competitiveness that were comparable to the wild target population. The Salaya1 strain was proven successful in a pilot-scale SIT project at suppressing the B. dorsalis wild population . Currently, the clean stream is being reared for further amplification in a recently established mass-rearing facility.
In this study, we report on the genetic stability, genetic variation, mating competitiveness, and fitness parameters (e.g., fertility, egg-to-adult recovery rate, and flight ability) of the B. dorsalis Salaya1 strain in the continuous FRS for production. The Salaya1 strain shows low genetic recombination, indicating high genetic stability and the efficiency of the FRS. In addition, its genetic background and biological properties are stable under the long-term mass-rearing conditions. This suggests that the quality control data could be used as baselines for the mass rearing of the Salaya1 strain.
The genetic stability of the Salaya1 strain from clean, initiation, injection, and release streams were measured. Recombination occurring in the different streams was characterised as the accumulation of females emerging from brown pupae (wp+) (WT females) and males emerging from white pupae (wp) (wp males) under semi-mass rearing conditions over 10 consecutive generations (Table 1 and Additional file 1: Table S1). There was a significant difference between two types of recombinants (t = 6.094, df = 78, P < 0.001): The occurrence of WT females ranged from 0.01 to 0.04% while that of wp males ranged from 0.21 to 0.43%. Within the same type, the accumulated numbers of recombinants were not significantly different among the different streams (WT female: F3,36 = 0.843, P = 0.479; wp male: F3,36 = 0.693, P = 0.563).
Genetic variation was evaluated in two generations (F40 and F108) of the clean stream and the wild population, using six inter simple sequence repeat (ISSR) markers. All 106 bands were used to estimate the number of alleles (Na), the number of effective alleles (Ne), expected heterozygosity (He), Shanon’s information index (I), and percentage of polymorphic loci (%P). The means of all the parameters estimated from the two generations of the clean stream, F40 (Na = 1.538; Ne = 1.427; He = 0.248; I = 0.371; %P = 70.75%) and F108 (Na = 1.387; Ne = 1.408; He = 0.236; I = 0.349; %P = 62.26%), showed slightly lower values than the wild population (Na = 1.726; Ne = 1.410; He = 0.258; I = 0.401; %P = 85.85%).
The Nei’s genetic identity of F40 and F108 were very close (0.960) as between those two generations and the wild population, 0.840 and 0.800, respectively. The Principle Coordinate Analysis (PCoA) plot using a genetic distance matrix illustrates the genetic divergence of fruit fly individuals from two generations of the clean stream and the wild population (Fig. 1). The first axis accounts for 21.93% of the total variation, which separates the two generations of the clean stream from the wild population while the second axis (7.21% of total variation) does not show any biological significance. The PCoA result is congruent with the genetic identity.
Mating competitiveness of the Salaya1 strain
The mating competition was observed between the sterile males of the Salaya1 release stream and wild males with wild females. Overall, the fruit flies used in mating performance evaluation were sexually mature (Proportion of mating (PM) ≥ 0.20). The mating analyses showed that the mean relative sterility index (RSI) was 0.80 ± 0.11 and 0.54 ± 0.10 when males of the Salaya1 strain were irradiated at 50 and 60 Gy, respectively (Fig. 2).
Flight ability testing of the Salaya1 strain
Flight ability was evaluated with a standard protocol using adult emergence (%), adult fliers (%), and rate of fliers . These values were compared between the clean and release streams for 10 consecutive generations (Additional file 2: Table S2). Although all of the observed values from the release stream were significantly lower than the clean stream (adult emergence: t = 4.521, df = 10.359, P = 0.001; adult fliers: t = 5.050, df = 10.127, P < 0.001; rate of fliers: t = 3.300, df = 18, P = 0.004), the averages of adult emergence and adult fliers were both well above 70%. The average rate of fliers was higher than 0.95 in both streams.
Traits related to fitness within the Salaya1 and wild-type strains
Fertility, pupal recovery, adult emergence, and sex ratio were estimated from the B. dorsalis wild-type Phayathai1 strain and the production lines (e.g., clean, initiation, and injection streams) of the Salaya1 strain through 10 consecutive generations (Additional file 3: Table S3). The fertility (%) was similar in the different production lines of the Salaya1 strain but was significantly lower than in the Phayathai1 strain (F3,36 = 484.459, P < 0.001). The mean fertility in the Salaya1 strain, ranging from 47.06 ± 1.24% to 49.41 ± 0.80% among different streams, was approximately half of the value observed in the Phayathai1 strain (88.63 ± 0.47%). The pupal recovery rate (%) also indicated a similar trend (F3,36 = 33.559, P < 0.001). The mean pupation value of the Salaya1 strain varied from 75.78 ± 0.98% to 77.20 ± 1.17% while that of the Phayathai1 strain was 89.25 ± 1.20%.
The rates of adult emergence (%) were similar between strains and among production streams (F3,36 = 2.272, P = 0.097). The mean sex ratio between males and females observed in the Phayathai1 strain was close to 1:1 but was significantly different from the clean, initiation, and injection streams of the Salaya1 strain (F3,36 = 41.461, P < 0.001), for which the sex ratios were close to 1.3:1 (Additional file 3: Table S3).
The mean of egg-to-adult recovery rate in the Salaya1 strain, ranging from 26.73 ± 3.11 to 27.48 ± 2.14 among different streams, was significantly lower than in the Phayathai1 strain (61.29 ± 1.24) (F3,36 = 321.662, P < 0.001) (Additional file 3: Table S3).
Production efficiency of the Salaya1 strain under semi-mass rearing conditions
Egg and pupae production, pupa weight, and egg-to-adult recovery rate were estimated from different streams of the Salaya1 strain over 10 consecutive generations (Additional file 3: Table S3). The clean stream had significantly higher mean egg production (eggs/female/day) than the initiation, injection, and release streams (F3,32 = 8.743, P < 0.001). The three latter production lines produced means for eggs/female/day of similar values, ranging from 23.18 ± 2.44 to 27.67 ± 3.34. However, the pupae production and egg-to-adult recovery rate were not significantly different among various streams (pupae production: F2,27 = 0.36, P = 0.791; egg-to-adult recovery rate: F2,27 = 0.237, P = 0.791). The pupa weight of each stream was approximately 11 to 12 mg for both brown and white pupae, which is comparable to the pupa weight of the wild-type Phayathai1 strain.
Genetic stability of the Salaya1 strain
During long-term strain management, the sexing system may be lost, which causes strain breakdown, mostly either due to the recombination in the chromosome region between the translocation break point and the selectable markers (e.g., wp) or the survival of adjacent-1 segregates [2, 19]. Accumulation and selection in favor of recombinants can occur rapidly and depend on the environmental stress conditions [4, 19]. Monitoring for stability of the sexing system is necessary to ensure their positive characteristics . Firstly, a GSS must present a relatively low level of recombination due to its short distance between the break point and the selectable markers of the T(Y:A). Secondly, a GSS should produce no or few progenies related to adjacent-1 segregation during male meiosis (Fig. 3b). The progenies with autosomal deficiency usually die at the embryonic stage while the triplication type adjacent-1 progenies may be viable through to the adult stage . Lastly, a GSS should have only one T(Y;A) break point. Higher complexity of chromosomal rearrangement including a potential multiple translocation system could lead to the generation of more unbalanced gametes (autosomal deficiency or duplication of chromosomal regions) of the translocation heterozygotes. This situation would increase the sterility level of the strain [19, 21, 22]. The genetic sexing Salaya1 strain for B. dorsalis seems to be stable. The number of accumulated recombinants was relatively low in both types of recombinants: females in brown pupae (WT female; 0.02%) and males in white pupae (wp male; 0.36%). This suggests that the T(Y:A) break point is closely linked to the wp locus. Different ratios between the two types of recombinants were also observed in the other GSSs (e.g., C. capitata  and Z. cucurbitae ). According to the potential outcome of segregation and crossing over during male meiosis (Fig. 3b and c), the two types of recombinants should occur at an equal ratio. The reliance on only pupae color sorting may sometimes cause results that are not as sufficiently accurate as the results found for other GSSs, such as in the case of C. capitata , Z. cucurbitae , and A. ludens [3, 24]. The puparium depigmentation in the wild-type pupae can be affected by nutritional deficiency [25, 26]. The sorted brown pupae females of the Salaya1 occur at very much lower rate (0.02%). Hence, they affect very little in the male-only release system for SIT. In addition, the numbers of wp male recombinants did not significantly accumulate in all production lines, even during the 10 generations of circulating the release stream. In this study, the low recombination rate suggested that the FRS is practical and effective for minimising the genetic instability. In addition, the egg-to-adult recovery rate of the Salaya1 strain (27.48%) was significantly different from its ancestral wild-type Phayathai1 strain (61.29%) by roughly 50%. This 50% reduction is consistent with the theoretical expectation that the Salaya1 possesses a simple T(Y;A) translocation which is suitable for mass-rearing production [3, 22, 24, 27]. When a simple Y-autosome translocation involving only one autosome was cytogenetically confirmed in the other GSS of B. dorsalis, the egg-to-adult recovery rate was also approximately 50% .
Standards for the maintenance of genetic integrity in the Salaya1 strain
A set of characteristics that are unique to the DNA of populations are the chief distinguishing aspects of the genetic profile of GSSs. The integrity of a shared genetic profile can be lost during strain improvement and colony management. This could affect the mass-rearing efficiency and field performance of the sterile males. The same GSS may manifest varied biological properties (e.g., fertility, eggs to pupae rate, and adult emergence) when reared in different places and/or conditions . This situation has been generally observed in several established GSSs such as the Vienna 7 and Vienna 8 strains of C. capitata [7, 11], the 49-wp strain of B. dorsalis , and the species Z. cucurbitae . In this study, small samples of fruit flies were evaluated under laboratory conditions. In the scaled-up rearing, it was found that the biological properties were similar in all production lines. The acquired data could be used as reference for the Salaya1 strain. The data will allow standard comparison among Salaya1 colonies from different time points, production levels, and mass-rearing conditions. In doing so, the standards of biological quality control should include the rate of egg hatch, pupation, adult emergence, egg-to-adult recovery rate, and recombinants.
Maintenance of the genetic integrity is a primary consideration during the colonisation stage and in mass-rearing management of the strain. The heterozygosity can be affected by the genetic diversity of founding individuals and the degree of inbreeding in production lines . A sufficient number of rearing individuals is therefore an important factor for sustaining genetic variation, mass rearing efficiency, and field performance of the sterile males. Under the current quality management system, the Salaya1 strain has retained its genetic integrity for 10 years when compared to the wild population. It has been maintained at the level of approximately 5000 individuals per generation. Isasawin and colleagues  and this study tested the genetic background of the different generations of the Salaya1 clean stream and wild population using microsatellite and ISSR markers, respectively. Even though the ISSR markers are dominant markers, a sufficient number of marker bands (depending on analysis method, e.g., a modest number of 30–50 bands for AMOVA and generally 90 or more for STRUCTURE) and individuals (more than 10) could provide satisfactory results . The results from 106 ISSR bands were similar to those from four polymorphic microsatellite loci. However, the ISSR marker system required less time and expense, as well as provided data as DNA fingerprinting. DNA banding patterns of ISSR_02 to 05 are examples of markers for primary screening of genetic background in the colony (Additional file 4: Fig. S1 to S6). Therefore, the ISSR system may be a potential tool for large-scale monitoring of genetic variation in a mass-rearing facility .
Potential for improvement of Salaya1 mass rearing
In this study, we present the minimum viable production for mass rearing of the Salaya1 strain, at least for maintaining the main production via amplification bridge (initiation and injection streams) at regular intervals . However, the strain fertility can be improved to increase egg-to-adult recovery rate at an industrial scale. Eggs can be initially incubated at a constant temperature using a bubbling system (e.g., an air pump for a fish tank) before seeding on a larval diet [32,33,34]. This step promotes embryonic development by providing higher oxygen levels. In order to enhance a high number of larvae production, a greater number of seeding eggs placed into the larval food can compensate for the low fertility rate [28, 35]. However, to offset production loss, a surplus of seed eggs should be avoided as this caused diminishing return. In addition, characterisation of the strain using cytogenetic analyses can provide more accurate information on the structure of the chromosome rearrangement and mapping of the chromosome break point [3, 28]. Introducing chromosomal inversion on the autosome carrying the selectable markers can cause a recombination-reducing effect [2, 27]. This is exemplified in the case of C. capitata, in which an autosomal inversion had been introduced to a translocation T(Y;5)101 line that was not very stable . When the autosomal inversion carrying the selectable markers (wp and tsl) was introduced, the stability was very high. This led to a new strain, called Vienna 8D53+ which was later adopted in several mass-rearing facilities [2, 7, 27].
To prevent possible loss of genetic variation during colonisation of GSSs, new genetic material can be sometimes introduced by outcrossing to wild flies from the target population  or by inter-crossing between independent domesticated lines . This could maintain the level of genetic variation, as well as improve the rearing efficiency and field performance of the genetically refreshed strain .
Mating competitiveness data of the Salaya1 strain suggested that it has equal mating performance to the wild males when they were sterilised with a sterilisation dose of 60 Gy. This was estimated by the RSI value of 0.54 ± 0.10. A much higher value of RSI (0.80 ± 0.11) demonstrated that the Salaya1 sterile males can succeed more often in the mating competition with the wild females when they were irradiated at 50 Gy. This RSI value is unusally high when compared with the values from GSSs Vienna 8 (0.22 to 0.34) [38, 39] and Tapachula-7 (0.4) . However, this is consistent with a previous finding. With a 50 Gy sterility dosage, the Salaya1 sterile males could generate 50% sterility in the wild females in a similar mating competition scheme . It was found that more wild females were attracted to the sterile males than the wild males of C. capitata in field cages . The sterile males showed more extensive calling behavior than the wild males . Comparative study of the mating behavior between the Salaya1 strain and the wild types should be carried out to improve the SIT. The results stemming from different sterilising doses (50 versus 60 Gy) suggest that the small radiation dosages utilised in this experiment can strongly relate to the mating competitiveness of the males. Radiation dosimetry and sterilising dose should be critically studied to maximise the mating competitiveness of the Salaya1 sterile males.
The Salaya1 strain, reared in both relaxed conditions and high-density rearing conditions at the mass-rearing facility, can retain its genetic stability, genetic variation, rearing efficiency, and field performance (e.g., competitive mating and flight ability). The stable parameters imply that the FRS implemented in this study is suitable. Incorporation of an autosomal inversion and additional selectable marker lines to the Salaya1 are desirable to refine the strain. The Salaya1 strain improvement can be evaluated based on the standard monitoring of the FRS system suggested in this work. According to the consistency of rearing profiles throughout different production lines, the Salaya1 strain is comparable to other pupae-color based GSSs and is ready for industrial-scale rearing and SIT programs.
Fruit fly strains
The Salaya1 strain is a GSS based on brown-white pupal colour dimorphisms for sex sorting . This strain had been developed and continuously maintained at a laboratory-scale and modular mass-rearing scale as a clean stream since 2009. In 2016, the standard mass-rearing facility named ‘Regional R&D Training Center for Insect Biotechnology-RCIB’ was launched. The initiation stream was initially colonised. Two years later, the injection and release streams were established and used in this study.
The Phayathai1 strain is a wild-type bisexual strain and also the original strain used for development of the Salaya1 strain . This strain was also used as a reference strain for studying genetic variation [42, 43] and has been continuously maintained under the laboratory conditions for more than 15 years.
Wild samples were obtained from larvae and eggs collected from infested organic orchards of mixed fruits (e.g., guava, mango, rose apple, and banana) in Nakhon Pathom Province, Thailand. These samples were reared on their host fruits in the quarantine section at the RCIB. The flies from the F1 generation were used for genetic variation tests while the rest of the samples were used to generate the further generations. Only the F2-F4 generations were used for the mating competitiveness experiments to ensure wild samples still maintained their wild performance equivalency .
Laboratory and semi-mass rearing conditions
The wild-type Phayathai1 strain has been continuously reared in a laboratory cage with the dimensions of 0.35 × 0.45 × 0.35 m (width x length x height). Every generation, approximately 600 adult flies are used to maintain their colony. On the other hand, at a semi-mass rearing scale, the Salaya1 strain has been continuously reared in a mass-rearing cage with the dimensions of 0.4 × 1.0 × 1.8 m (width x length x height). Every generation, approximately 6000, 17,000, and 17,000 adults are used to establish their clean, initiation, and injection streams, respectively (Fig. 4). Only one generation of the injection stream was used for colonisation of a release stream, and the release stream was maintained similar to the initiation and injection streams. All strains and streams have been maintained under a photoperiod of 13:11 [L:D] h at 25–28 °C with 70–80% RH. The adult diet contains a mixture of hydrolysate yeast and sugar in a ratio of 1:3 by weight. Water is supplied to the cages by polypropylene cup with a cotton wick, and a bottle placed over a paper towel.
After 14–15 days of adult emergence, small perforated plastic tubes and bottles painted inside with guava juice were placed into the laboratory and mass-rearing cages, respectively, for egg collection. These eggs were seeded on a standard larval semisolid diet (26% wheat bran, 12% sugar, 3.6% hydrolysate yeast, 0.1% sodium benzoate, 0.1% methyl-p-hydroxyl benzoate, 0.2% acetic acid and 58% water by weight) . Eight to ten days after egg collection, mature larvae popped out of the food and formed pupae in the sawdust. All strains were parallel-reared using the same protocol. The Phayathai1 strain and Salaya1 clean stream were maintained in the standard insectary room (96 m2) while the Salaya1 initiation, injection, and release streams were maintained at the standard mass-rearing facility (697 m2).
Both strains and each stream were maintained over 10 consecutive generations for this study. The life cycle of each generation was approximately 1 month.
The presence of recombinants was evaluated by separating white and brown pupae and adults according to their sex. Materials used were the same as those from the section of fitness measurement (mentioned later). The most frequent recombination of the genetic sexing strain takes place between the translocation break point and wp allele (type 1) . This phenomenon produces males emerging from white pupae (wp) and females emerging from brown pupae (wp+). Observed recombinants from the type 1 recombination was estimated in each stream at three to five replications per generation for 10 consecutive generations. One-way ANOVA was used for the comparisons among different streams. Tukey’s HSD test was performed as a post hoc test for all pairwise comparisons. The statistical program PASW statistics 18 was employed for the analysis .
In addition, the rate of recombinants was assessed in 10 consecutive release stream generations. One thousand brown and white pupae were randomly collected from each generation, and emerged adults were observed according to their sex.
Twelve fruit fly samples from the Salaya1 clean stream (F40 and F108) and the wild population were randomly selected. Genomic DNA was individually extracted from adult fruit flies as per Baruffi et al. 1995 . The ISSR markers were used to test genetic variation of all samples because ISSR can generate a multiple banding pattern by amplification of the regions between microsatellites [31, 48, 49]. Six primers containing core microsatellite and anchor sequences and their selected bands are described in Additional file 5: Table S4. Each PCR was carried in a 20 μl volume containing 10 ng of genomic DNA, 1x buffer, 1.5 mM MgCl2, 0.2 mM dNTPs, 1 U Taq polymerase (Vivantis), and 40 pmol primer. PCR conditions used are as follows: 94 °C for 5 min; 3 cycles of 94 °C for 30 s, 58 °C for 1 min, 72 °C for 1 min; 3 cycles of 94 °C for 30 s, 56 °C for 1 min, 72 °C for 1 min; 3 cycles of 94 °C for 30 s, 54 °C for 1 min, 72 °C for 1 min; 30 cycles of 94 °C for 30 s, 52 °C for 1 min, 72 °C for 1 min; 72 °C and one cycle of final extension at 72 °C for 5 min.
PCR products were evaluated using 2% agarose gel at 50–70 V for 3–5 h, compared to 100 bp plus DNA ladder (Thermo Fisher Scientific). Gel contained 0.5 μmol ethidium bromide while the solution in the chamber contained 1X TBE buffer and 0.5 μmol ethidium bromide.
The presence and absence of each ISSR-PCR amplicon were scored as ‘1’ and ‘0’, respectively. Multiple ISSR-PCR banding patterns derived from different ISSR primers could be used together to create a 1/0 binary genotypic profile for each sample (Additional file 6: Data sheet). These diploid binary profiles were used as input data for GenAlEx v.6.503 using the Hardy-Weinberg equilibrium . The number of alleles (Na), the number of effective alleles (Ne), expected and unbiased expected heterozygosity (He and uHe, respectively), Shanon’s information index (I), Percentage of polymorphic loci (%P), and Nei’s genetic identity were estimated. PCoA was analysed, and a two-dimensional graph was plotted using a subprogram in GenAlEx .
Measurements of strain performances
Mating competitiveness test
The experiments were carried out according to the standard protocol from FAO/IAEA/USDA . This was done in order to observe the ability of sterile males of the Salaya1 strain to compete with fertile wild males for mating with fertile wild females. The brown pupae from the release streams (F7-F9) were sterilised with a radiation dose of 50 and 60 Gy (MultiRad160, Faxitron) during their late pupal stage (2 days before emergence). The irradiation doses of 50 and 60 Gy were chosen because they provided 100% male sterility when approximately 2000 eggs were tested in three replications. Eye color change was used in parallel to assure the development of pupae. Emerged sterile males were maintained under a photoperiod of 13:11 [L:D] h at 25–28 °C with the adult diet and the displayed courtship behavior (e.g., male calling) was observed to ensure the sexual maturity. The tested wild flies were sex sorted within 1 day of emergence and similarly reared. All flies were maintained in the same manner as sterile males. At least 4 h before the mating experiment, the sterile and wild males were spot-marked with different water-based colors (white and gray, respectively) at their dorsal thoraxes. Rotation of the water-based dye marking between the sterile and wild males was not carried out due to limitations on the availability of the wild samples. One possible concern is that the colour could have an effect on the mating choices made by females, however, there have been no such observations from similar experiments [38, 39]. The study was conducted in three outdoor field cages 3 × 3.5 × 2.3 m (width x length x height) with a potted mango tree inside each. In each mating competitiveness test, 20 sterile males and 20 wild males were released 30 min earlier to allow lekking formation in the cages. Twenty wild females were subsequently released at 17:00. Only 20 wild males and 20 wild females were released in the control cages. Mating pairs were observed between 18.00–19.30. The mating couples of either sterile or wild males to the wild females were individually collected in plastic vials and counted.
Data from the control cage was used to calculate the proportion of mating (PM) . This parameter ensures the readiness of the flies and adequacy of environmental conditions for the mating test. The PM value was calculated by the number of mating pairs divided by the total number of released females. If the PM value was less than 0.2, the experiments were discarded. The relative sterility index (RSI) was calculated to determine the male sexual competitiveness. The values of RSI range from 0 (i.e., all of wild females mate with wild males) to 1 (i.e., all wild females mate with sterile males). A value of 0.5 indicates an equal mating performance ratio for sterile and wild males. The mating competitiveness experiment was repeated in triplicate, and the RSI values were compared using the Student’s t-test in PASW statistics 18 .
Flight ability testing
The experiments were carried out to estimate the ability of the brown pupae to emerge and develop into flying male adults, according to the standard protocol from FAO/IAEA/USDA . The brown pupae samples were taken from the clean and release streams. The experiment was performed in a flight testing apparatus; (i.e., black-painted plexiglass tube: outside diameter 8.9 cm with 3 mm wall thickness and 10.0 cm height). One hundred pupae were loaded into a closed flight testing apparatus approximately 2 days before their expected emergence. The pupae-loaded flight testing apparatus was placed in a transparent plastic container. The lid of the apparatus had been opened before the plastic container was closed to allow the adults to fly out. There are five outcome categories; non-emerged pupae, partially-emerged pupae, deformed adults, non-flier adults, and flier adults. The data from these categories were recorded daily and summarised after a period of 7 days.
The percentage of adult emergence was calculated by the total number of fully-emerged flies (i.e., deformed flies, non-fliers, and fliers) divided by the total number of pupae whereas the percentage of fliers was calculated from the number of fliers divided by the total number of pupae. The rate of fliers equals the percentage of fliers divided by the percentage of emergence . The flight ability test was done in five replicates. The Student’s t-test was used for the comparisons among different streams using PASW statistics 18 .
Measurements of traits related to fitness and production efficiency of the strain
Samples were randomly selected, observed, and recorded at all stages of development. The quality control parameters, including fertility (egg-hatching rate), egg to pupae recovery (pupal recovery rate), pupa weight, adult emergence, and sex ratio of the Salaya1 strain, were tested following the standard protocol .
Fertility was evaluated from 1400 eggs per replicate. The eggs were aligned on a piece of wet paper and placed on larval food in a small cup. Eggs were incubated at 25–28 °C for 3 days, and the egg-hatching rate was determined. This standard was followed for every generation of both strains and each stream collected in three to five replicates for 10 consecutive generations.
Pupal recovery was calculated by the number of pupae divided by the number of hatched eggs. Pupa weight was calculated by measuring the weight of 20 randomly selected pupae. Three replicates per generation were performed and adult emergence was evaluated by the number of normal adults divided by the number of pupae. The sex ratio of males to females was also calculated.
Egg-to-adult recovery rate was calculated to evaluate the production efficiency. This parameter was estimated by the number of normal adults divided by the number of reared eggs.
One-way ANOVA was used for the comparisons among different strains and streams. Tukey’s HSD test was performed as a post hoc test for all pairwise comparisons. The statistical program PASW statistics 18 was employed for the analysis .
Egg production was evaluated by the number of eggs produced per female per day. The number of eggs in 1.0 mL laid on the wet paper were counted three times, and the mean value was 17,500 ± 2% eggs/ml. This value was used to multiply with the total volume of eggs collected daily and divided by 74% of the number of white pupae for calculation of the egg production. This is because approximately 74% of pupae were emerged, and only females emerged from white pupae. Under the semi-mass rearing conditions, eggs were reared on larval food in a tray, and the fertility was not monitored. Pupae production was used instead of pupal recovery and was evaluated by the number of pupae divided by the number of reared eggs.
Egg-to-adult recovery rate was calculated to evaluate the production efficiency at the mass-rearing scale when all processes were done in the facility. This parameter was estimated by the number of normal adults divided by the number of reared eggs. The number for adult emergence was estimated from the number of pupae multiplied by 74%. One-way ANOVA was used for the comparisons among streams. Tukey’s HSD test was performed as a post hoc test for all pairwise comparisons. The statistical program PASW statistics 18 was employed for the analysis .
Availability of data and materials
All supporting data are included as additional information.
Analyses of Variance
Filter Rearing System
Genetic sexing strain
Shanon’s information index
Inter simple sequence repeat
Number of alleles
Number of effective alleles
Polymerase chain reaction
Principle Coordinate Analysis
Percentage of polymorphic loci
Proportion of mating
Relative sterility index
Sterile insect technique
Unbiased expected heterozygosity
Klassen W, Curtis CF. History of the sterile insect technique. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile insect technique: principles and practice in area-wide integrated Pest management. Dordrecht: Springer; 2005. p. 3–36.
Franz G. Genetic sexing strains in Mediterranean fruit fly, an example for other species amenable to large-scale rearing for the sterile insect technique. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile insect technique: principles and practice in area-wide integrated Pest management. Dordrecht: Springer; 2005. p. 427–51.
Zepeda-Cisneros CS, JSM H, Garcia-Martinez V, Ibañez-Palacios J, Zacharopoulou FG. Development, genetic and cytogenetic analyses of genetic sexing strains of the Mexican fruit fly, Anastrepha ludens Loew (Diptera: Tephriitidae). BMC Genetics. 2014;15(suppl 2):51.
Fisher K, Cáceres C. A filter rearing system for mass reared genetic sexing strains of Mediterranean fruit fly (Diptera: Tephritidae). In: Tan KH, editor. Area-Wide Control of Fruit Flies and other Insect Pests. Penang: Penerbit Universiti Sains Malaysia; 2000. p. 543–50.
Cáceres C. Mass rearing of temperature sensitive genetic sexing strains in the Mediterranean fruit fly (Ceratitis capitata). Genetica. 2002;116:107–16.
Cáceres C, Cayol JP, Enkerlin W, Franz G, Hendrichs J, Robinson AS. Comparison of Mediterranean Fruit Fly (Ceratitis capitata) (Tephritidae) Bisexual and Genetic Sexing Strains: Development, Evaluation and Economics. In: Braenes BN, editor. Proceedings of the 6th International Symposium on Fruit Flies of Economic Importance. Irene: Isteg Scientific Publications; 2004. p. 367–81.
Augustinos AA, Targovska A, Canino-Martinez E, Schorn E, Franz G, Cáceres C, et al. Ceratitis capitata genetic sexing strains: laboratory evaluation of strains from mass-rearing facilities worldwide. Entomol Exp Appl. 2017;164:305–17.
Krafsur ES. Role of population genetics in the sterile insect technique. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile insect technique: principles and practice in area-wide integrated Pest management. Dordrecht: Springer; 2005. p. 389–406.
Rempoulakis P, Tarat G, ul Haq I, Wornayoprn V, Ahmad S, Tomas US, et al. Evaluation of quality production parameters and mating behavior of novel genetic sexing strains of the Mediterranean fruit fly Ceratitis capitata (Wiedemann) (Diptera: Tephritidae). PLoS One. 2016;11(6):e0157679.
Meza JS, ul Haq I, MJB V, Bourtzis K, Kyritsis GA, Cáceres C. Comparison of classical and transgenic genetic sexing strains of Mediterranean fruit fly (Diptera: Tephritidae) for application of the sterile insect technique. PLoS One. 2018;13(12):e0208880.
Barnes B, Rosenberg S, Arnolds L, Johnson J. Production and quality assurance in the SIT Africa Mediterranean fruit fly (Diptera: Tephritidae) rearing facility in South Africa. Fla Entomol. 2007;90(1):41–52.
Clarke AR, Armstrong KF, Carmichael AE, Milne JR, Raghu S, Roderick GK, et al. Invasive phytophagous pests arising through a recent tropical evolutionary radiation: the Bactrocera dorsalis complex of fruit flies. Annu Rev Entomol. 2005;50(1):293–319.
Aketarawong N, Bonizzoni M, Thanaphum S, Gomulski LM, Gasperi G, Malacrida AR, et al. Inference on the population structure and colonization process of the invasive oriental fruit fly, Bactrocera dorsalis (Hendel). Mol Ecol. 2007;16:3522–32.
Aketarawong N, Guglielmino G, Karam N, Falchetto M, Manni M, Scolari F, et al. The oriental fruitfly Bactrocera dorsalis s.s. in East Asia: disentangling the different forces promoting the invasion and shaping the genetic make-up of populations. Genetica. 2014;142(3):201–13.
Qin YJ, Krosch MN, Schutze MK, Zhang Y, Wang XX, Prabhakar CS, et al. Population structure of a global agricultural invasive pest, Bactrocera dorsalis (Diptera: Tephritidae). Evol Appl. 2018;11:1990–2003.
Nugnes F, Russo E, Viggiani G, Bernardo U. First record of an invasive fruit fly belonging to Bactrocera dorsalis complex (Diptera: Tephritidae) in Europe. Insect. 2018;9(40):182.
Isasawin S, Aketarawong N, Thanaphum S. Characterization and evaluation of microsatellite markers in a strain of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae), with a genetic sexing character used in sterile insect population control. Eur J Entomol. 2012;109:331–8.
FAO/IAEA/USDA. Product Quality Control for Sterile Mass-Reared and Released Tephritid Fruit Flies. Version 7.0. Vienna: International Atomic Energy Agency; 2019.
Franz G. Recombination between homologous autosomes in medfly (Ceratitis capitata) males: type-1 recombination and the implications for the stability of genetic sexing strains. Genetica. 2002;116:73–84.
Kerremans P, Franz G. Isolation and cytogenetic analyses of genetic sexing strains for the medfly, Ceratitis capitata. Theor Appl Genet. 1995;91:225–61.
Zacharopoulou A, Franz G. Genetic and cytogenetic characterization of genetic sexing strains of Bactrocera dorsalis and Bactrocera cucurbitae (Diptera: Tephritidae). J Econ Entomol. 2013;106(2):995–1003.
Franz G. The “combi fly concept” revisited: how much radiation is required to sterilize males of a genetic sexing strain? In: Tan KH, editor. Area-Wide Control of Fruit Flies and other Insect Pests. Penang: Penerbit Universiti Sains Malaysia; 2000. p. 511–6.
McInnis DO, Tam S, Lim R, Komatsu J, Albrecht C. Development of a pupal-color based genetic sexing strain of the melon fly, Bactrocera cucurbitae (Diptera: Tephritidae). Ann Entomol Soc Am. 2004;97:1026–33.
Meza JS, Cáceres C, Bourtzis K. Slow Larvae mutant and its potential to improve the pupal color-based genetic sexing strain system in Mexican fruit fly, (Diptera: Tephritidae). J Econ Entomol. 2019;112(4):1604–10.
Lee KP, Simpson SJ, Wilson K. Dietary protein-quality influences melanization and immune function in an insect. Funct Ecol. 2008;22:1052–61.
Andersen SO. Insect cuticular sclerotization: a review. Insect Biochem Mol Biol. 2010;40:166–78.
Zacharopoulou A, Augustinos AA, Drosopoulou E, Tsoumani KT, Gariou-Papalexiou A, Franz G, et al. A review of more than 30 years of cytogenetic studies of Tephritidae in support of sterile insect technique and global trade. Entomol Exp Appl. 2017;164:204–25.
McCombs S, Saul SH. Translocation-based genetic sexing system for the oriental fruit fly (Diptera: Tephritidae) based on pupal color dimorphism. Ann Entomol Soc Am. 1995;88(5):695–8.
Parker AG. Mass-rearing for sterile insect release. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile insect technique: principles and practice in area-wide integrated Pest management. Dordrecht: Springer; 2005. p. 209–32.
Nelson MF, Anderson NO. How many marker loci are necessary? Analysis of dominant marker data sets using two popular population genetic algorithms. Ecol Evol. 2013;3(10):3455–70.
Abbot P. Individual and population variation in invertebrates revealed by inter-simple sequence repeats (ISSRs). J Insect Sci. 2001;1(8):1–3.
Moreno MP, Gomez Simuta Y, Zavila Lopez JL. Procedures for mass-rearing the west Indian fruit fly, Anastrepha obligua. https://inis.iaea.org/collection/NCLCollectionStore/_Public/30/007/30007844.pdf?r=1&r=1, Accessed 16 Oct 2019.
Orozco-Dávilla D, Quintero L, Hemández E, Solís E, Artiaga T, Hernández R, et al. Mass rearing and sterile insect releases for the control of Anastrepha spp. pests in Mexico – a review. Entomol Exp Appl. 2017;164:176–87.
Dominguez Gordillo JC. Mass rearing methods for fruit fly. https://inis.iaea.org/collection/NCLCollectionStore/_Public/30/007/30007845.pdf, Accessed 16 Oct 2019.
Fezza T, Geib SM, Shelly TE. Comparative rearing parameters for bisexual and genetic sexing strains of Zeugodacus cucurbitae and Bactrocera dorsalis (Diptera: Tephritidae) on an artificial diet. J Asia Pac Entomol. 2018;21:283–7.
Franz G, Kerremans P, Rendon P, Hendrichs J. Development and application of genetic sexing systems for the Mediterranean fruit fly based on a temperature sensitive lethal. In: McPheron BA, Steck GJ, editors. Fruit fly pests: a world assessment of their biology and management. United States: St. Lucie Press; 1996. p. 185–91.
Gilchrist AS, Cameron EC, Sved JA, Meats AW. Genetic consequences of domestication and mass rearing of pest fruit fly Bactrocera tryoni (Diptera: Tephritidae). J Econ Entomol. 2012;105(3):1051–6.
Paranhos BJ, McInnis D, Uramoto K, Damasceno I, Gonçalves N, Alves RM, et al. Sterile medfly males of the TSL Vienna 8 genetic sexing strain display improved mating performance with ginger root oil. In: Sugayama RL, Zucchi RA, Ovruski SM, Sivinski J, editors. Fruit flies of economic importance: from basic to applied knowledge. Salvador: SBPC; 2008. p. 313–8.
Virginio JF, Gómez M, Pinto AM, Aniely GG, Paranhos BJ, Gava CA, Cáceres C, Walder JM. Male sexual competitiveness of two Ceratitis capitata strains, tsl Vienna 8 and OX3864A transgenics, in field cage conditions. Entomol Exp Appl. 2017;164:318–26.
Orozco D, Meza JS, Zepeda S, Solís E, Quintero-Fong JL. Tapachula-7, a new genetic sexing strain of the Mexican fruit fly (Diptera: Tephritidae): sexual compatibility and competitiveness. J Econ Entomol. 2013;106(2):735–41.
Calkins CO, Parker AG. Sterile insect quality. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile insect technique: principles and practice in area-wide integrated Pest management. Dordrecht: Springer; 2005. p. 269–96.
Aketarawong N, Bonizzoni M, Malacrida AR, Gasperi G, Thanaphum S. Seventeen novel microsatellite markers from an enriched library of the pest species Bactrocera dorsalis sensu stricto. Mol Ecol Notes. 2006;6(4):1138–40.
Laohakieat K, Aketarawong N, Isasawin S, Thitamadee S, Thanaphum S. The study of the transformer gene from Bactrocera dorsalis and B correcta with putative core promoter regions. BMC Genetics. 2016;17:34.
Schutze MK, Dammalage T, Jessup A, Vreysen MJB, Wornoayporn V, Clarke AR. Effects of laboratory colonization on Bactrocera dorsalis (Diptera, Tephritidae) mating behavior: ‘what a difference a year makes’. ZooKeys. 2015;540:369–83.
Sutantawong M, Orankanok W, Enkerlin WR, Womoayporn V, Cáceres C. The sterile insect technique for control of the oriental fruit fly, Bactrocera dorsalis (Hendel), in mango orchards in Ratchaburi Province, Thailand. In: Braenes BN, editor. Proceeding of 6th International Symposium on Fruit Flies of Economic Importance. Irene: Isteg Scientific Publications; 2004. p. 223–32.
PASW. Statistics for windows, version 18.0. Chicago: SPSS Inc; 2009.
Baruffi L, Damini G, Guglielmino CR, Bandi C, Malacrida AR, Gasperi G. Polymorphism within and between populations of Ceratitis capitata: comparison between RAPD and multilocus enzyme electrophoresis data. Heredity. 1995;74:425–37.
Gupta M, Chyi Y-S, Romero-Severson J, Owen JL. Amplification of DNA markers from evolutionarily diverse genomes using single primers of simple-sequence repeats. Theor Appl Genet. 1994;89:998–1006.
Steffler LM, Dolabella SS, Ribolla PEM, Dreyer CS, Araújo ED, Oliveira RG, et al. Genetic variability and spatial distribution in small geographic scale of Aedes aegypti (Diptera: Culicidae) under different climatic conditions in northeastern Brazil. Parasite Vectors. 2016;9:530.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012;28:2537–9.
The authors gratefully acknowledge Mr. Robert Bachtell Eastland for his professional English editing service.
About this supplement
This article has been published as part of BMC Genetics Volume 21 Supplement 2, 2020: Comparing rearing efficiency and competitiveness of sterile male strains produced by genetic, transgenic or symbiont-based technologies. The full contents of the supplement are available online at https://bmcgenet.biomedcentral.com/articles/supplements/volume-21-supplement-2.
This research project was supported by International Atomic Energy Agency research contract no. 18795 as part of the Agency’s Coordinated Research Project “Comparing Rearing Efficiency and Competitiveness of Sterile Male Strains Produced by Genetics, Transgenic or Symbiont-based Technologies” to ST and NA. Additionally, S. Thanaphum also acquired joint research funding from the Faculty of Science and Mahidol University for ‘Setting up Regional R&D Training Center for Insect Biotechnology’. All of the funding bodies played no role in the design of the study and collection, analysis, and interpretation of data, and in writing the manuscript. Publication costs are funded by the Joint FAO/IAEA Division of Nuclear Techniques in Food and Agriculture, IAEA (CRP No.: D4.20.16) Vienna, Austria.
Ethics approval and consent to participate
This research project was a part of a project that obtained animal ethics approval from Mahidol University-Institute Animal Care and Use Committee (MU-IACUC) (CoA.No. MUSC 61–044-446).
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Genetic stability of the Salaya1 strain in the continuous filter rearing system under semi-mass rearing conditions.
Table S2. Flight ability (mean ± standard error) of the Salaya1 clean and release streams.
Table S3. Production efficiency of the Salaya1 strain in each production line under laboratory and semi-mass rearing conditions.
Figure S1. Results of ISSR_01 marker analysis. Figure S2. Results of ISSR_02 marker analysis. Figure S3. Results of ISSR_03 marker analysis. Figure S4. Results of ISSR_04 marker analysis. Figure S5. Results of ISSR_05 marker analysis. Figure S6. Results of ISSR_06 marker analysis.
Table S4. Six polymorphic ISSR primers and their matrices of obtained band sizes.
Figure S7. Raw gel images from ISSR_01 marker analysis. Figure S8. Raw gel images from ISSR_02 marker analysis. Figure S9. Raw gel images from ISSR_03 marker analysis. Figure S10. Raw gel images from ISSR_04 marker analysis. Figure S11. Raw gel images from ISSR_05 marker analysis. Figure S12. Raw gel images from ISSR_06 marker analysis.
About this article
Cite this article
Aketarawong, N., Isasawin, S., Laohakieat, K. et al. Genetic stability, genetic variation, and fitness performance of the genetic sexing Salaya1 strain for Bactrocera dorsalis, under long-term mass rearing conditions. BMC Genet 21 (Suppl 2), 131 (2020). https://doi.org/10.1186/s12863-020-00933-4
- Oriental fruit fly
- Male recombination
- Laboratory adaptation
- Filter rearing system (FRS)
- Routine monitoring