Phosphine resistance in the rust red flour beetle, Tribolium castaneum (Coleoptera: Tenebrionidae): inheritance, gene interactions and fitness costs

The recent emergence of heritable high level resistance to phosphine in stored grain pests is a serious concern among major grain growing countries around the world. Here we describe the genetics of phosphine resistance in the rust red flour beetle Tribolium castaneum (Herbst), a pest of stored grain as well as a genetic model organism. We investigated three field collected strains of T. castaneum viz., susceptible (QTC4), weakly resistant (QTC1012) and strongly resistant (QTC931) to phosphine. The dose-mortality responses of their test- and inter-cross progeny revealed that most resistance was conferred by a single major resistance gene in the weakly (3.2×) resistant strain. This gene was also found in the strongly resistant (431×) strain, together with a second major resistance gene and additional minor factors. The second major gene by itself confers only 12-20× resistance, suggesting that a strong synergistic epistatic interaction between the genes is responsible for the high level of resistance (431×) observed in the strongly resistant strain. Phosphine resistance is not sex linked and is inherited as an incompletely recessive, autosomal trait. The analysis of the phenotypic fitness response of a population derived from a single pair inter-strain cross between the susceptible and strongly resistant strains indicated the changes in the level of response in the strong resistance phenotype; however this effect was not consistent and apparently masked by the genetic background of the weakly resistant strain. The results from this work will inform phosphine resistance management strategies and provide a basis for the identification of the resistance genes.

The rph1 gene is a common contributor to the evolution of phosphine resistance in independent field isolates of Rhyzopertha dominica

Phosphine is the only economically viable fumigant for routine control of insect pests of stored food products, but its continued use is now threatened by the world-wide emergence of high-level resistance in key pest species. Phosphine has a unique mode of action relative to well-characterised contact pesticides. Similarly, the selective pressures that lead to resistance against field sprays differ dramatically from those encountered during fumigation. The consequences of these differences have not been investigated adequately. We determine the genetic basis of phosphine resistance in Rhyzopertha dominica strains collected from New South Wales and South Australia and compare this with resistance in a previously characterised strain from Queensland. The resistance levels range from 225 and 100 times the baseline response of a sensitive reference strain. Moreover, molecular and phenotypic data indicate that high-level resistance was derived independently in each of the three widely separated geographical regions. Despite the independent origins, resistance was due to two interacting genes in each instance. Furthermore, complementation analysis reveals that all three strains contain an incompletely recessive resistance allele of the autosomal rph1 resistance gene. This is particularly noteworthy as a resistance allele at rph1 was previously proposed to be a necessary first step in the evolution of high-level resistance. Despite the capacity of phosphine to disrupt a wide range of enzymes and biological processes, it is remarkable that the initial step in the selection of resistance is so similar in isolated outbreaks.

Two-tiered approach identifies a network of cancer and liver disease-related genes regulated by miR-122

MicroRNAs function as important regulators of gene expression and are commonly linked to development, differentiation, and diseases such as cancer. To better understand their roles in various biological processes, identification of genes targeted by microRNAs is necessary. Although prediction tools have significantly helped with this task, experimental approaches are ultimately required for extensive target search and validation. We employed two independent yet complementary high throughput approaches to map a large set of mRNAs regulated by miR-122, a liver-specific microRNA implicated in regulation of fatty acid and cholesterol metabolism, hepatitis C infection, and hepatocellular carcinoma. The combination of luciferase reporter-based screening and shotgun proteomics resulted in the identification of 260 proteins significantly down-regulated in response to miR-122 in at least one method, 113 of which contain predicted miR-122 target sites. These proteins are enriched for functions associated with the cell cycle, differentiation, proliferation, and apoptosis. Among these miR-122-sensitive proteins, we identified a large group with strong connections to liver metabolism, diseases, and hepatocellular carcinoma. Additional analyses, including examination of consensus binding motifs for both miR-122 and target sequences, provide further insight into miR-122 function.

The spatiotemporal dynamics of Tribolium castaneum (Herbst): adult flight and gene flow

