Extensive Genetic Differentiation between Homomorphic Sex Chromosomes in the Mosquito Vector, Aedes aegypti

Differential elimination of marked sex chromosomes enables production of nontransgenic male mosquitoes in a single strain

Diverse genetic strategies are being pursued to control mosquito-borne infectious diseases. These strategies often rely on the release of nonbiting males to either reduce the target mosquito population or render them resistant to pathogens. Male-only releases are important as any contaminating females can bite and potentially transmit pathogens. Despite significant efforts, it remains a major bottleneck to reliably and efficiently separate males from females, especially when nontransgenic males are preferred. In the yellow fever mosquito Aedes aegypti, sex is determined by a pair of homomorphic sex chromosomes, with the dominant male-determining locus (the M locus) and its counterpart (the m locus) embedded in an M-bearing and an m-bearing chromosome 1, respectively. We utilized both naturally occurring and engineered sex-linked recessive lethal alleles (RLAs) to create sex separation strains for Ae. aegypti on the basis of differential elimination of marked sex chromosomes (DeMark). DeMark strains are self-sustaining and produce nontransgenic males that are readily separated from individuals carrying RLA- and transgene-marked m chromosomes. For example, the marked m chromosome in the heterozygous mother in some strains was only inherited by her female progeny due to RLA-mediated incompatibility with the M-bearing chromosome in the father, producing nontransgenic males and transgenic females, generation after generation. We further explore strategies to conditionally eliminate females that contain marked sex chromosomes. We also discuss DeMark designs that are applicable for efficient sex separation in organisms with well-differentiated X and Y chromosomes, such as the Anopheles mosquitoes.

CRISPR-Cas9 knockout of membrane-bound alkaline phosphatase or cadherin does not confer resistance to Cry toxins in Aedes aegypti

The Aedes aegypti cadherin-like protein (Aae-Cad) and the membrane-bound alkaline phosphatase (Aae-mALP) are membrane proteins identified as putative receptors for the larvicidal Cry toxins produced by Bacillus thuringiensis subsp. israelensis bacteria. Cry toxins are the most used toxins in the control of different agricultural pest and mosquitos. Despite the relevance of Aae-Cad and Aae-mALP as possible toxin-receptors in mosquitoes, previous efforts to establish a clear functional connection among them and Cry toxins activity have been relatively limited. In this study, we used CRISPR-Cas9 to generate knockout (KO) mutations of Aae-Cad and Aae-mALP. The Aae-mALP KO was successfully generated, in contrast to the Aae-Cad KO which was obtained only in females. The female-linked genotype was due to the proximity of aae-cad gene to the sex-determining loci (M:m). Both A. aegypti KO mutant populations were viable and their insect-development was not affected, although a tendency on lower egg hatching rate was observed. Bioassays were performed to assess the effects of these KO mutations on the susceptibility of A. aegypti to Cry toxins, showing that the Aae-Cad female KO or Aae-mALP KO mutations did not significantly alter the susceptibility of A. aegypti larvae to the mosquitocidal Cry toxins, including Cry11Aa, Cry11Ba, Cry4Ba, and Cry4Aa. These findings suggest that besides the potential participation of Aae-Cad and Aae-mALP as Cry toxin receptors in A. aegypti, additional midgut membrane proteins are involved in the mode of action of these insecticidal toxins.

Finding divergent sequences of homomorphic sex chromosomes via diploidized nanopore-based assembly from a single male

Although homomorphic sex chromosomes can have non-recombining regions with elevated sequence divergence between its complements, such divergence signals can be difficult to detect bioinformatically. If found in genomes of e.g. insect pests, these sequences could be targeted by the engineered genetic sexing and control systems. Here, we report an approach that can leverage long-read nanopore sequencing of a single XY male to identify divergent regions of homomorphic sex chromosomes. Long-read data are used for de novo genome assembly that is diploidized in a way that maximizes sex-specific differences between its haploid complements. We show that the correct assembly phasing is supported by the mapping of nanopore reads from the male's haploid Y-bearing sperm cells. The approach revealed a highly divergent region (HDR) near the centromere of the homomorphic sex chromosome of Aedes aegypti, the most important arboviral vector, for which there is a great interest in creating new genetic control tools. HDR is located ~5Mb downstream of the known male-determining locus on chromosome 1 and is significantly enriched for ovary-biased genes. While recombination in HDR ceased relatively recently (~1.4 MYA), HDR gametologs have divergent exons and introns of protein coding genes, and most lncRNA genes became X-specific. Megabases of previously invisible sex-linked sequences provide new putative targets for engineering the genetic systems to control this deadly mosquito. Broadly, our approach expands the toolbox for studying cryptic structure of sex chromosomes.

