Linking molecular and population stress responses in Daphnia magna exposed to cadmium

DNA microarrays can be used to measure environmental stress responses. If they are to be predictive of environmental impact, we need to determine if altered gene expression translates into negative impacts on individuals and populations. A large cDNA microarray (14000 spots) was created to measure molecular stress responses to cadmium in Daphnia magna,the most widely used aquatic indicator species, and relate responses to population growth rate (pgr). We used the array to detect differences in the transcription of genes in juvenile D. magna (24 h old) after 24 h exposure to a control and three cadmium concentrations (6, 20, and 37 mu g Cd2+ L-1). Stress responses at the population level were estimated following a further 8 days exposure. Pgr was approximately linear negative with increasing cadmium concentration over this range. The microarray profile of gene expression in response to acute cadmium exposure begins to provide an overview of the molecular responses of D. magna, especially in relation to growth and development. Of the responding genes, 29% were involved with metabolism including carbohydrate, fat and peptide metabolism, and energy production, 31% were involved with transcription/translation, while 40% of responding genes were associated with cellular processes like growth and moulting, ion transport, and general stress responses (which included oxidative stress). Our production and application of a large Daphnia magna microarray has shown that measured gene responses can be logically linked to the impact of a toxicant such as cadmium on somatic growth and development, and consequently pgr.

Gene transcription in Daphnia magna: effects of acute exposure to a carbamate insecticide and an acetanilide herbicide.

