Transcription regulation of sex-biased genes during ontogeny in the malaria vector Anopheles gambiae

Comparative transcriptomics in East African cichlids reveals sex- and species-specific expression and new candidates for sex differentiation in fishes

Males and females of the same species differ largely in gene expression, which accounts for most of the morphological and physiological differences and sex-specific phenotypes. Here, we analyzed sex-specific gene expression in the brain and the gonads of cichlid fishes from Lake Tanganyika belonging to four different lineages, so-called tribes (Eretmodini, Ectodini, Haplochromini, and Lamprologini), using the outgroup Nile tilapia (Oreochromis niloticus) as reference. The comparison between male and female brains revealed few differences between the sexes, consistent in all investigated species. The gonads, on the other hand, showed a large fraction of differentially expressed transcripts with the majority of them showing the same direction of expression in all four species. All here-studied cichlids, especially the three investigated mouth-breeding species, showed a trend toward more male- than female biased transcripts. Transcripts, which were female-biased in expression in all four species, were overrepresented on linkage group (LG)1 in the reference genome and common male-biased transcripts showed accumulation on LG23, the presumable sex chromosomes of the Nile tilapia. Sex-specific transcripts contained candidate genes for sex determination and differentiation in fishes,especially members of the transforming growth factor-b-superfamily and the Wnt-pathway and also prominent members of the sox-, dm-domain-, and high mobility group-box families. We further confirmed our previous finding on species/lineage-specific gene expression shifts in the sex steroid pathway, including synthesizing enzymes as the aromatase cyp19a1 and estrogen and androgen receptors.

Multi-character approach reveals a discordant pattern of phenotypic variation during ontogeny in Culex pipiens biotypes (Diptera: Culicidae)

Culex (Culex) pipiens s.l. (Diptera: Culicidae) comprises two distinct biotypes, pipiens ('rural') and molestus ('urban'), both of which are thought to have differing capacities due to different host preferences. To better understand West Nile encephalitis epidemiology and improve risk assessment, local distinction between these forms is essential. This study assesses phenotypic variation at larval and adult stages of 'urban' and 'rural' biotypes of the species by complementary use of meristic, univariate and multivariate traits analyzed by traditional and geometric morphometrics. Third- and fourth-instar larvae from a broad area of the city of Novi Sad (Serbia) were collected and reared in the laboratory. After adult eclosion, the sex of each larva was recorded based on the sex of the corresponding adult. Examination of the association between variations of larval traits revealed contrasting variations regarding pecten spines vs. siphonal size and siphonal shape in the 'rural' biotype. Siphons of larvae collected in marshes and forest ecosystems outside urban areas were found to be the largest, but possessed the smallest number of pecten spines. In addition, statistically significant female-biased sexual dimorphism was observed in siphonal size, wing size and wing shape. Finally, we propose that an integrative approach is essential in delimitation of Cx. pipiens s.l. biotypes, since their differentiation was not possible based solely on larval and adult traits. Our findings shed light on the phenotypic plasticity important for population persistence in the changing environment of these medically important taxa.

The molecular and immunochemical expression of innexins in the yellow fever mosquito, Aedes aegypti: insights into putative life stage- and tissue-specific functions of gap junctions

Gap junctions (GJ) mediate direct intercellular communication by forming channels through which certain small molecules and/or ions can pass. Connexins, the proteins that form vertebrate GJ, are well studied and known to contribute to neuronal, muscular and epithelial physiology. Innexins, the GJ proteins of insects, have only recently received much investigative attention and many of their physiological roles remain to be determined. Here we characterize the molecular expression of six innexin (Inx) genes in the yellow fever mosquito Aedes aegypti (AeInx1, AeInx2, AeInx3, AeInx4, AeInx7, and AeInx8) and the immunochemical expression of one innexin protein, AeInx3, in the alimentary canal. We detected the expression of no less than four innexin genes in each mosquito life stage (larva, pupa, adult) and tissue/body region from adult males and females (midgut, Malpighian tubules, hindgut, head, carcass, gonads), suggesting a remarkable potential molecular diversity of GJ in mosquitoes. Moreover, the expression patterns of some innexins were life stage and/or tissue specific, suggestive of potential functional specializations. Cloning of the four full-length cDNAs expressed in the Malpighian tubules of adult females (AeInx1, AeInx2, AeInx3, and AeInx7) revealed evidence for 1) alternative splicing of AeInx1 and AeInx3 transcripts, and 2) putative N-glycosylation of AeInx3 and AeInx7. Finally, immunohistochemistry of AeInx3 in the alimentary canal of larval and adult female mosquitoes confirmed localization of this innexin to the intercellular regions of Malpighian tubule and hindgut epithelial cells, suggesting that it is an important component of GJ in these tissues.