Heteromorphic ZZ/ZW sex chromosomes sharing gene content with mammalian XX/XY are conserved in Madagascan chameleons of the genus Furcifer

Chameleons are well-known lizards with unique morphology and physiology, but their sex determination has remained poorly studied. Madagascan chameleons of the genus Furcifer have cytogenetically distinct Z and W sex chromosomes and occasionally Z1Z1Z2Z2/Z1Z2W multiple neo-sex chromosomes. To identify the gene content of their sex chromosomes, we microdissected and sequenced the sex chromosomes of F. oustaleti (ZZ/ZW) and F. pardalis (Z1Z1Z2Z2/Z1Z2W). In addition, we sequenced the genomes of a male and a female of F. lateralis (ZZ/ZW) and F. pardalis and performed a comparative coverage analysis between the sexes. Despite the notable heteromorphy and distinctiveness in heterochromatin content, the Z and W sex chromosomes share approximately 90% of their gene content. This finding demonstrates poor correlation of the degree of differentiation of sex chromosomes at the cytogenetic and gene level. The test of homology based on the comparison of gene copy number variation revealed that female heterogamety with differentiated sex chromosomes remained stable in the genus Furcifer for at least 20 million years. These chameleons co-opted for the role of sex chromosomes the same genomic region as viviparous mammals, lacertids and geckos of the genus Paroedura, which makes these groups excellent model for studies of convergent and divergent evolution of sex chromosomes.

On the Origin and Evolution of the Mosquito Male-determining Factor Nix

The mosquito family Culicidae is divided into 2 subfamilies named the Culicinae and Anophelinae. Nix, the dominant male-determining factor, has only been found in the culicines Aedes aegypti and Aedes albopictus, 2 important arboviral vectors that belong to the subgenus Stegomyia. Here we performed sex-specific whole-genome sequencing and RNAseq of divergent mosquito species and explored additional male-inclusive datasets to investigate the distribution of Nix. Except for the Culex genus, Nix homologs were found in all species surveyed from the Culicinae subfamily, including 12 additional species from 3 highly divergent tribes comprising 4 genera, suggesting Nix originated at least 133 to 165 million years ago (MYA). Heterologous expression of 1 of 3 divergent Nix open reading frames (ORFs) in Ae. aegypti resulted in partial masculinization of genetic females as evidenced by morphology and doublesex splicing. Phylogenetic analysis suggests Nix is related to femaleless (fle), a recently described intermediate sex-determining factor found exclusively in anopheline mosquitoes. Nix from all species has a conserved structure, including 3 RNA-recognition motifs (RRMs), as does fle. However, Nix has evolved at a much faster rate than fle. The RRM3 of both Nix and fle are distantly related to the single RRM of a widely distributed and conserved splicing factor transformer-2 (tra2). The RRM3-based phylogenetic analysis suggests this domain in Nix and fle may have evolved from tra2 or a tra2-related gene in a common ancestor of mosquitoes. Our results provide insights into the evolution of sex determination in mosquitoes and will inform broad applications of mosquito-control strategies based on manipulating sex ratios toward nonbiting males.

Combining two genetic sexing strains allows sorting of non-transgenic males for Aedes genetic control

Chemical control of disease vectoring mosquitoes Aedes albopictus and Aedes aegypti is costly, unsustainable, and increasingly ineffective due to the spread of insecticide resistance. The Sterile Insect Technique is a valuable alternative but is limited by slow, error-prone, and wasteful sex-separation methods. Here, we present four Genetic Sexing Strains (two for each Aedes species) based on fluorescence markers linked to the m and M sex loci, allowing for the isolation of transgenic males. Furthermore, we demonstrate how combining these sexing strains enables the production of non-transgenic males. In a mass-rearing facility, 100,000 first instar male larvae could be sorted in under 1.5 h with an estimated 0.01-0.1% female contamination on a single machine. Cost-efficiency analyses revealed that using these strains could result in important savings while setting up and running a mass-rearing facility. Altogether, these Genetic Sexing Strains should enable a major upscaling in control programmes against these important vectors.

Strong Positive Selection in Aedes aegypti and the Rapid Evolution of Insecticide Resistance

Aedes aegypti vectors the pathogens that cause dengue, yellow fever, Zika virus, and chikungunya and is a serious threat to public health in tropical regions. Decades of work has illuminated many aspects of Ae. aegypti's biology and global population structure and has identified insecticide resistance genes; however, the size and repetitive nature of the Ae. aegypti genome have limited our ability to detect positive selection in this mosquito. Combining new whole genome sequences from Colombia with publicly available data from Africa and the Americas, we identify multiple strong candidate selective sweeps in Ae. aegypti, many of which overlap genes linked to or implicated in insecticide resistance. We examine the voltage-gated sodium channel gene in three American cohorts and find evidence for successive selective sweeps in Colombia. The most recent sweep encompasses an intermediate-frequency haplotype containing four candidate insecticide resistance mutations that are in near-perfect linkage disequilibrium with one another in the Colombian sample. We hypothesize that this haplotype may continue to rapidly increase in frequency and perhaps spread geographically in the coming years. These results extend our knowledge of how insecticide resistance has evolved in this species and add to a growing body of evidence suggesting that Ae. aegypti has an extensive genomic capacity to rapidly adapt to insecticide-based vector control.

