Bradley J. Main

Adaptive introgression in an African malaria mosquito coincident with the increased usage of insecticide-treated bed nets

Significance We report that during a recent period of hybridization between two major African malaria mosquitoes, Anopheles gambiae and Anopheles coluzzii, an island of divergence on chromosome 2 introgressed from the A. gambiae into the A. coluzzii genome and its frequency subsequently increased. This introgression was coincident with the start of a major insecticide-treated bed net campaigns in Mali. These observations suggest that increased insecticide exposure acted as a selective force sufficient to drive introgression of an entire genomic island of divergence across the reproductive barrier separating these two species. This study provides a rare example of adaptive introgression in an animal species and elucidates the dynamics of how insecticide resistance evolved in A. coluzzii. Animal species adapt to changes in their environment, including man-made changes such as the introduction of insecticides, through selection for advantageous genes already present in populations or newly arisen through mutation. A possible alternative mechanism is the acquisition of adaptive genes from related species via a process known as adaptive introgression. Differing levels of insecticide resistance between two African malaria vectors, Anopheles coluzzii and Anopheles gambiae, have been attributed to assortative mating between the two species. In a previous study, we reported two bouts of hybridization observed in the town of Selinkenyi, Mali in 2002 and 2006. These hybridization events did not appear to be directly associated with insecticide-resistance genes. We demonstrate that during a brief breakdown in assortative mating in 2006, A. coluzzii inherited the entire A. gambiae-associated 2L divergence island, which includes a suite of insecticide-resistance alleles. In this case, introgression was coincident with the start of a major insecticide-treated bed net distribution campaign in Mali. This suggests that insecticide exposure altered the fitness landscape, favoring the survival of A. coluzzii/A. gambiae hybrids, and provided selection pressure that swept the 2L divergence island through A. coluzzii populations in Mali. We propose that the work described herein presents a unique description of the temporal dynamics of adaptive introgression in an animal species and represents a mechanism for the rapid evolution of insecticide resistance in this important vector of human malaria in Africa.

Regulatory Divergence in Drosophila melanogaster and D. simulans, a Genomewide Analysis of Allele-Specific Expression

Species-specific regulation of gene expression contributes to the development and maintenance of reproductive isolation and to species differences in ecologically important traits. A better understanding of the evolutionary forces that shape regulatory variation and divergence can be developed by comparing expression differences among species and interspecific hybrids. Once expression differences are identified, the underlying genetics of regulatory variation or divergence can be explored. With the goal of associating cis and/or trans components of regulatory divergence with differences in gene expression, overall and allele-specific expression levels were assayed genomewide in female adult heads of Drosophila melanogaster, D. simulans, and their F1 hybrids. A greater proportion of cis differences than trans differences were identified for genes expressed in heads and, in accordance with previous studies, cis differences also explained a larger number of species differences in overall expression level. Regulatory divergence was found to be prevalent among genes associated with defense, olfaction, and among genes downstream of the Drosophila sex determination hierarchy. In addition, two genes, with critical roles in sex determination and micro RNA processing, Sxl and loqs, were identified as misexpressed in hybrid female heads, potentially contributing to hybrid incompatibility.

Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, Anopheles gambiae

Significance Populations of the African malaria vector, Anopheles gambiae, are structured into M and S forms. All current work assumes the two rarely hybridize. Here we show this assumption is false. We demonstrate (i) significant exchange of genes between the two forms, even though (ii) hybrids have reduced fitness and (iii) the gene exchange process is spatially and temporally dynamic. For malaria, it is important to determine if genes for traits like insecticide resistance are shared between forms. For evolutionary biologists, this work confirms that this mosquito is a good model for studying if and how species may evolve in cases where there is ongoing gene flow. The M and S forms of Anopheles gambiae have been the focus of intense study by malaria researchers and evolutionary biologists interested in ecological speciation. Divergence occurs at three discrete islands in genomes that are otherwise nearly identical. An “islands of speciation” model proposes that diverged regions contain genes that are maintained by selection in the face of gene flow. An alternative “incidental island” model maintains that gene flow between M and S is effectively zero and that divergence islands are unrelated to speciation. A “divergence island SNP” assay was used to explore the spatial and temporal distributions of hybrid genotypes. Results revealed that hybrid individuals occur at frequencies ranging between 5% and 97% in every population examined. A temporal analysis revealed that assortative mating is unstable and periodically breaks down, resulting in extensive hybridization. Results suggest that hybrids suffer a fitness disadvantage, but at least some hybrid genotypes are viable. Stable introgression of the 2L speciation island occurred at one site following a hybridization event.

Allele-specific expression assays using Solexa

