The duplicated P450s CYP6P9a/b drive carbamates and pyrethroids cross-resistance in the major African malaria vector Anopheles funestus.

Cross-resistance to insecticides in multiple resistant malaria vectors is hampering resistance management. Understanding its underlying molecular basis is critical to implementation of suitable insecticide-based interventions. Here, we established that the tandemly duplicated cytochrome P450s, CYP6P9a/b are driving carbamate and pyrethroid cross-resistance in Southern African populations of the major malaria vector Anopheles funestus. Transcriptome sequencing revealed that cytochrome P450s are the most over-expressed genes in bendiocarb and permethrin-resistant An. funestus. The CYP6P9a and CYP6P9b genes are overexpressed in resistant An. funestus from Southern Africa (Malawi) versus susceptible An. funestus (Fold change (FC) is 53.4 and 17 respectively), while the CYP6P4a and CYP6P4b genes are overexpressed in resistant An. funestus in Ghana, West Africa, (FC is 41.1 and 17.2 respectively). Other up-regulated genes in resistant An. funestus include several additional cytochrome P450s (e.g. CYP9J5, CYP6P2, CYP6P5), glutathione-S transferases, ATP-binding cassette transporters, digestive enzymes, microRNA and transcription factors (FC<7). Targeted enrichment sequencing strongly linked a known major pyrethroid resistance locus (rp1) to carbamate resistance centering around CYP6P9a/b. In bendiocarb resistant An. funestus, this locus exhibits a reduced nucleotide diversity, significant p-values when comparing allele frequencies, and the most non-synonymous substitutions. Recombinant enzyme metabolism assays showed that both CYP6P9a/b metabolize carbamates. Transgenic expression of CYP6P9a/b in Drosophila melanogaster revealed that flies expressing both genes were significantly more resistant to carbamates than controls. Furthermore, a strong correlation was observed between carbamate resistance and CYP6P9a genotypes with homozygote resistant An. funestus (CYP6P9a and the 6.5kb enhancer structural variant) exhibiting a greater ability to withstand bendiocarb/propoxur exposure than homozygote CYP6P9a_susceptible (e.g Odds ratio = 20.8, P<0.0001 for bendiocarb) and heterozygotes (OR = 9.7, P<0.0001). Double homozygote resistant genotype (RR/RR) were even more able to survive than any other genotype combination showing an additive effect. This study highlights the risk that pyrethroid resistance escalation poses to the efficacy of other classes of insecticides. Available metabolic resistance DNA-based diagnostic assays should be used by control programs to monitor cross-resistance between insecticides before implementing new interventions.

Overexpression and nonsynonymous mutations of UDP-glycosyltransferases are potentially associated with pyrethroid resistance in Anopheles funestus.

UDP-glycosyltransferases (UGTs) enzymes are pivotal in insecticide resistance by transforming hydrophobic substrates into more hydrophilic forms for efficient cell elimination. This study provides the first comprehensive investigation of Anopheles funestus UGT genes, their evolution, and their association with pyrethroid resistance. We employed a genome-wide association study using pooled sequencing (GWAS-PoolSeq) and transcriptomics on pyrethroid-resistant An. funestus, along with deep-targeted sequencing of UGTs in 80 mosquitoes Africa-wide. UGT310B2 was consistently overexpressed Africa-wide and significant gene-wise Fst differentiation was observed between resistant and susceptible populations: UGT301C2 and UGT302A3 in Malawi, and UGT306C2 in Uganda. Additionally, nonsynonymous mutations in UGT genes were identified. Gene-wise Tajima's D density curves provide insights into population structures within populations across these countries, supporting previous observations. These findings have important implications for current An. funestus control strategies facilitating the prediction of cross-resistance to other UGT-metabolised polar insecticides, thereby guiding more effective and targeted insecticide resistance management efforts.

Molecular drivers of insecticide resistance in the Sahelo-Sudanian populations of a major malaria vector Anopheles coluzzii.