Genomic basis of Y-linked dwarfism in cichlids pursuing alternative reproductive tactics

Sexually antagonistic selection, which favours different optima in males and females, is predicted to play an important role in the evolution of sex chromosomes. Body size is a sexually antagonistic trait in the shell-brooding cichlid fish Lamprologous callipterus, as "bourgeois" males must be large enough to carry empty snail shells to build nests whereas females must be small enough to fit into shells for breeding. In this species, there is also a second male morph: smaller "dwarf" males employ an alternative reproductive strategy by wriggling past spawning females into shells to fertilize eggs. L. callipterus male morphology is passed strictly from father to son, suggesting Y-linkage. However, sex chromosomes had not been previously identified in this species, and the genomic basis of size dimorphism was unknown. Here we used whole-genome sequencing to identify a 2.4-Mb sex-linked region on scaffold_23 with reduced coverage and single nucleotide polymorphism density in both male morphs compared to females. Within this sex region, distinct Y-haplotypes delineate the two male morphs, and candidate genes for body size (GHRHR, a known dwarfism gene) and sex determination (ADCYAP1R1) are in high linkage disequilibrium. Because differences in body size between females and males are under strong selection in L. callipterus, we hypothesize that sexual antagonism over body size initiated early events in sex chromosome evolution, followed by Y divergence to give rise to bourgeois and dwarf male reproductive strategies. Our results are consistent with the hypothesis that sexually antagonistic traits should be linked to young sex chromosomes.

Marker-assisted mapping enables forward genetic analysis in Aedes aegypti, an arboviral vector with vast recombination deserts

Aedes aegypti is a major vector of arboviruses that cause dengue, chikungunya, yellow fever, and Zika. Although recent success in reverse genetics has facilitated rapid progress in basic and applied research, integration of forward genetics with modern technologies remains challenging in this important species, as up to 47% of its chromosome is refractory to genetic mapping due to extremely low rate of recombination. Here, we report the development of a marker-assisted mapping strategy to readily screen for and genotype only the rare but informative recombinants, drastically increasing both the resolution and signal-to-noise ratio. Using marker-assisted mapping, we mapped a transgene that was inserted in a >100-Mb recombination desert and a sex-linked spontaneous red-eye (re) mutation just outside the region. We subsequently determined, by CRISPR/Cas9-mediated knockout, that cardinal is the causal gene of re, which is the first forward genetic identification of a causal gene in Ae. aegypti. The identification of the causal gene of the sex-linked re mutation provides the molecular foundation for using gene editing to develop versatile and stable genetic sexing methods. To facilitate genome-wide forward genetics in Ae. aegypti, we generated and compiled a number of lines with markers throughout the genome. Thus, by overcoming the challenges presented by the vast recombination deserts and the scarcity of markers, we have shown that effective forward genetic analysis is increasingly feasible in this important arboviral vector species.

Selective targeting of biting females to control mosquito-borne infectious diseases

Mosquitoes are vectors for a number of infectious diseases. Only females feed on blood to provision for their embryos and, in doing so, transmit pathogens to the associated vertebrate hosts. Therefore, sex is an important phenotype in the context of genetic control programs, both for sex separation in the rearing facilities to avoid releasing biting females and for ways to distort the sex ratio towards nonbiting males. We review recent progress in the fundamental knowledge of sex determination and sex chromosomes in mosquitoes and discuss new methods to achieve sex separation and sex ratio distortion to help control mosquito-borne infectious diseases. We conclude by suggesting a few critical areas for future research.

Retrogene Duplication and Expression Patterns Shaped by the Evolution of Sex Chromosomes in Malaria Mosquitoes