Tribolium castaneum (Herbst) has been used as a model organism to develop and test important ecological and evolutionary concepts and is also a major pest of grain and grain products globally. This beetle species is assumed to be a good colonizer of grain storages through anthropogenic movement of grain, and active dispersal by flight is considered unlikely. Studies using T. castaneum have therefore used confined or walking insects. We combine an ecological study of dispersal with an analysis of gene flow using microsatellites to investigate the spatiotemporal dynamics and adult flight of T. castaneum in an ecological landscape in eastern Australia. Flying beetles were caught in traps at grain storages and in fields at least 1 km from the nearest stored grain at regular intervals for an entire year. Significantly more beetles were trapped at storages than in fields, and almost no beetles were caught in native vegetation reserves many kilometres from the nearest stored grain. Genetic differentiation between beetles caught at storages and in fields was low, indicating that although T. castaneum is predominantly aggregated around grain storages, active dispersal takes place to the extent that significant gene flow occurs between them, mitigating founder effects and genetic drift. By combining ecological and molecular techniques, we reveal much higher levels of active dispersal through adult flight in T. castaneum than previously thought. We conclude that the implications of adult flight to previous and future studies on this model organism warrant consideration.

Identification of a G2-like porcine rotavirus bearing a novel VP4 type, P[32]

A porcine group A rotavirus (GARV) strain, 61/07/Ire, was isolated from a 4–5 week asymptomatic piglet, during an epidemiological survey of porcine herds in Southern Ireland, in 2007. The nucleotide (nt) and amino acid (aa) sequence of the full-length VP4 protein of the PoRV strain 61/07/Ire was determined. Based on the entire VP4 open reading frame (nt), strain 61/07/Ire displayed ≤ 76.5% identity to representatives of the established 31 P-types, a value far lower than the percentage identity cutoff value (80%) established by the Rotavirus Classification Working Group (RCWG) to define a novel P genotype. Strain 61/07/Ire revealed low aa identity, ranging from 57.1% to 83.6%, to the cognate sequences of representatives of the various P genotypes. The aa identity was lower in the VP8* trypsin-cleavage fragment of the VP4, which encompasses the VP4 hypervariable region, ranging from 36.9% to 75.3%. Sequence analyses of the VP7, VP6, and NSP4 genes revealed that the GARV strain 61/07/Ire possessed a G2-like VP7, an E9 NSP4 genotype and an I5 VP6 genotype. Altogether, these results indicate that the GARV strain 61/07/Ire should be considered as a prototype of a new VP4 genotype, P[32], and provide further evidence for the vast heterogeneity of group A rotaviruses.

[Impact of avian influenza virus H5N1 neuraminidase mutations on the activity of neuraminidase and the sensibility to neuraminidase inhibitors]

OBJECTIVE

To study the impact of avian influenza virus H5N1 neuraminidase mutations I117V, I314V and I117V + I314V on the sensibility of neuraminidase inhibitors (NAIs) and the activity of neuraminidase (NA).

METHODS

The mutations were introduced into NA genes of virus strain A/Vietnam/1203/04 (H5N1) by site-directed mutagenesis. With the A/WSN/33 (H1N1) background, recombinant influenza viruses containing NA mutations were rescued by reverse genetics. After viral propagation in chicken embryos, fluorimetric assays were conducted to assess the sensibility to NAIs and NA activity (IC(50), Km & Ki).

RESULTS

Compared to the wild-type virus VN1203, the mutation I117V decreased the susceptibility to oseltamivir (17-fold increment of IC(50) value, 20-fold increment of Ki value) and the NA activity (23-fold increment of Km value) while there was little impact on zanamivir sensitivity (2-fold increment of IC(50) value, 3-fold increment of Ki value). The mutation I314V had no marked influence on either the NA activity or the NAIs susceptibility.

CONCLUSION

It appears that the NA mutations of I117V and I314V can not cause NAIs resistance. Oseltamivir or zanamivir may still be prescribed for anti-viral treatment.

Detection and characterisation of group A rotavirus in asymptomatic piglets in southern Ireland

Porcine group A rotaviruses (GARV) are causative agents of enteritis in piglets and are a large reservoir of genetic material for the diversification of human GARVs. Accumulation of information on the genetic heterogeneity of porcine viruses is pivotal for readily characterising unusual human strains. Screening of 292 fecal samples, collected from 4-5- to 8-9-week-old asymptomatic pigs from four herds in Ireland between 2005 and 2007 resulted in 19 (6.5%) samples testing positive by reverse-transcription PCR (RT-PCR) for GARV. The strains were molecularly characterized to collate data on the VP7 and partial VP4 outer capsid genes. By sequence analysis of the VP7 gene, the Irish strains were identified as G2, G4, G5, G9 and G11 viruses. The G11 strains were closely related to other human and porcine G11 strains, while the G2 strains resembled porcine G2 viruses detected recently in Europe and southern Asia. The G4 strains were distantly related to other G4 human and animal strains, constituting a separate G4 VP7 lineage. Analysis of the G5 strains revealed that they were similar to a selection of G5 human and porcine strains, while the G9 strains resembled other porcine G9 viruses. By sequence analysis of the VP8* fragment of the VP4, the Irish viruses were characterised as P[6], P[7], P[13], P[13]/[22], P[26] and P[32].