BACKGROUND: Information on common markers of metabolic resistance in malaria vectors from countries sharing similar eco-climatic characteristics can facilitate coordination of malaria control. Here, we characterized populations of the major malaria vector Anopheles coluzzii from Sahel region, spanning four sub-Saharan African countries: Nigeria, Niger, Chad and Cameroon. RESULTS: Genome-wide transcriptional analysis identified major genes previously implicated in pyrethroid and/or cross-resistance to other insecticides, overexpressed across the Sahel, including CYP450s, glutathione S-transferases, carboxylesterases and cuticular proteins. Several, well-known markers of insecticide resistance were found in high frequencies-including in the voltage-gated sodium channel (V402L, I940T, L995F, I1527T and N1570Y), the acetylcholinesterase-1 gene (G280S) and the CYP4J5-L43F (which is fixed). High frequencies of the epidemiologically important chromosomal inversion polymorphisms, 2La, 2Rb and 2Rc, were observed (~80% for 2Rb and 2Rc). The 2La alternative arrangement is fixed across the Sahel. Low frequencies of these inversions (<10%) were observed in the fully insecticide susceptible laboratory colony of An. coluzzii (Ngoussou). Several of the most commonly overexpressed metabolic resistance genes sit in these three inversions. Two commonly overexpressed genes, GSTe2 and CYP6Z2, were functionally validated. Transgenic Drosophila melanogaster flies expressing GSTe2 exhibited extremely high DDT and permethrin resistance (mortalities <10% in 24h). Serial deletion of the 5' intergenic region, to identify putative nucleotide(s) associated with GSTe2 overexpression, revealed that simultaneous insertion of adenine nucleotide and a transition (T->C), between Forkhead box L1 and c-EST putative binding sites, were responsible for the high overexpression of GSTe2 in the resistant mosquitoes. Transgenic flies expressing CYP6Z2 exhibited marginal resistance towards 3-phenoxybenzylalcohol (a primary product of pyrethroid hydrolysis by carboxylesterases) and a type II pyrethroid, α-cypermethrin. However, significantly higher mortalities were observed in CYP6Z2 transgenic flies compared with controls, on exposure to the neonicotinoid, clothianidin. This suggests a possible bioactivation of clothianidin into a toxic intermediate, which may make it an ideal insecticide against populations of An. coluzzii overexpressing this P450. CONCLUSIONS: These findings will facilitate regional collaborations within the Sahel region and refine implementation strategies through re-focusing interventions, improving evidence-based, cross-border policies towards local and regional malaria pre-elimination.

An Africa-wide genomic evolution of insecticide resistance in the malaria vector Anopheles funestus involves selective sweeps, copy number variations, gene conversion and transposons.

Insecticide resistance in malaria vectors threatens to reverse recent gains in malaria control. Deciphering patterns of gene flow and resistance evolution in malaria vectors is crucial to improving control strategies and preventing malaria resurgence. A genome-wide survey of Anopheles funestus genetic diversity Africa-wide revealed evidences of a major division between southern Africa and elsewhere, associated with different population histories. Three genomic regions exhibited strong signatures of selective sweeps, each spanning major resistance loci (CYP6P9a/b, GSTe2 and CYP9K1). However, a sharp regional contrast was observed between populations correlating with gene flow barriers. Signatures of complex molecular evolution of resistance were detected with evidence of copy number variation, transposon insertion and a gene conversion between CYP6P9a/b paralog genes. Temporal analyses of samples before and after bed net scale up suggest that these genomic changes are driven by this control intervention. Multiple independent selective sweeps at the same locus in different parts of Africa suggests that local evolution of resistance in malaria vectors may be a greater threat than trans-regional spread of resistance haplotypes.

Daphnia magna egg piRNA cluster expression profiles change as mothers age.