Genes that originate during evolution are an important source of novel biological functions. Retrogenes are functional copies of genes produced by retroduplication and as such are located in different genomic positions. To investigate retroposition patterns and retrogene expression, we computationally identified interchromosomal retroduplication events in nine portions of the phylogenetic history of malaria mosquitoes, making use of species that do or do not have classical sex chromosomes to test the roles of sex-linkage. We found 40 interchromosomal events and a significant excess of retroduplications from the X chromosome to autosomes among a set of young retrogenes. These young retroposition events occurred within the last 100 million years in lineages where all species possessed differentiated sex chromosomes. An analysis of available microarray and RNA-seq expression data for Anopheles gambiae showed that many of the young retrogenes evolved male-biased expression in the reproductive organs. Young autosomal retrogenes with increased meiotic or postmeiotic expression in the testes tend to be male biased. In contrast, older retrogenes, i.e., in lineages with undifferentiated sex chromosomes, do not show this particular chromosomal bias and are enriched for female-biased expression in reproductive organs. Our reverse-transcription PCR data indicates that most of the youngest retrogenes, which originated within the last 47.6 million years in the subgenus Cellia, evolved non-uniform expression patterns across body parts in the males and females of An. coluzzii. Finally, gene annotation revealed that mitochondrial function is a prominent feature of the young autosomal retrogenes. We conclude that mRNA-mediated gene duplication has produced a set of genes that contribute to mosquito reproductive functions and that different biases are revealed after the sex chromosomes evolve. Overall, these results suggest potential roles for the evolution of meiotic sex chromosome inactivation in males and of sexually antagonistic conflict related to mitochondrial energy function as the main selective pressures for X-to-autosome gene reduplication and testis-biased expression in these mosquito lineages.

Sex-specific distribution and classification of Wolbachia infections and mitochondrial DNA haplogroups in Aedes albopictus from the Indo-Pacific

The arbovirus vector Aedes albopictus (Asian tiger mosquito) is common throughout the Indo-Pacific region, where most global dengue transmission occurs. We analysed population genomic data and tested for cryptic species in 160 Ae. albopictus sampled from 16 locations across this region. We found no evidence of cryptic Ae. albopictus but found multiple intraspecific COI haplotypes partitioned into groups representing three Asian lineages: East Asia, Southeast Asia and Indonesia. Papua New Guinea (PNG), Vanuatu and Christmas Island shared recent coancestry, and Indonesia and Timor-Leste were likely invaded from East Asia. We used a machine learning trained on morphologically sexed samples to classify sexes using multiple genetic features and then characterized the wAlbA and wAlbB Wolbachia infections in 664 other samples. The wAlbA and wAlbB infections as detected by qPCR showed markedly different patterns in the sexes. For females, most populations had a very high double infection incidence, with 67% being the lowest value (from Timor-Leste). For males, the incidence of double infections ranged from 100% (PNG) to 0% (Vanuatu). Only 6 females were infected solely by the wAlbA infection, while rare uninfected mosquitoes were found in both sexes. The wAlbA and wAlbB densities varied significantly among populations. For mosquitoes from Torres Strait and Vietnam, the wAlbB density was similar in single-infected and superinfected (wAlbA and wAlbB) mosquitoes. There was a positive association between wAlbA and wAlbB infection densities in superinfected Ae. albopictus. Our findings provide no evidence of cryptic species of Ae. albopictus in the region and suggest site-specific factors influencing the incidence of Wolbachia infections and their densities. We also demonstrate the usefulness of ddRAD tag depths as sex-specific mosquito markers. The results provide baseline data for the exploitation of Wolbachia-induced cytoplasmic incompatibility (CI) in dengue control.

The mosquito Aedes albopictus transmits dengue and other arboviruses. This study investigates the genetics of these mosquitoes and their endosymbiont Wolbachia in the Indo-Pacific region, where 70% of global dengue transmission occurs. The analysis of mitochondrial DNA sequences showed no evidence of cryptic Ae. albopictus but suggested three Asian lineages: East Asia, Southeast Asia and Indonesia. Papua New Guinea, Vanuatu and Christmas Island shared recent coancestry, and Indonesia and Timor-Leste were likely invaded from East Asia. We used bioinformatics to classify sexes and then characterized the wAlbA and wAlbB Wolbachia infections via both bioinformatics and quantitative PCR. We found markedly different patterns of wAlbA and wAlbB infections in the sexes. The wAlbA and wAlbB densities varied significantly among populations, suggesting site-specific factors influencing the incidence of Wolbachia infections and their densities. We also demonstrate the usefulness of next generation sequencing data in developing molecular markers that can be repeatedly reanalysed to investigate new issues as these arise. These results provide baseline data for the exploitation of Wolbachia-induced cytoplasmic incompatibility in dengue control.

Transgenic expression of Nix converts genetic females into males and allows automated sex sorting in Aedes albopictus

Aedes albopictus is a major vector of arboviruses. Better understanding of its sex determination is crucial for developing mosquito control tools, especially genetic sexing strains. In Aedes aegypti, Nix is the primary gene responsible for masculinization and Nix-expressing genetic females develop into fertile, albeit flightless, males. In Ae. albopictus, Nix has also been implicated in masculinization but its role remains to be further characterized. In this work, we establish Ae. albopictus transgenic lines ectopically expressing Nix. Several are composed exclusively of genetic females, with transgenic individuals being phenotypic and functional males due to the expression of the Nix transgene. Their reproductive fitness is marginally impaired, while their flight performance is similar to controls. Overall, our results show that Nix is sufficient for full masculinization in Ae. albopictus. Moreover, the transgene construct contains a fluorescence marker allowing efficient automated sex sorting. Consequently, such strains constitute valuable sexing strains for genetic control.

