A conserved female-specific larval requirement for MtnB function facilitates sex separation in multiple species of disease vector mosquitoes
Parasites & Vectors volume 14, Article number: 338 (2021)
Clusters of sex-specific loci are predicted to shape the boundaries of the M/m sex-determination locus of the dengue vector mosquito Aedes aegypti, but the identities of these genes are not known. Identification and characterization of these loci could promote a better understanding of mosquito sex chromosome evolution and lead to the elucidation of new strategies for male mosquito sex separation, a requirement for several emerging mosquito population control strategies that are dependent on the mass rearing and release of male mosquitoes. This investigation revealed that the methylthioribulose-1-phosphate dehydratase (MtnB) gene, which resides adjacent to the M/m locus and encodes an evolutionarily conserved component of the methionine salvage pathway, is required for survival of female larvae.
Larval consumption of Saccharomyces cerevisiae (yeast) strains engineered to express interfering RNA corresponding to MtnB resulted in target gene silencing and significant female death, yet had no impact on A. aegypti male survival or fitness. Integration of the yeast larvicides into mass culturing protocols permitted scaled production of fit adult male mosquitoes. Moreover, silencing MtnB orthologs in Aedes albopictus, Anopheles gambiae, and Culex quinquefasciatus revealed a conserved female-specific larval requirement for MtnB among different species of mosquitoes.
The results of this investigation, which may have important implications for the study of mosquito sex chromosome evolution, indicate that silencing MtnB can facilitate sex separation in multiple species of disease vector insects.
Blood-feeding vector mosquitoes transmit disease-causing pathogens that lead to hundreds of thousands of human deaths each year. For example, Zika, yellow fever, chikungunya, and dengue result from infections with arboviruses transmitted through the bites of Aedes mosquitoes, including Aedes aegypti and Aedes albopictus . Given poor progress in vaccine development and distribution, mosquito control is the primary mechanism for prevention of these and other mosquito-borne illnesses. Unfortunately, the emergence of insecticide resistance, concerns for the negative impacts of pesticides on the environment, and a lack of support for mosquito control programs compromise current strategies for managing mosquitoes .
Female mosquitoes possess sex-specific morphological, physiological, and behavioral traits, such as blood-feeding behavior, that promote the spread of disease-causing pathogens. Although genes that regulate female-specific mosquito development may represent novel targets for vector control, a majority of these genes have not yet been identified . Discovery and functional characterization of female-specific genes in multiple mosquito species could reveal mechanisms that regulate female mosquito development, provide insight into the evolution of these mechanisms, and promote the development of novel mosquito control interventions. For example, gene drive strategies that involve the manipulation of sex-specific genes are being developed and characterized [3,4,5]. Moreover, renewed interest in the sterile insect technique (SIT), which involves the release of sterile male insects to reduce insect pest populations , has recently emerged in the mosquito control community [7, 8]. Broad use of SIT for mosquito control has been hampered by the need to remove female mosquitoes, which is required to reduce health and nuisance biting risks posed by the release of female mosquitoes. Removal of females is also believed to increase the efficacy of this intervention, while also reducing costs [7, 8]. Likewise, the Wolbachia-based incompatible insect technique (IIT), which is proving useful for mosquito population suppression  and can be paired with SIT [10, 11], also requires mass rearing and release of adult male mosquitoes [7, 8]. Transgenic-based population suppression approaches involving the release of insects carrying a dominant lethal (RIDL) [12, 13], as well as various emerging gene drive technologies to block pathogen transmission or to suppress mosquito populations , also involve the release of male mosquitoes . Thus, the development and implementation of effective and affordable sex-separation technologies that can deployed worldwide are rate-limiting steps in the global introduction of numerous population-based mosquito control methods . The identification of genes with sex-specific functions could permit the elucidation of novel male sex-sorting technologies that could support these interventions.
Aedes aegypti sex determination is regulated by a non-recombining Y-chromosome-like male-determining region, the M locus, which is present on chromosome 1 [15, 16] and contains the male-determining factor Nix . Males, which possess one copy of the chromosome bearing the M locus and one which lacks it, have an M/m genotype, while A. aegypti females, which lack the male-determining locus, have an m/m genotype . Clusters of loci that cause sex-specific lethal effects are believed to reside in the vicinity of the sex-determining region [19, 20]. However, the identities of these genes are not known. Characterization and sequencing of the M/m locus region had been thwarted by the repetitive nature of DNA located in this centromeric region of chromosome 1. However, recent innovations in sequencing technology generated an improved and re-annotated genome assembly that facilitated better estimation of the M/m locus, providing opportunities to study the evolution of the sex-determining chromosomes and prospects for evaluating the functional contributions of genes linked to the M/m locus . In this investigation, methylthioribulose-1-phosphate dehydratase (MtnB), which resides adjacent to the M/m locus, was identified as a female-specific larval lethal locus.
MtnB, an evolutionarily conserved component of the methionine salvage pathway, functions to convert 5-methylthioadenosine (MTA) to methionine. This process occurs through six enzymatic reactions, the third of which, the dehydration of 5-methylthioribulose-1-phosphate (MTRu-1-P) to 2,3-diketo-5-methylthiopentyl-1-phosphate (DK-MTP-1-P), is catalyzed by MtnB [22,23,24,25]. This pathway facilitates recycling of MTA, a by-product of polyamine synthesis, back to methionine, an essential amino acid, while also serving as a useful means of recycling the sulfur present in MTA . The methionine salvage pathway has also been implicated in cell death and inflammation [26,27,28,29]. In this investigation, RNAi-based silencing of MtnB in A. aegypti revealed a female-specific requirement for this gene in larvae that was exploited for the development of a scalable yeast-based oral RNAi method for production of fit adult males. Characterization of MtnB in several additional species of disease vector mosquitoes, including Aedes albopictus (vector of multiple arboviruses), Anopheles gambiae (malaria parasite vector), and Culex quinquefasciatus (lymphatic filariasis and West Nile viral vector) , suggested that targeting MtnB during larval development could facilitate male sex sorting in multiple species of disease vector mosquitoes.
The A. aegypti LVP-IB12, A. albopictus Gainesville, A. gambiae G3, and C. quinquefasciatus JHB (BEI Resources, NIAID, NIH) strains were reared as described , except that sheep blood (purchased from HemoStat Laboratories, Dixon, CA, USA) was provided with an artificial membrane feeding system (Hemotek Limited, Blackburn, UK).
Custom small interfering RNAs (siRNAs) corresponding to mosquito MtnB genes (Table 1) were selected using the Integrated DNA Technologies (IDT) Custom Dicer-Substrate siRNA (DsiRNA) tool . Custom siRNAs (see Table 1 for further information), as well as a control siRNA lacking a known target in mosquitoes , were purchased from IDT (Coralville, IA, USA) for use in larval soaking screening experiments that were performed in duplicate using the method of Singh et al. . For each soaking experiment, 20 first instar larvae (L1) were soaked in 0.5 µg/µl control or experimental siRNA (Table 1) for 4 h and then reared as described . Soaking data were statistically evaluated using the chi-squared test to identify statistically significant differences between the observed and expected survival of female and male mosquitoes, assuming expected 1 male:1 female ratios on the basis of laboratory observations of the mosquitostrains when reared in the laboratory insectary.
Production of yeast interfering RNA larvicides
Yeast strains were prepared as described  through cloning of custom shRNA expression cassettes corresponding to target sequences of interest (Table 1; synthesized by Invitrogen, Carlsbad, CA, USA) into the pRS426 GPD shuttle vector . After confirming the clones through restriction digestion and sequencing, the resulting plasmids were used to transform S. cerevisiae CEN.PK yeast (genotype MATa/α ura3-52/ura3-52 trp1-289/trp1-289 leu2-3_112/leu2-3_112 his3 Δ1/his3 Δ1 MAL2-8C/MAL2-8C SUC2/SUC2 . Transformants were selected based on growth on minimal media lacking uracil. shRNA expression was confirmed through PCR conducted as described  using the following primers: forward 5′-CAGGATATTCTTTCATAGCATGTTC-3′ (specific to the 3′ end of the MtnB hairpin) and reverse primer 5′-TCCTTCCTTTTCGGTTAGAGC-3′ (which corresponds to the terminator sequence). An image of the PCR products visualized with ethidium bromide staining (Additional file 1: Fig. S1) was edited using Adobe Photoshop 2021 software (to remove non-relevant lanes and to label the figure). Following strain confirmation, dried inactivated yeast interfering RNA was prepared for use in larvicide assays as described .
Yeast interfering RNA larvicide screening
Yeast interfering RNA larvicides were evaluated through larvicide assays that conformed to the World Health Organization (WHO)  guidelines, which were performed as previously described . In summary, 20 L1 larvae were placed in 500 ml plastic cups containing 50 ml of distilled water and 50 mg of dried inactivated yeast (either experimental larvicidal yeast or yeast prepared from a control interfering RNA strain ). The emergence rates and sexes of the resulting adults were recorded in each of nine replicate trials that were performed for each control and larvicidal yeast treatment. Data were evaluated using the chi-squared test to identify statistically significant differences between the observed and expected 1 male:1 female ratios, female survival, and male survival for each treatment.
MtnB silencing was confirmed through qRT-PCR assays that were performed as described . In summary, TRIzol Reagent (Invitrogen, Carlsbad, CA, USA) was used according to the manufacturer’s instructions to prepare total RNA from 20 pooled third instar larvae that had fed on either control or MtnB.496 yeast throughout larval development. cDNA was prepared from total RNA extracts using the High Capacity RNA to cDNA Kit (Applied Biosystems, Foster City, CA, USA) according to the manufacturer’s instructions. Real-time PCR quantification assays were performed using an Applied Biosystems Step One Plus Real-Time PCR System in conjunction with the Power SYBR Green PCR Master Mix as described by the manufacturer (Applied Biosystems, Foster City, CA, USA). The following primer sets were used in these reactions: forward 5′-GGC AAT GGC TGA ATG TTA CG-3′ and reverse: 5′-GTG ATT GGA ATC AGA ATT GAC TTA-3′. Amplification of rpS17, an internal standard for data normalization, was performed as described  with the following primers: forward 5′ AGA CAA CTA CGT GCC GGA AG 3′ and reverse 5′ TTG GTG ACC TGG ACA ACG ATG 3′. PCR amplifications were performed in six replicate wells in each of two separate biological replicate trials, and results were quantified by standardizing reactions to rpS17 transcript levels using the ΔΔCt method as described . Data were statistically evaluated with Student’s t-test.
MtnB developmental transcript expression was quantified in two separate pools of ten untreated individuals for L1 larvae, L2 larvae, L3 larvae, L4 larvae, male pupae, and female pupae (note that the pupal time point is the first stage in which A. aegypti sexes are visibly discernable). Total RNA from the separate pools was prepared and analyzed as described above. Results were quantified by standardizing reactions to rpS17 transcript levels using the ΔΔCt method as described . Setting the expression level of L4 larvae to one permitted standardization and comparison of qRT-PCR results across the different stages. Transcript expression data were analyzed using analysis of variance (ANOVA).
Male life history trait assessments
Life history traits were evaluated as described by Hill et al. . In these assays, males that had been reared on control interfering RNA or MtnB.496 yeast (n = 52 combined from two separate biological replicate trials) were mated as individuals to wild-type females, which had been reared as described . Mating success was measured by confirming that females had laid eggs; if females failed to lay eggs, mating success/failure was further assessed through dissection of the spermatheca to confirm the presence or absence of sperm. Mating data were analyzed with the chi-squared test. For successful matings (n = 51 females mated to control males and n = 49 females mated to Mtn.496-treated males), the number of eggs laid per individual female (fecundity) and number of larvae that hatched (fertility) were recorded. Fertility and fecundity data were analyzed with Student’s t-test.
Two hundred L1 LVP-IB12 larvae were placed in 34 × 25 × 4 cm trays (1426B, Bioquip, Rancho Dominguez, CA, USA) bearing 1 L of distilled water. Per Zhang et al. , larvae were cultured using a mass-rearing diet that consisted of a slurry of 250 mg bovine liver powder (MP Biomedicals, Santa Ana, CA, USA) and 150 mg of ground shrimp (Tetra GMBH, Melle, Germany) mixed with 10 ml water, but with the brewer’s yeast component of the diet  replaced with yeast interfering RNA. The yeast interfering RNA component (control or Mtn.496) was combined with the liver powder and shrimp slurry and fed to larvae as follows: 100 mg at L1, 100 mg at L2, 100 mg at L3, and 200 mg at L4. Each day, larval trays were examined for pupae, which were removed, sorted by sex, and counted. Mortality and sex ratio distortions were evaluated in three biological replicate trials, each with two replicate trays of larvae fed using diets prepared with control interfering RNA or MtnB.496 yeast. Statistically significant differences between the observed and 1:1 expected survival levels for males and females were evaluated using the chi-squared test. The fitness of surviving males was ascertained through comparisons of wing lengths, which were measured as described . Wing length data were evaluated with Student’s t-test.
Silencing MtnB during A. aegypti larval development selectively kills female mosquitoes
Given that the M/m locus is believed to be tightly linked to developmental genes that confer sex-specific effects in A. aegypti, MtnB, which is located on chromosome 1 and flanks the M/m locus [21, 44], was hypothesized to function as a potential sex-specific lethal gene. MtnB was silenced through larval soaking with siRNA #496, which corresponds to the MtnB transcript. Treatment with siRNA #496 resulted in significant female mortality yet had no significant impact on male survival (further details, including P values, are provided in Table 1). Given the results of this brief siRNA soaking-mediated transient gene silencing assay, S. cerevisiae, which was recently developed as a system for expression and sustained delivery of interfering RNA to mosquito larvae, was used to express a short hairpin RNA (shRNA) corresponding to siRNA #496. shRNA expression in this yeast strain, hereafter referred to as MtnB.496, was confirmed through PCR amplification of cDNA corresponding to the shRNA transcript (Additional file 1: Fig. S1). Dried inactivated MtnB.496 yeast was fed to larvae throughout larval development. Larvae were also reared on yeast that expresses an shRNA with no known target in mosquitoes  for comparison. MtnB transcripts, which are detected during larval and pupal development (Fig. 1a), are silenced in larvae that consume MtnB.496 yeast (Fig. 1b; 88 ± 4% reduction in transcript levels in comparison to larvae reared on control interfering RNA yeast; t-test, P < 0.001). Larval consumption of MtnB.496 yeast, but not control interfering RNA yeast, resulted in significantly higher male:female ratios among the surviving mosquitoes (Fig. 1c; χ2 = 31.387, P < 0.0001). Although no significant impact on male survival was observed following consumption of MtnB.496 yeast (χ2 = 3.214, P > 0.05), only 13 ± 2% of expected adult female mosquitoes survived, with a majority of animals dying as late larvae or early pupae (χ2 = 128.571, P < 0.0001).
Scaled production of fit adult A. aegypti male mosquitoes
To confirm that the impacts of MtnB.496 larvicide are specific to female larvae, life history traits were assessed in adult male mosquitoes that had consumed the yeast throughout larval development. The percentage of treated males that mated with wild-type females was not significantly different than adult males that were reared on control yeast interfering RNA larvicides (Fig. 2a, χ2 = 0.230, P > 0.05). When the MtnB.496-treated males mated with wild-type female mosquitoes, no significant differences in female fertility (Fig. 2b, t = −1.925, P > 0.05) or fecundity were observed (Fig. 2b, t = 0.698, P > 0.05). These findings indicated that MtnB.496 yeast larvicide treatment did not produce any detectable effects on male fitness in the assays conducted.
To assess whether the RNAi-based yeast larvicides could also facilitate scaled production of males, a larval diet and feeding regimen employed at Aedes mass-rearing facilities  was modified through replacement of the nutritional yeast component of the diet with dried inactivated yeast larvicide MtnB.496. The modified diet was tested on A. aegypti larvae grown in mass-rearing trays. Under these conditions, the larvicides continued to induce significant A. aegypti female mortality (Fig. 2c, χ2 = 1275.764, P < 0.0001) with no impacts on male survival (Fig. 2c, χ2 = 0.5777, P > 0.05). The fitness of male survivors, which was estimated through measurement of wing lengths, a correlate of body size and fitness, was not significantly different than that of males reared on the standard mass-rearing diet (Fig. 2d, t = 1.031, P > 0.05).
MtnB orthologs are female-specific lethal genes in multiple species of disease vector mosquitoes
MtnB orthologs have been identified in A. gambiae and C. quinquefasciatus, and BLAST searches identified an ortholog in A. albopictus (Table 1). RNAi-mediated targeting of these genes through brief soaking treatments in siRNAs #529 (corresponds to A. albopictus MtnB), #523 (corresponds to A. gambiae MtnB), and #534 (corresponds to C. quinquefasciatus MtnB) resulted in significant female-specific larval death, yet did not impact male survival through adulthood (Table 1, statistics provided in table). Yeast interfering RNA larvicide strains that express shRNAs corresponding to siRNAs #529, #523, and #534 were therefore constructed and are hereafter referred to as strains MtnB.529 (targets A. albopictus), MtnB.523 (targets A. gambiae), and MtnB.534 (targets C. quinquefasciatus).
Larval consumption of MtnB.529 by A. albopictus larvae resulted in significantly higher than expected male:female ratios (Fig. 3a; χ2 = 24.202, P < 0.0001). Only 24 ± 1% of expected A. albopictus female mosquitoes survived (Fig. 3a, χ2 = 101.829, P < 0.00001), but no significant impact on male survival was observed (Fig. 3a, χ2 = 2.067, P > 0.05). Likewise, MtnB.523 resulted in significantly higher male:female ratios among surviving A. gambiae adults (Fig. 3b; χ2 = 7.453; P < 0.01). Although MtnB.523 had no significant impact on the survival of males through adulthood (χ2 = 0.004; P > 0.05), only 43 ± 2% of expected A. gambiae female mosquitoes survived (Fig. 3b, χ2 = 37.551, P < 0.0001). Finally, significantly higher than expected male:female ratios were also observed among C. quinquefasciatus adult survivors following larval consumption of MtnB.534 (Fig. 3c; χ2 = 28.751, P < 0.0001), a difference that coincided with only 21 ± 1% of expected C. quinquefasciatus adult female emergence (Fig. 3c; χ2 = 95.890; P < 0.00001), with no significant corresponding changes in male survival (Fig. 3c; χ2 = 3.051, P > 0.05).
Implications for understanding the evolution of sex chromosomes in mosquitoes
Sex chromosome evolution is believed to occur in several stages [45,46,47], starting with a homologous pair of autosomes that acquire sex-determining loci. A proto-Y chromosome, which has a male fertility locus (M) and a dominant female suppressor locus (SuF), as well as a proto-X chromosome carrying a female fertility locus (F) and a male sterility locus (m) form. Suppressed recombination in the sex-determining region evolves, eventually spreading along the proto-sex chromosomes, which accumulate transposable elements and non-coding sequences, undergo chromosomal rearrangements, and accumulate sexually antagonistic genes with different alleles that differentially benefit either males or females. The proto-sex chromosomes eventually evolve into heteromorphic X and Y chromosomes, with the Y chromosome eventually degenerating and reducing in size [45,46,47].
Sex chromosomes in different mosquito species have progressed to varying stages of this evolution. Culicines, including Aedes and Culex mosquitoes, possess homomorphic sex chromosomes in which the sex-determination alleles are linked to the smallest homomorphic first chromosome . In A. aegypti, the homomorphic sex chromosomes have evolved into proto-sex chromosomes bearing a sex determining M/m locus [15, 16]. The M locus at the proto-Y chromosome contains the male-determining factor, Nix , which directs male-specific splicing of the sex chromosome gene doublesex (dsx). In females, which have two copies of the proto-X (m) chromosome, the female-specific DsxF form regulates ovary development and fertility . Although a sex-differentiated region of suppressed recombination extends ~ 100 Mb beyond the M/m locus [21, 49, 50], heteromorphic X and Y chromosomes have not yet evolved in A. aegypti. Likewise, A. albopictus also possesses homomorphic sex chromosomes, and the A. albopictus Nix gene was identified as a strong candidate for a male-determining factor in this species [17, 51, 52]. In contrast, A. gambiae has already evolved separate X and Y chromosomes, with the Yob gene functioning as a male-determining factor on the Y chromosome . Despite different mosquito species having advanced to various stages of sex chromosome evolution, this investigation functionally verified that MtnB is required for survival of female larvae in A. aegypti (Figs. 1c, 2c), A. albopictus (Fig. 3a), A. gambiae (Fig. 3b), and C. quinquefasciatus (Fig. 3c) mosquitoes. These findings suggest that loci which cause sex-specific lethal effects and shape the boundaries of the sex-determination locus were accumulating on the sex chromosomes prior to separation of the culicine and anopheline mosquitoes. Future studies will be directed toward understanding the nature of this female-specific requirement during larval development.
Krzywinska et al.  discuss possible mechanisms underlying the sex-specific lethality they observed in A. aegypti that had inherited sex chromosomes produced through male meiotic recombination events that occurred between the sex-determination locus and a linked enhanced green fluorescent protein (EGFP)-marked transgene. An explanation which is consistent with the MtnB silencing results presented here proposed that such recombination events could lead to the loss of a haploinsufficient sex chromosome gene for which two copies are required in m/m females. If MtnB is haploinsufficient in females, any female that inherits a recombinant m chromosome lacking a copy of the MtnB gene will die, an event that is sometimes referred to as “suppression of recombination” because resulting female offspring are not observed in adult populations. Likewise, as observed in this investigation, silencing MtnB using RNAi-based yeast larvicides can also result in female death. These silencing results suggest that the presence of haploinsufficient genes adjacent to the sex-determination locus may have helped to shape the boundaries of the sex-determination region during mosquito sex chromosome evolution, as previously suggested .
Papathanos et al.  recommend that sex-sorting methodology should be simple, stable, and produce competitive males, ideally through female elimination during development. The yeast female-specific larvicide system described herein meets these criteria (Fig. 2). It is also recommended that the system be cost-effective, scalable, and suitable for implementation in remote or resource-limited regions of the world , criteria that would be met through the use of yeast, which is routinely produced at commercial scale by the food industry for worldwide distribution. Moreover, the demonstrated ability to seamlessly incorporate yeast into a standard mass-rearing diet (Fig. 2c) makes this system simple to implement. Sex-sorting technology should meet regulatory standards, ideally requiring no modification of existing regulatory protocols that may have already been acquired for the mosquito technology under development. It is also important that the technology for elimination of females be conditional, so as not to interfere with strain maintenance . The yeast sex-sorting system meets these criteria, as it induces no permanent changes to the mosquito genome, but acts through conditional gene silencing during female development. The use of yeast RNAi-based larvicides would therefore circumvent the need to further genetically manipulate existing mosquito strains developed for population control strategies. Furthermore, the yeast would be used only during the larval stages and would be heat-killed prior to use indoors at mass-rearing facilities. This may mitigate a need for modification of existing regulatory permits previously obtained for male mosquito releases, a prospect that can be further discussed with the appropriate regulatory bodies. Thus, sex-specific yeast interfering RNA technology meets many of the criteria for a useful sex-separation system . Future studies will be directed toward scaled industrial-sized production of yeast interfering RNA larvicides for distribution of yeast to mass-rearing mosquito facilities worldwide.
In this investigation, replacement of the nutritional yeast component of a larval mass-rearing diet  resulted in female-specific larval killing, facilitating production of fit male mosquitoes at a ratio of 5 male:1 female mosquitoes (Fig. 2d). Silencing MtnB did not eliminate all females (Fig. 2b, c) and this yeast strain could not likely be used as a stand-alone strategy for male mosquito production. It would be interesting to explore methods for further improving the female-killing capacities of the yeast strains (i.e. increased shRNA expression levels, stability, or uptake). We hope to explore this in future studies, as a stand-alone strategy that generates > 99% male separation, as observed through use of a sophisticated camera system recently developed for A. aegypti , would be preferable. Although the automated camera system is highly effective, it is not yet clear if such technology can be deployed worldwide or adapted for use is all mosquito species. In such instances, and particularly when resources are limited, replacing nutritional yeast with larvicidal yeast during the mass-rearing process could be used in combination with existing sex-separation technologies, particularly manual separation strategies, which are highly labor-intensive .
Female-specific yeast larvicide technology could particularly benefit the development and implementation of population-based control strategies for C. quinquefasciatus (Fig. 3c), in which sex-sorting methods other than physical separation have yet to be established [7, 8]. The female-specific mortality levels observed in the current investigation were relatively lower in A. gambiae with respect to the other mosquitoes evaluated (Fig. 3). However, given the conservation of MtnB in many species of Anopheles mosquitoes , it may nevertheless still be worth exploring the potential for using MtnB silencing to sex other Anopheles species. For example, the importance of developing sex-separation strategies that could be implemented in malaria vectors such as Anopheles arabiensis  and Anopheles stephensi  has been discussed , and MtnB is conserved in both species, as well as several other anophelines. It may therefore be helpful to pair a yeast larvicide diet with the manual or visual sorting strategies currently employed in these Anopheles mosquitoes . One could also envision pairing the RNAi yeast diet sex-separation approach with various transgenic sex-sorting strains that have been generated for several species, including female lethal strains or transgenic fluorescent marker strains [57, 58] that facilitate automated sorting. Likewise, yeast feeding could be combined with imaging-based sorting strategies  and automated pupal size estimators , or paired with feeding of toxicant-infused blood meals to females . Such combined approaches may better permit the > 99% removal of females that is generally accepted as a minimal requirement for male-only release control strategies .
These studies have uncovered a conserved female-specific larval requirement for MtnB, which encodes a component of the methionine salvage pathway, in distantly related mosquitoes. The results of this investigation may have important implications for the study of mosquito sex chromosome evolution, and it will be interesting to further investigate this topic in future studies. It is anticipated that female-specific yeast interfering RNA technologies that capitalize on the conserved female-specific requirement for MtnB function in mosquito larvae could permit mass production of fit males, thereby facilitating emerging population-based strategies for control of multiple disease vector mosquito species.
Availability of data and materials
All data generated or analyzed during this investigation are included in this manuscript. Yeast strains and corresponding plasmids generated during the course of this study will be made available upon satisfactory completion of a material transfer agreement with Indiana University and following procurement of required import permits by the requesting party.
- Dsx :
First instar larvae
Integrated DNA Technologies
Short hairpin RNA
Small interfering RNA
Centers for Disease Control CDC. Mosquitoes and diseases: A-Z. https://www.cdc.gov/mosquitoes/about/diseases.html. Accessed Jan 2020.
Duman-Scheel M, Syed Z. Developmental neurogenetics of sexual dimorphism in Aedes aegypti. Front Ecol Evol. 2015;3:61.
Kyrou K, Hammond AM, Galizi R, Kranjc N, Burt A, Beaghton AK, et al. A CRISPR-Cas9 gene drive targeting doublesex causes complete population suppression in caged Anopheles gambiae mosquitoes. Nat Biotechnol. 2018;36(11):1062–6.
Simoni A, Hammond AM, Beaghton AK, Galizi R, Taxiarchi C, Kyrou K, et al. A male-biased sex-distorter gene drive for the human malaria vector Anopheles gambiae. Nat Biotechnol. 2020;38(9):1054–60.
Adelman ZN, Tu Z. Control of mosquito-borne infectious diseases: sex and gene drive. Trends Parasitol. 2016;32(3):219–29.
Knipling EF. Possibilities of insect control or eradication through the use of sexually sterile males. J Econ Entomol. 1955;48:459–62.
Lutrat C, Giesbrecht D, Marois E, Whyard S, Baldet T, Bouyer J. Sex sorting for pest control: It’s raining men! Trends Parasitol. 2019;35(8):649–62.
Papathanos PA, Bourtzis K, Tripet F, Bossin H, Virginio JF, Capurro ML, et al. A perspective on the need and current status of efficient sex separation methods for mosquito genetic control. Parasit Vectors. 2018;11(Suppl 2):654.
Crawford JE, Clarke DW, Criswell V, Desnoyer M, Cornel D, Deegan B, et al. Efficient production of male Wolbachia-infected Aedes aegypti mosquitoes enables large-scale suppression of wild populations. Nat Biotechnol. 2020;38(4):482–92.
Zheng X, Zhang D, Li Y, Yang C, Wu Y, Liang X, et al. Incompatible and sterile insect techniques combined eliminate mosquitoes. Nature. 2019;572(7767):56–61.
Kittayapong P, Ninphanomchai S, Limohpasmanee W, Chansang C, Chansang U, Mongkalangoon P. Combined sterile insect technique and incompatible insect technique: The first proof-of-concept to suppress Aedes aegypti vector populations in semi-rural settings in Thailand. PLoS Negl Trop Dis. 2019;13(10):e0007771.
Carvalho DO, McKemey AR, Garziera L, Lacroix R, Donnelly CA, Alphey L, et al. Suppression of a field population of Aedes aegypti in Brazil by sustained release of transgenic male mosquitoes. PLoS Negl Trop Dis. 2015;9(7):e0003864.
Wise de Valdez MR, Nimmo D, Betz J, Gong HF, James AA, Alphey L, et al. Genetic elimination of dengue vector mosquitoes. Proc Natl Acad Sci USA. 2011;108(12):4772–5.
Macias VM, Ohm JR, Rasgon JL. Gene drive for mosquito control: where did it come from and where are we headed? Int J Environ Res Public Health. 2017;14(9):1006.
McClelland GAH. Sex-linkage in Aedes aegypti. Trans R Soc Trop Med Hyg. 1962;56(4).
Newton ME, Wood RJ, Southern DI. Cytological mapping of the M and D loci in the mosquito, Aedes aegypti (L.). Genetica. 1978;48(2):137–43.
Hall AB, Basu S, Jiang X, Qi Y, Timoshevskiy VA, Biedler JK, et al. Sex determination. A male-determining factor in the mosquito Aedes aegypti. Science. 2015;348(6240):1268–70.
Motara MA, Rai KS. Giemsa C-banding patterns in Aedes (Stegomyia) mosquitoes. Chromosoma. 1978;70(1):51–8.
Wood RJ. Lethal genes on the sex chromosomes concealed in a population of the mosquito Aedes aegypti L. Genetica. 1990;46(1):49–66.
Krzywinska E, Kokoza V, Morris M, de la Casa-Esperon E, Raikhel AS, Krzywinski J. The sex locus is tightly linked to factors conferring sex-specific lethal effects in the mosquito Aedes aegypti. Heredity. 2016;117(6):408–16.
Matthews BJ, Dudchenko O, Kingan SB, Koren S, Antoshechkin I, Crawford JE, et al. Improved reference genome of Aedes aegypti informs arbovirus vector control. Nature. 2018;563(7732):501–7.
Sekowska A, Denervaud V, Ashida H, Michoud K, Haas D, Yokota A, et al. Bacterial variations on the methionine salvage pathway. BMC Microbiol. 2004;4:9.
Ashida H, Saito Y, Kojima C, Yokota A. Enzymatic characterization of 5-methylthioribulose-1-phosphate dehydratase of the methionine salvage pathway in Bacillus subtilis. Biosci Biotechnol Biochem. 2008;72(4):959–67.
Pirkov I, Norbeck J, Gustafsson L, Albers E. A complete inventory of all enzymes in the eukaryotic methionine salvage pathway. FEBS J. 2008;275(16):4111–20.
Mary C, Duek P, Salleron L, Tienz P, Bumann D, Bairoch A, et al. Functional identification of APIP as human mtnB, a key enzyme in the methionine salvage pathway. PLoS One. 2012;7(12):e52877.
Albers E. Metabolic characteristics and importance of the universal methionine salvage pathway recycling methionine from 5’-methylthioadenosine. IUBMB Life. 2009;61(12):1132–42.
Avila MA, Garcia-Trevijano ER, Lu SC, Corrales FJ, Mato JM. Methylthioadenosine. Int J Biochem Cell Biol. 2004;36(11):2125–30.
Tang B, Kadariya Y, Murphy ME, Kruger WD. The methionine salvage pathway compound 4-methylthio-2-oxobutanate causes apoptosis independent of down-regulation of ornithine decarboxylase. Biochem Pharmacol. 2006;72(7):806–15.
Li C, Li Z, Sletten E, Arnesano F, Losacco M, Natile G, et al. Methionine can favor DNA platination by trans-coordinated platinum antitumor drugs. Angew Chem Int Ed Engl. 2009;48(45):8497–500.
Clemons A, Mori A, Haugen M, Severson DW, Duman-Scheel M. Culturing and egg collection of Aedes aegypti. Cold Spring Harb Protoc. 2010;2010(10):pdbprot5507.
IDT. Custom Dicer-Substrate siRNA (DsiRNA). https://www.idtdna.com/site/order/designtool/index/DSIRNA_CUSTOM (2017). Accessed Apr 2017.
Tomchaney M, Mysore K, Sun L, Li P, Emrich SJ, Severson DW, et al. Examination of the genetic basis for sexual dimorphism in the Aedes aegypti (dengue vector mosquito) pupal brain. Biol Sex Differ. 2014;5:10.
Singh AD, Wong S, Ryan CP, Whyard S. Oral delivery of double-stranded RNA in larvae of the yellow fever mosquito, Aedes aegypti: implications for pest mosquito control. J Insect Sci. 2013;13:69.
Hapairai LK, Mysore K, Chen Y, Harper EI, Scheel MP, Lesnik AM, et al. Lure-and-kill yeast interfering RNA larvicides targeting neural genes in the human disease \vector mosquito Aedes aegypti. Sci Rep. 2017;7(1):13223.
Mysore K, Hapairai LK, Wei N, Realey JS, Scheel ND, Severson DW, et al. Preparation and use of a yeast shRNA delivery system for gene silencing in mosquito larvae. Methods Mol Biol. 2019;1858:213–31.
Mumberg D, Muller R, Funk M. Yeast vectors for the controlled expression of heterologous proteins in different genetic backgrounds. Gene. 1995;156(1):119–22.
van Dijken JP, Bauer J, Brambilla L, Duboc P, Francois JM, Gancedo C, et al. An interlaboratory comparison of physiological and genetic properties of four Saccharomyces cerevisiae strains. Enzyme Microb Technol. 2000;26(9–10):706–14.
Mysore K, Hapairai LK, Sun L, Li P, Wang CW, Scheel ND, et al. Characterization of a dual-action adulticidal and larvicidal interfering RNA pesticide targeting the Shaker gene of multiple disease vector mosquitoes. PLoS Negl Trop Dis. 2020;14(7):e0008479.
WHO. Guidelines for laboratory and field testing of mosquito larvicides. Geneva: World Health Organization; 2005.
Morlais I, Mori A, Schneider JR, Severson DW. A targeted approach to the identification of candidate genes determining susceptibility to Plasmodium gallinaceum in Aedes aegypti. Mol Genet Genomics. 2003;269(6):753–64.
Hill CL, Sharma A, Shouche Y, Severson DW. Dynamics of midgut microflora and dengue virus impact on life history traits in Aedes aegypti. Acta Trop. 2014;140:151–7.
Zhang D, Zhang M, Wu Y, Gilles JRL, Yamada H, Wu Z, et al. Establishment of a medium-scale mosquito facility: optimization of the larval mass-rearing unit for Aedes albopictus (Diptera: Culicidae). Parasit Vectors. 2017;10(1):569.
Mysore K, Sun L, Tomchaney M, Sullivan G, Adams H, Piscoya AS, et al. siRNA-mediated silencing of doublesex during female development of the dengue vector mosquito Aedes aegypti. PLoS Negl Trop Dis. 2015;9(11):e0004213.
Giraldo-Calderon GI, Emrich SJ, MacCallum RM, Maslen G, Dialynas E, Topalis P, et al. VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. Nucleic Acids Res. 2015;43:D707-13.
Charlesworth B. Model for evolution of Y chromosomes and dosage compensation. Proc Natl Acad Sci USA. 1978;75(11):5618–22.
Charlesworth B. The evolution of sex chromosomes. Science. 1991;251(4997):1030–3.
Charlesworth D, Charlesworth B, Marais G. Steps in the evolution of heteromorphic sex chromosomes. Heredity. 2005;95(2):118–28.
Sinkins SP, Gould F. Gene drive systems for insect disease vectors. Nat Rev Genet. 2006;7(6):427–35.
Juneja P, Osei-Poku J, Ho YS, Ariani CV, Palmer WJ, Pain A, et al. Assembly of the genome of the disease vector Aedes aegypti onto a genetic linkage map allows mapping of genes affecting disease transmission. PLoS Negl Trop Dis. 2014;8(1):e2652.
Fontaine A, Filipovic I, Fansiri T, Hoffmann AA, Cheng C, Kirkpatrick M, et al. Extensive genetic differentiation between homomorphic sex chromosomes in the mosquito vector, Aedes aegypti. Genome Biol Evol. 2017;9(9):2322–35.
Miller JR, Koren S, Dilley KA, Puri V, Brown DM, Harkins DM, et al. Analysis of the Aedes albopictus C6/36 genome provides insight into cell line utility for viral propagation. Gigascience. 2018;7(3):1–13.
Gomulski LM, Mariconti M, Di Cosimo A, Scolari F, Manni M, Savini G, et al. The Nix locus on the male-specific homologue of chromosome 1 in Aedes albopictus is a strong candidate for a male-determining factor. Parasit Vectors. 2018;11(Suppl 2):647.
Krzywinska E, Krzywinski J. Effects of stable ectopic expression of the primary sex determination gene Yob in the mosquito Anopheles gambiae. Parasit Vectors. 2018;11(Suppl 2):648.
Papathanos PA, Bossin HC, Benedict MQ, Catteruccia F, Malcolm CA, Alphey L, et al. Sex separation strategies: past experience and new approaches. Malar J. 2009;8(Suppl 2):S5.
Mashatola T, Ndo C, Koekemoer LL, Dandalo LC, Wood OR, Malakoane L, et al. A review on the progress of sex-separation techniques for sterile insect technique applications against Anopheles arabiensis. Parasit Vectors. 2018;11(Suppl 2):646.
Wu Y, Hu W, Biedler JK, Chen XG, Tu ZJ. Pure early zygotic genes in the Asian malaria mosquito Anopheles stephensi. Parasit Vectors. 2018;11(Suppl 2):652.
Catteruccia F, Benton JP, Crisanti A. An Anopheles transgenic sexing strain for vector control. Nat Biotechnol. 2005;23(11):1414–7.
Magnusson K, Mendes AM, Windbichler N, Papathanos PA, Nolan T, Dottorini T, et al. Transcription regulation of sex-biased genes during ontogeny in the malaria vector Anopheles gambiae. PLoS ONE. 2011;6(6):e21572.
Zacares M, Salvador-Herranz G, Almenar D, Tur C, Argiles R, Bourtzis K, et al. Exploring the potential of computer vision analysis of pupae size dimorphism for adaptive sex sorting systems of various vector mosquito species. Parasit Vectors. 2018;11(Suppl 2):656.
Mukiama TK. Y-autosome genetic sexing strain of Anopheles albimanus (Diptera: Culicidae). Int J Trop Insect Sci. 1985;6:649–52.
We thank Na Wei and Chien-Wei Wang for their suggestions on yeast strain production and Zhiyong Xi for his recommendations on preparation of the larval mass-rearing diet. We thank members of the lab for useful discussions.
This work was funded by NIH-NIAID Awards 1 R21 AI144256-01 and R21AI156170-01 to MDS. The funders were not involved in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics approval and consent to participate
Not applicable. No vertebrate animals or human subjects were used in this research investigation.
Consent for publication
MDS was named as the inventor on U.S. patent application number 62/751,052. The pending application did not impact her data interpretation or decision to publish this work. All other authors declare that they lack competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Confirmation of shRNA expression in recombinant yeast strain MtnB.496. PCR reactions performed with primers corresponding to the MtnB.496 shRNA transcript generated a ~ 100 bp amplicon (see DNA standard at far left; cDNA template was prepared from MtnB.496 yeast total RNA). Negative control PCR reactions included an amplification with cDNA prepared from non-transformed yeast (marked by – sign) and a reaction with no cDNA added (marked by knot symbol). A representative ethidium bromide-stained agarose gel from one of two comparable biological replicate experiments is shown.
About this article
Cite this article
Mysore, K., Sun, L., Roethele, J.B. et al. A conserved female-specific larval requirement for MtnB function facilitates sex separation in multiple species of disease vector mosquitoes. Parasites Vectors 14, 338 (2021). https://doi.org/10.1186/s13071-021-04844-w
- Aedes aegypti
- Aedes albopictus
- Anopheles gambiae
- Culex quinquefasciatus
- Saccharomyces cerevisiae