Efficient synthesis of highly active phospha-isosteres of the influenza neuraminidase inhibitor oseltamivir

With a Hunsdiecker-Barton iododecarboxylation strategy, we converted the carboxylate group of the oseltamivir precursor into exemplary phosphonate monoesters. In all cases, K(i) values towards influenza virus sialidase remained in the sub-nanomolar range. We have thus made valuable structural space available for the design of novel oseltamivir-based tools for influenza virus research.

Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs

Accurate and complete measurement of single nucleotide (SNP) and copy number (CNV) variants, both common and rare, will be required to understand the role of genetic variation in disease. We present Birdsuite, a four-stage analytical framework instantiated in software for deriving integrated and mutually consistent copy number and SNP genotypes. The method sequentially assigns copy number across regions of common copy number polymorphisms (CNPs), calls genotypes of SNPs, identifies rare CNVs via a hidden Markov model (HMM), and generates an integrated sequence and copy number genotype at every locus (for example, including genotypes such as A-null, AAB and BBB in addition to AA, AB and BB calls). Such genotypes more accurately depict the underlying sequence of each individual, reducing the rate of apparent mendelian inconsistencies. The Birdsuite software is applied here to data from the Affymetrix SNP 6.0 array. Additionally, we describe a method, implemented in PLINK, to utilize these combined SNP and CNV genotypes for association testing with a phenotype.

Integrated detection and population-genetic analysis of SNPs and copy number variation

Dissecting the genetic basis of disease risk requires measuring all forms of genetic variation, including SNPs and copy number variants (CNVs), and is enabled by accurate maps of their locations, frequencies and population-genetic properties. We designed a hybrid genotyping array (Affymetrix SNP 6.0) to simultaneously measure 906,600 SNPs and copy number at 1.8 million genomic locations. By characterizing 270 HapMap samples, we developed a map of human CNV (at 2-kb breakpoint resolution) informed by integer genotypes for 1,320 copy number polymorphisms (CNPs) that segregate at an allele frequency >1%. More than 80% of the sequence in previously reported CNV regions fell outside our estimated CNV boundaries, indicating that large (>100 kb) CNVs affect much less of the genome than initially reported. Approximately 80% of observed copy number differences between pairs of individuals were due to common CNPs with an allele frequency >5%, and more than 99% derived from inheritance rather than new mutation. Most common, diallelic CNPs were in strong linkage disequilibrium with SNPs, and most low-frequency CNVs segregated on specific SNP haplotypes.

Influence of concentration, temperature and humidity on the toxicity of phosphine to the strongly phosphine-resistant psocid Liposcelis bostrychophila Badonnel (Psocoptera: Liposcelididae)

BACKGROUND

The psocid Liposcelis bostrychophila Badonnel, is a widespread, significant pest of stored commodities, has developed strong resistance to phosphine, the major grain disinfestant. The aim was to develop effective fumigation protocols to control this resistant pest.

RESULTS

Time to population extinction of all life stages (TPE) in days was evaluated at a series of phosphine concentrations and temperatures at two relative humidities. Regression analysis showed that temperature, concentration and relative humidity all contributed significantly to describing TPE (P<0.001, R(2)=0.95), with temperature being the dominant variable, accounting for 74.4% of the variation. Irrespective of phosphine concentration, TPE was longer at lower temperatures and high humidity (70% RH) and shorter at higher temperatures and low humidity (55% RH). At any concentration of phosphine, a combination of higher temperature and lower humidity provides the shortest fumigation period to control resistant L. bostrychophila. For example, 19 and 11 days of fumigation are required at 15 degrees C and 70% RH at 0.1 and 1.0 mg L(-1) of phosphine respectively, whereas only 4 and 2 days are required at 35 degrees C and 55% RH for the same respective concentrations.

CONCLUSIONS

The developed fumigation protocols will provide industry with flexibility in application of phosphine.

The balance of reproducibility, sensitivity, and specificity of lists of differentially expressed genes in microarray studies

BACKGROUND