Nix expression with a fluorescent marker in genetically female Ae. albopictus causes masculinization with minimal effects to reproductive fitness and flight performance.

A Conserved Female-Specific Requirement for the GGT Gene in Mosquito Larvae Facilitates RNAi-Mediated Sex Separation in Multiple Species of Disease Vector Mosquitoes

Although several emerging mosquito control technologies are dependent on mass releases of adult males, methods of sex-sorting that can be implemented globally have not yet been established. RNAi screens led to the discovery of siRNA, which targets gamma-glutamyl transpeptidase (GGT), a gene which is well conserved in multiple species of mosquitoes and located at the sex-determining M locus region in Aedes aegypti. Silencing the A. aegypti, Aedes albopictus, Anopheles gambiae, Culex pipiens, and Culex quinquefasciatus GGT genes resulted in female larval death, with no significant impact on male survival. Generation of yeast strains that permitted affordable expression and oral delivery of shRNA corresponding to mosquito GGT genes facilitated larval target gene silencing and generated significantly increased 5 males:1 female adult ratios in each species. Yeast targeting a conserved sequence in Culex GGT genes was incorporated into a larval mass-rearing diet, permitting the generation of fit adult male C. pipiens and C. quinquefasciatus, two species for which labor-intensive manual sex separation had previously been utilized. The results of this study indicate that female-specific yeast-based RNAi larvicides may facilitate global implementation of population-based control strategies that require releases of sterile or genetically modified adult males, and that yeast RNAi strategies can be utilized in various species of mosquitoes that have progressed to different stages of sex chromosome evolution.

A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata)

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7–10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

A conserved female-specific larval requirement for MtnB function facilitates sex separation in multiple species of disease vector mosquitoes

Background

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.

Results

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.

Conclusions

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.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04844-w.

A functional requirement for sex-determination M/m locus region lncRNA genes in Aedes aegypti female larvae

Although many putative long non-coding RNA (lncRNA) genes have been identified in insect genomes, few of these genes have been functionally validated. A screen for female-specific larvicides that facilitate Aedes aegypti male sex separation uncovered multiple interfering RNAs with target sites in lncRNA genes located in the M/m locus region, including loci within or tightly linked to the sex determination locus. Larval consumption of a Saccharomyces cerevisiae (yeast) strain engineered to express interfering RNA corresponding to lncRNA transcripts resulted in significant female death, yet had no impact on male survival or fitness. Incorporation of the yeast larvicides into mass culturing protocols facilitated scaled production and separation of fit adult males, indicating that yeast larvicides could benefit mosquito population control strategies that rely on mass releases of male mosquitoes. These studies functionally verified a female-specific developmental requirement for M/m locus region lncRNA genes, suggesting that sexually antagonistic lncRNA genes found within this highly repetitive pericentromeric DNA sequence may be contributing to the evolution of A. aegypti sex chromosomes.

Disagreement in FST estimators: A case study from sex chromosomes

Sewall Wright developed FST for describing population differentiation and it has since been extended to many novel applications, including the detection of homomorphic sex chromosomes. However, there has been confusion regarding the expected estimate of FST for a fixed difference between the X- and Y-chromosome when comparing males and females. Here, we attempt to resolve this confusion by contrasting two common FST estimators and explain why they yield different estimates when applied to the case of sex chromosomes. We show that this difference is true for many allele frequencies, but the situation characterized by fixed differences between the X- and Y-chromosome is among the most extreme. To avoid additional confusion, we recommend that all authors using FST clearly state which estimator of FST their work uses.

Exome-wide association study reveals largely distinct gene sets underlying specific resistance to dengue virus types 1 and 3 in Aedes aegypti