Daphnia magna is a key invertebrate in the freshwater environment and is used widely as a model in ecotoxicological measurements and risk assessment. Understanding the genomic responses of D. magna to chemical challenges will be of value to regulatory authorities worldwide. Here we exposed D. magna to the insecticide methomyl and the herbicide propanil to compare phenotypic effects with changes in mRNA expression levels. Both pesticides are found in drainage ditches and surface water bodies standing adjacent to crops. Methomyl, a carbamate insecticide widely used in agriculture, inhibits acetylcholinesterase, a key enzyme in nerve transmission. Propanil, an acetanilide herbicide, is used to control grass and broad-leaf weeds. The phenotypic effects of single doses of each chemical were evaluated using a standard immobilisation assay. Immobilisation was linked to global mRNA expression levels using the previously estimated 48h-EC(1)s, followed by hybridization to a cDNA microarray with more than 13,000 redundant cDNA clones representing >5000 unique genes. Following exposure to methomyl and propanil, differential expression was found for 624 and 551 cDNAs, respectively (one-way ANOVA with Bonferroni correction, P

A specific group of genes respond to cold dehydration stress in cut Alstroemeria flowers whereas ambient dehydration stress accelerates developmental senescence expression patterns

Petal development and senescence entails a normally irreversible process. It starts with petal expansion and pigment production, and ends with nutrient remobilization and ultimately cell death. In many species this is accompanied by petal abscission. Post-harvest stress is an important factor in limiting petal longevity in cut flowers and accelerates some of the processes of senescence such as petal wilting and abscission. However, some of the effects of moderate stress in young flowers are reversible with appropriate treatments. Transcriptomic studies have shown that distinct gene sets are expressed during petal development and senescence. Despite this, the overlap in gene expression between developmental and stress-induced senescence in petals has not been fully investigated in any species. Here a custom-made cDNA microarray from Alstroemeria petals was used to investigate the overlap in gene expression between developmental changes (bud to first sign of senescence) and typical post-harvest stress treatments. Young flowers were stressed by cold or ambient temperatures without water followed by a recovery and rehydration period. Stressed flowers were still at the bud stage after stress treatments. Microarray analysis showed that ambient dehydration stress accelerates many of the changes in gene expression patterns that would normally occur during developmental senescence. However, a higher proportion of gene expression changes in response to cold stress were specific to this stimulus and not senescence related. The expression of 21 transcription factors was characterized, showing that overlapping sets of regulatory genes are activated during developmental senescence and by different stresses.

Display of a Maize cDNA library on baculovirus infected insect cells

Background: Maize is a good model system for cereal crop genetics and development because of its rich genetic heritage and well-characterized morphology. The sequencing of its genome is well advanced, and new technologies for efficient proteomic analysis are needed. Baculovirus expression systems have been used for the last twenty years to express in insect cells a wide variety of eukaryotic proteins that require complex folding or extensive posttranslational modification. More recently, baculovirus display technologies based on the expression of foreign sequences on the surface of Autographa californica ( AcMNPV) have been developed. We investigated the potential of a display methodology for a cDNA library of maize young seedlings. Results: We constructed a full-length cDNA library of young maize etiolated seedlings in the transfer vector pAcTMVSVG. The library contained a total of 2.5 x 10(5) independent clones. Expression of two known maize proteins, calreticulin and auxin binding protein (ABPI), was shown by western blot analysis of protein extracts from insect cells infected with the cDNA library. Display of the two proteins in infected insect cells was shown by selective biopanning using magnetic cell sorting and demonstrated proof of concept that the baculovirus maize cDNA display library could be used to identify and isolate proteins. Conclusion: The maize cDNA library constructed in this study relies on the novel technology of baculovirus display and is unique in currently published cDNA libraries. Produced to demonstrate proof of principle, it opens the way for the development of a eukaryotic in vivo display tool which would be ideally suited for rapid screening of the maize proteome for binding partners, such as proteins involved in hormone regulation or defence.

The tagged microarray marker (TAM) method for high-throughput detection of single nucleotide and indel polymorphisms

The tagged microarray marker (TAM) method allows high-throughput differentiation between predicted alternative PCR products. Typically, the method is used as a molecular marker approach to determining the allelic states of single nucleotide polymorphisms (SNPs) or insertion-deletion (indel) alleles at genomic loci in multiple individuals. Biotin-labeled PCR products are spotted, unpurified, onto a streptavidin-coated glass slide and the alternative products are differentiated by hybridization to fluorescent detector oligonucleotides that recognize corresponding allele-specific tags on the PCR primers. The main attractions of this method are its high throughput (thousands of PCRs are analyzed per slide), flexibility of scoring (any combination, from a single marker in thousands of samples to thousands of markers in a single sample, can be analyzed) and flexibility of scale (any experimental scale, from a small lab setting up to a large project). This protocol describes an experiment involving 3,072 PCRs scored on a slide. The whole process from the start of PCR setup to receiving the data spreadsheet takes 2 d.

Expression of target and reference genes in Daphnia magna exposed to ibuprofen

Background: Transcriptomic techniques are now being applied in ecotoxicology and toxicology to measure the impact of stressors and develop understanding of mechanisms of toxicity. Microarray technology in particular offers the potential to measure thousands of gene responses simultaneously. However, it is important that microarrays responses should be validated, at least initially, using real-time quantitative polymerase chain reaction (QPCR). The accurate measurement of target gene expression requires normalisation to an invariant internal control e. g., total RNA or reference genes. Reference genes are preferable, as they control for variation inherent in the cDNA synthesis and PCR. However, reference gene expression can vary between tissues and experimental conditions, which makes it crucial to validate them prior to application. Results: We evaluated 10 candidate reference genes for QPCR in Daphnia magna following a 24 h exposure to the non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) at 0, 20, 40 and 80 mg IB l(-1). Six of the 10 candidates appeared suitable for use as reference genes. As a robust approach, we used a combination normalisation factor (NF), calculated using the geNorm application, based on the geometric mean of three selected reference genes: glyceraldehyde-3-phosphate dehydrogenase, ubiquitin conjugating enzyme and actin. The effects of normalisation are illustrated using as target gene leukotriene B4 12-hydroxydehydrogenase (Ltb4dh), which was upregulated following 24 h exposure to 63-81 mg IB l(-1). Conclusions: As anticipated, use of the NF clarified the response of Ltb4dh in daphnids exposed to sublethal levels of ibuprofen. Our findings emphasise the importance in toxicogenomics of finding and applying invariant internal QPCR control(s) relevant to the study conditions.

GeneMCL in microarray analysis

