The availability of genome sequences from 16 anopheline species provides unprecedented opportunities to study the evolution of reproductive traits relevant for malaria transmission. In Anopheles gambiae, a likely candidate for sexual selection is male 20-hydroxyecdysone (20E). Sexual transfer of this steroid hormone as part of a mating plug dramatically changes female physiological processes intimately tied to vectorial capacity. By combining phenotypic studies with ancestral state reconstructions and phylogenetic analyses, we show that mating plug transfer and male 20E synthesis are both derived characters that have coevolved in anophelines, driving the adaptation of a female 20E-interacting protein that promotes oogenesis via mechanisms also favoring Plasmodium survival. Our data reveal coevolutionary dynamics of reproductive traits between the sexes likely to have shaped the ability of anophelines to transmit malaria.
Fig S1 Species phylogeny based on single copy orthologs_protein sequence alignment
The concatenated protein sequence alignments of nine anopheline and two outgroup dipteran species (created using MUSCLE and trimmed with trimAl) of 4,829 relaxed single-copy orthologs (a maximum of three paralogs allowed in no more than two species, longest protein selected) determined by OrthoDB orthology delineation (available at http://cegg.unige.ch/orthodbmoz2), which resulted in 3,400,861 amino acid columns (with 906,126 distinct alignment patterns).
Fig 4 Phylogenetic analysis of MISO and AGAP002621_protein sequence alignment
MISO and AGAP002621 alignment of orthologous protein sequences from nine anopheline and two outgroup dipteran species retrieved from VectorBase (https://www.vectorbase.org). For each species, MISO andAGAP002621 orthologs were identified via BLAST/Translated blat search in VectorBase. However in some instances the gene models of MISO and AGAP002621 were fused and as a result were manually resolved using Apollo Genome Annotation and Curation Tool, Version 1.11.8. Gene identifiers for MISO and 2621 orthologs respectively are listed below with single identifiers corresponding to fused gene models: A. arabiensis (AARA004231); A. stephensi (ASTE005249); A. funestus (AFUN011077); A. dirus (ADIR010237); A. farauti (AFAF007180, AFAF018121); A. sinensis (ASIS013552, ASIS009594); A. atroparvus (AATE003240, AATE008584); A. albimanus (AALB003668); A. aegypti (AAEL003936, AAEL017464); D. melanogaster (FBgn0003065, no AGAP002621 ortholog present).