Positional cloning of rp2 QTL associates the P450 genes CYP6Z1, CYP6Z3 and CYP6M7 with pyrethroid resistance in the malaria vector Anopheles funestus

Pyrethroid resistance in Anopheles funestus is threatening malaria control in Africa. Elucidation of underlying resistance mechanisms is crucial to improve the success of future control programs. A positional cloning approach was used to identify genes conferring resistance in the uncharacterised rp2 quantitative trait locus (QTL) previously detected in this vector using F6 advanced intercross lines (AIL). A 113 kb BAC clone spanning rp2 was identified and sequenced revealing a cluster of 15 P450 genes and one salivary protein gene (SG7-2). Contrary to A. gambiae, AfCYP6M1 is triplicated in A. funestus, while AgCYP6Z2 orthologue is absent. Five hundred and sixty-five new single nucleotide polymorphisms (SNPs)were identified for genetic mapping from rp2 P450s and other genes revealing high genetic polymorphisms with one SNP every 36 bp. A significant genotype/phenotype association was detected for rp2 P450s but not for a cluster of cuticular protein genes previously associated with resistance in A. gambiae. QTL mapping using F6 AIL confirms the rp2 QTL with an increase logarithm of odds score of 5. Multiplex gene expression profiling of 15 P450s and other genes around rp2 followed by individual validation using qRT–PCR indicated a significant overexpression in the resistant FUMOZ-R strain of the P450s AfCYP6Z1, AfCYP6Z3, AfCYP6M7 and the glutathione-s-transferase GSTe2 with respective fold change of 11.2,6.3, 5.5 and 2.8. Polymorphisms analysis of AfCYP6Z1 and AfCYP6Z3 identified amino acid changes potentially associated with resistance further indicating that these genes are controlling the pyrethroid resistance explained by the rp2 QTL. The characterisation of this rp2 QTL significantly improves our understanding of resistance mechanisms in A. funestus

Cross-protection against European swine influenza viruses in the context of infection immunity against the 2009 pandemic H1N1 virus : studies in the pig model of influenza

Pigs are natural hosts for the same influenza virus subtypes as humans and are a valuable model for cross-protection studies with influenza. In this study, we have used the pig model to examine the extent of virological protection between a) the 2009 pandemic H1N1 (pH1N1) virus and three different European H1 swine influenza virus (SIV) lineages, and b) these H1 viruses and a European H3N2 SIV. Pigs were inoculated intranasally with representative strains of each virus lineage with 6- and 17-week intervals between H1 inoculations and between H1 and H3 inoculations, respectively. Virus titers in nasal swabs and/or tissues of the respiratory tract were determined after each inoculation. There was substantial though differing cross-protection between pH1N1 and other H1 viruses, which was directly correlated with the relatedness in the viral hemagglutinin (HA) and neuraminidase (NA) proteins. Cross-protection against H3N2 was almost complete in pigs with immunity against H1N2, but was weak in H1N1/pH1N1-immune pigs. In conclusion, infection with a live, wild type influenza virus may offer substantial cross-lineage protection against viruses of the same HA and/or NA subtype. True heterosubtypic protection, in contrast, appears to be minimal in natural influenza virus hosts. We discuss our findings in the light of the zoonotic and pandemic risks of SIVs

Subcellular location prediction of proteins using support vector machines with alignment of block sequences utilizing amino acid composition

Background: Subcellular location prediction of proteins is an important and well-studied problem in bioinformatics. This is a problem of predicting which part in a cell a given protein is transported to, where an amino acid sequence of the protein is given as an input. This problem is becoming more important since information on subcellular location is helpful for annotation of proteins and genes and the number of complete genomes is rapidly increasing. Since existing predictors are based on various heuristics, it is important to develop a simple method with high prediction accuracies. Results: In this paper, we propose a novel and general predicting method by combining techniques for sequence alignment and feature vectors based on amino acid composition. We implemented this method with support vector machines on plant data sets extracted from the TargetP database. Through fivefold cross validation tests, the obtained overall accuracies and average MCC were 0.9096 and 0.8655 respectively. We also applied our method to other datasets including that of WoLF PSORT. Conclusion: Although there is a predictor which uses the information of gene ontology and yields higher accuracy than ours, our accuracies are higher than existing predictors which use only sequence information. Since such information as gene ontology can be obtained only for known proteins, our predictor is considered to be useful for subcellular location prediction of newly-discovered proteins. Furthermore, the idea of combination of alignment and amino acid frequency is novel and general so that it may be applied to other problems in bioinformatics. Our method for plant is also implemented as a web-system and available on http://sunflower.kuicr.kyoto-u.ac.jp/~tamura/slpfa.html webcite

Systematic search for putative new domain families in Mycoplasma gallisepticum genome