Maternal germline-specific genes in the Asian malaria mosquito Anopheles stephensi: characterization and application for disease control

Anopheles stephensi is a principal vector of urban malaria on the Indian subcontinent and an emerging model for molecular and genetic studies of mosquito biology. To enhance our understanding of female mosquito reproduction, and to develop new tools for basic research and for genetic strategies to control mosquito-borne infectious diseases, we identified 79 genes that displayed previtellogenic germline-specific expression based on RNA-Seq data generated from 11 life stage-specific and sex-specific samples. Analysis of this gene set provided insights into the biology and evolution of female reproduction. Promoters from two of these candidates, vitellogenin receptor and nanos, were used in independent transgenic cassettes for the expression of artificial microRNAs against suspected mosquito maternal-effect genes, discontinuous actin hexagon and myd88. We show these promoters have early germline-specific expression and demonstrate 73% and 42% knockdown of myd88 and discontinuous actin hexagon mRNA in ovaries 48 hr after blood meal, respectively. Additionally, we demonstrate maternal-specific delivery of mRNA and protein to progeny embryos. We discuss the application of this system of maternal delivery of mRNA/miRNA/protein in research on mosquito reproduction and embryonic development, and for the development of a gene drive system based on maternal-effect dominant embryonic arrest.

Sex-biased gene expression and sexual conflict throughout development

Sex-biased gene expression is likely to account for most sexually dimorphic traits because males and females share much of their genome. When fitness optima differ between sexes for a shared trait, sexual dimorphism can allow each sex to express their optimum trait phenotype, and in this way, the evolution of sex-biased gene expression is one mechanism that could help to resolve intralocus sexual conflict. Genome-wide patterns of sex-biased gene expression have been identified in a number of studies, which we review here. However, very little is known about how sex-biased gene expression relates to sex-specific fitness and about how sex-biased gene expression and conflict vary throughout development or across different genotypes, populations, and environments. We discuss the importance of these neglected areas of research and use data from a small-scale experiment on sex-specific expression of genes throughout development to highlight potentially interesting avenues for future research.

Innexin AGAP001476 is critical for mediating anti-Plasmodium responses in Anopheles mosquitoes

The Toll and IMD pathways are known to be induced upon Plasmodium berghei and Plasmodium falciparum infection, respectively. It is unclear how Plasmodium or other pathogens in the blood meal and their invasion of the midgut epithelium would trigger the innate immune responses in immune cells, in particular hemocytes. Gap junctions, which can mediate both cell-to-cell and cell-to-extracellular communication, may participate in this signal transduction. This study examined whether innexins, gap junction proteins in insects, are involved in anti-Plasmodium responses in Anopheles gambiae. Inhibitor studies using carbenoxolone indicated that blocking innexons resulted in an increase in Plasmodium oocyst number and infection prevalence. This was accompanied by a decline in TEP1 levels in carbenoxolone-treated mosquitoes. Innexin AGAP001476 mRNA levels in midguts were induced during Plasmodium infection and a knockdown of AGAP001476, but not AGAP006241, caused an induction in oocyst number. Silencing AGAP001476 caused a concurrent increase in vitellogenin levels, a TEP1 inhibitor, in addition to a reduced level of TEP1-LRIM1-APL1C complex in hemolymph. Both vitellogenin and TEP1 are regulated by Cactus under the Toll pathway. Simultaneous knockdown of cactus and AGAP001476 failed to reverse the near refractoriness induced by the knockdown of cactus, suggesting that the AGAP001476-mediated anti-Plasmodium response is Cactus-dependent. These data demonstrate a critical role for innexin AGAP001476 in mediating innate immune responses against Plasmodium through Toll pathway in mosquitoes.

Analysis of the genitalia rotation in the male Anopheles funestus (Diptera: Culicidae)