Accurately and reliably identifying the actual number of clusters present with a dataset of gene expression profiles, when no additional information on cluster structure is available, is a problem addressed by few algorithms. GeneMCL transforms microarray analysis data into a graph consisting of nodes connected by edges, where the nodes represent genes, and the edges represent the similarity in expression of those genes, as given by a proximity measurement. This measurement is taken to be the Pearson correlation coefficient combined with a local non-linear rescaling step. The resulting graph is input to the Markov Cluster (MCL) algorithm, which is an elegant, deterministic, non-specific and scalable method, which models stochastic flow through the graph. The algorithm is inherently affected by any cluster structure present, and rapidly decomposes a graph into cohesive clusters. The potential of the GeneMCL algorithm is demonstrated with a 5730 gene subset (IGS) of the Van't Veer breast cancer database, for which the clusterings are shown to reflect underlying biological mechanisms. (c) 2005 Elsevier Ltd. All rights reserved.

Neoglycolipid probes prepared via oxime ligation for microarray analysis of oligosaccharide-protein interactions

Neoglycolipid technology is the basis of a microarray platform for assigning oligosaccharide ligands for carbohydrate-binding proteins. The strategy for generating the neoglycolipid probes by reductive amination results in ring opening of the core monosaccharides. This often limits applicability to short-chain saccharides, although the majority of recognition motifs are satisfactorily presented with neoglycolipids of longer oligosaccharides. Here, we describe neoglycolipids prepared by oxime ligation. We provide evidence from NMR studies that a significant proportion of the oxime-linked core monosaccharide is in the ring-closed form, and this form selectively interacts with a carbohydrate-binding protein. By microarray analyses we demonstrate the effective presentation with oxime-linked neoglycolipids of (1) Lewis(x) trisaccharide to antibodies to Lewisx, (2) sialyllactose analogs to the sialic acid-binding receptors, siglecs, and (3) N-glycans to a plant lectin that requires an intact N-acetylglucosamine core.

The transcriptome of Populus in elevated CO2

The consequences of increasing atmospheric carbon dioxide for long-term adaptation of forest ecosystems remain uncertain, with virtually no studies undertaken at the genetic level. A global analysis using cDNA microarrays was conducted following 6 yr exposure of Populus × euramericana (clone I-214) to elevated [CO2] in a FACE (free-air CO2 enrichment) experiment.• Gene expression was sensitive to elevated [CO2] but the response depended on the developmental age of the leaves, and < 50 transcripts differed significantly between different CO2 environments. For young leaves most differentially expressed genes were upregulated in elevated [CO2], while in semimature leaves most were downregulated in elevated [CO2].• For transcripts related only to the small subunit of Rubisco, upregulation in LPI 3 and downregulation in LPI 6 leaves in elevated CO2 was confirmed by anova. Similar patterns of gene expression for young leaves were also confirmed independently across year 3 and year 6 microarray data, and using real-time RT–PCR.• This study provides the first clues to the long-term genetic expression changes that may occur during long-term plant response to elevated CO2.

Amphioxus type I keratin cDNA and the evolution of intermediate filament genes

We report the cloning of an intermediate filament (IF) cDNA from the cephalochordate amphioxus that encodes a protein assignable to the type I keratin group. This is the first type I keratin reported from an invertebrate. Molecular phylogenetic analyses reveal that amphioxus also possesses a type II keratin, and that the genes encoding short-rod IF proteins underwent different patterns of duplication in vertebrates and their closest relatives, the cephalochordates. Extensive IF gene duplication and divergence may have facilitated the origin of new specialised cell types in vertebrates.

Phenotype MicroArray (TM) in the metabolic characterisation of Salmonella serotypes Agona, Enteritidis, Give, Hvittingfoss, Infantis, Newport and Typhimurium

The Phenotype MicroArray (TM) (PM) technology was used to study the metabolic characteristics of 29 Salmonella strains belonging to seven serotypes of S. enterica spp. enterica. Strains of serotypes Typhimurium (six strains among definite phage types DTs 1, 40 and 104) and Agona (two strains) were tested for 949 substrates, Enteritidis (six strains of phage type PT1), Give, Hvittingfoss, Infantis and Newport strains (two of each) were tested for 190 substrates and seven other Agona strains for 95 substrates. The strains represented 18 genotypes in pulsed-field gel electrophoresis (PFGE). Among 949 substrates, 18 were identified that could be used to differentiate between the strains of those seven serotypes or within a single serotype. Unique metabolic differences between the Finnish endemic Typhimurium DT1 and Agona strains were detected, for example, in the metabolism of d-tagatose, d-galactonic acid gamma-lactone and l-proline as a carbon source. Thus, the PM technique is a useful tool for identifying potential differential markers on a metabolic basis that could be used for epidemiological surveillance.