BACKGROUND: PiRNAs prevent transposable elements wreaking havoc on the germline genome. Changes in piRNA expression over the lifetime of an individual may impact on ageing through continued suppression, or release, of transposable element expression. We identified piRNA producing clusters in the genome of Daphnia magna by a combination of bioinformatic methods, and then contrasted their expression between parthenogenetically produced eggs representing maternally-deposited germline piRNAs of young (having their 1st clutch) and old (having their 5th clutch) mothers. Results from eggs were compared to cluster expression in three generations of adults. RESULTS: As for other arthropods, D. magna encodes long uni-directionally transcribed non-coding RNAs consisting of fragmented transposable elements which account for most piRNAs expressed. Egg tissues showed extensive differences between clutches from young mothers and those from old mothers, with 578 and 686 piRNA clusters upregulated, respectively. Most log fold-change differences for significant clusters were modest, however. When considering only highly expressed clusters, there was a bias towards 1st clutch eggs at 41 upregulated versus eight clusters in the eggs from older mothers. F0 generation differences between young and old mothers were fewer than eggs, as 179 clusters were up-regulated in young versus 170 old mothers. This dropped to 31 versus 22 piRNA clusters when comparing adults in the F1 generation, and no differences were detected in the F3 generation. Inter-generational losses of differential piRNA cluster were similar to that observed for D. magna micro-RNA expression. CONCLUSIONS: Little overlap in differentially expressed clusters was found between adults containing mixed somatic and germline (ovary) tissues and germ-line representing eggs. A cluster encompassing a Tudor domain containing gene important in the piRNA pathway was upregulated in the eggs from old mothers. We hypothesise that regulation of this gene could form part of a feedback loop that reduces piRNA pathway activity explaining the reduced number of highly-expressed clusters in eggs from old mothers.

Multi-omics analysis identifies a CYP9K1 haplotype conferring pyrethroid resistance in the malaria vector Anopheles funestus in East Africa.

Metabolic resistance to pyrethroids is a menace to the continued effectiveness of malaria vector controls. Its molecular basis is complex and varies geographically across Africa. Here, we used a multi-omics approach, followed-up with functional validation to show that a directionally selected haplotype of a cytochrome P450, CYP9K1 is a major driver of resistance in Anopheles funestus. A PoolSeq GWAS using mosquitoes alive and dead after permethrin exposure, from Malawi and Cameroon, detected candidate genomic regions, but lacked consistency across replicates. Targeted sequencing of candidate resistance genes detected several SNPs associated with known pyrethroid resistance QTLs. The most significant SNPs were in the cytochrome P450 CYP304B1 (Cameroon), CYP315A1 (Uganda) and the ABC transporter gene ABCG4 (Malawi). However, when comparing field resistant mosquitoes to laboratory susceptible, the pyrethroid resistance locus rp1 and SNPs around the ABC transporter ABCG4 were consistently significant, except for Uganda where SNPs in the P450 CYP9K1 was markedly significant. In vitro heterologous metabolism assays with recombinant CYP9K1 revealed that it metabolises type II pyrethroid (deltamethrin; 64% depletion) but not type I (permethrin; 0%), while moderately metabolising DDT (17%). CYP9K1 exhibited reduced genetic diversity in Uganda underlying an extensive selective sweep. Furthermore, a glycine to alanine (G454A) amino acid change in CYP9K1 was fixed in Ugandan mosquitoes but not in other An. funestus populations. This study sheds further light on the evolution of metabolic resistance in a major malaria vector by implicating more genes and variants that can be used to design field-applicable markers to better track resistance Africa-wide.

Identification of Parachlamydiaceae DNA in nasal and rectal passages of healthy dairy cattle.

AIMS: The order Chlamydiales comprises a broad range of bacterial pathogens and endosymbionts, which infect a wide variety of host species. Within this order, members of the family Parachlamydiaceae, which includes Parachlamydia and Neochlamydia species, have been particularly associated with infections in both humans and cattle, including having a potential pathogenic role in cases of bovine abortion. While the route of transmission has yet to be defined, it has been hypothesised that asymptomatic carriage and contamination of the immediate environment may be a route of inter-animal transmission. We investigated the asymptomatic carriage of Chlamydia-related organisms in healthy cattle. METHODS & RESULTS: DNA was isolated from nasal and rectal swabs obtained from 38 healthy dairy heifers. A Chlamydiales sp. 16S rRNA qPCR was performed on each sample. A total of 18/38 nasal samples and all 38/38 rectal samples were identified as positive for Chlamydiales sp. Each positive sample was sequenced confirming the presence of DNA belonging to the Parachlamydiaceae. CONCLUSIONS: The presence of Parachlamydiaceae DNA in nasal and rectal swab samples of healthy cattle provides evidence for the asymptomatic carriage of parachlamydial organisms within cattle. SIGNIFICANCE & IMPACT OF THE STUDY: The study provides evidence of potential routes of environmental contamination that could provide a route for inter-animal and animal transmission of Parachlamydiaceae.