Although specific interactions between host and pathogen genotypes have been well documented in invertebrates, the identification of host genes involved in discriminating pathogen genotypes remains a challenge. In the mosquito Aedes aegypti, the main dengue virus (DENV) vector worldwide, statistical associations between host genetic markers and DENV types or strains were previously detected, but the host genes underlying this genetic specificity have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets. Here, we investigated the genetic architecture of DENV resistance in a population of Ae. aegypti from Bakoumba, Gabon, which displays a stronger resistance phenotype to DENV type 1 (DENV-1) than to DENV type 3 (DENV-3) infection. Following experimental exposure to either DENV-1 or DENV-3, we sequenced the exomes of large phenotypic pools of mosquitoes that are either resistant or susceptible to each DENV type. Using variation in single-nucleotide polymorphism (SNP) frequencies among the pools, we computed empirical p values based on average gene scores adjusted for the differences in SNP counts, to identify genes associated with infection in a DENV type-specific manner. Among the top 5% most significant genes, 263 genes were significantly associated with resistance to both DENV-1 and DENV-3, 287 genes were only associated with DENV-1 resistance and 290 were only associated with DENV-3 resistance. The shared significant genes were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity, whereas the genes uniquely associated with DENV-3 resistance were enriched in genes with zinc ion binding activity. Together, these results indicate that specific resistance to different DENV types relies on largely non-overlapping sets of genes in this Ae. aegypti population and pave the way for further mechanistic studies.

Compatibility between hosts and pathogens is often genetically specific in invertebrates but host genes underlying this genetic specificity have not been elucidated. We investigated the genetic architecture of dengue virus type-specific resistance in the mosquito vector Aedes aegypti. We used a natural population of Ae. aegypti from Bakoumba, Gabon, which is differentially resistant to dengue virus type 1 and dengue virus type 3. We surveyed genetic variation in protein-coding regions of the mosquito genome and compared the frequency of genetic polymorphisms between groups of mosquitoes that are either resistant or susceptible to each dengue virus type. We found that the Ae. aegypti genes associated with resistance to dengue virus type 1 or dengue virus type 3 were largely non-overlapping. This finding indicates that different sets of host genes, rather than different variants of the same genes, confer pathogen-specific resistance in this population. This study is an important step towards identification of mechanisms underlying the genetic specificity of invertebrate host-pathogen interactions.

Development of a confinable gene drive system in the human disease vector Aedes aegypti

Aedes aegypti is the principal mosquito vector for many arboviruses that increasingly infect millions of people every year. With an escalating burden of infections and the relative failure of traditional control methods, the development of innovative control measures has become of paramount importance. The use of gene drives has sparked significant enthusiasm for genetic control of mosquitoes; however, no such system has been developed in Ae. aegypti. To fill this void, here we develop several CRISPR-based split gene drives for use in this vector. With cleavage rates up to 100% and transmission rates as high as 94%, mathematical models predict that these systems could spread anti-pathogen effector genes into wild populations in a safe, confinable and reversible manner appropriate for field trials and effective for controlling disease. These findings could expedite the development of effector-linked gene drives that could safely control wild populations of Ae. aegypti to combat local pathogen transmission.

Proteins, Transcripts, and Genetic Architecture of Seminal Fluid and Sperm in the Mosquito Aedes aegypti

The yellow fever mosquito, Aedes aegypti,, transmits several viruses causative of serious diseases, including dengue, Zika, and chikungunya. Some proposed efforts to control this vector involve manipulating reproduction to suppress wild populations or to replace them with disease-resistant mosquitoes. The design of such strategies requires an intimate knowledge of reproductive processes, yet our basic understanding of reproductive genetics in this vector remains largely incomplete. To accelerate future investigations, we have comprehensively catalogued sperm and seminal fluid proteins (SFPs) transferred to females in the ejaculate using tandem mass spectrometry. By excluding female-derived proteins using an isotopic labeling approach, we identified 870 sperm proteins and 280 SFPs. Functional composition analysis revealed parallels with known aspects of sperm biology and SFP function in other insects. To corroborate our proteome characterization, we also generated transcriptomes for testes and the male accessory glands-the primary contributors to Ae. aegypti, sperm and seminal fluid, respectively. Differential gene expression of accessory glands from virgin and mated males suggests that transcripts encoding proteins involved in protein translation are upregulated post-mating. Several SFP transcripts were also modulated after mating, but >90% remained unchanged. Finally, a significant enrichment of SFPs was observed on chromosome 1, which harbors the male sex determining locus in this species. Our study provides a comprehensive proteomic and transcriptomic characterization of ejaculate production and composition and thus provides a foundation for future investigations of Ae. aegypti, reproductive biology, from functional analysis of individual proteins to broader examination of reproductive processes.

Repeated Evolution Versus Common Ancestry: Sex Chromosome Evolution in the Haplochromine Cichlid Pseudocrenilabrus philander

