Genome-wide interaction studies identify sex-specific risk alleles for nonsyndromic orofacial clefts.

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is the most common craniofacial birth defect in humans and is notable for its apparent sexual dimorphism where approximately twice as many males are affected as females. The sources of this disparity are largely unknown, but interactions between genetic and sex effects are likely contributors. We examined gene-by-sex (G × S) interactions in a worldwide sample of 2,142 NSCL/P cases and 1,700 controls recruited from 13 countries. First, we performed genome-wide joint tests of the genetic (G) and G × S effects genome-wide using logistic regression assuming an additive genetic model and adjusting for 18 principal components of ancestry. We further interrogated loci with suggestive results from the joint test ( p < 1.00 × 10 -5 ) by examining the G × S effects from the same model. Out of the 133 loci with suggestive results ( p < 1.00 × 10 -5 ) for the joint test, we observed one genome-wide significant G × S effect in the 10q21 locus (rs72804706; p = 6.69 × 10 -9 ; OR = 2.62 CI [1.89, 3.62]) and 16 suggestive G × S effects. At the intergenic 10q21 locus, the risk of NSCL/P is estimated to increase with additional copies of the minor allele for females, but the opposite effect for males. Our observation that the impact of genetic variants on NSCL/P risk differs for males and females may further our understanding of the genetic architecture of NSCL/P and the sex differences underlying clefts and other birth defects.

A Genome-wide Association Study of Nonsyndromic Cleft Palate Identifies an Etiologic Missense Variant in GRHL3.

Cleft palate (CP) is a common birth defect occurring in 1 in 2,500 live births. Approximately half of infants with CP have a syndromic form, exhibiting other physical and cognitive disabilities. The other half have nonsyndromic CP, and to date, few genes associated with risk for nonsyndromic CP have been characterized. To identify such risk factors, we performed a genome-wide association study of this disorder. We discovered a genome-wide significant association with a missense variant in GRHL3 (p.Thr454Met [c.1361C>T]; rs41268753; p = 4.08 × 10(-9)) and replicated the result in an independent sample of case and control subjects. In both the discovery and replication samples, rs41268753 conferred increased risk for CP (OR = 8.3, 95% CI 4.1-16.8; OR = 2.16, 95% CI 1.43-3.27, respectively). In luciferase transactivation assays, p.Thr454Met had about one-third of the activity of wild-type GRHL3, and in zebrafish embryos, perturbed periderm development. We conclude that this mutation is an etiologic variant for nonsyndromic CP and is one of few functional variants identified to date for nonsyndromic orofacial clefting. This finding advances our understanding of the genetic basis of craniofacial development and might ultimately lead to improvements in recurrence risk prediction, treatment, and prognosis.

Structural and biochemical evaluation of Ceratitis capitata odorant-binding protein 22 affinity for odorants involved in intersex communication.

In insects, odorant-binding proteins (OBPs) connect the peripheral sensory system to receptors of olfactory organs. Medfly Ceratitis capitata CcapObp22 shows 37% identity and close phylogenetic affinities with Drosophila melanogaster OBP69a/pheromone-binding protein related protein 1. The CcapObp22 gene is transcribed in the antennae and maxillary palps, suggesting an active role in olfaction. Here, we recombinantly produced CcapObp22, obtaining a 13.5 kDa protein capable of binding multiple strongly hydrophobic terpene compounds, including medfly male pheromone components. The highest binding affinity [half maximal effective concentration (EC50) = 0.48 µM] was to (E,E)-α-farnesene, one of the most abundant compounds in the male pheromone blend. This odorant was used in cocrystallization experiments, yielding the structure of CcapOBP22. The monomeric structure shows the typical OBP folding, constituted by six α-helical elements interconnected by three disulphide bridges. A C-terminal seventh α-helix constitutes the wall of a deep, L-shaped hydrophobic cavity. Analysis of the electron density in this cavity suggested trapping of farnesene in the crystal structure, although with partial occupancy. Superposition of the CcapOBP22 structure with related seven-helical OBPs highlights striking similarity in the organization of the C-terminal segment of these proteins. Collectively, our molecular and physiological data on medfly CcapOBP22 suggest its involvement in intersex olfactory communication.

Association of low-frequency genetic variants in regulatory regions with nonsyndromic orofacial clefts.

