Wolbachia Blocks Viral Genome Replication Early in Infection without a Transcriptional Response by the Endosymbiont or Host Small RNA Pathways.

The intracellular endosymbiotic bacterium Wolbachia can protect insects against viral infection, and is being introduced into mosquito populations in the wild to block the transmission of arboviruses that infect humans and are a major public health concern. To investigate the mechanisms underlying this antiviral protection, we have developed a new model system combining Wolbachia-infected Drosophila melanogaster cell culture with the model mosquito-borne Semliki Forest virus (SFV; Togaviridae, Alphavirus). Wolbachia provides strong antiviral protection rapidly after infection, suggesting that an early stage post-infection is being blocked. Wolbachia does appear to have major effects on events distinct from entry, assembly or exit as it inhibits the replication of an SFV replicon transfected into the cells. Furthermore, it causes a far greater reduction in the expression of proteins from the 3' open reading frame than the 5' non-structural protein open reading frame, indicating that it is blocking the replication of viral RNA. Further to this separation of the replicase proteins and viral RNA in transreplication assays shows that uncoupling of viral RNA and replicase proteins does not overcome Wolbachia's antiviral activity. This further suggests that replicative processes are disrupted, such as translation or replication, by Wolbachia infection. This may occur by Wolbachia mounting an active antiviral response, but the virus did not cause any transcriptional response by the bacterium, suggesting that this is not the case. Host microRNAs (miRNAs) have been implicated in protection, but again we found that host cell miRNA expression was unaffected by the bacterium and neither do our findings suggest any involvement of the antiviral siRNA pathway. We conclude that Wolbachia may directly interfere with early events in virus replication such as translation of incoming viral RNA or RNA transcription, and this likely involves an intrinsic (as opposed to an induced) mechanism.

Alternative patterns of sex chromosome differentiation in Aedes aegypti (L).

BACKGROUND: Some populations of West African Aedes aegypti, the dengue and zika vector, are reproductively incompatible; our earlier study showed that divergence and rearrangements of genes on chromosome 1, which bears the sex locus (M), may be involved. We also previously described a proposed cryptic subspecies SenAae (PK10, Senegal) that had many more high inter-sex FST genes on chromosome 1 than did Ae.aegypti aegypti (Aaa, Pai Lom, Thailand). The current work more thoroughly explores the significance of those findings. RESULTS: Intersex standardized variance (FST) of single nucleotide polymorphisms (SNPs) was characterized from genomic exome capture libraries of both sexes in representative natural populations of Aaa and SenAae. Our goal was to identify SNPs that varied in frequency between males and females, and most were expected to occur on chromosome 1. Use of the assembled AaegL4 reference alleviated the previous problem of unmapped genes. Because the M locus gene nix was not captured and not present in AaegL4, the male-determining locus, per se, was not explored. Sex-associated genes were those with FST values ≥ 0.100 and/or with increased expected heterozygosity (H exp , one-sided T-test, p < 0.05) in males. There were 85 genes common to both collections with high inter-sex FST values; all genes but one were located on chromosome 1. Aaa showed the expected cluster of high inter-sex FST genes proximal to the M locus, whereas SenAae had inter-sex FST genes along the length of chromosome 1. In addition, the Aaa M-locus proximal region showed increased H exp levels in males, whereas SenAae did not. In SenAae, chromosomal rearrangements and subsequent suppressed recombination may have accelerated X-Y differentiation. CONCLUSIONS: The evidence presented here is consistent with differential evolution of proto-Y chromosomes in Aaa and SenAae.

Exome and transcriptome sequencing of Aedes aegypti identifies a locus that confers resistance to Brugia malayi and alters the immune response.

Many mosquito species are naturally polymorphic for their abilities to transmit parasites, a feature which is of great interest for controlling vector-borne disease. Aedes aegypti, the primary vector of dengue and yellow fever and a laboratory model for studying lymphatic filariasis, is genetically variable for its capacity to harbor the filarial nematode Brugia malayi. The genome of Ae. aegypti is large and repetitive, making genome resequencing difficult and expensive. We designed exome captures to target protein-coding regions of the genome, and used association mapping in a wild Kenyan population to identify a single, dominant, sex-linked locus underlying resistance. This falls in a region of the genome where a resistance locus was previously mapped in a line established in 1936, suggesting that this polymorphism has been maintained in the wild for the at least 80 years. We then crossed resistant and susceptible mosquitoes to place both alleles of the gene into a common genetic background, and used RNA-seq to measure the effect of this locus on gene expression. We found evidence for Toll, IMD, and JAK-STAT pathway activity in response to early stages of B. malayi infection when the parasites are beginning to die in the resistant genotype. We also found that resistant mosquitoes express anti-microbial peptides at the time of parasite-killing, and that this expression is suppressed in susceptible mosquitoes. Together, we have found that a single resistance locus leads to a higher immune response in resistant mosquitoes, and we identify genes in this region that may be responsible for this trait.