Reproducibility is a fundamental requirement in scientific experiments. Some recent publications have claimed that microarrays are unreliable because lists of differentially expressed genes (DEGs) are not reproducible in similar experiments. Meanwhile, new statistical methods for identifying DEGs continue to appear in the scientific literature. The resultant variety of existing and emerging methods exacerbates confusion and continuing debate in the microarray community on the appropriate choice of methods for identifying reliable DEG lists.

RESULTS

Using the data sets generated by the MicroArray Quality Control (MAQC) project, we investigated the impact on the reproducibility of DEG lists of a few widely used gene selection procedures. We present comprehensive results from inter-site comparisons using the same microarray platform, cross-platform comparisons using multiple microarray platforms, and comparisons between microarray results and those from TaqMan - the widely regarded "standard" gene expression platform. Our results demonstrate that (1) previously reported discordance between DEG lists could simply result from ranking and selecting DEGs solely by statistical significance (P) derived from widely used simple t-tests; (2) when fold change (FC) is used as the ranking criterion with a non-stringent P-value cutoff filtering, the DEG lists become much more reproducible, especially when fewer genes are selected as differentially expressed, as is the case in most microarray studies; and (3) the instability of short DEG lists solely based on P-value ranking is an expected mathematical consequence of the high variability of the t-values; the more stringent the P-value threshold, the less reproducible the DEG list is. These observations are also consistent with results from extensive simulation calculations.

CONCLUSION

We recommend the use of FC-ranking plus a non-stringent P cutoff as a straightforward and baseline practice in order to generate more reproducible DEG lists. Specifically, the P-value cutoff should not be stringent (too small) and FC should be as large as possible. Our results provide practical guidance to choose the appropriate FC and P-value cutoffs when selecting a given number of DEGs. The FC criterion enhances reproducibility, whereas the P criterion balances sensitivity and specificity.

The ets-related transcription factor GABP directs bidirectional transcription

Approximately 10% of genes in the human genome are distributed such that their transcription start sites are located less than 1 kb apart on opposite strands. These divergent gene pairs have a single intergenic segment of DNA, which in some cases appears to share regulatory elements, but it is unclear whether these regions represent functional bidirectional promoters or two overlapping promoters. A recent study showed that divergent promoters are enriched for consensus binding sequences of a small group of transcription factors, including the ubiquitous ets-family transcription factor GA-binding protein (GABP). Here we show that GABP binds to more than 80% of divergent promoters in at least one cell type. Furthermore, we demonstrate that GABP binding is correlated and associated with bidirectional transcriptional activity in a luciferase transfection assay. In addition, we find that the addition of a strict consensus GABP site into a set of promoters that normally function in only one direction significantly increases activity in the opposite direction in 67% of cases. Our findings demonstrate that GABP regulates the majority of divergent promoters and suggest that bidirectional transcriptional activity is mediated through GABP binding and transactivation at both divergent and nondivergent promoters.

Surveys of the locations of genes in the human genome have revealed that a surprising number of genes, greater than 10%, have transcription start sites within 1 kb of one another on opposite strands. These divergent gene pairs, sometimes referred to as bidirectional genes, are common in organisms such as bacteria and yeast, but it is unknown why such an arrangement exists in large, mammalian genomes. Recently, it has become apparent that the promoters of these divergent genes are regulated by a subset of transcription factors, and we have focused on one of these, GA-binding protein (GABP). We find that it regulates a large number of human genes, including the majority of divergent genes, and that its binding is associated with, correlated with, and sufficient for bidirectional transcriptional activity. Although clearly GABP is a major regulator of divergent genes, which carry out a variety of roles critical for the function and survival of the cell, these data also propose novel roles for GABP as a transcription factor. For example, the ability of GABP to promote bidirectional transcription may prove to be biologically relevant in generating many of the transcripts that have been observed outside of protein coding genes.

Integrated analysis of experimental data sets reveals many novel promoters in 1% of the human genome