Anopheles funestus is a major malaria vector in Africa. Insecticide resistance has developed in populations of this species in several African countries, prompting the need to develop additional vector control methods such as the sterile insect technique (SIT). This technique requires an understanding of those underlying physiological events that lead to sexual maturity of An. funestus males, the rotation of their genitalia in particular. The aim of this study was to qualitatively and quantitatively describe genital rotation in An. funestus males as it is an essential function of sexual maturation. Genital rotation of all the males reached its final rotation stage (135-180° rotation) 36 h post emergence at 23 ± 1 °C in laboratory colonised An. funestus males. These males had a comparable rotation rate to wild caught An. funestus at the same temperature setting. A temperature change (either 18 ± 1 °C or 29 ± 1 °C versus 23 ± 1 °C) significantly influenced the genital rotation rate such that this rate increased with increasing temperature. This information enhances our knowledge of the An. funestus male biology. This is important in terms of applying the sterile insect technique as the understanding and manipulation of the rate of sexual maturation in males has implications for the timing of sterile male release.

The developmental transcriptome of the mosquito Aedes aegypti, an invasive species and major arbovirus vector

Mosquitoes are vectors of a number of important human and animal diseases. The development of novel vector control strategies requires a thorough understanding of mosquito biology. To facilitate this, we used RNA-seq to identify novel genes and provide the first high-resolution view of the transcriptome throughout development and in response to blood feeding in a mosquito vector of human disease, Aedes aegypti, the primary vector for Dengue and yellow fever. We characterized mRNA expression at 34 distinct time points throughout Aedes development, including adult somatic and germline tissues, by using polyA+ RNA-seq. We identify a total of 14,238 novel new transcribed regions corresponding to 12,597 new loci, as well as many novel transcript isoforms of previously annotated genes. Altogether these results increase the annotated fraction of the transcribed genome into long polyA+ RNAs by more than twofold. We also identified a number of patterns of shared gene expression, as well as genes and/or exons expressed sex-specifically or sex-differentially. Expression profiles of small RNAs in ovaries, early embryos, testes, and adult male and female somatic tissues also were determined, resulting in the identification of 38 new Aedes-specific miRNAs, and ~291,000 small RNA new transcribed regions, many of which are likely to be endogenous small-interfering RNAs and Piwi-interacting RNAs. Genes of potential interest for transgene-based vector control strategies also are highlighted. Our data have been incorporated into a user-friendly genome browser located at www.Aedes.caltech.edu, with relevant links to Vectorbase (www.vectorbase.org)

Masculinization of gene expression is associated with exaggeration of male sexual dimorphism

Gene expression differences between the sexes account for the majority of sexually dimorphic phenotypes, and the study of sex-biased gene expression is important for understanding the genetic basis of complex sexual dimorphisms. However, it has been difficult to test the nature of this relationship due to the fact that sexual dimorphism has traditionally been conceptualized as a dichotomy between males and females, rather than an axis with individuals distributed at intermediate points. The wild turkey (Meleagris gallopavo) exhibits just this sort of continuum, with dominant and subordinate males forming a gradient in male secondary sexual characteristics. This makes it possible for the first time to test the correlation between sex-biased gene expression and sexually dimorphic phenotypes, a relationship crucial to molecular studies of sexual selection and sexual conflict. Here, we show that subordinate male transcriptomes show striking multiple concordances with their relative phenotypic sexual dimorphism. Subordinate males were clearly male rather than intersex, and when compared to dominant males, their transcriptomes were simultaneously demasculinized for male-biased genes and feminized for female-biased genes across the majority of the transcriptome. These results provide the first evidence linking sexually dimorphic transcription and sexually dimorphic phenotypes. More importantly, they indicate that evolutionary changes in sexual dimorphism can be achieved by varying the magnitude of sex-bias in expression across a large proportion of the coding content of a genome.

Males and females exhibit many differences in morphology, behavior and physiology, yet they share the vast majority of their genomes. Most differences between the sexes are therefore thought to be the product of gene expression differences between females and males. Studies of sex differences in expression assume that genes expressed more in males encode male traits, and genes expressed more in females encode female traits, and this assumption is a key foundation to genetic studies of sexual dimorphism and sexual conflict. Despite this key assumption, this relationship has yet to be empirically tested, as the main model organisms for studies of sex-biased gene expression lack multiple male and female morphs. Here, we use the two male morphs in the wild turkey to show that the magnitude of male-biased gene expression correlates with the manifestation of sexually dimorphic traits. Males with less manifestation of sexual dimorphism in phenotype were both demasculinized for male-biased genes, as well as feminized for female-biased genes. This pattern encompassed the majority of expressed loci, suggesting that evolutionary changes in the magnitude of sexual dimorphism may be achieved by small changes in the magnitude of sex-biased transcription across thousands of genes.