A process for analysis of microarray comparative genomics hybridisation studies for bacterial genomes

Background: Microarray based comparative genomic hybridisation (CGH) experiments have been used to study numerous biological problems including understanding genome plasticity in pathogenic bacteria. Typically such experiments produce large data sets that are difficult for biologists to handle. Although there are some programmes available for interpretation of bacterial transcriptomics data and CGH microarray data for looking at genetic stability in oncogenes, there are none specifically to understand the mosaic nature of bacterial genomes. Consequently a bottle neck still persists in accurate processing and mathematical analysis of these data. To address this shortfall we have produced a simple and robust CGH microarray data analysis process that may be automated in the future to understand bacterial genomic diversity. Results: The process involves five steps: cleaning, normalisation, estimating gene presence and absence or divergence, validation, and analysis of data from test against three reference strains simultaneously. Each stage of the process is described and we have compared a number of methods available for characterising bacterial genomic diversity, for calculating the cut-off between gene presence and absence or divergence, and shown that a simple dynamic approach using a kernel density estimator performed better than both established, as well as a more sophisticated mixture modelling technique. We have also shown that current methods commonly used for CGH microarray analysis in tumour and cancer cell lines are not appropriate for analysing our data. Conclusion: After carrying out the analysis and validation for three sequenced Escherichia coli strains, CGH microarray data from 19 E. coli O157 pathogenic test strains were used to demonstrate the benefits of applying this simple and robust process to CGH microarray studies using bacterial genomes.

Development of a miniaturised microarray-based assay for the rapid identification of antimicrobial resistance genes in Gram-negative bacteria

We describe the development of a miniaturised microarray for the detection of antimicrobial resistance genes in Gram-negative bacteria. Included on the array are genes encoding resistance to aminoglycosides, trimethoprim, sulphonamides, tetracyclines and beta-lactams, including extended-spectrum beta-lactamases. Validation of the array with control strains demonstrated a 99% correlation between polymerase chain reaction and array results. There was also good correlation between phenotypic and genotypic results for a large panel of Escherichia coli and Salmonella isolates. Some differences were also seen in the number and type of resistance genes harboured by E. coli and Salmonella strains. The array provides an effective, fast and simple method for detection of resistance genes in clinical isolates suitable for use in diagnostic laboratories, which in future will help to understand the epidemiology of isolates and to detect gene linkage in bacterial populations. (C) 2008 Published by Elsevier B.V. and the International Society of Chemotherapy.

Identification of core and variable components of the Salmonella enterica subspecies I genome by microarray

We have performed microarray hybridization studies on 40 clinical isolates from 12 common serovars within Salmonella enterica subspecies I to identify the conserved chromosomal gene pool. We were able to separate the core invariant portion of the genome by a novel mathematical approach using a decision tree based on genes ranked by increasing variance. All genes within the core component were confirmed using available sequence and microarray information for S. enterica subspecies I strains. The majority of genes within the core component had conserved homologues in Escherichia coli K-12 strain MG1655. However, many genes present in the conserved set which were absent or highly divergent in K-12 had close homologues in pathogenic bacteria such as Shigella flexneri and Pseudomonas aeruginosa. Genes within previously established virulence determinants such as SPI1 to SPI5 were conserved. In addition several genes within SPI6, all of SPI9, and three fimbrial operons (fim, bcf, and stb) were conserved within all S. enterica strains included in this study. Although many phage and insertion sequence elements were missing from the core component, approximately half the pseudogenes present in S. enterica serovar Typhi were conserved. Furthermore, approximately half the genes conserved in the core set encoded hypothetical proteins. Separation of the core and variant gene sets within S. enterica subspecies I has offered fundamental biological insight into the genetic basis of phenotypic similarity and diversity across S. enterica subspecies I and shown how the core genome of these pathogens differs from the closely related E. coli K-12 laboratory strain.