The regulation of transcriptional initiation in the human genome is a critical component of global gene regulation, but a complete catalog of human promoters currently does not exist. In order to identify regulatory regions, we developed four computational methods to integrate 129 sets of ENCODE-wide chromatin immunoprecipitation data. They collectively predicted 1393 regions. Roughly 47% of the regions were unique to one method, as each method makes different assumptions about the data. Overall, predicted regions tend to localize to highly conserved, DNase I hypersensitive, and actively transcribed regions in the genome. Interestingly, a significant portion of the regions overlaps with annotated 3'-UTRs, suggesting that some of them might regulate anti-sense transcription. The majority of the predicted regions are >2 kb away from the 5'-ends of previously annotated human cDNAs and hence are novel. These novel regions may regulate unannotated transcripts or may represent new alternative transcription start sites of known genes. We tested 163 such regions for promoter activity in four cell lines using transient transfection assays, and 25% of them showed transcriptional activity above background in at least one cell line. We also performed 5'-RACE experiments on 62 novel regions, and 76% of the regions were associated with the 5'-ends of at least two RACE products. Our results suggest that there are at least 35% more functional promoters in the human genome than currently annotated.

Transcription factor binding and modified histones in human bidirectional promoters

Bidirectional promoters have received considerable attention because of their ability to regulate two downstream genes (divergent genes). They are also highly abundant, directing the transcription of approximately 11% of genes in the human genome. We categorized the presence of DNA sequence motifs, binding of transcription factors, and modified histones as overrepresented, shared, or underrepresented in bidirectional promoters with respect to unidirectional promoters. We found that a small set of motifs, including GABPA, MYC, E2F1, E2F4, NRF-1, CCAAT, YY1, and ACTACAnnTCC are overrepresented in bidirectional promoters, while the majority (73%) of known vertebrate motifs are underrepresented. We performed chromatin-immunoprecipitation (ChIP), followed by quantitative PCR for GABPA, on 118 regions in the human genome and showed that it binds to bidirectional promoters more frequently than unidirectional promoters, and its position-specific scoring matrix is highly predictive of binding. Signatures of active transcription, such as occupancy of RNA polymerase II and the modified histones H3K4me2, H3K4me3, and H3ac, are overrepresented in regions around bidirectional promoters, suggesting that a higher fraction of divergent genes are transcribed in a given cell than the fraction of other genes. Accordingly, analysis of whole-genome microarray data indicates that 68% of divergent genes are transcribed compared with 44% of all human genes. By combining the analysis of publicly available ENCODE data and a detailed study of GABPA, we survey bidirectional promoters with breadth and depth, leading to biological insights concerning their motif composition and bidirectional regulatory mode.

Detection and characterisation of strong resistance to phosphine in Brazilian Rhyzopertha dominica (F.) (Coleoptera: Bostrychidae)

As failure to control Rhyzopertha dominica (F.) with phosphine is a common problem in the grain-growing regions of Brazil, a study was undertaken to investigate the frequency, distribution and strength of phosphine resistance in R. dominica in Brazil. Nineteen samples of R. dominica were collected between 1991 and 2003 from central storages where phosphine fumigation had failed to control this species. Insects were cultured without selection until testing in 2005. Each sample was tested for resistance to phosphine on the basis of the response of adults to discriminating concentrations of phosphine (20 and 48 h exposures) and full dose-response assays (48 h exposure). Responses of the Brazilian R. dominica samples were compared with reference susceptible, weak-resistance and strong-resistance strains from Australia in parallel assays. All Brazilian population samples showed resistance to phosphine: five were diagnosed with weak resistance and 14 with strong resistance. Five samples showed levels of resistance similar to the reference strong-resistance strain. A representative highly resistant sample was characterised by exposing mixed-age cultures to a range of constant concentrations of phosphine for various exposure periods. Time to population extinction (TPE) and time to 99.9% suppression of population (LT(99.9)) values of this sample were generally similar to those of the reference strong-resistance strain. For example, at 0.1, 0.5 and 1.0 mg L(-1), LT(99.9) values for BR33 and the reference strong-resistance strain were respectively 21, 6.4 and 3.7 days and 17, 6.2 and 3.8 days. With both strains, doubling phosphine concentrations to 2 mg L(-1) resulted in increased LT(99.9) and TPE. High level and frequency of resistance in all population samples, some of which had been cultured without selection for up to 12 years, suggest little or no fitness deficit associated with phosphine resistance. The present research indicates that widespread phosphine resistance may be developing in Brazil. Fumigation practices should be monitored and resistance management plans implemented to alleviate further resistance development.

Disruption of Iron Homeostasis Increases Phosphine Toxicity in Caenorhabditis elegans