A novel, magnetic bead-based extraction method for the isolation of antimicrobial resistance genes with a case study in river water in Malawi.

AIMS: The environment is increasingly recognized as an important reservoir of antimicrobial resistance genes (ARGs), which can be identified using molecular platforms. Yet, environmental surveillance remains an underutilised tool as there is no agreement on the best strategy for sample processing. We aim to develop a low-cost extraction method independent to commercial kits or reagents. METHODS AND RESULTS: We present a novel, magnetic bead-based method for the isolation of ARGs from river water named MagnaExtract. We present this with analytic limit of detection as well as a case study in Southern Malawi. Here we compare the DNA yield from MagnaExtract with commercially available QIAGEN kits and the crude boil and spin method, using a high-resolution melt analysis PCR panel designed for the detection of third-generation cephalosporin and carbapenem-resistant genes from 98 water samples. CONCLUSION: The MagnaExtract method is comparable, and in some instance's superior to commercially available kits for the isolation of ARGs from river water samples. SIGNIFICANCE AND IMPACT OF THE STUDY: The MagnaExtract approach offers a simple, affordable, high yielding extraction method that could be used for the detection of ARGs from river water samples in surveillance campaigns in East Africa.

Genomic dissection of an extended phenotype: Oak galling by a cynipid gall wasp.

Galls are plant tissues whose development is induced by another organism for the inducer's benefit. 30,000 arthropod species induce galls, and in most cases the inducing effectors and target plant systems are unknown. Cynipid gall wasps are a speciose monophyletic radiation that induce structurally complex galls on oaks and other plants. We used a model system comprising the gall wasp Biorhiza pallida and the oak Quercus robur to characterise inducer and host plant gene expression at defined stages through the development of galled and ungalled plant tissues, and tested alternative hypotheses for the origin and type of galling effectors and plant metabolic pathways involved. Oak gene expression patterns diverged markedly during development of galled and normal buds. Young galls showed elevated expression of oak genes similar to legume root nodule Nod factor-induced early nodulin (ENOD) genes and developmental parallels with oak buds. In contrast, mature galls showed substantially different patterns of gene expression to mature leaves. While most oak transcripts could be functionally annotated, many gall wasp transcripts of interest were novel. We found no evidence in the gall wasp for involvement of third-party symbionts in gall induction, for effector delivery using virus-like-particles, or for gallwasp expression of genes coding for plant hormones. Many differentially and highly expressed genes in young larvae encoded secretory peptides, which we hypothesise are effector proteins exported to plant tissues. Specifically, we propose that host arabinogalactan proteins and gall wasp chitinases interact in young galls to generate a somatic embryogenesis-like process in oak tissues surrounding the gall wasp larvae. Gall wasp larvae also expressed genes encoding multiple plant cell wall degrading enzymes (PCWDEs). These have functional orthologues in other gall inducing cynipids but not in figitid parasitoid sister groups, suggesting that they may be evolutionary innovations associated with cynipid gall induction.

Pyrethroid resistance in the New World malaria vector Anopheles albimanus is mediated by cytochrome P450 CYP6P5.

Pyrethroid resistance in the malaria vector Anopheles albimanus presents an obstacle to malaria elimination in the Americas. Here, An. albimanus CYP6P5 (the most overexpressed P450 in a Peruvian population) was functionally characterized. Recombinant CYP6P5 metabolized the type II pyrethroids, deltamethrin and α-cypermethrin with comparable affinities (KM of 3.3 µM ± 0.4 and 3.6 µM ± 0.5, respectively), but exhibited a 2.7-fold higher catalytic rate for α-cypermethrin (kcat of 6.02 min-1 ± 0.2) versus deltamethrin (2.68 min-1 ± 0.09). Time-course assays revealed progressive depletion of the above pyrethroids with production of four HPLC-detectable metabolites. Low depletion was obtained with type I pyrethroid, permethrin. Transgenic expression in Drosophila melanogaster demonstrated that overexpression of CYP6P5 alone conferred type II pyrethroid resistance, with only 16% and 55.3% mortalities in flies exposed to 0.25% α-cypermethrin and 0.15% deltamethrin, respectively. Synergist bioassays using P450 inhibitor piperonylbutoxide significantly recovered susceptibility (mortality = 73.6%, p < 0.001) in synergized flies exposed to 4% piperonylbutoxide, plus 0.25% α-cypermethrin, compared to non-synergized flies (mortality = 4.9%). Moderate resistance was also observed towards 4% DDT. These findings established the preeminent role of CYP6P5 in metabolic resistance in An. albimanus, highlighting challenges associated with deployment of insecticide-based control tools in the Americas.