Genome-wide scans have shown that common risk alleles for orofacial clefts (OFC) tend to be located in noncoding regulatory elements and cumulatively explain only part of the heritability of OFCs. Low-frequency variants may account for some of the "missing" heritability. Therefore, we scanned low-frequency variants located within putative craniofacial enhancers to identify novel OFC risk variants and implicate new regulatory elements in OFC pathogenesis. Analyses were performed in a multiethnic sample of 1,995 cases of cleft lip with or without cleft palate (CL/P), 221 cases with cleft palate (CP) only, and 1,576 unaffected controls. One hundred and nineteen putative craniofacial enhancers identified from ChIP-Seq studies in craniofacial tissues or cell lines contained multiple low-frequency (0.01-1%) variants, which we genotyped in participants using a custom Illumina panel. Two complementary statistical approaches, sequence kernel association test and combined multivariate and collapsing, were used to test association of the aggregated low-frequency variants across each enhancer region with CL/P and CP. We discovered a significant association between CP and a branchial arch enhancer near FOXP1 (mm60; p-value = .0002). Additionally, we observed a suggestive association between CL/P and a forebrain enhancer near FOXE1 (hs1717; p-value = .001). These findings suggest that low-frequency variants in craniofacial enhancer regions contribute to the complex etiology of nonsyndromic OFCs.

Identification of 16q21 as a modifier of nonsyndromic orofacial cleft phenotypes.

Orofacial clefts (OFCs) are common, complex birth defects with extremely heterogeneous phenotypic presentations. Two common subtypes-cleft lip alone (CL) and CL plus cleft palate (CLP)-are typically grouped into a single phenotype for genetic analysis (i.e., CL with or without cleft palate, CL/P). However, mounting evidence suggests there may be unique underlying pathophysiology and/or genetic modifiers influencing expression of these two phenotypes. To this end, we performed a genome-wide scan for genetic modifiers by directly comparing 450 CL cases with 1,692 CLP cases from 18 recruitment sites across 13 countries from North America, Central or South America, Asia, Europe, and Africa. We identified a region on 16q21 that is strongly associated with different cleft type (P = 5.611 × 10-8 ). We also identified significant evidence of gene-gene interactions between this modifier locus and two recognized CL/P risk loci: 8q21 and 9q22 (FOXE1) (P = 0.012 and 0.023, respectively). Single nucleotide polymorphism (SNPs) in the 16q21 modifier locus demonstrated significant association with CL over CLP. The marker alleles on 16q21 that increased risk for CL were found at highest frequencies among individuals with a family history of CL (P = 0.003). Our results demonstrate the existence of modifiers for which type of OFC develops and suggest plausible elements responsible for phenotypic heterogeneity, further elucidating the complex genetic architecture of OFCs.

A multi-ethnic genome-wide association study identifies novel loci for non-syndromic cleft lip with or without cleft palate on 2p24.2, 17q23 and 19q13.

Orofacial clefts (OFCs), which include non-syndromic cleft lip with or without cleft palate (CL/P), are among the most common birth defects in humans, affecting approximately 1 in 700 newborns. CL/P is phenotypically heterogeneous and has a complex etiology caused by genetic and environmental factors. Previous genome-wide association studies (GWASs) have identified at least 15 risk loci for CL/P. As these loci do not account for all of the genetic variance of CL/P, we hypothesized the existence of additional risk loci. We conducted a multiethnic GWAS in 6480 participants (823 unrelated cases, 1700 unrelated controls and 1319 case-parent trios) with European, Asian, African and Central and South American ancestry. Our GWAS revealed novel associations on 2p24 near FAM49A, a gene of unknown function (P = 4.22 × 10-8), and 19q13 near RHPN2, a gene involved in organizing the actin cytoskeleton (P = 4.17 × 10-8). Other regions reaching genome-wide significance were 1p36 (PAX7), 1p22 (ARHGAP29), 1q32 (IRF6), 8q24 and 17p13 (NTN1), all reported in previous GWASs. Stratification by ancestry group revealed a novel association with a region on 17q23 (P = 2.92 × 10-8) among individuals with European ancestry. This region included several promising candidates including TANC2, an oncogene required for development, and DCAF7, a scaffolding protein required for craniofacial development. In the Central and South American ancestry group, significant associations with loci previously identified in Asian or European ancestry groups reflected their admixed ancestry. In summary, we have identified novel CL/P risk loci and suggest new genes involved in craniofacial development, confirming the highly heterogeneous etiology of OFCs.