The aim of this study is to identify the biochemical mechanism of phosphine toxicity and resistance, using Caenorhabditis elegans as a model organism. To date, the precise mode of phosphine action is unclear. In this report, we demonstrate the following dose-dependent actions of phosphine, in vitro: (1) reduction of ferric iron (Fe3+) to ferrous iron (Fe2+), (2) release of iron from horse ferritin, (3) and the peroxidation of lipid as a result of iron release from ferritin. Using in situ hybridization, we show that the ferritin genes of C. elegans, both ferritin-1 and ferritin-2, are expressed along the digestive tract with greatest expression at the proximal and distal ends. Basal expression of the ferritin-2 gene, as determined by quantitative PCR, is approximately 80 times that of ferritin-1. However, transcript levels of ferritin-1 are induced at least 20-fold in response to phosphine, whereas there is no change in the level of ferritin-2. This resembles the reported pattern of ferritin gene regulation by iron, suggesting that phosphine toxicity may be related to an increase in the level of free iron. Indeed, iron overload increases phosphine toxicity in C. elegans at least threefold. Moreover, we demonstrate that suppression of ferritin-2 gene expression by RNAi, significantly increases sensitivity to phosphine. This study identifies similarities between phosphine toxicity and iron overload and demonstrates that phosphine can trigger iron release from storage proteins, increasing lipid peroxidation, leading to cell injury and/or cell death.

Evaluation of external RNA controls for the assessment of microarray performance

External RNA controls (ERCs), although important for microarray assay performance assessment, have yet to be fully implemented in the research community. As part of the MicroArray Quality Control (MAQC) study, two types of ERCs were implemented and evaluated; one was added to the total RNA in the samples before amplification and labeling; the other was added to the copyRNAs (cRNAs) before hybridization. ERC concentration-response curves were used across multiple commercial microarray platforms to identify problematic assays and potential sources of variation in the analytical process. In addition, the behavior of different ERC types was investigated, resulting in several important observations, such as the sample-dependent attributes of performance and the potential of using these control RNAs in a combinatorial fashion. This multiplatform investigation of the behavior and utility of ERCs provides a basis for articulating specific recommendations for their future use in evaluating assay performance across multiple platforms.

Performance comparison of one-color and two-color platforms within the MicroArray Quality Control (MAQC) project

Microarray-based expression profiling experiments typically use either a one-color or a two-color design to measure mRNA abundance. The validity of each approach has been amply demonstrated. Here we provide a simultaneous comparison of results from one- and two-color labeling designs, using two independent RNA samples from the Microarray Quality Control (MAQC) project, tested on each of three different microarray platforms. The data were evaluated in terms of reproducibility, specificity, sensitivity and accuracy to determine if the two approaches provide comparable results. For each of the three microarray platforms tested, the results show good agreement with high correlation coefficients and high concordance of differentially expressed gene lists within each platform. Cumulatively, these comparisons indicate that data quality is essentially equivalent between the one- and two-color approaches and strongly suggest that this variable need not be a primary factor in decisions regarding experimental microarray design.

The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements

Over the last decade, the introduction of microarray technology has had a profound impact on gene expression research. The publication of studies with dissimilar or altogether contradictory results, obtained using different microarray platforms to analyze identical RNA samples, has raised concerns about the reliability of this technology. The MicroArray Quality Control (MAQC) project was initiated to address these concerns, as well as other performance and data analysis issues. Expression data on four titration pools from two distinct reference RNA samples were generated at multiple test sites using a variety of microarray-based and alternative technology platforms. Here we describe the experimental design and probe mapping efforts behind the MAQC project. We show intraplatform consistency across test sites as well as a high level of interplatform concordance in terms of genes identified as differentially expressed. This study provides a resource that represents an important first step toward establishing a framework for the use of microarrays in clinical and regulatory settings.

The structure of H5N1 avian influenza neuraminidase suggests new opportunities for drug design

The worldwide spread of H5N1 avian influenza has raised concerns that this virus might acquire the ability to pass readily among humans and cause a pandemic. Two anti-influenza drugs currently being used to treat infected patients are oseltamivir (Tamiflu) and zanamivir (Relenza), both of which target the neuraminidase enzyme of the virus. Reports of the emergence of drug resistance make the development of new anti-influenza molecules a priority. Neuraminidases from influenza type A viruses form two genetically distinct groups: group-1 contains the N1 neuraminidase of the H5N1 avian virus and group-2 contains the N2 and N9 enzymes used for the structure-based design of current drugs. Here we show by X-ray crystallography that these two groups are structurally distinct. Group-1 neuraminidases contain a cavity adjacent to their active sites that closes on ligand binding. Our analysis suggests that it may be possible to exploit the size and location of the group-1 cavity to develop new anti-influenza drugs.