Latitudinal clines in gene expression and cis-regulatory element variation in Drosophila melanogaster.

BACKGROUND: Organisms can rapidly adapt to their environment when colonizing a new habitat, and this could occur by changing protein sequences or by altering patterns of gene expression. The importance of gene expression in driving local adaptation is increasingly being appreciated, and cis-regulatory elements (CREs), which control and modify the expression of the nearby genes, are predicted to play an important role. Here we investigate genetic variation in gene expression in immune-challenged Drosophila melanogaster from temperate and tropical or sub-tropical populations in Australia and United States. RESULTS: We find parallel latitudinal changes in gene expression, with genes involved in immunity, insecticide resistance, reproduction, and the response to the environment being especially likely to differ between latitudes. By measuring allele-specific gene expression (ASE), we show that cis-regulatory variation also shows parallel latitudinal differences between the two continents and contributes to the latitudinal differences in gene expression. CONCLUSIONS: Both Australia and United States were relatively recently colonized by D. melanogaster, and it was recently shown that introductions of both African and European flies occurred, with African genotypes contributing disproportionately to tropical populations. Therefore, both the demographic history of the populations and local adaptation may be causing the patterns that we see.

Estimates of allele-specific expression in Drosophila with a single genome sequence and RNA-seq data.

MOTIVATION: Genetic variation in cis-regulatory elements is an important cause of variation in gene expression. Cis-regulatory variation can be detected by using high-throughput RNA sequencing (RNA-seq) to identify differences in the expression of the two alleles of a gene. This requires that reads from the two alleles are equally likely to map to a reference genome(s), and that single-nucleotide polymorphisms (SNPs) are accurately called, so that reads derived from the different alleles can be identified. Both of these prerequisites can be achieved by sequencing the genomes of the parents of the individual being studied, but this is often prohibitively costly. RESULTS: In Drosophila, we demonstrate that biases during read mapping can be avoided by mapping reads to two alternative genomes that incorporate SNPs called from the RNA-seq data. The SNPs can be reliably called from the RNA-seq data itself, provided any variants not found in high-quality SNP databases are filtered out. Finally, we suggest a way of measuring allele-specific expression (ASE) by crossing the line of interest to a reference line with a high-quality genome sequence. Combined with our bioinformatic methods, this approach minimizes mapping biases, allows poor-quality data to be identified and removed and aides in the biological interpretation of the data as the parent of origin of each allele is known. In conclusion, our results suggest that accurate estimates of ASE do not require the parental genomes of the individual being studied to be sequenced. AVAILABILITY AND IMPLEMENTATION: Scripts used to perform our analysis are available at https://github.com/d-quinn/bio_quinn2013.

Haplotype structure and expression divergence at the Drosophila cellular immune gene eater.

The protein Eater plays an important role in microbial recognition and defensive phagocytosis in Drosophila melanogaster. We sequenced multiple alleles of the eater gene from an African and a North American population of D. melanogaster and found signatures of a partial selective sweep in North America that is localized around the second intron. This pattern is consistent with local adaptation to novel selective pressures during range expansion out of Africa. The North American sample is divided into two predominant haplotype groups, and the putatively selected haplotype is associated with a significantly higher gene expression level, suggesting that gene regulation is a possible target of selection. The eater alleles contain from 22 to 40 repeat units that are characterized by the presence of a cysteine-rich NIM motif. NIM repeats in the structural stalk of the protein exhibit concerted evolution as a function of physical location in the repeat array. Several NIM repeats within eater have previously been implicated in binding to microbial ligands, a function which in principle might subject them to special evolutionary pressures. However, we find no evidence of elevated positive selection on these pathogen-interacting units. Our study presents an instance where gene expression rather than protein structure is thought to drive the adaptive evolution of a pathogen recognition molecule in the immune system.

A robust automated system elucidates mouse home cage behavioral structure.

Patterns of behavior exhibited by mice in their home cages reflect the function and interaction of numerous behavioral and physiological systems. Detailed assessment of these patterns thus has the potential to provide a powerful tool for understanding basic aspects of behavioral regulation and their perturbation by disease processes. However, the capacity to identify and examine these patterns in terms of their discrete levels of organization across diverse behaviors has been difficult to achieve and automate. Here, we describe an automated approach for the quantitative characterization of fundamental behavioral elements and their patterns in the freely behaving mouse. We demonstrate the utility of this approach by identifying unique features of home cage behavioral structure and changes in distinct levels of behavioral organization in mice with single gene mutations altering energy balance. The robust, automated, reproducible quantification of mouse home cage behavioral structure detailed here should have wide applicability for the study of mammalian physiology, behavior, and disease.