CluGene: A Bioinformatics Framework for the Identification of Co-Localized, Co-Expressed and Co-Regulated Genes Aimed at the Investigation of Transcriptional Regulatory Networks from High-Throughput Expression Data

The full understanding of the mechanisms underlying transcriptional regulatory networks requires unravelling of complex causal relationships. Genome high-throughput technologies produce a huge amount of information pertaining gene expression and regulation; however, the complexity of the available data is often overwhelming and tools are needed to extract and organize the relevant information. This work starts from the assumption that the observation of co-occurrent events (in particular co-localization, co-expression and co-regulation) may provide a powerful starting point to begin unravelling transcriptional regulatory networks. Co-expressed genes often imply shared functional pathways; co-expressed and functionally related genes are often co-localized, too; moreover, co-expressed and co-localized genes are also potential targets for co-regulation; finally, co-regulation seems more frequent for genes mapped to proximal chromosome regions. Despite the recognized importance of analysing co-occurrent events, no bioinformatics solution allowing the simultaneous analysis of co-expression, co-localization and co-regulation is currently available. Our work resulted in developing and valuating CluGene, a software providing tools to analyze multiple types of co-occurrences within a single interactive environment allowing the interactive investigation of combined co-expression, co-localization and co-regulation of genes. The use of CluGene will enhance the power of testing hypothesis and experimental approaches aimed at unravelling transcriptional regulatory networks. The software is freely available at http://bioinfolab.unipg.it/.

Regulation of Anopheles gambiae male accessory gland genes influences postmating response in female

In Drosophila, the accessory gland proteins (Acps) secreted from the male accessory glands (MAGs) and transferred along with sperm into the female reproductive tract have been implicated in triggering postmating behavioral changes, including refractoriness to subsequent mating and propensity to egg laying. Recently, Acps have been found also in Anopheles, suggesting similar functions. Understanding the mechanisms underlying transcriptional regulation of Acps and their functional role in modulating Anopheles postmating behavior may lead to the identification of novel vector control strategies to reduce mosquito populations. We identified heat-shock factor (HSF) binding sites within the Acp promoters of male Anopheles gambiae and discovered three distinct Hsf isoforms; one being significantly up-regulated in the MAGs after mating. Through genome-wide transcription analysis of Hsf-silenced males, we observed significant down-regulation in 50% of the Acp genes if compared to control males treated with a construct directed against an unrelated bacterial sequence. Treated males retained normal life span and reproductive behavior compared to control males. However, mated wild-type females showed a ∼46% reduction of egg deposition rate and a ∼23% reduction of hatching rate (∼58% combined reduction of progeny). Our results highlight an unsuspected role of HSF in regulating Acp transcription in A. gambiae and provide evidence that Acp down-regulation in males leads a significant reduction of progeny, thus opening new avenues toward the development of novel vector control strategies.

The genetics of obligate parthenogenesis in an aphid species and its consequences for the maintenance of alternative reproductive modes

Although loss of sex is widespread among metazoans, the genetic mechanisms underlying the transition to asexuality are poorly understood. Aphids are good models to address this issue because they frequently show reproductive-mode variation at the species level, involving cyclical parthenogens (CP) that reproduce sexually once a year and obligate parthenogens (OP) that reproduce asexually all year round. Here, we explore the genetic basis of OP in the cereal aphid Sitobion avenae by crossing several genotypes with contrasting reproductive modes and then characterising the reproductive phenotypes of F1 and F2 offspring. The analysis of phenotypic variation in F1 and F2 progenies suggests that at least two autosomal loci control OP in S. avenae. First, the transition to asexuality seems to depend on a single recessive locus, because the offspring from self-crossed cyclical parthenogenetic genotypes contain either 0 or 25% OP. Second, as we observed OP in the F1 progenies from crosses between CP and OP, and some CP in the offspring from outcrossed OP, a dominant 'suppressor' gene may also be involved, being inactive when in a recessive homozygous state in CP; this is the most parsimonious explanation for these results. This oligogenic inheritance of OP in S. avenae appears to be an efficient genetic system to generate new OP genotypes continually. It also allows asexuality-inducing alleles to be protected locally during harsh winters when extreme frost kills most OP, and then to spread very quickly after winter.

