Monitoring of malaria parasite resistance to chloroquine and sulphadoxine-pyrimethamine in the Solomon Islands by DNA microarray technology.

BACKGROUND: Little information is available on resistance to anti-malarial drugs in the Solomon Islands (SI). The analysis of single nucleotide polymorphisms (SNPs) in drug resistance associated parasite genes is a potential alternative to classical time- and resource-consuming in vivo studies to monitor drug resistance. Mutations in pfmdr1 and pfcrt were shown to indicate chloroquine (CQ) resistance, mutations in pfdhfr and pfdhps indicate sulphadoxine-pyrimethamine (SP) resistance, and mutations in pfATPase6 indicate resistance to artemisinin derivatives. METHODS: The relationship between the rate of treatment failure among 25 symptomatic Plasmodium falciparum-infected patients presenting at the clinic and the pattern of resistance-associated SNPs in P. falciparum infecting 76 asymptomatic individuals from the surrounding population was investigated. The study was conducted in the SI in 2004. Patients presenting at a local clinic with microscopically confirmed P. falciparum malaria were recruited and treated with CQ+SP. Rates of treatment failure were estimated during a 28-day follow-up period. In parallel, a DNA microarray technology was used to analyse mutations associated with CQ, SP, and artemisinin derivative resistance among samples from the asymptomatic community. Mutation and haplotype frequencies were determined, as well as the multiplicity of infection. RESULTS: The in vivo study showed an efficacy of 88% for CQ+SP to treat P. falciparum infections. DNA microarray analyses indicated a low diversity in the parasite population with one major haplotype present in 98.7% of the cases. It was composed of fixed mutations at position 86 in pfmdr1, positions 72, 75, 76, 220, 326 and 356 in pfcrt, and positions 59 and 108 in pfdhfr. No mutation was observed in pfdhps or in pfATPase6. The mean multiplicity of infection was 1.39. CONCLUSION: This work provides the first insight into drug resistance markers of P. falciparum in the SI. The obtained results indicated the presence of a very homogenous P. falciparum population circulating in the community. Although CQ+SP could still clear most infections, seven fixed mutations associated with CQ resistance and two fixed mutations related to SP resistance were observed. Whether the absence of mutations in pfATPase6 indicates the efficacy of artemisinin derivatives remains to be proven.

A microarray-based system for the simultaneous analysis of single nucleotide polymorphisms in human genes involved in the metabolism of anti-malarial drugs