<p>Abstract</p> <p>Background</p> <p>Protein domains are the fundamental units of protein structure, function and evolution. The delineation of different domains in proteins is important for classification, understanding of structure, function and evolution. The delineation of protein domains within a polypeptide chain, namely at the genome scale, can be achieved in several ways but may remain problematic in many instances. Difficulties in identifying the domain content of a given sequence arise when the query sequence has no homologues with experimentally determined structure and searching against sequence domain databases also results in insignificant matches. Identification of domains under low sequence identity conditions and lack of structural homologues acquire a crucial importance especially at the genomic scale.</p> <p>Findings</p> <p>We have developed a new method for the identification of domains in unassigned regions through indirect connections and scaled up its application to the analysis of 434 unassigned regions in 726 protein sequences of <it>Mycoplasma gallisepticum </it>genome. We could establish 71 new domain relationships and probable 63 putative new domain families through intermediate sequences in the unassigned regions, which importantly represent an overall 10% increase in PfamA domain annotation over the direct assignment in this genome.</p> <p>Conclusions</p> <p>The systematic analysis of the unassigned regions in the <it>Mycoplasma gallisepticum </it>genome has provided some insight into the possible new domain relationships and putative new domain families. Further investigation of these predicted new domains may prove beneficial in improving the existing domain prediction algorithms.</p

Genomic Footprints of Selective Sweeps from Metabolic Resistance to Pyrethroids in African Malaria Vectors Are Driven by Scale up of Insecticide-Based Vector Control