Genome-wide meta-analyses of nonsyndromic orofacial clefts identify novel associations between FOXE1 and all orofacial clefts, and TP63 and cleft lip with or without cleft palate.

Nonsyndromic orofacial clefts (OFCs) are a heterogeneous group of common craniofacial birth defects with complex etiologies that include genetic and environmental risk factors. OFCs are commonly categorized as cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which have historically been analyzed as distinct entities. Genes for both CL/P and CP have been identified via multiple genome-wide linkage and association studies (GWAS); however, altogether, known variants account for a minority of the estimated heritability in risk to these craniofacial birth defects. We performed genome-wide meta-analyses of CL/P, CP, and all OFCs across two large, multiethnic studies. We then performed population-specific meta-analyses in sub-samples of Asian and European ancestry. In addition to observing associations with known variants, we identified a novel genome-wide significant association between SNPs located in an intronic TP63 enhancer and CL/P (p = 1.16 × 10-8). Several novel loci with compelling candidate genes approached genome-wide significance on 4q21.1 (SHROOM3), 12q13.13 (KRT18), and 8p21 (NRG1). In the analysis of all OFCs combined, SNPs near FOXE1 reached genome-wide significance (p = 1.33 × 10-9). Our results support the highly heterogeneous nature of OFCs and illustrate the utility of meta-analysis for discovering new genetic risk factors.

Use of the synergist piperonyl butoxide can slow the development of alpha-cypermethrin resistance in the whitefly Bemisia tabaci.

The development of insecticide resistance in insect pests of crops is a growing threat to sustainable food production, and strategies that slow the development of resistance are therefore urgently required. The insecticide synergist piperonyl butoxide (PBO) inhibits certain insect detoxification systems and so may delay the evolution of metabolic resistance. In the current study we characterized resistance development in the silverleaf whitefly, Bemisia tabaci, after selection with either a neonicotinoid (thiacloprid) or pyrethroid (alpha-cypermethrin) insecticide alone or in combination with PBO. Resistance development was significantly suppressed (> 60%) in the line selected with alpha-cypermethrin + PBO compared to the line selected with alpha-cypermethrin alone. RNA sequencing (RNAseq) analyses revealed an increase in frequency of a knock-down resistance mutation but no differentially expressed genes were identified that could explain the sensitivity shift. No significant difference was observed in the level of resistance between the thiacloprid and thiacloprid + PBO selected lines, and RNA sequencing (RNAseq) analyses revealed that the cytochrome P450 monooxygenase CYP6CM1, known to metabolize neonicotinoids, was significantly upregulated (>10-fold) in both lines. The findings of this study demonstrate that PBO used in combination with certain insecticides can suppress the development of resistance in a laboratory setting; however, the mechanism by which PBO supresses resistance development remains unclear.

Induced thiacloprid insensitivity in honeybees (Apis mellifera L.) is associated with up-regulation of detoxification genes.

Honey bees, Apis mellifera, are markedly less sensitive to neonicotinoid insecticides containing a cyanoimino pharmacophore than to those with a nitroimino group. Although previous work has suggested that this results from enhanced metabolism of the former by detoxification enzymes, the specific enzyme(s) involved remain to be characterized. In this work, a pretreatment of honey bees with a sublethal dose of thiacloprid resulted in induced insensitivity to the same compound immediately following thiacloprid feeding. A longer pretreatment time resulted in no, or increased, sensitivity. Transcriptome profiling, using microarrays, identified a number of genes encoding detoxification enzymes that were over-expressed significantly in insecticide-treated bees compared with untreated controls. These included five P450s, CYP6BE1, CYP305D1, CYP6AS5, CYP315A1, CYP301A1, and a carboxyl/cholinesterase (CCE) CCE8. Four of these P450s were functionally expressed in Escherichia coli and their ability to metabolize thiacloprid examined by liquid chromatography-mass spectrometry (LC-MS) analysis.

A locus on chromosome 15q26 (IDDM3) produces susceptibility to insulin-dependent diabetes mellitus.