Background: In order to provide a cost-effective tool to analyse pharmacogenetic markers in malaria treatment, DNA microarray technology was compared with sequencing of polymerase chain reaction (PCR) fragments to detect single nucleotide polymorphisms (SNPs) in a larger number of samples. Methods: The microarray was developed to affordably generate SNP data of genes encoding the human cytochrome P450 enzyme family (CYP) and N-acetyltransferase-2 (NAT2) involved in antimalarial drug metabolisms and with known polymorphisms, i.e. CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, and NAT2. Results: For some SNPs, i.e. CYP2A6*2, CYP2B6*5, CYP2C8*3, CYP2C9*3/*5, CYP2C19*3, CYP2D6*4 and NAT2*6/*7/*14, agreement between both techniques ranged from substantial to almost perfect (kappa index between 0.61 and 1.00), whilst for other SNPs a large variability from slight to substantial agreement (kappa index between 0.39 and 1.00) was found, e. g. CYP2D6*17 (2850C>T), CYP3A4*1B and CYP3A5*3. Conclusion: The major limit of the microarray technology for this purpose was lack of robustness and with a large number of missing data or with incorrect specificity.

Gene expression databases and data mining.

The DNA microarray technology has arguably caught the attention of the worldwide life science community and is now systematically supporting major discoveries in many fields of study. The majority of the initial technical challenges of conducting experiments are being resolved, only to be replaced with new informatics hurdles, including statistical analysis, data visualization, interpretation, and storage. Two systems of databases, one containing expression data and one containing annotation data are quickly becoming essential knowledge repositories of the research community. This present paper surveys several databases, which are considered "pillars" of research and important nodes in the network. This paper focuses on a generalized workflow scheme typical for microarray experiments using two examples related to cancer research. The workflow is used to reference appropriate databases and tools for each step in the process of array experimentation. Additionally, benefits and drawbacks of current array databases are addressed, and suggestions are made for their improvement.

Cross-species and cross-platform gene expression studies with the Bioconductor-compliant R package 'annotationTools'.

Background: The variety of DNA microarray formats and datasets presently available offers an unprecedented opportunity to perform insightful comparisons of heterogeneous data. Cross-species studies, in particular, have the power of identifying conserved, functionally important molecular processes. Validation of discoveries can now often be performed in readily available public data which frequently requires cross-platform studies.Cross-platform and cross-species analyses require matching probes on different microarray formats. This can be achieved using the information in microarray annotations and additional molecular biology databases, such as orthology databases. Although annotations and other biological information are stored using modern database models ( e. g. relational), they are very often distributed and shared as tables in text files, i.e. flat file databases. This common flat database format thus provides a simple and robust solution to flexibly integrate various sources of information and a basis for the combined analysis of heterogeneous gene expression profiles.Results: We provide annotationTools, a Bioconductor-compliant R package to annotate microarray experiments and integrate heterogeneous gene expression profiles using annotation and other molecular biology information available as flat file databases. First, annotationTools contains a specialized set of functions for mining this widely used database format in a systematic manner. It thus offers a straightforward solution for annotating microarray experiments. Second, building on these basic functions and relying on the combination of information from several databases, it provides tools to easily perform cross-species analyses of gene expression data.Here, we present two example applications of annotationTools that are of direct relevance for the analysis of heterogeneous gene expression profiles, namely a cross-platform mapping of probes and a cross-species mapping of orthologous probes using different orthology databases. We also show how to perform an explorative comparison of disease-related transcriptional changes in human patients and in a genetic mouse model.Conclusion: The R package annotationTools provides a simple solution to handle microarray annotation and orthology tables, as well as other flat molecular biology databases. Thereby, it allows easy integration and analysis of heterogeneous microarray experiments across different technological platforms or species.

The cooperative induction of hypoxia-inducible factor-1 alpha and STAT3 during hypoxia induced an impairment of tumor susceptibility to CTL-mediated cell lysis.

Hypoxia is an essential component of tumor microenvironment. In this study, we investigated the influence of hypoxia (1% PO(2)) on CTL-mediated tumor cell lysis. We demonstrate that exposure of target tumor cells to hypoxia has an inhibitory effect on the CTL clone (Heu171)-induced autologous target cell lysis. Such inhibition correlates with hypoxia-inducible factor-1alpha (HIF-1alpha) induction but is not associated with an alteration of CTL reactivity as revealed by granzyme B polarization or morphological change. Western blot analysis indicates that although hypoxia had no effect on p53 accumulation, it induced the phosphorylation of STAT3 in tumor cells by a mechanism at least in part involving vascular endothelial growth factor secretion. We additionally show that a simultaneous nuclear translocation of HIF-1alpha and phospho-STAT3 was observed. Interestingly, gene silencing of STAT3 by small interfering RNA resulted in HIF-1alpha inhibition and a significant restoration of target cell susceptibility to CTL-induced killing under hypoxic conditions by a mechanism involving at least in part down-regulation of AKT phosphorylation. Moreover, knockdown of HIF-1alpha resulted in the restoration of target cell lysis under hypoxic conditions. This was further supported by DNA microarray analysis where STAT3 inhibition resulted in a partly reversal of the hypoxia-induced gene expression profile. The present study demonstrates that the concomitant hypoxic induction of phospho-STAT3 and HIF-1alpha are functionally linked to the alteration of non-small cell lung carcinoma target susceptibility to CTL-mediated killing. Considering the eminent functions of STAT3 and HIF-1alpha in the tumor microenvironment, their targeting may represent novel strategies for immunotherapeutic intervention.

A conserved transcript pattern in response to a specialist and a generalist herbivore.

Transcript patterns elicited in response to attack reveal, at the molecular level, how plants respond to aggressors. These patterns are fashioned both by inflicted physical damage as well as by biological components displayed or released by the attacker. Different types of attacking organisms might therefore be expected to elicit different transcription programs in the host. Using a large-scale DNA microarray, we characterized gene expression in damaged as well as in distal Arabidopsis thaliana leaves in response to the specialist insect, Pieris rapae. More than 100 insect-responsive genes potentially involved in defense were identified, including genes involved in pathogenesis, indole glucosinolate metabolism, detoxification and cell survival, and signal transduction. Of these 114 genes, 111 were induced in Pieris feeding, and only three were repressed. Expression patterns in distal leaves were markedly similar to those of local leaves. Analysis of wild-type and jasmonate mutant plants, coupled with jasmonate treatment, showed that between 67 and 84% of Pieris-regulated gene expression was controlled, totally or in part, by the jasmonate pathway. This was correlated with increased larval performance on the coronatine insensitive1 glabrous1 (coi1-1 gl1) mutant. Independent mutations in COI1 and GL1 led to a faster larval weight gain, but the gl1 mutation had relatively little effect on the expression of the insect-responsive genes examined. Finally, we compared transcript patterns in Arabidopis in response to larvae of the specialist P. rapae and to a generalist insect, Spodoptera littoralis. Surprisingly, given the complex nature of insect salivary components and reported differences between species, almost identical transcript profiles were observed. This study also provides a robustly characterized gene set for the further investigation of plant-insect interaction.

Proof of concept for microarray-based detection of DNA-binding oncogenes in cell extracts.

The function of DNA-binding proteins is controlled not just by their abundance, but mainly at the level of their activity in terms of their interactions with DNA and protein targets. Moreover, the affinity of such transcription factors to their target sequences is often controlled by co-factors and/or modifications that are not easily assessed from biological samples. Here, we describe a scalable method for monitoring protein-DNA interactions on a microarray surface. This approach was designed to determine the DNA-binding activity of proteins in crude cell extracts, complementing conventional expression profiling arrays. Enzymatic labeling of DNA enables direct normalization of the protein binding to the microarray, allowing the estimation of relative binding affinities. Using DNA sequences covering a range of affinities, we show that the new microarray-based method yields binding strength estimates similar to low-throughput gel mobility-shift assays. The microarray is also of high sensitivity, as it allows the detection of a rare DNA-binding protein from breast cancer cells, the human tumor suppressor AP-2. This approach thus mediates precise and robust assessment of the activity of DNA-binding proteins and takes present DNA-binding assays to a high throughput level.

Protein-Binding Microarray Analysis of Tumor Suppressor AP2α Target Gene Specificity.

Cheap and massively parallel methods to assess the DNA-binding specificity of transcription factors are actively sought, given their prominent regulatory role in cellular processes and diseases. Here we evaluated the use of protein-binding microarrays (PBM) to probe the association of the tumor suppressor AP2α with 6000 human genomic DNA regulatory sequences. We show that the PBM provides accurate relative binding affinities when compared to quantitative surface plasmon resonance assays. A PBM-based study of human healthy and breast tumor tissue extracts allowed the identification of previously unknown AP2α target genes and it revealed genes whose direct or indirect interactions with AP2α are affected in the diseased tissues. AP2α binding and regulation was confirmed experimentally in human carcinoma cells for novel target genes involved in tumor progression and resistance to chemotherapeutics, providing a molecular interpretation of AP2α role in cancer chemoresistance. Overall, we conclude that this approach provides quantitative and accurate assays of the specificity and activity of tumor suppressor and oncogenic proteins in clinical samples, interfacing genomic and proteomic assays.

Genotyping microarray for CSNB-associated genes.

PURPOSE: Congenital stationary night blindness (CSNB) is a clinically and genetically heterogeneous retinal disease. Although electroretinographic (ERG) measurements can discriminate clinical subgroups, the identification of the underlying genetic defects has been complicated for CSNB because of genetic heterogeneity, the uncertainty about the mode of inheritance, and time-consuming and costly mutation scanning and direct sequencing approaches. METHODS: To overcome these challenges and to generate a time- and cost-efficient mutation screening tool, the authors developed a CSNB genotyping microarray with arrayed primer extension (APEX) technology. To cover as many mutations as possible, a comprehensive literature search was performed, and DNA samples from a cohort of patients with CSNB were first sequenced directly in known CSNB genes. Subsequently, oligonucleotides were designed representing 126 sequence variations in RHO, CABP4, CACNA1F, CACNA2D4, GNAT1, GRM6, NYX, PDE6B, and SAG and spotted on the chip. RESULTS: Direct sequencing of genes known to be associated with CSNB in the study cohort revealed 21 mutations (12 novel and 9 previously reported). The resultant microarray containing oligonucleotides, which allow to detect 126 known and novel mutations, was 100% effective in determining the expected sequence changes in all known samples assessed. In addition, investigation of 34 patients with CSNB who were previously not genotyped revealed sequence variants in 18%, of which 15% are thought to be disease-causing mutations. CONCLUSIONS: This relatively inexpensive first-pass genetic testing device for patients with a diagnosis of CSNB will improve molecular diagnostics and genetic counseling of patients and their families and gives the opportunity to analyze whether, for example, more progressive disorders such as cone or cone-rod dystrophies underlie the same gene defects.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

A glycine-arginine domain in control of the human MRE11 DNA repair protein.

Human MRE11 is a key enzyme in DNA double-strand break repair and genome stability. Human MRE11 bears a glycine-arginine-rich (GAR) motif that is conserved among multicellular eukaryotic species. We investigated how this motif influences MRE11 function. Human MRE11 alone or a complex of MRE11, RAD50, and NBS1 (MRN) was methylated in insect cells, suggesting that this modification is conserved during evolution. We demonstrate that PRMT1 interacts with MRE11 but not with the MRN complex, suggesting that MRE11 arginine methylation occurs prior to the binding of NBS1 and RAD50. Moreover, the first six methylated arginines are essential for the regulation of MRE11 DNA binding and nuclease activity. The inhibition of arginine methylation leads to a reduction in MRE11 and RAD51 focus formation on a unique double-strand break in vivo. Furthermore, the MRE11-methylated GAR domain is sufficient for its targeting to DNA damage foci and colocalization with gamma-H2AX. These studies highlight an important role for the GAR domain in regulating MRE11 function at the biochemical and cellular levels during DNA double-strand break repair.

Results of a collaborative study of the EDNAP group regarding mitochondrial DNA heteroplasmy and segregation in hair shafts.

A collaborative exercise was carried out by the European DNA Profiling Group (EDNAP) in order to evaluate the distribution of mitochondrial DNA (mtDNA) heteroplasmy amongst the hairs of an individual who displays point heteroplasmy in blood and buccal cells. A second aim of the exercise was to study reproducibility of mtDNA sequencing of hairs between laboratories using differing chemistries, further to the first mtDNA reproducibility study carried out by the EDNAP group. Laboratories were asked to type 2 sections from each of 10 hairs, such that each hair was typed by at least two laboratories. Ten laboratories participated in the study, and a total of 55 hairs were typed. The results showed that the C/T point heteroplasmy observed in blood and buccal cells at position 16234 segregated differentially between hairs, such that some hairs showed only C, others only T and the remainder, C/T heteroplasmy at varying ratios. Additionally, differential segregation of heteroplasmic variants was confirmed in independent extracts at positions 16093 and the poly(C) tract at 302-309, whilst a complete A-G transition was confirmed at position 16129 in one hair. Heteroplasmy was observed at position 16195 on both strands of a single extract from one hair segment, but was not observed in the extracts from any other segment of the same hair. Similarly, heteroplasmy at position 16304 was observed on both strands of a single extract from one hair. Additional variants at positions 73, 249 and the HVII poly(C) region were reported by one laboratory; as these were not confirmed in independent extracts, the possibility of contamination cannot be excluded. Additionally, the electrophoresis and detection equipment used by this laboratory was different to those of the other laboratories, and the discrepancies at position 249 and the HVII poly(C) region appear to be due to reading errors that may be associated with this technology. The results, and their implications for forensic mtDNA typing, are discussed in the light of the biology of hair formation.

Single-cell gene-expression profiling reveals qualitatively distinct CD8 T cells elicited by different gene-based vaccines.

CD8 T cells play a key role in mediating protective immunity against selected pathogens after vaccination. Understanding the mechanism of this protection is dependent upon definition of the heterogeneity and complexity of cellular immune responses generated by different vaccines. Here, we identify previously unrecognized subsets of CD8 T cells based upon analysis of gene-expression patterns within single cells and show that they are differentially induced by different vaccines. Three prime-boost vector combinations encoding HIV Env stimulated antigen-specific CD8 T-cell populations of similar magnitude, phenotype, and functionality. Remarkably, however, analysis of single-cell gene-expression profiles enabled discrimination of a majority of central memory (CM) and effector memory (EM) CD8 T cells elicited by the three vaccines. Subsets of T cells could be defined based on their expression of Eomes, Cxcr3, and Ccr7, or Klrk1, Klrg1, and Ccr5 in CM and EM cells, respectively. Of CM cells elicited by DNA prime-recombinant adenoviral (rAd) boost vectors, 67% were Eomes(-) Ccr7(+) Cxcr3(-), in contrast to only 7% and 2% stimulated by rAd5-rAd5 or rAd-LCMV, respectively. Of EM cells elicited by DNA-rAd, 74% were Klrk1(-) Klrg1(-)Ccr5(-) compared with only 26% and 20% for rAd5-rAd5 or rAd5-LCMV. Definition by single-cell gene profiling of specific CM and EM CD8 T-cell subsets that are differentially induced by different gene-based vaccines will facilitate the design and evaluation of vaccines, as well as enable our understanding of mechanisms of protective immunity.

Use of DNA profiles for investigation using a simulated national DNA database: Part I. Partial SGM Plus profiles.

In traditional criminal investigation, uncertainties are often dealt with using a combination of common sense, practical considerations and experience, but rarely with tailored statistical models. For example, in some countries, in order to search for a given profile in the national DNA database, it must have allelic information for six or more of the ten SGM Plus loci for a simple trace. If the profile does not have this amount of information then it cannot be searched in the national DNA database (NDNAD). This requirement (of a result at six or more loci) is not based on a statistical approach, but rather on the feeling that six or more would be sufficient. A statistical approach, however, could be more rigorous and objective and would take into consideration factors such as the probability of adventitious matches relative to the actual database size and/or investigator's requirements in a sensible way. Therefore, this research was undertaken to establish scientific foundations pertaining to the use of partial SGM Plus loci profiles (or similar) for investigation.