Irradiation induced inversions suppress recombination between the M locus and morphological markers in Aedes aegypti

Background Aedes aegypti is the primary vector of arthropod-borne viruses and one of the most widespread and invasive mosquito species. Due to the lack of efficient specific drugs or vaccination strategies, vector population control methods, such as the sterile insect technique, are receiving renewed interest. However, availability of a reliable genetic sexing strategy is crucial, since there is almost zero tolerance for accidentally released females. Development of genetic sexing strains through classical genetics is hindered by genetic recombination that is not suppressed in males as is the case in many Diptera. Isolation of naturally-occurring or irradiation-induced inversions can enhance the genetic stability of genetic sexing strains developed through genetically linking desirable phenotypes with the male determining region. Results For the induction and isolation of inversions through irradiation, 200 male pupae of the ‘BRA’ wild type strain were irradiated at 30 Gy and 100 isomale lines were set up by crossing with homozygous ‘red-eye’ (re) mutant females. Recombination between re and the M locus and the white (w) gene (causing a recessive white eye phenotype when mutated) and the M locus was tested in 45 and 32 lines, respectively. One inversion (Inv35) reduced recombination between both re and the M locus, and wand the M locus, consistent with the presence of a rather extended inversion between the two morphological mutations, that includes the M locus. Another inversion (Inv5) reduced recombination only between w and the M locus. In search of naturally-occurring, recombination-suppressing inversions, homozygous females from the red eye and the white eye strains were crossed with seventeen and fourteen wild type strains collected worldwide, representing either recently colonized or long-established laboratory populations. Despite evidence of varying frequencies of recombination, no combination led to the elimination or substantial reduction of recombination. Conclusion Inducing inversions through irradiation is a feasible strategy to isolate recombination suppressors either on the M or the m chromosome for Aedes aegypti. Such inversions can be incorporated in genetic sexing strains developed through classical genetics to enhance their genetic stability and support SIT or other approaches that aim to population suppression through male-delivered sterility. Supplementary Information The online version contains supplementary material available at 10.1186/s12863-020-00949-w.