Synergistic impairment of glucose homeostasis in ob/ob mice lacking functional serotonin 2C receptors.

To investigate how serotonin and leptin interact in the regulation of energy balance and glucose homeostasis, we generated a genetic mouse model, the OB2C mouse, which lacks functional serotonin 2C receptors and the adipocyte hormone leptin. The OB2C mice exhibited a dramatic diabetes phenotype, evidenced by a synergistic increase in serum glucose levels and water intake. The severity of the animals' diabetes phenotype would not have been predicted from the phenotypic characterization of mice bearing mutations of either the leptin (OB mutant mice) or the serotonin 2C receptor gene (2C mutant mice). The synergistic impairment in glucose homeostasis developed at an age when OB2C mice did not differ in body weight from OB mice, suggesting that this impairment was not an indirect consequence of increased adiposity. We also demonstrated that the improvement in glucose tolerance in wild-type mice treated with the serotonin releaser and reuptake inhibitor fenfluramine was blunted in 2C mutant mice. These pharmacological and genetic findings provide evidence that the serotonin 2C receptor has direct effects on glucose homeostasis.

Exon-Enriched Libraries Reveal Large Genic Differences Between Aedes aegypti from Senegal, West Africa, and Populations Outside Africa.

Aedes aegypti is one of the most studied mosquito species, and the principal vector of several arboviruses pathogenic to humans. Recently failure to oviposit, low fecundity, and poor egg-to-adult survival were observed when Ae. aegypti from Senegal (SenAae) West Africa were crossed with Ae. aegypti (Aaa) from outside of Africa, and in SenAae intercrosses. Fluorescent in situ hybridization analyses indicated rearrangements on chromosome 1, and pericentric inversions on chromosomes 2 and 3. Herein, high throughput sequencing (HTS) of exon-enriched libraries was used to compare chromosome-wide genetic diversity among Aaa collections from rural Thailand and Mexico, a sylvatic collection from southeastern Senegal (PK10), and an urban collection from western Senegal (Kaolack). Sex-specific polymorphisms were analyzed in Thailand and PK10 to assess genetic differences between sexes. Expected heterozygosity was greatest in SenAae FST distributions of 15,735 genes among all six pairwise comparisons of the four collections indicated that Mexican and Thailand collections are genetically similar, while FST distributions between PK10 and Kaolack were distinct. All four comparisons of SenAae with Aaa indicated extreme differentiation. FST was uniform between sexes across all chromosomes in Thailand, but were different, especially on the sex autosome 1, in PK10. These patterns correlate with the reproductive isolation noted earlier. We hypothesize that cryptic Ae. aegypti taxa may exist in West Africa, and the large genic differences between Aaa and SenAae detected in the present study have accumulated over a long period following the evolution of chromosome rearrangements in allopatric populations that subsequently cause reproductive isolation when these populations became sympatric.

Vector competence of Aedes aegypti mosquitoes for filarial nematodes is affected by age and nutrient limitation.

Mosquitoes are one of the most important vectors of human disease. The ability of mosquitoes to transmit disease is dependent on the age structure of the population, as mosquitoes must survive long enough for the parasites to complete their development and infect another human. Age could have additional effects due to mortality rates and vector competence changing as mosquitoes senesce, but these are comparatively poorly understood. We have investigated these factors using the mosquito Aedes aegypti and the filarial nematode Brugia malayi. Rather than observing any effects of immune senescence, we found that older mosquitoes were more resistant, but this only occurred if they had previously been maintained on a nutrient-poor diet of fructose. Constant blood feeding reversed this decline in vector competence, meaning that the number of parasites remained relatively unchanged as mosquitoes aged. Old females that had been maintained on fructose also experienced a sharp spike in mortality after an infected blood meal ("refeeding syndrome") and few survived long enough for the parasite to develop. Again, this effect was prevented by frequent blood meals. Our results indicate that old mosquitoes may be inefficient vectors due to low vector competence and high mortality, but that frequent blood meals can prevent these effects of age.

Environmental and genetic factors determine whether the mosquito Aedes aegypti lays eggs without a blood meal.

Some mosquito strains or species are able to lay eggs without taking a blood meal, a trait named autogeny. This may allow populations to persist through times or places where vertebrate hosts are scarce. Autogenous egg production is highly dependent on the environment in some species, but the ideal conditions for its expression in Aedes aegypti mosquitoes are unknown. We found that 3.2% of females in a population of Ae. aegypti from Kenya were autogenous. Autogeny was strongly influenced by temperature, with many more eggs laid at 28°C compared with 22°C. Good nutrition in larval stages and feeding on higher concentrations of sugar solution during the adult stage both result in more autogenous eggs being produced. The trait also has a genetic basis, as not all Ae. aegypti genotypes can lay autogenously. We conclude that Ae. aegypti requires a favorable environment and a suitable genotype to be able to lay eggs without a blood meal.