Whole-genome data reveal the complex history of a diverse ecological community.

How widespread ecological communities assemble remains a key question in ecology. Trophic interactions between widespread species may reflect a shared population history or ecological fitting of local pools of species with very different population histories. Which scenario applies is central to the stability of trophic associations and the potential for coevolution between species. Here we show how alternative community assembly hypotheses can be discriminated using whole-genome data for component species and provide a likelihood framework that overcomes current limitations in formal comparison of multispecies histories. We illustrate our approach by inferring the assembly history of a Western Palearctic community of insect herbivores and parasitoid natural enemies, trophic groups that together comprise 50% of terrestrial species. We reject models of codispersal from a shared origin and of delayed enemy pursuit of their herbivore hosts, arguing against herbivore attainment of "enemy-free space." The community-wide distribution of species expansion times is also incompatible with a random, neutral model of assembly. Instead, we reveal a complex assembly history of single- and multispecies range expansions through the Pleistocene from different directions and over a range of timescales. Our results suggest substantial turnover in species associations and argue against tight coevolution in this system. The approach we illustrate is widely applicable to natural communities of nonmodel species and makes it possible to reveal the historical backdrop against which natural selection acts.

Daphnia magna modifies its gene expression extensively in response to caloric restriction revealing a novel effect on haemoglobin isoform preference.

Caloric restriction (CR) produces clear phenotypic effects within and between generations of the model crustacean Daphnia magna. We have previously established that micro-RNAs and cytosine methylation change in response to CR in this organism, and we demonstrate here that CR has a dramatic effect on gene expression. Over 6,000 genes were differentially expressed between CR and well-fed D. magna, with a bias towards up-regulation of genes under caloric restriction. We identified a highly expressed haemoglobin gene that responds to CR by changing isoform proportions. Specifically, a transcript containing three haem-binding erythrocruorin domains was strongly down-regulated under CR in favour of transcripts containing fewer or no such domains. This change in the haemoglobin mix is similar to the response to hypoxia in Daphnia, which is mediated through the transcription factor hypoxia-inducible factor 1, and ultimately the mTOR signalling pathway. This is the first report of a role for haemoglobin in the response to CR. We also observed high absolute expression of superoxide dismutase (SOD) in normally fed individuals, which contrasts with observations of high SOD levels under CR in other taxa. However, key differentially expressed genes, like SOD, were not targeted by differentially expressed micro-RNAs. Whether the link between haemoglobin and CR occurs in other organisms, or is related to the aquatic lifestyle, remains to be tested. It suggests that one response to CR may be to simply transport less oxygen and lower respiration.

Low-coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp.

Population divergence and gene flow are key processes in evolution and ecology. Model-based analysis of genome-wide data sets allows discrimination between alternative scenarios for these processes even in nonmodel taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum [bSFS], the second on the pairwise sequentially Markovian coalescent [PSMC]) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia and are separated into The Northern population is to the north of the Southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of postdivergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca. 196 kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca. 57 kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome-level data for small population samples of tiny, nonmodel organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.

A 6.5-kb intergenic structural variation enhances P450-mediated resistance to pyrethroids in malaria vectors lowering bed net efficacy.