Why sex chromosomes turn over and remain undifferentiated in some taxa, whereas they degenerate in others, is still an area of ongoing research. The recurrent occurrence of homologous and homomorphic sex chromosomes in distantly related taxa suggests their independent evolution or continued recombination since their first emergence. Fishes display a great diversity of sex-determining systems. Here, we focus on sex chromosome evolution in haplochromines, the most species-rich lineage of cichlid fishes. We investigate sex-specific signatures in the Pseudocrenilabrus philander species complex, which belongs to a haplochromine genus found in many river systems and ichthyogeographic regions in northern, eastern, central, and southern Africa. Using whole-genome sequencing and population genetic, phylogenetic, and read-coverage analyses, we show that one population of P. philander has an XX-XY sex-determining system on LG7 with a large region of suppressed recombination. However, in a second bottlenecked population, we did not find any sign of a sex chromosome. Interestingly, LG7 also carries an XX-XY system in the phylogenetically more derived Lake Malawi haplochromine cichlids. Although the genomic regions determining sex are the same in Lake Malawi cichlids and P. philander, we did not find evidence for shared ancestry, suggesting that LG7 evolved as sex chromosome at least twice in haplochromine cichlids. Hence, our work provides further evidence for the labile nature of sex determination in fishes and supports the hypothesis that the same genomic regions can repeatedly and rapidly be recruited as sex chromosomes in more distantly related lineages.

Vgsc-interacting proteins are genetically associated with pyrethroid resistance in Aedes aegypti

Association mapping of factors that condition pyrethroid resistance in Aedes aegypti has consistently identified genes in multiple functional groups. Toward better understanding of the mechanisms involved, we examined high throughput sequencing data (HTS) from two Aedes aegypti aegypti collections from Merida, Yucatan, Mexico treated with either permethrin or deltamethrin. Exome capture enrichment for coding regions and the AaegL5 annotation were used to identify genes statistically associated with resistance. The frequencies of single nucleotide polymorphisms (SNPs) were compared between resistant and susceptible mosquito pools using a contingency χ² analysis. The -log₁₀(χ²p value) was calculated at each SNP site, with a weighted average determined from all sites in each gene. Genes with -log₁₀(χ²p value) ≥ 4.0 and present among all 3 treatment groups were subjected to gene set enrichment analysis (GSEA). We found that several functional groups were enriched compared to all coding genes. These categories were transport, signal transduction and metabolism, in order from highest to lowest statistical significance. Strikingly, 21 genes with demonstrated association to synaptic function were identified. In the high association group (n = 1,053 genes), several genes were identified that also genetically or physically interact with the voltage-gated sodium channel (VGSC). These genes were eg., CHARLATAN (CHL), a transcriptional regulator, several ankyrin-domain proteins, PUMILIO (PUM), a translational repressor, and NEDD4 (E3 ubiquitin-protein ligase). There were 13 genes that ranked among the top 10%: these included VGSC; CINGULIN, a predicted neuronal gap junction protein, and the aedine ortholog of NERVY (NVY), a transcriptional regulator. Silencing of CHL and NVY followed by standard permethrin bottle bioassays validated their association with permethrin resistance. Importantly, VGSC levels were also reduced about 50% in chl- or nvy-dsRNA treated mosquitoes. These results are consistent with the contribution of a variety of neuronal pathways to pyrethroid resistance in Ae. aegypti.

The Nix locus on the male-specific homologue of chromosome 1 in Aedes albopictus is a strong candidate for a male-determining factor

Background

Global concern over the rapid expansion of the Asian tiger mosquito, Aedes albopictus, and its vector competence has highlighted an urgent need to improve currently available population control methods, like the Sterile Insect Technique. Knowledge of the sex determination cascade is a prerequisite for the development of early-stage sexing systems. To this end, we have characterised the putative sex determination gene, Nix, in this species. In Aedes species the chromosome complement consists of three pairs of chromosomes. The sex determination alleles are linked to the smallest homomorphic chromosome.

Results

We identified the male-specific chromosome 1 of Ae. albopictus that carries the putative male-determining gene Nix. We have also characterised the complete genomic sequence of the Nix gene which is composed of two exons and a short intron. The gene displays different levels of intron retention during development. Comparison of DNA sequences covering most of the Nix gene from individuals across the species range revealed no polymorphism.

Conclusions

Our characterisation of the Nix gene in Ae. albopictus represents an initial step in the analysis of the sex determination cascade in this species. We found evidence of intron retention (IR) in Nix. IR might play a role in regulating the expression of Nix during development. Our results provide the basis for the development of new genetic control strategies.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3215-8) contains supplementary material, which is available to authorized users.

Improved reference genome of Aedes aegypti informs arbovirus vector control

Female Aedes aegypti mosquitoes infect more than 400 million people each year with dangerous viral pathogens including dengue, yellow fever, Zika and chikungunya. Progress in understanding the biology of mosquitoes and developing the tools to fight them has been slowed by the lack of a high-quality genome assembly. Here we combine diverse technologies to produce the markedly improved, fully re-annotated AaegL5 genome assembly, and demonstrate how it accelerates mosquito science. We anchored physical and cytogenetic maps, doubled the number of known chemosensory ionotropic receptors that guide mosquitoes to human hosts and egg-laying sites, provided further insight into the size and composition of the sex-determining M locus, and revealed copy-number variation among glutathione S-transferase genes that are important for insecticide resistance. Using high-resolution quantitative trait locus and population genomic analyses, we mapped new candidates for dengue vector competence and insecticide resistance. AaegL5 will catalyse new biological insights and intervention strategies to fight this deadly disease vector.