To identify new loci predisposing to insulin-dependent diabetes mellitus (IDDM), we have investigated 250 families with more than one diabetic child. Affected sibling pair linkage analysis revealed strong evidence for an IDDM susceptibility locus near D15S107 on chromosome 15q26 (P = 0.0010) termed IDDM3. Families less predisposed through genes in the HLA region provided most of the evidence for linkage. In these families, discordant sibling pairs also showed linkage (P = 0.0052), and sibling pair disease concordance or discordance was strongly related to the proportion of genes the pair shared at D15S107 (P = 0.0003). Our study also revealed evidence for an IDDM locus on chromosome 11q13 (IDDM4) using affected siblings (P = 0.0043), but no evidence using discordant siblings.

A specific requirement for PDGF-C in palate formation and PDGFR-alpha signaling.

PDGF-C is a member of the platelet-derived growth factor (PDGF) family, which signals through PDGF receptor (PDGFR) alphaalpha and alphabeta dimers. Here we show that Pdgfc(-/-) mice die in the perinatal period owing to feeding and respiratory difficulties associated with a complete cleft of the secondary palate. This phenotype was less severe than that of Pdgfra(-/-) embryos. Pdgfc(-/-) Pdgfa(-/-) embryos developed a cleft face, subepidermal blistering, deficiency of renal cortex mesenchyme, spina bifida and skeletal and vascular defects. Complete loss of function of both ligands, therefore, phenocopied the loss of PDGFR-alpha function, suggesting that both PDGF-A and PDGF-C signal through PDGFR-alpha to regulate the development of craniofacial structures, the neural tube and mesodermal organs. Our results also show that PDGF-C signaling is a new pathway in palatogenesis, different from, and independent of, those previously implicated.

A family-based association analysis and meta-analysis of the reading disabilities candidate gene DYX1C1.

Reading disabilities (RD) have a significant genetic basis and have shown linkage to multiple regions including chromosome 15q. Dyslexia susceptibility 1 candidate gene 1 (DYX1C1) on chromosome 15q21 was originally proposed as a candidate gene with two potentially functional polymorphisms at the -3G/A and 1249G/T positions showing association with RD. However, subsequent studies have yielded mixed results. We performed a literature review and meta-analysis of the -3G/A and 1249G/T polymorphisms, including new unpublished data from two family-based samples. Ten markers in DYX1C1 were genotyped in the two independently ascertained samples. Single marker and -3G/A:1249G/T haplotype analyses were performed for RD in both samples, and quantitative trait analyses using standardized reading-related measures was performed in one of the samples. For the meta-analysis, we used a random-effects model to summarize studies that tested for association between -3G/A or 1249G/T and RD. No significant association was found between the DYX1C1 SNPs and RD or any of the reading-related measures tested after correction for the number of tests performed. The previously reported risk haplotype (-3A:1249T) was not biased in transmission. A total of 9 and 10 study samples were included in the meta-analysis of the -3G/A and 1249G/T polymorphisms, respectively. Neither polymorphism reached statistical significance, but the heterogeneity for the 1249G/T polymorphism was high. The results of this study do not provide evidence for association between the putatively functional SNPs -3G/A and 1249G/T and RD.

Expression of chemosensory proteins in the tsetse fly Glossina morsitans morsitans is related to female host-seeking behaviour.

Chemosensory proteins (CSPs) are a class of soluble proteins present in high concentrations in the sensilla of insect antennae. It has been proposed that they play an important role in insect olfaction by mediating interactions between odorants and odorant receptors. Here we report, for the first time, the presence of five CSP genes in the tsetse fly Glossina morsitans morsitans, a major vector transmitting nagana in livestock. Real-time quantitative reverse transcription PCR showed that three of the CSPs are expressed in antennae. One of them, GmmCSP2, is transcribed at a very high level and could be involved in olfaction. We also determined expression in the antennae of both males and females at different life stages and with different blood feeding regimes. The transcription of GmmCSP2 was lower in male antennae than in females, with a sharp increase in 10-week-old flies, 48 h after a bloodmeal. Thus there is a clear relationship between CSP gene transcription and host searching behaviour. Genome annotation and phylogenetic analyses comparing G. morsitans morsitans CSPs with those of other Diptera showed rapid evolution after speciation of mosquitoes.

Induction and activation of antiviral enzyme 2',5'-oligoadenylate synthetase by in vitro transcribed insulin mRNA and other cellular RNAs.