Elucidating the complex evolutionary armory that mosquitoes deploy against insecticides is crucial to maintain the effectiveness of insecticide-based interventions. Here, we deciphered the role of a 6.5-kb structural variation (SV) in driving cytochrome P450-mediated pyrethroid resistance in the malaria vector, Anopheles funestus. Whole-genome pooled sequencing detected an intergenic 6.5-kb SV between duplicated CYP6P9a/b P450s in pyrethroid-resistant mosquitoes through a translocation event. Promoter analysis revealed a 17.5-fold higher activity (p < .0001) for the SV- carrying fragment than the SV- free one. Quantitative real-time PCR expression profiling of CYP6P9a/b for each SV genotype supported its role as an enhancer because SV+/SV+ homozygote mosquitoes had a significantly greater expression for both genes than heterozygotes SV+/SV- (1.7- to 2-fold) and homozygotes SV-/SV- (4-to 5-fold). Designing a PCR assay revealed a strong association between this SV and pyrethroid resistance (SV+/SV+ vs. SV-/SV-; odds ratio [OR] = 2,079.4, p < .001). The 6.5-kb SV is present at high frequency in southern Africa (80%-100%) but absent in East/Central/West Africa. Experimental hut trials revealed that homozygote SV mosquitoes had a significantly greater chance to survive exposure to pyrethroid-treated nets (OR 27.7; p < .0001) and to blood feed than susceptible mosquitoes. Furthermore, mosquitoes homozygote-resistant at the three loci (SV+/CYP6P9a_R/CYP6P9b_R) exhibited a higher resistance level, leading to a far superior ability to survive exposure to nets than those homozygotes susceptible at the three loci, revealing a strong additive effect. This study highlights the important role of structural variations in the development of insecticide resistance in malaria vectors and their detrimental impact on the effectiveness of pyrethroid-based nets.

Genome-wide methylation is modified by caloric restriction in Daphnia magna.

BACKGROUND: The degradation of epigenetic control with age is associated with progressive diseases of ageing, including cancers, immunodeficiency and diabetes. Reduced caloric intake slows the effects of ageing and age-related disease in vertebrates and invertebrates, a process potentially mediated by the impact of caloric restriction on epigenetic factors such as DNA methylation. We used whole genome bisulphite sequencing to study how DNA methylation patterns change with diet in a small invertebrate, the crustacean Daphnia magna. Daphnia show the classic response of longer life under caloric restriction (CR), and they reproduce clonally, which permits the study of epigenetic changes in the absence of genetic variation. RESULTS: Global cytosine followed by guanine (CpG) methylation was 0.7-0.9%, and there was no difference in overall methylation levels between normal and calorie restricted replicates. However, 333 differentially methylated regions (DMRs) were evident between the normally fed and CR replicates post-filtering. Of these 65% were hypomethylated in the CR group, and 35% were hypermethylated in the CR group. CONCLUSIONS: Our results demonstrate an effect of CR on the genome-wide methylation profile. This adds to a growing body of research in Daphnia magna that demonstrate an epigenomic response to environmental stimuli. Specifically, gene Ontology (GO) term enrichment of genes associated with hyper and hypo-methylated DMRs showed significant enrichment for methylation and acyl-CoA dehydrogenase activity, which are linked to current understanding of their roles in CR in invertebrate model organisms.

Daphnia magna microRNAs respond to nutritional stress and ageing but are not transgenerational.

Maternal effects, where the performance of offspring is determined by the condition of their mother, are widespread and may in some cases be adaptive. The crustacean Daphnia magna shows strong maternal effects: offspring size at birth and other proxies for fitness are altered when their mothers are older or when mothers have experienced dietary restriction. The mechanisms for this transgenerational transmission of maternal experience are unknown, but could include changes in epigenetic patterning. MicroRNAs (miRNAs) are regulators of gene expression that have been shown to play roles in intergenerational information transfer, and here, we test whether miRNAs are involved in D. magna maternal effects. We found that miRNAs were differentially expressed in mothers of different ages or nutritional state. We then examined miRNA expression in their eggs, their adult daughters and great granddaughters, which did not experience any treatments. The maternal (treatment) generation exhibited differential expression of miRNAs, as did their eggs, but this was reduced in adult daughters and lost by great granddaughters. Thus, miRNAs are a component of maternal provisioning, but do not appear to be the cause of transgenerational responses under these experimental conditions. MicroRNAs may act in tandem with egg provisioning (e.g., with carbohydrates or fats), and possibly other small RNAs or epigenetic modifications.

Likelihood-based inference of population history from low-coverage de novo genome assemblies.