The sequence of a male-specific genome region containing the sex determination switch in Aedes aegypti

Background

Aedes aegypti is the principal vector of several important arboviruses. Among the methods of vector control to limit transmission of disease are genetic strategies that involve the release of sterile or genetically modified non-biting males, which has generated interest in manipulating mosquito sex ratios. Sex determination in Ae. aegypti is controlled by a non-recombining Y chromosome-like region called the M locus, yet characterisation of this locus has been thwarted by the repetitive nature of the genome. In 2015, an M locus gene named Nix was identified that displays the qualities of a sex determination switch.

Results

With the use of a whole-genome bacterial artificial chromosome (BAC) library, we amplified and sequenced a ~200 kb region containing the male-determining gene Nix. In this study, we show that Nix is comprised of two exons separated by a 99 kb intron primarily composed of repetitive DNA, especially transposable elements.

Conclusions

Nix, an unusually large and highly repetitive gene, exhibits features in common with Y chromosome genes in other organisms. We speculate that the lack of recombination at the M locus has allowed the expansion of repeats in a manner characteristic of a sex-limited chromosome, in accordance with proposed models of sex chromosome evolution in insects.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3090-3) contains supplementary material, which is available to authorized users.

Shared and Species-Specific Patterns of Nascent Y Chromosome Evolution in Two Guppy Species

Sex chromosomes form once recombination is halted around the sex-determining locus between a homologous pair of chromosomes, resulting in a male-limited Y chromosome. We recently characterized the nascent sex chromosome system in the Trinidadian guppy (Poeciliareticulata). The guppy Y is one of the youngest animal sex chromosomes yet identified, and therefore offers a unique window into the early evolutionary forces shaping sex chromosome formation, particularly the rate of accumulation of repetitive elements and Y-specific sequence. We used comparisons between male and female genomes in P. reticulata and its sister species, Endler’s guppy (P. wingei), which share an ancestral sex chromosome, to identify male-specific sequences and to characterize the degree of differentiation between the X and Y chromosomes. We identified male-specific sequence shared between P. reticulata and P. wingei consistent with a small ancestral non-recombining region. Our assembly of this Y-specific sequence shows substantial homology to the X chromosome, and appears to be significantly enriched for genes implicated in pigmentation. We also found two plausible candidates that may be involved in sex determination. Furthermore, we found that the P. wingei Y chromosome exhibits a greater signature of repetitive element accumulation than the P. reticulata Y chromosome. This suggests that Y chromosome divergence does not necessarily correlate with the time since recombination suppression. Overall, our results reveal the early stages of Y chromosome divergence in the guppy.

Alternative patterns of sex chromosome differentiation in Aedes aegypti (L)

Background

Some populations of West African Aedes aegypti, the dengue and zika vector, are reproductively incompatible; our earlier study showed that divergence and rearrangements of genes on chromosome 1, which bears the sex locus (M), may be involved. We also previously described a proposed cryptic subspecies SenAae (PK10, Senegal) that had many more high inter-sex F_ST genes on chromosome 1 than did Ae.aegypti aegypti (Aaa, Pai Lom, Thailand). The current work more thoroughly explores the significance of those findings.

Results

Intersex standardized variance (F_ST) of single nucleotide polymorphisms (SNPs) was characterized from genomic exome capture libraries of both sexes in representative natural populations of Aaa and SenAae. Our goal was to identify SNPs that varied in frequency between males and females, and most were expected to occur on chromosome 1. Use of the assembled AaegL4 reference alleviated the previous problem of unmapped genes. Because the M locus gene nix was not captured and not present in AaegL4, the male-determining locus, per se, was not explored. Sex-associated genes were those with F_ST values ≥ 0.100 and/or with increased expected heterozygosity (H_exp, one-sided T-test, p < 0.05) in males. There were 85 genes common to both collections with high inter-sex F_ST values; all genes but one were located on chromosome 1. Aaa showed the expected cluster of high inter-sex F_ST genes proximal to the M locus, whereas SenAae had inter-sex F_ST genes along the length of chromosome 1. In addition, the Aaa M-locus proximal region showed increased H_exp levels in males, whereas SenAae did not. In SenAae, chromosomal rearrangements and subsequent suppressed recombination may have accelerated X-Y differentiation.

Conclusions

The evidence presented here is consistent with differential evolution of proto-Y chromosomes in Aaa and SenAae.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4348-4) contains supplementary material, which is available to authorized users.