Background
Aedes aegypti is the primary vector of arthropod-borne viruses (arboviruses) such as dengue [1,2], chikugunya [3], Zika [4,5], and the yellow fever [6,7]. It is one of the most widespread and invasive mosquito species globally, originated from Western Africa and spread worldwide in the past 70 years through human trading and travelling activities [8,9].
Due to the lack of efficient specific drugs or vaccination strategies (except yellow fever) to impede disease transmission [10,11], vector population control methods are receiving renewed interest. Main strategies so far mostly rely on the extensive use of insecticides and the community engagement for habitat management. However, both tools have been proven inefficient due to emerging insecticide resistance in Ae. aegypti populations, negative environmental and ecological impact of pesticide use, and difficulty in identifying and destroying mosquito breeding sites, particularly the cryptic ones [12][13][14][15][16][17]. Therefore, more effective, sustainable, and environmentally friendly control approaches are needed, including genetically based population suppression methods, such as the sterile insect technique (SIT) and other related methods, all of which rely mainly on the induction of sterility in natural populations.
Sterile males can be produced by irradiation-based (SIT), symbiont-based (incompatible insect technique, IIT), combined SIT/IIT or transgenic approaches [18][19][20][21][22][23][24]. An important requirement for any of them is the availability of an efficient and robust sex separation system that will result to the accurate separation of males from females [21]. Sex separation is currently possible in Ae. aegypti using mechanical tools based on pupal sexual dimorphism but their efficacy depends heavily on the rearing conditions [25][26][27]. Furthermore, these methods are currently appropriate only for small scale manipulations, hence it is essential to develop advanced sexing strategies based on genetic and molecular approaches for large scale mosquito male releases, including genetic sexing strains (GSS) [21,26]. The development of a GSS through classical genetics requires two basic components: a) a visible or conditional recessive mutation that can be used as a selectable marker (e.g. eye or pupal color, insecticide resistance, etc.) and b) the linkage of the phenotypic marker to the sex-determining genetic locus [28]. There is a renewed interest in revival and application of SIT for mosquitos [24,[29][30][31], as evidenced also by recent advances in refinement of irradiation doses and drone-mediated releases in the field as well as the recently published guidance framework and phased conditional approach for testing the SIT against Aedes mosquitoes [32][33][34]. Such strategies will be greatly facilitated by the availability of GSS developed through classical genetics, since there is documented efficiency through time and wide public acceptance [35].
In addition, recent studies suggest that the impact of the irradiation on the biological quality of the males does not severely compromises their performance [31,36,37].
In Aedes mosquitoes, male development depends on a dominant male-determining locus (M-locus) that resides on a homomorphic sex-determining chromosome [38][39][40]. The M-locus has been mapped to chromosome 1, band 1q21 and the Ae. aegypti males are heterogametic (Mm) while the females are homogametic (mm) [41]. Recent studies suggest that the Nix gene has the properties of being the M-locus since it is both 'required and sufficient to initiate male development' [40,42]. Nix resides in a genomic region that is protected from recombination and other functional genes have been identified tightly linked with Nix, with myo-sex being almost exclusively found in males and only sporadically in females, due to recombination [43]. An ideal selectable marker for Ae. aegypti would be located on chromosome 1, closely linked to M-locus. In order to construct a GSS, the wild-type allele of the selectable marker should be physically linked to the maledetermining factor which in the Diptera model of GSS, the Mediterranean fruit fly, Ceratitis capitata (medfly), has been accomplished by Y-autosome translocation(s). In the resulting strain, the males are heterozygous with a normal "wild-type" phenotype. On the other hand, the females are homozygous for the recessive alleles of the selectable marker thus exhibiting the mutant phenotype and can be separated from males [28].
There is a wide range of potential markers known from previous studies that could be used for GSS development in Ae. aegypti. These are related to eye color [44][45][46][47], insecticide resistance [48,49], and body color [50,51]. Some of the eye color genes, such as the red eye (re) and the white (w), have already been mapped to the sex determining chromosome, which makes them suitable selectable markers since there is no need of inducing chromosome 1-autosome translocations [44,45] (Fig. 1a). However, the presence of genetic recombination in males of Aedes species (unlike many other Diptera) would reduce the genetic stability of such a GSS.
Chromosomal inversions are quite common in Aedes species and are considered as an important factor contributing to speciation [52,53]. Inversions can also be detected within species or populations undergoing speciation (incipient speciation) and small inversions can be widely spread in the genome. Intraspecies inversions can have multiple effects in respect to recombination frequencies and fitness [54][55][56]. In the case of homozygous inversions, where genomic regions of both homologous chromosomes are inverted compared to the 'standard' orientation, changes in recombination frequencies are expected, consistent with the effect of the inversion on the chromosomal localization of the markers studied. At the same time, no major negative effects on fitness are expected, unless inversion interferes with gene expression of important genes. On the other hand, heterozygous inversions (where a genomic region of one of the two homologues has been inverted) are considered as recombination suppressors, since the recombination is eliminated within the inversion, a phenomenon that can be extended also outside the inversion [54,55]. Recombination suppression in this case mainly happens because recombination within the inversion leads to the production of imbalanced gametes, which are eliminated during development. Therefore, longer heterozygous inversions may have a cost on the fitness of a strain (in respect to productivity).
The construction of GSSs in the medfly has been facilitated by the fact that male recombination is highly suppressed [28]. The development of a filtering system for the mass rearing facilities further supported the adaptation of GSSs in large scale rearing [28,57]. The incorporation of additional tools, such as the D53 inversion (Inv D53) further enhanced the genetic stability of the GSS [28,58,59]. However, in Aedes species, recombination frequency in males is not suppressed. Therefore, GSS developed through classical genetics may be unstable and additional elements that suppress recombination are needed. Recombination suppressing inversions have been developed in the past for both M and m chromosomes [60], leading to reduced recombination among the M locus and morphological markers of chromosome I.
The aim of the present study was to induce through irradiation (or identify naturally-occurring) inversions on the M chromosome of Ae. aegypti that would suppress recombination between re and M locus and/or w and M locus. These two morphological markers are being considered for the development of GSSs, therefore reducing recombination with M locus is critical for the enhancement of the genetic stability of resulting strains.

Establishment of iso-male lines following irradiation
Approximately 300 male pupae of the 'BRA' strain were irradiated in batches of~100 pupae each. Following the crossing scheme described in Methods section, one hundred iso-male lines were set up (Fig. 1b). Eggs from three consecutive gonotrophic cycles were collected and iso-male lines produced approximately 53 eggs per line (Additional file 1 Table S1). The first gonotrophic cycle produced most of the eggs (34.33 eggs per line), followed by the second (14.04 eggs per line) and the third (4.19 eggs per line). Only 57 of the one hundred lines produced eggs in at least one of the gonotrophic cycles. After excluding lines that produced zero eggs, the average number of eggs was 92.21 per line.

Isolation of recombination suppressors (re -M)
Larvae hatched efficiently from 45 of the 57 lines that produced eggs (Table 1), allowing estimation of recombination frequency in F2 and, if needed, F3. As a rule of thumb and considering that recombination between the re and the M locus has been reported to vary between 1 and 7%, all lines that showed less than 2% percentage of recombinants in F2 were upscaled and evaluated again in F3. Therefore, thirteen lines were upscaled, as shown

Isolation of recombination suppressors (w -M)
A total of 32 iso-male lines were successfully crossed and backcrossed with homozygous w mutant virgin females thus allowing the estimation of recombination between w and the M locus (Fig. 1b). Since recombination between these two loci has been reported to usually vary between 8 and 16%, lines that showed less than 6% of recombinants in F2 were evaluated again in F3. Only three lines (5, 35, and 67) gave less than 6% recombination in F2 ( Table 2). One of the lines (67) showed increased recombination frequency in F3 and was not further followed up. Two of the lines showed consistent recombination suppressing properties. Line 5 showed recombination frequencies of approximately 1.4% in generations F2 and F3, with a rather large sample size tested (more than 1000 individuals per generation) (χ 2 = 250.25; df = 1; p < < 0.001). Even more promising, line 35 exhibited stable recombination suppressing properties down to 1% in F2 and F3 (χ 2 = 904.14; df = 1; p < < 0.001), with a sample size of more than 4000 and 5000 mosquitos screened, respectively.

Searching of naturally-occurring recombination suppressors (re -M and w -M)
To identify naturally-occurring recombination suppressors among re and the M locus, males from sixteen wild type strains were crossed with homozygous re mutant virgin females 'en masse'. All F1 males exhibited the wild type phenotype and were backcrossed to virgin re females. Recombination frequencies were measured for one (F2) or two (F2 and F3) generations. Despite the presence of varying levels of recombination, there was no combination leading to an important recombination suppression (Additional file 2 Table S2). Following the same approach, we searched for naturally-occurring inversions that could act as recombination suppressors by crossing homozygous w mutant females with males from fourteen different genomic backgrounds. Again, although varying recombination frequencies were observed between the two loci (w -M) no combination could severely reduce recombination (Additional file 3 Table S3).

Main findings
Induction of inversions was attempted through irradiation at 30 Gy. Analysis of recombination frequencies revealed at least one line that suppressed recombination between re and the M locus (line 35) and at least two lines that suppressed recombination between w and the M locus (lines 5 and 35). Searching for naturallyoccurring mutations in many different genomic backgrounds did not result in any combination that efficiently suppressed recombination either between re and M or w and M. Recombination in the re-M-w region: Different studies up to now have assessed recombination between these markers and all agree that the M locus resides between re and w (Additional file 4 Table S4). Recombination in the re-M region has been shown to vary between 1 and 12 cM, but usually in the more narrow range of 2.5-8.5 cM [45,47,[61][62][63][64][65][66]. Data from our group generated from recombination frequencies from consecutive generations using the same strains (BRA and Red-eye) suggest that recombination ranges between 1 and 2.5% [36]. At the same time, most studies suggest that recombination in the w-M  [44,[64][65][66][67][68] while our baseline recombination data collected from multiple generations using the same strains (BRA and Higgs Whiteeye) indicate that recombination ranges between 9 and 13% [36]. Variation in recombination frequencies have been attributed to factors such as the age, the sex, and the temperature, with the presence of widespread (small or extended) chromosomal rearrangements being the most probable explanation [45,47,61,62,68]. Following irradiation, most of the iso-male families screened exhibited recombination frequencies within the expected range. Line 35 significantly suppressed recombination both in the re -M and the w-M regions, whereas line 5 significantly suppressed recombination only in the w-M region, pointing to the presence of different inversions. Such low recombination frequencies (less than 0.2% for the re-M and less than 2% for the w-M region) have not been described before pointing to these lines as irradiation-induced inversions. Additional cytogenetic and/or genomic work can shed light to the genomic organization of chromosome I in these lines. It is encouraging that these two lines have been maintained without problems since they were isolated, suggesting that the chromosomal rearrangements involved do not have a severe negative effect on fitness.

Conclusions
SIT, IIT, and combined approaches, are receiving renewed interest for vector control (mosquitoes). Male recombination is restricting the development of GSS through classical genetics, since the desirable mutations must be very closely linked to the M locus. Inducing inversions through irradiation is a feasible strategy to isolate recombination suppressors either on the M or the m chromosome for Aedes species. Such inversions can be incorporated in genetic sexing strains to enhance their genetic stability and support vector control strategies that aim to population suppression through maledelivered sterility.

Strains description and rearing conditions
The Red Eye, Higgs White Eye (HWE), and a Brazilian wildtype strain (BRA) were used in the present study. The "wildtype" color of the eye in Ae. aegypti is dark brown/black and stable during all developmental stages. were inserted into the cages 48 h after the last blood feeding round in order to collect the mosquito eggs.

Irradiation
Wild type male pupae of the 'BRA' strain were irradiated in batches of 100 at 30 Gy, 30 to 36 h post pupation, using the method described as a stackable petri dish canister, in a Gammacell 220, self-shielded, gamma-ray Cobalt 60 irradiator [69]. Dosimetry was performed according to standard operating procedures regarding dosimetry systems for SIT [70] and all readings were within the 95% confidence intervals. The irradiated male pupae emerged in BugDorm-1 cages (30x30x30 cm), with access to 10% sucrose solution.

Crosses
Crosses are presented in Fig. 1 and followed the experimental set up described in the past [60] for the isolation of recombination suppressors on the M chromosome.

Isolation of recombination suppressors (re -M)
Parental cross was performed 'en masse'. Virgin homozygous re mutant females were transferred in the cages with the irradiated 'BRA' males, in a ratio of 3:1 (approximately 150 females for 50 males). F1 males were separated and screened in respect to the eye color at the pupal stage. All F1 progeny exhibited the wild type eye phenotype. One hundred F1 males were individually placed in BugDorm cages (BD4S1515: 15 cm × 15 cm × 15 cm) and backcrossed to five homozygous re mutant virgin females. F2 progeny were sorted based on sex and eye color at the pupal stage. Parental (males with wild type eyes and females with red eyes) and recombinant (males with red eye and females with black eye) genotypes were recorded. All families were kept through inbreeding after the removal of the recombinant genotypes, at the pupal stage (Additional file 5 Fig. S1). Promising families (showing evidence of reduced recombination) were upscaled and recombination frequencies were recorded in subsequent generations as described for F2 generation. χ 2 statistics were used to compare the observed recombination frequencies of the F2 and F3 of the promising families against the recombination frequency observed in the BRA genetic background (Additional file 2 Table S2). Calculations were performed using the Microsoft Excel 2016 formulas.

Isolation of recombination suppressors between (w -M)
Since resulting families represented iso-male lines, five F2 males from each family were transferred to new cages and crossed with ten homozygous w mutant virgin females. This is considered as the parental cross for this experiment. We acknowledge that genetic recombination may influence the genomic background of these males preventing them from being 'genetically identical'. However, the fact that these two mutations cannot be simultaneously tested due to the inconclusive phenotypes of double mutants led us to the compromise of sequentially testing them. F1 progeny exhibited the expected, wild type phenotype, and 5-10 F1 males were backcrossed with approximately 25 homozygous w mutant virgin females. F2 progeny were sorted based on sex and eye color at the pupal stage. Sex was also verified at the adult stage (Additional file 5 Fig. S1). Parental (males with wild type eyes and females with white eyes) and recombinant (males with white eye and females with black eye) genotypes were recorded. All families were kept through inbreeding after the removal of the recombinant genotypes at the pupal stage. Promising families (showing evidence of reduced recombination) were upscaled and recombination frequencies were recorded in F3. χ 2 statistics were used to compare the observed recombination frequencies of the F2 and F3 of the promising families against the recombination frequency observed in the BRA genetic background (Additional file 3 Table S3). Calculations were performed using the Microsoft Excel 2016 formulas.

Isolation of naturally-occurring recombination suppressors (re -M and w -M)
Since naturally-occurring mutations have been reported in different mosquito species, males from sixteen wild type colonized populations were crossed with virgin females of the 'Red Eye' and the 'Higgs White Eye' strains, respectively. Parental crosses were performed 'en masse' with 5-10 males and approximately 20-25 females. F1 progeny exhibited the expected, wild type phenotype, and at least ten F1 males were backcrossed with 25 virgin females of the respective homozygous mutant strain (either re or w). F2 progeny were sorted based on sex and eye color at the pupal stage. Parental (males with wild type eyes and females with either white or red eyes) and recombinant (males with either white or red eye and females with black eye) genotypes were recorded. All strains were kept through inbreeding after the removal of the recombinant genotypes at the pupal stage. Promising strains (showing evidence of reduced recombination) were upscaled and recombination frequencies were recorded in the following generation (F3).