Assembly of the genome of the disease vector Aedes aegypti onto a genetic linkage map allows mapping of genes affecting disease transmission.

The mosquito Aedes aegypti transmits some of the most important human arboviruses, including dengue, yellow fever and chikungunya viruses. It has a large genome containing many repetitive sequences, which has resulted in the genome being poorly assembled - there are 4,758 scaffolds, few of which have been assigned to a chromosome. To allow the mapping of genes affecting disease transmission, we have improved the genome assembly by scoring a large number of SNPs in recombinant progeny from a cross between two strains of Ae. aegypti, and used these to generate a genetic map. This revealed a high rate of misassemblies in the current genome, where, for example, sequences from different chromosomes were found on the same scaffold. Once these were corrected, we were able to assign 60% of the genome sequence to chromosomes and approximately order the scaffolds along the chromosome. We found that there are very large regions of suppressed recombination around the centromeres, which can extend to as much as 47% of the chromosome. To illustrate the utility of this new genome assembly, we mapped a gene that makes Ae. aegypti resistant to the human parasite Brugia malayi, and generated a list of candidate genes that could be affecting the trait.

Providencia sneebia sp. nov. and Providencia burhodogranariea sp. nov., isolated from wild Drosophila melanogaster.

Multiple isolates of the genus Providencia were obtained from the haemolymph of wild-caught Drosophila melanogaster fruit flies. Sixteen isolates were distinguished from the six previously described species based on 16S rRNA gene sequences. These isolates belonged to two distinct groups, which we propose each comprise previously undescribed species. Two isolates, designated A(T) and B(T), were characterized by DNA sequences of the fusA, lepA, leuS, gyrB and ileS housekeeping genes, whole-genome DNA-DNA hybridizations with their nearest relatives and utilization of substrates for metabolism. The closest phylogenetic relatives of strain A(T) are strain B(T) (86.9 % identity for the housekeeping genes) and Providencia stuartii DSM 4539(T) (86.0 % identity). The closest phylogenetic relatives of strain B(T) are strain A(T) (86.9 % identity) and P. stuartii DSM 4539(T) (86.6 % identity). The type strains of described species in this genus shared between 84.1 and 90.1 % identity for these sequences. DNA-DNA hybridization between the strain pairs A(T)-B(T), A(T)-P. stuartii DSM 4539(T) and B(T)-P. stuartii DSM 4539(T) all resulted in less than 25 % relatedness. In addition, patterns of utilization of amygdalin, arbutin, aesculin, salicin, d-sorbitol, trehalose, inositol, d-adonitol and d-galactose distinguish strains A(T) and B(T) from other members of this genus. Strains A(T) and B(T) therefore represent novel species, for which the names Providencia sneebia sp. nov. (type strain A(T) =DSM 19967(T) =ATCC BAA-1589(T)) and Providencia burhodogranariea sp. nov. (type strain B(T) =DSM 19968(T) =ATCC BAA-1590(T)) are proposed.

Detection and quantification of Entomophaga maimaiga resting spores in forest soil using real-time PCR.

Environmental sampling to monitor entomopathogen titre in forest soil, a known reservoir of insect pathogens such as fungi and viruses, is important in the evaluation of conditions that could trigger epizootics and in the development of strategies for insect pest management. Molecular or PCR-based analysis of environmental samples provides a sensitive method for strain- or species-based detection, and real-time PCR, in particular, allows quantification of the organism of interest. In this study we developed a DNA extraction method and a real-time PCR assay for detection and quantification of Entomophaga maimaiga (Zygomycetes: Entomophthorales), a fungal pathogen of the gypsy moth, in the organic layer of forest soil. DNA from fungal resting spores (azygospores) in soil was extracted using a detergent and bead mill homogenization treatment followed by purification of the crude DNA extract using Sephadex-polyvinylpolypyrrolidone microcolumns. The purification step eliminated most of the environmental contaminants commonly co-extracted with genomic DNA from soil samples but detection assays still required the addition of bovine serum albumin to relieve PCR inhibition. The real-time PCR assay used primers and probe based on sequence analysis of the nuclear ribosomal ITS region of several E. maimaiga and two E. aulicae strains. Comparison of threshold cycle values from different soil samples spiked with E. maimaiga DNA showed that soil background DNA and remaining co-extracted contaminants are critical factors determining detection sensitivity. Based on our results from comparisons of resting spore titres among different forest soils, estimates were best for organic soils with comparatively high densities of resting spores.