BackgroundAllele-specific expression (ASE) assays can be used to identify cis, trans, and cis-by-trans regulatory variation. Understanding the source of expression variation has important implications for disease susceptibility, phenotypic diversity, and adaptation. While ASE is commonly measured via relative fluorescence at a SNP, next generation sequencing provides an opportunity to measure ASE in an accurate and high-throughput manner using read counts.ResultsWe introduce a Solexa-based method to perform large numbers of ASE assays using only a single lane of a Solexa flowcell. In brief, transcripts of interest, which contain a known SNP, are PCR enriched and barcoded to enable multiplexing. Then high-throughput sequencing is used to estimate allele-specific expression using sequencing counts. To validate this method, we measured the allelic bias in a dilution series and found high correlations between measured and expected values (r>0.9, p < 0.001). We applied this method to a set of 5 genes in a Drosophila simulans parental mix, F1 and introgression and found that for these genes the majority of expression divergence can be explained by cis-regulatory variation.ConclusionWe present a new method with the capacity to measure ASE for large numbers of assays using as little as one lane of a Solexa flowcell. This will be a valuable technique for molecular and population genetic studies, as well as for verification of genome-wide data sets.

GBshape: a genome browser database for DNA shape annotations

Many regulatory mechanisms require a high degree of specificity in protein-DNA binding. Nucleotide sequence does not provide an answer to the question of why a protein binds only to a small subset of the many putative binding sites in the genome that share the same core motif. Whereas higher-order effects, such as chromatin accessibility, cooperativity and cofactors, have been described, DNA shape recently gained attention as another feature that fine-tunes the DNA binding specificities of some transcription factor families. Our Genome Browser for DNA shape annotations (GBshape; freely available at http://rohslab.cmb.usc.edu/GBshape/) provides minor groove width, propeller twist, roll, helix twist and hydroxyl radical cleavage predictions for the entire genomes of 94 organisms. Additional genomes can easily be added using the GBshape framework. GBshape can be used to visualize DNA shape annotations qualitatively in a genome browser track format, and to download quantitative values of DNA shape features as a function of genomic position at nucleotide resolution. As biological applications, we illustrate the periodicity of DNA shape features that are present in nucleosome-occupied sequences from human, fly and worm, and we demonstrate structural similarities between transcription start sites in the genomes of four Drosophila species.

The Genetic Basis of Host Preference and Resting Behavior in the Major African Malaria Vector, Anopheles arabiensis

Malaria transmission is dependent on the propensity of Anopheles mosquitoes to bite humans (anthropophily) instead of other dead end hosts. Recent increases in the usage of Long Lasting Insecticide Treated Nets (LLINs) in Africa have been associated with reductions in highly anthropophilic and endophilic vectors such as Anopheles gambiae s.s., leaving species with a broader host range, such as Anopheles arabiensis, as the most prominent remaining source of transmission in many settings. An. arabiensis appears to be more of a generalist in terms of its host choice and resting behavior, which may be due to phenotypic plasticity and/or segregating allelic variation. To investigate the genetic basis of host choice and resting behavior in An. arabiensis we sequenced the genomes of 23 human-fed and 25 cattle-fed mosquitoes collected both in-doors and out-doors in the Kilombero Valley, Tanzania. We identified a total of 4,820,851 SNPs, which were used to conduct the first genome-wide estimates of “SNP heritability” for host choice and resting behavior in this species. A genetic component was detected for host choice (human vs cow fed; permuted P = 0.002), but there was no evidence of a genetic component for resting behavior (indoors versus outside; permuted P = 0.465). A principal component analysis (PCA) segregated individuals based on genomic variation into three groups which were characterized by differences at the 2Rb and/or 3Ra paracentromeric chromosome inversions. There was a non-random distribution of cattle-fed mosquitoes between the PCA clusters, suggesting that alleles linked to the 2Rb and/or 3Ra inversions may influence host choice. Using a novel inversion genotyping assay, we detected a significant enrichment of the standard arrangement (non-inverted) of 3Ra among cattle-fed mosquitoes (N = 129) versus all non-cattle-fed individuals (N = 234; χ2, p = 0.007). Thus, tracking the frequency of the 3Ra in An. arabiensis populations may be of use to infer selection on host choice behavior within these vector populations; possibly in response to vector control. Controlled host-choice assays are needed to discern whether the observed genetic component has a direct relationship with innate host preference. A better understanding of the genetic basis for host feeding behavior in An. arabiensis may also open avenues for novel vector control strategies based on driving genes for zoophily into wild mosquito populations.

Complex genome evolution in Anopheles coluzzii associated with increased insecticide usage in Mali.

In certain cases, a species may have access to important genetic variation present in a related species via adaptive introgression. These novel alleles may interact with their new genetic background, resulting in unexpected phenotypes. In this study, we describe a selective sweep on standing variation on the X chromosome in the mosquito Anopheles coluzzii, a principal malaria vector in West Africa. This event may have been influenced by the recent adaptive introgression of the insecticide resistance gene known as kdr from the sister species Anopheles gambiae. Individuals carrying both kdr and a nearly fixed X‐linked haplotype, encompassing at least four genes including the P450 gene CYP9K1 and the cuticular protein CPR125, have rapidly increased in relative frequency. In parallel, a reproductively isolated insecticide‐susceptible A. gambiae population (Bamako form) has been driven to local extinction, likely due to strong selection from increased insecticide‐treated bed net usage.

Transcription Start Site Evolution in Drosophila

Transcription start site (TSS) evolution remains largely undescribed in Drosophila, likely due to limited annotations in non-melanogaster species. In this study, we introduce a concise new method that selectively sequences from the 5′-end of mRNA and used it to identify TSS in four Drosophila species, including Drosophila melanogaster, D. simulans, D. sechellia, and D. pseudoobscura. For verification, we compared our results in D. melanogaster with known annotations, published 5′-rapid amplification of cDNA ends data, and with RNAseq from the same mRNA pool. Then, we paired 2,849 D. melanogaster TSS with its closest equivalent TSS in each species (likely to be its true ortholog) using the available multiple sequence alignments. Most of the D. melanogaster TSSs were successfully paired with an ortholog in each species (83%, 86%, and 55% for D. simulans, D. sechellia, and D. pseudoobscura, respectively). On the basis of the number and distribution of reads mapped at each TSS, we also estimated promoter-specific expression (PSE) and TSS peak shape, respectively. Among paired TSS orthologs, the location and promoter activity were largely conserved. TSS location appears important as PSE, and TSS peak shape was more frequently divergent among TSS that had moved. Unpaired TSS were surprisingly common in D. pseudoobscura. An increased mutation rate upstream of TSS might explain this pattern. We found an enrichment of ribosomal protein genes among diverged TSS, suggesting that TSS evolution is not uniform across the genome.

Salinity adaptation and the contribution of parental environmental effects in medicago truncatula

High soil salinity negatively influences plant growth and yield. Some taxa have evolved mechanisms for avoiding or tolerating elevated soil salinity, which can be modulated by the environment experienced by parents or offspring. We tested the contribution of the parental and offspring environments on salinity adaptation and their potential underlying mechanisms. In a two-generation greenhouse experiment, we factorially manipulated salinity concentrations for genotypes of Medicago truncatula that were originally collected from natural populations that differed in soil salinity. To compare population level adaptation to soil salinity and to test the potential mechanisms involved we measured two aspects of plant performance, reproduction and vegetative biomass, and phenological and physiological traits associated with salinity avoidance and tolerance. Saline-origin populations had greater biomass and reproduction under saline conditions than non-saline populations, consistent with local adaptation to saline soils. Additionally, parental environmental exposure to salt increased this difference in performance. In terms of environmental effects on mechanisms of salinity adaptation, parental exposure to salt spurred phenological differences that facilitated salt avoidance, while offspring exposure to salt resulted in traits associated with greater salt tolerance. Non-saline origin populations expressed traits associated with greater growth in the absence of salt while, for saline adapted populations, the ability to maintain greater performance in saline environments was also associated with lower growth potential in the absence of salt. Plastic responses induced by parental and offspring environments in phenology, leaf traits, and gas exchange contribute to salinity adaptation in M. truncatula. The ability of plants to tolerate environmental stress, such as high soil salinity, is likely modulated by a combination of parental effects and within-generation phenotypic plasticity, which are likely to vary in populations from contrasting environments.

Genetic variation associated with increased insecticide resistance in the malaria mosquito, Anopheles coluzzii.

BackgroundMalaria mortality rates in sub-Saharan Africa have declined significantly in recent years as a result of increased insecticide-treated bed net (ITN) usage. A major challenge to further progress is the emergence and spread of insecticide resistance alleles in the Anopheles mosquito vectors, like An. coluzzii. A non-synonymous mutation in the para voltage-gated sodium channel gene reduces pyrethroid-binding affinity, resulting in knockdown resistance (kdr). Metabolic mechanisms of insecticide resistance involving detoxification genes like cytochrome P450 genes, carboxylesterases, and glutathione S-transferases are also important. As some gene activity is tissue-specific and/or environmentally induced, gene regulatory variation may be overlooked when comparing expression from whole mosquito bodies under standard rearing conditions.ResultsWe detected complex insecticide resistance in a 2014 An. coluzzii colony from southern Mali using bottle bioassays. Additional bioassays involving recombinant genotypes from a cross with a relatively susceptible 1995 An. coluzzii colony from Mali confirmed the importance of kdr and associated increased permethrin resistance to the CYP9K1 locus on the X chromosome. Significant differential expression of CYP9K1 was not observed among these colonies in Malpighian tubules. However, the P450 gene CYP6Z1 was overexpressed in resistant individuals following sublethal permethrin exposure and the carboxylesterase gene COEAE5G was constitutively overexpressed.ConclusionsThe significant P450-related insecticide resistance observed in the 2014 An. coluzzii colony indicates that ITNs treated with the P450 inhibitor piperonyl butoxide (PBO) would be more effective in this region. The known insecticide resistance gene CYP6Z1 was differentially expressed exclusively in the context of sublethal permethrin exposure, highlighting the importance of tissue-specificity and environmental conditions in gene expression studies. The increased activity of the carboxylesterase COEAE5G in the resistant An. coluzzii colony suggests resistance to other insecticides like organophosphates. Additional gene expression studies involving other tissues (e.g. fat body) would provide a more comprehensive view of genes underlying metabolic insecticide resistance in An. coluzzii from Mali. Identifying genetic markers linked to these regulatory alleles is an important next step that would substantially improve insecticide resistance surveillance and population genetic studies in this important vector species.