Genetic basis for reproductive diapause is correlated with life history traits within the Culex pipiens complex

Random Forest Analysis of Impact of Abiotic Factors on Culex pipiens and Culex quinquefasciatus Occurrence

The Culex pipiens complex of mosquitoes are significant vectors of several pathogens resulting in infectious human diseases in North America, including but not limited to West Nile encephalitis, Rift Valley Fever, and Lymphatic filariasis. Among this complex are C. pipiens form pipiens and Culex quinquefasciatus. While morphologically similar, the mosquitoes exhibit unique life histories that suit them uniquely to divergent niches, wherein C. pipiens can thrive despite the cold winters of the northern United States and C. quinquefasciatus is able to survive periods of drought typical in the southern states. Here, Random Forests machine-learning algorithms were employed to model and explore which environmental parameters best explain mosquito occurrence in historical trapping data across the continental United States of America, and test correlation with abundance data. The models explained between 71 and 97% of the presence or absence of the two mosquitoes based on historical climatic data. The results of this study will improve vector management programs by explaining which environmental variables will provide the most accurate predictions of mosquito presence at a given site.

Differential Gene Expression in the Heads of Behaviorally Divergent Culex pipiens Mosquitoes

Host preferences of Cx. pipiens, a bridge vector for West Nile virus to humans, have the potential to drive pathogen transmission dynamics. Yet much remains unknown about the extent of variation in these preferences and their molecular basis. We conducted host choice assays in a laboratory setting to quantify multi-day human and avian landing rates for Cx. pipiens females. Assayed populations originated from five above-ground and three below-ground breeding and overwintering habitats. All three below-ground populations were biased toward human landings, with rates of human landing ranging from 69-85%. Of the five above-ground populations, four had avian landing rates of >80%, while one landed on the avian host only 44% of the time. Overall response rates and willingness to alternate landing on the human and avian hosts across multiple days of testing also varied by population. For one human- and one avian-preferring population, we examined patterns of differential expression and splice site variation at genes expressed in female heads. We also compared gene expression and splice site variation within human-seeking females in either gravid or host-seeking physiological states to identify genes that may regulate blood feeding behaviors. Overall, we identified genes with metabolic and regulatory function that were differentially expressed in our comparison of gravid and host-seeking females. Differentially expressed genes in our comparison of avian- and human-seeking females were enriched for those involved in sensory perception. We conclude with a discussion of specific sensory genes and their potential influence on the divergent behaviors of avian- and human-seeking Cx. pipiens.

Unique genetic signatures of local adaptation over space and time for diapause, an ecologically relevant complex trait, in Drosophila melanogaster

Organisms living in seasonally variable environments utilize cues such as light and temperature to induce plastic responses, enabling them to exploit favorable seasons and avoid unfavorable ones. Local adapation can result in variation in seasonal responses, but the genetic basis and evolutionary history of this variation remains elusive. Many insects, including Drosophila melanogaster, are able to undergo an arrest of reproductive development (diapause) in response to unfavorable conditions. In D. melanogaster, the ability to diapause is more common in high latitude populations, where flies endure harsher winters, and in the spring, reflecting differential survivorship of overwintering populations. Using a novel hybrid swarm-based genome wide association study, we examined the genetic basis and evolutionary history of ovarian diapause. We exposed outbred females to different temperatures and day lengths, characterized ovarian development for over 2800 flies, and reconstructed their complete, phased genomes. We found that diapause, scored at two different developmental cutoffs, has modest heritability, and we identified hundreds of SNPs associated with each of the two phenotypes. Alleles associated with one of the diapause phenotypes tend to be more common at higher latitudes, but these alleles do not show predictable seasonal variation. The collective signal of many small-effect, clinally varying SNPs can plausibly explain latitudinal variation in diapause seen in North America. Alleles associated with diapause are segregating in Zambia, suggesting that variation in diapause relies on ancestral polymorphisms, and both pro- and anti-diapause alleles have experienced selection in North America. Finally, we utilized outdoor mesocosms to track diapause under natural conditions. We found that hybrid swarms reared outdoors evolved increased propensity for diapause in late fall, whereas indoor control populations experienced no such change. Our results indicate that diapause is a complex, quantitative trait with different evolutionary patterns across time and space.

Using targeted next-generation sequencing to characterize genetic differences associated with insecticide resistance in Culex quinquefasciatus populations from the southern U.S