Short-read sequencing technologies have in principle made it feasible to draw detailed inferences about the recent history of any organism. In practice, however, this remains challenging due to the difficulty of genome assembly in most organisms and the lack of statistical methods powerful enough to discriminate between recent, nonequilibrium histories. We address both the assembly and inference challenges. We develop a bioinformatic pipeline for generating outgroup-rooted alignments of orthologous sequence blocks from de novo low-coverage short-read data for a small number of genomes, and show how such sequence blocks can be used to fit explicit models of population divergence and admixture in a likelihood framework. To illustrate our approach, we reconstruct the Pleistocene history of an oak-feeding insect (the oak gallwasp Biorhiza pallida), which, in common with many other taxa, was restricted during Pleistocene ice ages to a longitudinal series of southern refugia spanning the Western Palaearctic. Our analysis of sequence blocks sampled from a single genome from each of three major glacial refugia reveals support for an unexpected history dominated by recent admixture. Despite the fact that 80% of the genome is affected by admixture during the last glacial cycle, we are able to infer the deeper divergence history of these populations. These inferences are robust to variation in block length, mutation model and the sampling location of individual genomes within refugia. This combination of de novo assembly and numerical likelihood calculation provides a powerful framework for estimating recent population history that can be applied to any organism without the need for prior genetic resources.

ABC inference of multi-population divergence with admixture from unphased population genomic data.

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome-scale data from population-level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three-population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation-based cross-validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.

Exploring the molecular mechanisms of increased intensity of pyrethroid resistance in Central African population of a major malaria vector Anopheles coluzzii.

Molecular mechanisms driving the escalation of pyrethroid resistance in the major malaria mosquitoes of Central Africa remain largely uncharacterized, hindering effective management strategies. Here, resistance intensity and the molecular mechanisms driving it were investigated in a population of Anopheles coluzzii from northern Cameroon. High levels of pyrethroid and organochloride resistance were observed in An. coluzzii population, with no mortality for 1× permethrin; only 11% and 33% mortalities for 5× and 10× permethrin diagnostic concentrations, and <2% mortalities for deltamethrin and DDT, respectively. Moderate bendiocarb resistance (88% mortality) and full susceptibility to malathion were observed. Synergist bioassays with piperonyl butoxide recovered permethrin susceptibility, with mortalities increasing to 53.39%, and 87.30% for 5× and 10× permethrin, respectively, implicating P450 monooxygenases. Synergist bioassays with diethyl maleate (DEM) recovered permethrin and DDT susceptibilities (mortalities increasing to 34.75% and 14.88%, respectively), implicating glutathione S-transferases. RNA-seq-based genome-wide transcriptional analyses supported by quantitative PCR identified glutathione S-transferase, GSTe2 (RNA-seqFC = 2.93 and qRT-PCRFC = 8.4, p < 0.0043) and CYP450, CYP6Z2 (RNA-seqFC = 2.39 and qRT-PCRFC = 11.7, p < 0.0177) as the most overexpressed detoxification genes in the pyrethroid-resistant mosquitoes, compared to mosquitoes of the susceptible Ngousso colony. Other overexpressed genes include P450s, CYP6M2 (FC = 1.68, p < 0.0114), CYP4G16 (FC = 2.02, p < 0.0005), and CYP4G17 (FC = 1.86, p < 0.0276). While high frequency of the 1014F kdr mutation (50%) and low frequencies of 1014S (6.61%) and 1575Y (10.29%) were observed, no ace-1 mutation was detected in bendiocarb-resistant populations, suggesting the preeminent role of metabolic mechanism. Overexpression of metabolic resistance genes (including GSTe2 and CYP6Z2 known to confer resistance to multiple insecticides) in An. coluzzii from the Sudan Savannah of Cameroon highlights the need for alternative management strategies to reduce malaria burden in northern Cameroon.

A new genus of oak gallwasp, Cyclocynips Melika, Tang & Sinclair (Hymenoptera: Cynipidae: Cynipini), with descriptions of two new species from Taiwan.

A new genus of cynipid oak gallwasp-Cyclocynips Melika, Tang, & Sinclair (Hymenoptera: Cynipidae: Cynipini), with two new species--C. uberis and C. tumorvirgae--reared from galls on oaks of the Quercus subgenus Cyclobalanopsis is described from Taiwan. Descriptions of asexual generation adults and their diagnostic characters are presented. The likelihood of yet undiscovered sexual generations and the evolution of host-plant associations in these species are discussed.