Insecticide resistance in mosquito populations threatens recent successes in malaria prevention. Elucidating patterns of genetic structure in malaria vectors to predict the speed and direction of the spread of resistance is essential to get ahead of the `resistance curve' and to avert a public health catastrophe. Here, applying a combination of microsatellite analysis, whole genome sequencing and targeted sequencing of a resistance locus, we elucidated the continent-wide population structure of a major African malaria vector, Anopheles funestus. We identified a major selective sweep in a genomic region controlling cytochrome P450-based metabolic resistance conferring high resistance to pyrethroids. This selective sweep occurred since 2002, likely as a direct consequence of scaled up vector control as revealed by whole genome and fine-scale sequencing of pre- and post-intervention populations. Fine-scaled analysis of the pyrethroid resistance locus revealed that a resistanceassociated allele of the cytochrome P450 monooxygenase CYP6P9a has swept through southern Africa to near fixation, in contrast to high polymorphism levels before interventions, conferring high levels of pyrethroid resistance linked to control failure. Population structure analysis revealed a barrier to gene flow between southern Africa and other areas, which may prevent or slow the spread of the southern mechanism of pyrethroid resistance to other regions. By identifying a genetic signature of pyrethroid-based interventions, we have demonstrated the intense selective pressure that control interventions exert on mosquito populations. If this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised

Detection and quantitation of copy number variation in the voltage-gated sodium channel gene of the mosquito Culex quinquefasciatus

Insecticide resistance is typically associated with alterations to the insecticidal target-site or with gene expression variation at loci involved in insecticide detoxification. In some species copy number variation (CNV) of target site loci (e.g. the Ace-1 target site of carbamate insecticides) or detoxification genes has been implicated in the resistance phenotype. We show that field-collected Ugandan Culex quinquefasciatus display CNV for the voltage-gated sodium channel gene (Vgsc), target-site of pyrethroid and organochlorine insecticides. In order to develop field-applicable diagnostics for Vgsc CN, and as a prelude to investigating the possible association of CN with insecticide resistance, three assays were compared for their accuracy in CN estimation in this species. The gold standard method is droplet digital PCR (ddPCR), however, the hardware is prohibitively expensive for widespread utility. Here, ddPCR was compared to quantitative PCR (qPCR) and pyrosequencing. Across all platforms, CNV was detected in ≈10% of mosquitoes, corresponding to three or four copies (per diploid genome). ddPCR and qPCR-Std-curve yielded similar predictions for Vgsc CN, indicating that the qPCR protocol developed here can be applied as a diagnostic assay, facilitating monitoring of Vgsc CN in wild populations and the elucidation of association between the Vgsc CN and insecticide resistance

Large introns in relation to alternative splicing and gene evolution: a case study of Drosophila bruno-3

Background: Alternative splicing (AS) of maturing mRNA can generate structurally and functionally distinct transcripts from the same gene. Recent bioinformatic analyses of available genome databases inferred a positive correlation between intron length and AS. To study the interplay between intron length and AS empirically and in more detail, we analyzed the diversity of alternatively spliced transcripts (ASTs) in the Drosophila RNA-binding Bruno-3 (Bru-3) gene. This gene was known to encode thirteen exons separated by introns of diverse sizes, ranging from 71 to 41,973 nucleotides in D. melanogaster. Although Bru-3's structure is expected to be conducive to AS, only two ASTs of this gene were previously described. Results: Cloning of RT-PCR products of the entire ORF from four species representing three diverged Drosophila lineages provided an evolutionary perspective, high sensitivity, and long-range contiguity of splice choices currently unattainable by high-throughput methods. Consequently, we identified three new exons, a new exon fragment and thirty-three previously unknown ASTs of Bru-3. All exon-skipping events in the gene were mapped to the exons surrounded by introns of at least 800 nucleotides, whereas exons split by introns of less than 250 nucleotides were always spliced contiguously in mRNA. Cases of exon loss and creation during Bru-3 evolution in Drosophila were also localized within large introns. Notably, we identified a true de novo exon gain: exon 8 was created along the lineage of the obscura group from intronic sequence between cryptic splice sites conserved among all Drosophila species surveyed. Exon 8 was included in mature mRNA by the species representing all the major branches of the obscura group. To our knowledge, the origin of exon 8 is the first documented case of exonization of intronic sequence outside vertebrates. Conclusion: We found that large introns can promote AS via exon-skipping and exon turnover during evolution likely due to frequent errors in their removal from maturing mRNA. Large introns could be a reservoir of genetic diversity, because they have a greater number of mutable sites than short introns. Taken together, gene structure can constrain and/or promote gene evolution

Genomic perspectives on the evolution of fungal entomopathogenicity in Beauveria bassiana

The ascomycete fungus Beauveria bassiana is a pathogen of hundreds of insect species and is commercially produced as an environmentally friendly mycoinsecticide. We sequenced the genome of B. bassiana and a phylogenomic analysis confirmed that ascomycete entomopathogenicity is polyphyletic, but also revealed convergent evolution to insect pathogenicity. We also found many species-specific virulence genes and gene family expansions and contractions that correlate with host ranges and pathogenic strategies. These include B. bassiana having many more bacterial-like toxins (suggesting an unsuspected potential for oral toxicity) and effector-type proteins. The genome also revealed that B. bassiana resembles the closely related Cordyceps militaris in being heterothallic, although its sexual stage is rarely observed. A high throughput RNA-seq transcriptomic analysis revealed that B. bassiana could sense and adapt to different environmental niches by activating well-defined gene sets. The information from this study will facilitate further development of B. bassiana as a cost-effective mycoinsecticide

PCR-based methods for the detection of L1014 kdr mutation in Anopheles culicifacies sensu lato

<p>Abstract</p> <p>Background</p> <p><it>Anopheles culicifacies s.l</it>., a major malaria vector in India, has developed widespread resistance to DDT and is becoming resistant to pyrethroids–the only insecticide class recommended for the impregnation of bed nets. Knock-down resistance due to a point mutation in the voltage gated sodium channel at L1014 residue (<it>kdr</it>) is a common mechanism of resistance to DDT and pyrethroids. The selection of this resistance may pose a serious threat to the success of the pyrethroid-impregnated bed net programme. This study reports the presence of <it>kdr </it>mutation (L1014F) in a field population of <it>An. culicifacies s.l</it>. and three new PCR-based methods for <it>kdr </it>genotyping.</p> <p>Methods</p> <p>The IIS4-IIS5 linker to IIS6 segments of the para type voltage gated sodium channel gene of DDT and pyrethroid resistant <it>An. culicifacies s.l</it>. population from the Surat district of India was sequenced. This revealed the presence of an A-to-T substitution at position 1014 leading to a leucine-phenylalanine mutation (L1014F) in a few individuals. Three molecular methods viz. Allele Specific PCR (AS-PCR), an Amplification Refractory Mutation System (ARMS) and Primer Introduced Restriction Analysis-PCR (PIRA-PCR) were developed and tested for <it>kdr </it>genotyping. The specificity of the three assays was validated following DNA sequencing of the samples genotyped.</p> <p>Results</p> <p>The genotyping of this <it>An. culicifacies s.l</it>. population by the three PCR based assays provided consistent result and were in agreement with DNA sequencing result. A low frequency of the <it>kdr </it>allele mostly in heterozygous condition was observed in the resistant population. Frequencies of the different genotypes were in Hardy-Weinberg equilibrium.</p> <p>Conclusion</p> <p>The Leu-Phe mutation, which generates the <it>kdr </it>phenotype in many insects, was detected in a pyrethroid and DDT resistant <it>An. culicifacies s.l</it>. population. Three PCR-based methods were developed for <it>kdr </it>genotyping. All the three assays were specific. The ARMS method was refractory to non-specific amplification in non-stringent amplification conditions. The PIRA-PCR assay is able to detect both the codons for the phenylalanine mutation at <it>kdr </it>locus, i.e., TTT and TTC, in a single assay, although the latter codon was not found in the population genotyped.</p