Resistance to insecticides can hamper the control of mosquitoes such as Culex quinquefasciatus, known to vector arboviruses such as West Nile virus and others. The strong selective pressure exerted on a mosquito population by the use of insecticides can result in heritable genetic changes associated with resistance. We sought to characterize genetic differences between insecticide resistant and susceptible Culex quinquefasciatus mosquitoes using targeted DNA sequencing. To that end, we developed a panel of 122 genes known or hypothesized to be involved in insecticide resistance, and used an Ion Torrent PGM sequencer to sequence 125 unrelated individuals from seven populations in the southern U.S. whose resistance phenotypes to permethrin and malathion were known from previous CDC bottle bioassay testing. Data analysis consisted of discovering SNPs (Single Nucleotide Polymorphism) and genes with evidence of copy number variants (CNVs) statistically associated with resistance. Ten of the seventeen genes found to be present in higher copy numbers were experimentally validated with real-time PCR. Of those, six, including the gene with the knock-down resistance (kdr) mutation, showed evidence of a ≥ 1.5 fold increase compared to control DNA. The SNP analysis revealed 228 unique SNPs that had significant p-values for both a Fisher’s Exact Test and the Cochran-Armitage Test for Trend. We calculated the population frequency for each of the 64 nonsynonymous SNPs in this group. Several genes not previously well characterized represent potential candidates for diagnostic assays when further validation is conducted.

Diapause and quiescence: dormancy mechanisms that contribute to the geographical expansion of mosquitoes and their evolutionary success

Mosquitoes are insects belonging to the order Diptera and family Culicidae. They are distributed worldwide and include approximately 3500 species, of which about 300 have medical and veterinary importance. The evolutionary success of mosquitoes, in both tropical and temperate regions, is due to the various survival strategies these insects have developed throughout their life histories. Of the many adaptive mechanisms, diapause and quiescence, two different types of dormancy, likely contribute to the establishment, maintenance and spread of natural mosquito populations. This review seeks to objectively and coherently describe the terms diapause and quiescence, which can be confused in the literature because the phenotypic effects of these mechanisms are often similar.

Comparative Transcriptomics Reveals Key Gene Expression Differences between Diapausing and Non-Diapausing Adults of Culex pipiens

Diapause is a critical eco-physiological adaptation for winter survival in the West Nile Virus vector, Culex pipiens, but little is known about the molecular mechanisms that distinguish diapause from non-diapause in this important mosquito species. We used Illumina RNA-seq to simultaneously identify and quantify relative transcript levels in diapausing and non-diapausing adult females. Among 65,623,095 read pairs, we identified 41 genes with significantly different transcript abundances between these two groups. Transcriptome divergences between these two phenotypes include genes related to juvenile hormone synthesis, anaerobic metabolism, innate immunity and cold tolerance.

Mitotic-chromosome-based physical mapping of the Culex quinquefasciatus genome

The genome assembly of southern house mosquito Cx. quinquefasciatus is represented by a high number of supercontigs with no order or orientation on the chromosomes. Although cytogenetic maps for the polytene chromosomes of this mosquito have been developed, their utilization for the genome mapping remains difficult because of the low number of high-quality spreads in chromosome preparations. Therefore, a simple and robust mitotic-chromosome-based approach for the genome mapping of Cx. quinquefasciatus still needs to be developed. In this study, we performed physical mapping of 37 genomic supercontigs using fluorescent in situ hybridization on mitotic chromosomes from imaginal discs of 4th instar larvae. The genetic linkage map nomenclature was adopted for the chromosome numbering based on the direct positioning of 58 markers that were previously genetically mapped. The smallest, largest, and intermediate chromosomes were numbered as 1, 2, and 3, respectively. For idiogram development, we analyzed and described in detail the morphology and proportions of the mitotic chromosomes. Chromosomes were subdivided into 19 divisions and 72 bands of four different intensities. These idiograms were used for mapping the genomic supercontigs/genetic markers. We also determined the presence of length polymorphism in the q arm of sex-determining chromosome 1 in Cx. quinquefasciatus related to the size of ribosomal locus. Our physical mapping and previous genetic linkage mapping resulted in the chromosomal assignment of 13% of the total genome assembly to the chromosome bands. We provided the first detailed description, nomenclature, and idiograms for the mitotic chromosomes of Cx. quinquefasciatus. Further application of the approach developed in this study will help to improve the quality of the southern house mosquito genome.

Composite linkage map and enhanced genome map for Culex pipiens complex mosquitoes

We report here the development of 65 novel microsatellite loci and construction of a composite genetic linkage map for Culex pipiens complex mosquitoes. Microsatellites were identified by in silico screening of the Culex quinquefasciatus genome assembly. Cross-species utility of 73 microsatellites for population studies in C. pipiens sensu stricto and C. quinquefasciatus was evaluated by genotyping a subset of samples collected in Indiana, United States, and Point Fortin, Trinidad. Allele frequencies of 67 microsatellites were within Hardy-Weinberg expectations in both population subsets. A composite linkage map was constructed based on restriction fragment length polymorphism and microsatellite polymorphisms in 12 independent F1 intercross mapping populations. The composite map consists of 61 marker loci totaling 183.9 cM distributed across the 3 linkage groups. These loci cover 29.5, 88.8, and 65.6 cM on chromosomes I-III, respectively, and allow for assignment of 10.4% of the genome assembly and 13.5% of the protein coding genes to chromosome position. Our results suggest that these microsatellites will be useful for mapping and population studies of 2 pervasive species in the C. pipiens complex. Moreover, the composite map presented here will serve as a basis for the construction of high-resolution genetic and physical maps, as well as detection of quantitative trait loci to aid in the investigation of complex genetic traits influencing phenotypes of interest.