Double-stranded RNA (dsRNA) can induce antiviral enzyme 2',5'-oligoadenylate synthetase (2'5'AS) expression and activate latent 2'5'AS. Our previous data have shown pancreatic ß cells are sensitive to dsRNA-induced 2'5'AS expression, and constitutive high basal 2'5'AS expression is associated with susceptibility to developing type 1 diabetes (T1D), a disease due to pancreatic ß cell loss. Here we report that in vitro transcribed human insulin mRNA induces the activation of human OAS gene promoter sequences, and specifically and dose-dependently induces 2'5'AS expression in murine pancreatic ßTC3 cells. Over-expression of dsRNA receptor retinoic acid-inducible gene-1 enhances insulin mRNA-induced 2'5'AS expression. In vitro transcribed insulin and other mRNAs, as well as total cellular RNAs, activate latent 2'5'AS in vitro with activation ability likely associated with the sequence and length of individual mRNAs or the sample source of total cellular RNA. Insulin mRNA does not show any specificity to activate 2'5'AS, but total cellular RNA from ßTC3 cells has high activation ability. Constitutive 2'5'AS expression in ßTC3 cells leads to cell proliferation inhibition and apoptosis. Our study suggests the possibility of cellular RNA-regulated 2'5'AS expression and activation, and the potential risk of high insulin gene transcription in pancreatic ß cells, and may help explain genetic predisposition to T1D associated with INS VNTR class I alleles.

The re-emergence of the bed bug as a nuisance pest: implications of resistance to the pyrethroid insecticides.

A global resurgence of bed bugs (Hemiptera: Cimicidae) has led to renewed scientific interest in these insects. The current bed bug upsurge appears to have started almost synchronously in the late 1990 s in Europe, the U.S.A. and Australia. Several factors have led to this situation, with resistance to applied insecticides making a significant contribution. With a growing number of insecticides (DDT, carbamates, organophosphates etc.) being no longer available as a result of regulatory restrictions, the mainstay chemistry used for bed bug control over the past few decades has been the pyrethroid insecticides. With reports of increasing tolerance to pyrethroids leading to control failures on the rise, containing and eradicating bed bugs is proving to be a difficult task. Consequently, several recent studies have focused on determining the mode of action of pyrethroid resistance in bed bug populations sourced from different locations. Correct identification of the factor(s) responsible for the increasing resistance is critical to the development of effective management strategies, which need to be based, wherever possible, on firm scientific evidence. Here we review the literature on this topic, highlighting the mechanisms thought to be involved and the problems currently faced by pest control professionals in dealing with a developing pandemic.

Overexpression of a cytochrome P450 monooxygenase, CYP6ER1, is associated with resistance to imidacloprid in the brown planthopper, Nilaparvata lugens.

The brown planthopper, Nilaparvata lugens, is an economically significant pest of rice throughout Asia and has evolved resistance to many insecticides including the neonicotinoid imidacloprid. The resistance of field populations of N. lugens to imidacloprid has been attributed to enhanced detoxification by cytochrome P450 monooxygenases (P450s), although, to date, the causative P450(s) has (have) not been identified. In the present study, biochemical assays using the model substrate 7-ethoxycoumarin showed enhanced P450 activity in several resistant N. lugens field strains when compared with a susceptible reference strain. Thirty three cDNA sequences encoding tentative unique P450s were identified from two recent sequencing projects and by degenerate PCR. The mRNA expression level of 32 of these was examined in susceptible, moderately resistant and highly resistant N. lugens strains using quantitative real-time PCR. A single P450 gene (CYP6ER1) was highly overexpressed in all resistant strains (up to 40-fold) and the level of expression observed in the different N. lugens strains was significantly correlated with the resistance phenotype. These results provide strong evidence for a role of CYP6ER1 in the resistance of N. lugens to imidacloprid.

Baseline ecological risk assessment of the Calcasieu Estuary, Louisiana: Part 3. An evaluation of the risks to benthic invertebrates associated with exposure to contaminated sediments.