High-throughput sorting of mosquito larvae for laboratory studies and for future vector control interventions

Background

Mosquito transgenesis offers new promises for the genetic control of vector-borne infectious diseases such as malaria and dengue fever. Genetic control strategies require the release of large number of male mosquitoes into field populations, whether they are based on the use of sterile males (sterile insect technique, SIT) or on introducing genetic traits conferring refractoriness to disease transmission (population replacement). However, the current absence of high-throughput techniques for sorting different mosquito populations impairs the application of these control measures.

Methods

A method was developed to generate large mosquito populations of the desired sex and genotype. This method combines flow cytometry and the use of Anopheles gambiae transgenic lines that differentially express fluorescent markers in males and females.

Results

Fluorescence-assisted sorting allowed single-step isolation of homozygous transgenic mosquitoes from a mixed population. This method was also used to select wild-type males only with high efficiency and accuracy, a highly desirable tool for genetic control strategies where the release of transgenic individuals may be problematic. Importantly, sorted males showed normal mating ability compared to their unsorted brothers.

Conclusions

The developed method will greatly facilitate both laboratory studies of mosquito vectorial capacity requiring high-throughput approaches and future field interventions in the fight against infectious disease vectors.

Demasculinization of the Anopheles gambiae X chromosome

Background

In a number of organisms sex-biased genes are non-randomly distributed between autosomes and the shared sex chromosome X (or Z). Studies on Anopheles gambiae have produced conflicting results regarding the underrepresentation of male-biased genes on the X chromosome and it is unclear to what extent sexual antagonism, dosage compensation or X-inactivation in the male germline, the evolutionary forces that have been suggested to affect the chromosomal distribution of sex-biased genes, are operational in Anopheles.

Results

We performed a meta-analysis of sex-biased gene expression in Anopheles gambiae which provides evidence for a general underrepresentation of male-biased genes on the X-chromosome that increased in significance with the observed degree of sex-bias. A phylogenomic comparison between Drosophila melanogaster, Aedes aegypti and Culex quinquefasciatus also indicates that the Anopheles X chromosome strongly disfavours the evolutionary conservation of male-biased expression and that novel male-biased genes are more likely to arise on autosomes. Finally, we demonstrate experimentally that transgenes situated on the Anopheles gambiae X chromosome are transcriptionally silenced in the male germline.

Conclusion

The data presented here support the hypothesis that the observed demasculinization of the Anopheles X chromosome is driven by X-chromosome inactivation in the male germline and by sexual antagonism. The demasculinization appears to be the consequence of a loss of male-biased expression, rather than a failure in the establishment or the extinction of male-biased genes.

Spermless males elicit large-scale female responses to mating in the malaria mosquito Anopheles gambiae

Anopheles gambiae sensu stricto is the major vector of malaria, a disease with devastating consequences for human health. Given the constant spread of the disease, alternative approaches to the use of insecticides are urgently needed to control vector populations. Females of this species undergo large behavioral changes after mating, which include a life-long refractoriness to further insemination and the induction of egg laying in blood-fed individuals. Genetic control strategies aimed at impacting Anopheles fertility through the release of sterile males are being advocated to reduce the size of mosquito field populations. Such strategies depend on the ability of the released sterile males to mate successfully with wild females and to switch off the female receptivity to further copulation. Here we evaluate the role of sperm in regulating female behavioral responses after mating in An. gambiae. We developed spermless males by RNAi silencing of a germ cell differentiation gene. These males mated successfully and preserved standard accessory gland functions. Females mated to spermless males exhibited normal postcopulatory responses, which included laying large numbers of eggs upon blood feeding and becoming refractory to subsequent insemination. Moreover, spermless males induced transcriptional changes in female reproductive genes comparable to those elicited by fertile males. Our data demonstrate that, in contrast to Drosophila, targeting sperm in An. gambiae preserves normal male and female reproductive behavior for the traits and time frame analyzed and validate the use of approaches based on incapacitation or elimination of sperm for genetic control of vector populations to block malaria transmission.