Genome-based microsatellite development in the Culex pipiens complex and comparative microsatellite frequency with Aedes aegypti and Anopheles gambiae

Background

Mosquitoes in the Culex pipiens complex are among the most medically important vectors for human disease worldwide and include major vectors for lymphatic filariasis and West Nile virus transmission. However, detailed genetic studies in the complex are limited by the number of genetic markers available. Here, we describe methods for the rapid and efficient identification and development of single locus, highly polymorphic microsatellite markers for Cx. pipiens complex mosquitoes via in silico screening of the Cx. quinquefasciatus genome sequence.

Methodology/Principal Findings

Six lab colonies representing four Cx. pipiens and two Cx. quinquefasciatus populations were utilized for preliminary assessment of 38 putative loci identified within 16 Cx. quinquefasciatus supercontig assemblies (CpipJ1) containing previously mapped genetic marker sequences. We identified and validated 12 new microsatellite markers distributed across all three linkage groups that amplify consistently among strains representing the complex. We also developed four groups of 3–5 microsatellite loci each for multiplex-ready PCR. Field collections from three cities in Indiana were used to assess the multiplex groups for their application to natural populations. All were highly polymorphic (Mean  = 13.0 alleles) per locus and reflected high polymorphism information content (PIC) (Mean  = 0.701). Pairwise F_ST indicated population structuring between Terre Haute and Fort Wayne and between Terre Haute and Indianapolis, but not between Fort Wayne and Indianapolis. In addition, we performed whole genome comparisons of microsatellite motifs and abundance between Cx. quinquefasciatus and the primary vectors for dengue virus and malaria parasites, Aedes aegypti and Anopheles gambiae, respectively.

Conclusions/Significance

We demonstrate a systematic approach for isolation and validation of microsatellites for the Cx. pipiens complex by direct screen of the Cx. quinquefasciatus genome supercontig assemblies. The genome density of microsatellites is greater in Cx. quinquefasciatus (0.26%) than in Ae. aegypti (0.14%), but considerably lower than in An. gambiae (0.77%).

Analysis of expression in the Anopheles gambiae developing testes reveals rapidly evolving lineage-specific genes in mosquitoes

Background

Male mosquitoes do not feed on blood and are not involved in delivery of pathogens to humans. Consequently, they are seldom the subjects of research, which results in a very poor understanding of their biology. To gain insights into male developmental processes we sought to identify genes transcribed exclusively in the reproductive tissues of male Anopheles gambiae pupae.

Results

Using a cDNA subtraction strategy, five male-specifically or highly male-biased expressed genes were isolated, four of which remain unannotated in the An. gambiae genome. Spatial and temporal expression patterns suggest that each of these genes is involved in the mid-late stages of spermatogenesis. Their sequences are rapidly evolving; however, two genes possess clear homologs in a wide range of taxa and one of these probably acts in a sperm motility control mechanism conserved in many organisms, including humans. The other three genes have no match to sequences from non-mosquito taxa, thus can be regarded as orphans. RNA in situ hybridization demonstrated that one of the orphans is transcribed in spermatids, which suggests its involvement in sperm maturation. Two other orphans have unknown functions. Expression analysis of orthologs of all five genes indicated that male-biased transcription was not conserved in the majority of cases in Aedes and Culex.

Conclusion

Discovery of testis-expressed orphan genes in mosquitoes opens new prospects for the development of innovative control methods. The orphan encoded proteins may represent unique targets of selective anti-mosquito sterilizing agents that will not affect non-target organisms.

Microsatellite characterization of subspecies and their hybrids in Culex pipiens complex (Diptera: Culicidae) mosquitoes along a north-south transect in the central United States

Mosquitoes in the Culex pipiens complex, Cx. p. pipiens L. and Cx. p. quinquefasciatus Say are morphologically similar and important vectors of West Nile and St. Louis Encephalitis viruses in the United States. The subspecies differ with respect to overwintering strategies, with Cx. p. pipiens entering diapause in response to winter conditions and Cx. p. quinquefasciatus lacking this ability, yet they hybridize when found in sympatry. Specimens (n = 646) were collected using gravid traps set along a transect from New Orleans, LA, to Chicago, IL. Microsatellite markers were used to: genetically characterize subspecies and hybrids, determine the degree and extent of hybridization to better define the hybrid zone, and examine the efficacy of hybrid detection between microsatellites and a single-gene assay based on the acetylcholinesterase.2 gene (HotAce.2). The results support the presence of two distinct genetic entities, with abroad, stable hybrid zone in between. Admixture analyses classified >40% of individuals as hybrids. Allelic richness was markedly different at the northern and southern ends of the transect, and there was a significant isolation by distance effect. The hybrid zone seems to be wider and extends further to the south than previous work indicated, and as a result, we propose new boundaries compared with those indicated by a previous study. Microsatellites detected more hybrids than the HotAce.2 assay, although the latter assay remains useful as a screening tool. We suggest that the parental subspecies and the hybrid zone are maintained by heterosis combined with selection for diapause at northern latitudes.