The sediments in the Calcasieu Estuary are contaminated with a wide variety of chemicals of potential concern (COPCs), including heavy metals, polycyclic aromatic hydrocarbons, polychlorinated biphenyls, phthalates, chlorinated benzenes, and polychlorinated dibenzo-p-dioxins and dibenzofurans. The sources of these COPCs include both point and non-point source discharges. As part of a baseline ecological risk assessment, the risks to benthic invertebrates posed by exposure to sediment-associated COPCs were assessed using five lines of evidence, including whole-sediment chemistry, pore-water chemistry, whole-sediment toxicity, pore-water toxicity, and benthic invertebrate community structure. The results of this assessment indicated that exposure to whole sediments and/or pore water from the Calcasieu Estuary generally posed low risks to benthic invertebrate communities (i.e., risks were classified as low for 68% of the sampling locations investigated). However, incremental risks to benthic invertebrates (i.e., compared with those associated with exposure to conditions in reference areas) were indicated for 32% of the sampling locations within the estuary. Of the three areas of concern (AOCs) investigated, the risks to benthic invertebrates were highest in the Bayou d'Inde AOC; risks were generally lower in the Upper Calcasieu River AOC and Middle Calcasieu River AOC. The areas showing the highest risks to sediment-dwelling organisms were generally located in the vicinity of point source discharges of COPCs. These results provided risk managers with the information required to make decisions regarding the need for remedial actions at the site.

Review and optimisation of foot radiography technique.

INTRODUCTION: Radiography is an important tool in the diagnosis of trauma and diseases of the foot. Consistent image quality and anatomical appearances are essential to ensure that images are interpreted correctly and this study aims to evaluate current technique and use anthropomorphical phantoms in order to investigate the most optimal radiographic technique. METHODS: A two part study was conducted. An electronic survey was administered in order to ascertain the current radiographic technique utilised within the United Kingdom (UK) National Health Service (NHS) and compare this with techniques published within radiographic positioning textbooks in common use within the UK. A phantom experiment also sought to identify the most appropriate angulation of the central beam to employ in order to optimally visualise the anatomy for standard non-trauma radiography of the foot. A relative visual grading analysis technique was used to assess radiographic quality. RESULTS: The electronic survey was completed by 38% of departments (n = 69/182) and demonstrated wide inconsistencies in projections, centring points and beam angulations employed across the NHS. The most commonly employed techniques no longer match the majority of positions described in textbooks. The phantom experiment modelled the ability to accurately visualise anatomy with a variation of the central beam angulation. The results demonstrated that the most commonly described textbook techniques remain the most optimal for accurately visualising the anatomy of the foot. CONCLUSION: This study has validated the radiographic techniques of the foot described in many of the key texts, however self-reported current practice has diverged from this textbook description. IMPLICATIONS FOR PRACTICE: Reinforcement of fundamental principles will benefit the patient with increased image quality, improved diagnosis and improved consistency of images, particularly when patient care is transferred between centres.

Genome annotation and comparative analyses of the odorant-binding proteins and chemosensory proteins in the pea aphid Acyrthosiphon pisum.

Odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) are two families of small water-soluble proteins, abundant in the aqueous fluid surrounding olfactory receptor neurons in insect antennae. OBPs are involved in the first step of olfactory signal transduction, carrying airborne semiochemicals to the odorant receptors and can be classified into three groups: Classic OBPs, Plus-C OBPs and Atypical OBPs. Here, we identified and annotated genes encoding putative OBPs and CSPs in the pea aphid Acyrthosiphon pisum using bioinformatics. This identified genes encoding 13 Classic and two Plus-C OBPs and 13 CSPs. Homologous OBP sequences were also identified in nine other aphid species, allowing us to compare OBPs across several aphid and non-aphid species. We show that, although OBP sequences are divergent within a species and between different orders, there is a high similarity between orthologs within a range of aphid species. Furthermore, the phylogenetic relationships between OBP orthologs reflect the divergence of aphid evolution lineages. Our results support the 'birth-and-death' model as the major mechanism explaining aphid OBP sequence evolution, with the main force acting on the evolution being purifying selection.

Identification of ion channel genes in the Acyrthosiphon pisum genome.

Aphids are major pests of crops, causing hundreds of millions of dollars worth of damage annually. Ion channel proteins are often the targets of modern insecticides and mutations in ion channel genes can lead to resistance to many leading classes of insecticides. The sequencing of the pea aphid, Acyrthosiphon pisum, genome has now allowed detailed in silico analysis of the aphid ion channels. The study has revealed significant differences in the composition of the ion channel families between the aphid and other insects. For example A. pisum does not appear to contain a homologue of the nACh receptor alpha 5 gene whilst the calcium channel beta subunit has been duplicated. These variations could result in differences in function or sensitivity to insecticides. The genome sequence will allow the study of aphid ion channels to be accelerated, leading to a better understanding of the function of these economically important channels. The potential for identifying novel insecticide targets within the aphid is now a step closer.