Jeanine S. Morey

Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR

Quantitative real-time PCR (qPCR) is a commonly used validation tool for confirming gene expression results obtained from microarray analysis; however, microarray and qPCR data often result in disagreement. The current study assesses factors contributing to the correlation between these methods in five separate experiments employing two-color 60-mer oligonucleotide microarrays and qPCR using SYBR green. Overall, significant correlation was observed between microarray and qPCR results (ρ=0.708, p<0.0001, n=277) using these platforms. The contribution of factors including up — vs. down-regulation, spot intensity, ρ-value, fold-change, cycle threshold (Ct), array averaging, tissue type, and tissue preparation was assessed. Filtering of microarray data for measures of quality (fold-change and ρ-value) proves to be the most critical factor, with significant correlations of ρτ;0.80 consistently observed when quality scores are applied.

MICROARRAY ANALYSIS OF DIURNAL- AND CIRCADIAN-REGULATED GENES IN THE FLORIDA RED-TIDE DINOFLAGELLATE KARENIA BREVIS (DINOPHYCEAE)1

The photoperiod plays a central role in regulating the physiology and behavior of photosynthetic phytoplankton, and many of these processes are controlled by an underlying circadian rhythm. In dinoflagellates, circadian rhythms have been shown to depend largely on posttranscriptional regulation. However, the extent to which dinoflagellates modulate transcript levels to regulate gene expression in response to light and dark has not been addressed. Here we utilized an oligonucleotide microarray containing probes to 4629 unique genes from the Florida red‐tide dinoflagellate Karenia brevis (C. C. Davis) G. Hansen et Moestrup to characterize global gene expression patterns over a 24 h day in cells exposed to a 16:8 light:dark cycle (LD treatment) and cells under constant light (LL treatment). We determined that 9.8% of the queried genes were differentially expressed during LD, while 3% of genes were differentially expressed over the 24 h day when exposed to LL. Most genes exhibited either peak or minimum expression in early dark phase. Of the 453 differentially expressed genes in LD, 104 were assigned predicted annotations based on BLASTX searches. Few of the identifiable genes that were differentially expressed in LD appear to be under circadian control, as their cyclical expression did not persist in LL. Analysis of coordinately expressed genes revealed several novel insights into dinoflagellate gene expression, including an unusually high representation of genes involved in post‐transcriptional processing of RNA and protein turnover, the regulation of PSII genes in response to light and dark, and the absence of transcriptional regulation of cell‐cycle genes.

Acute phase gene expression in mice exposed to the marine neurotoxin domoic acid

Domoic acid is a rigid analog of the neurotransmitter glutamate and a potent agonist of kainate subtype glutamate receptors. Persistent activation of these receptor subtypes results in rapid excitotoxicity, calcium dependent cell death and neuronal lesions in areas of the brain where kainate pathways are concentrated. To better understand responses to domoic acid induced excitotoxicity, microarrays were used to profile gene expression in mouse brain following domoic acid exposure. Adult female mice were subjected intraperitoneally to domoic acid at the lethal dose 50, killed and dissected at 30, 60 and 240 min post-injection. Total brain RNA from treated mice was compared with time-matched controls on Agilent 22K feature microarrays. Real-time PCR was performed on selected genes. For the 30, 60 and 240 min time points, 3.96%, 3.94% and 4.36% of the genes interrogated were differentially expressed (P-value < or = 0.01), respectively. Rigorous filtering of the data resulted in a set of 56 genes used for trending analysis and K-medians and agglomerative clustering. The earliest genes induced consisted primarily of early response gene families (Jun, Fos, Ier, Egr, growth arrest and DNA damage 45) and the inflammatory response element cyclooxygenase 2. Some later responding genes involved glucocorticoid responses (Gilz, Sgk), cold inducible proteins (Cirbp, Rbm3), Map kinases (Map3k6) and NF-kappaB inhibition. Real-time PCR in male mice from an additional study confirmed the expression of several of these genes across gender. The transcriptional profile induced by domoic acid shared similarity with expression profiles of brain ischemia and other excitotoxins, suggesting a common transcriptional response.

Transcriptional profiling and inhibition of cholesterol biosynthesis in human T lymphocyte cells by the marine toxin azaspiracid.

Azaspiracid-1 (AZA-1) is a marine biotoxin reported to accumulate in shellfish from several countries, including eastern Canada, Morocco, and much of western Europe, and is frequently associated with severe gastrointestinal human intoxication. As the mechanism of action of AZA-1 is currently unknown, human DNA microarrays and qPCR were used to profile gene expression patterns in human T lymphocyte cells following AZA-1 exposure. Some of the early (1 h) responding genes consisted of transcription factors, membrane proteins, receptors, and inflammatory genes. Four- and 24-h responding genes were dominated by genes involved in de novo lipid biosynthesis of which 17 of 18 involved in cholesterol biosynthesis were significantly up regulated. The up regulation of synthesis genes was likely in response to the ca. 50% reduction in cellular cholesterol, which correlated with up regulated protein expression levels of the low-density lipoprotein receptor. These data collectively detail the inhibition of de novo cholesterol synthesis, which is the likely cause of cytotoxicity and potentially a target pathway of the toxin.

Liver Genomic Responses to Ciguatoxin: Evidence for Activation of Phase I and Phase II Detoxification Pathways following an Acute Hypothermic Response in Mice

Ciguatoxins (CTX) are polyether neurotoxins that target voltage-gated sodium channels and are responsible for ciguatera, the most common fish-borne food poisoning in humans. This study characterizes the global transcriptional response of mouse liver to a symptomatic dose (0.26 ng/g) of the highly potent Pacific ciguatoxin-1 (P-CTX-1). At 1 h post-exposure 2.4% of features on a 44K whole genome array were differentially expressed (p < or = 0.0001), increasing to 5.2% at 4 h and decreasing to 1.4% by 24 h post-CTX exposure. Data were filtered (/fold change/ > or = 1.5 and p < or = 0.0001 in at least one time point) and a trend set of 1550 genes were used for further analysis. Early gene expression was likely influenced prominently by an acute 4 degrees C decline in core body temperature by 1 h, which resolved by 8 h following exposure. An initial downregulation of 32 different solute carriers, many involved in sodium transport, was observed. Differential gene expression in pathways involving eicosanoid biosynthesis and cholesterol homeostasis was also noted. Cytochrome P450s (Cyps) were of particular interest due to their role in xenobiotic metabolism. Twenty-seven genes, mostly members of Cyp2 and Cyp4 families, showed significant changes in expression. Many Cyps underwent an initial downregulation at 1 h but were quickly and strongly upregulated at 4 and 24 h post-exposure. In addition to Cyps, increases in several glutathione S-transferases were observed, an indication that both phase I and phase II metabolic reactions are involved in the hepatic response to CTX in mice.

Gene expression profiling in brain of mice exposed to the marine neurotoxin ciguatoxin reveals an acute anti-inflammatory, neuroprotective response

BackgroundCiguatoxins (CTXs) are polyether marine neurotoxins and potent activators of voltage-gated sodium channels. This toxin is carried by multiple reef-fish species and human consumption of ciguatoxins can result in an explosive gastrointestinal/neurologic illness. This study characterizes the global transcriptional response in mouse brain to a symptomatic dose of the highly toxic Pacific ciguatoxin P-CTX-1 and additionally compares this data to transcriptional profiles from liver and whole blood examined previously. Adult male C57/BL6 mice were injected with 0.26 ng/g P-CTX-1 while controls received only vehicle. Animals were sacrificed at 1, 4 and 24 hrs and transcriptional profiling was performed on brain RNA with Agilent whole genome microarrays. RT-PCR was used to independently validate gene expression and the web tool DAVID was used to analyze gene ontology (GO) and molecular pathway enrichment of the gene expression data.ResultsA pronounced 4°C hypothermic response was recorded in these mice, reaching a minimum at 1 hr and lasting for 8 hrs post toxin exposure. Ratio expression data were filtered by intensity, fold change and p-value, with the resulting data used for time course analysis, K-means clustering, ontology classification and KEGG pathway enrichment. Top GO hits for this gene set included acute phase response and mono-oxygenase activity. Molecular pathway analysis showed enrichment for complement/coagulation cascades and metabolism of xenobiotics. Many immediate early genes such as Fos, Jun and Early Growth Response isoforms were down-regulated although others associated with stress such as glucocorticoid responsive genes were up-regulated. Real time PCR confirmation was performed on 22 differentially expressed genes with a correlation of 0.9 (Spearmans Rho, p < 0.0001) with microarray results.ConclusionsMany of the genes differentially expressed in this study, in parallel with the hypothermia, figure prominently in protection against neuroinflammation. Pathologic activity of the complement/coagulation cascade has been shown in patients suffering from a chronic form of ciguatera poisoning and is of particular interest in this model. Anti-inflammatory processes were at work not only in the brain but were also seen in whole blood and liver of these animals, creating a systemic anti-inflammatory environment to protect against the initial cellular damage caused by the toxin.

Osmotic stress does not trigger brevetoxin production in the dinoflagellate Karenia brevis

With the global proliferation of toxic harmful algal bloom species, there is a need to identify the environmental and biological factors that regulate toxin production. One such species, Karenia brevis, forms nearly annual blooms that threaten coastal regions throughout the Gulf of Mexico. This dinoflagellate produces brevetoxins, which are potent neurotoxins that cause neurotoxic shellfish poisoning and respiratory illness in humans, as well as massive fish kills. A recent publication reported that a rapid decrease in salinity increased cellular toxin quotas in K. brevis and hypothesized that brevetoxins serve a role in osmoregulation. This finding implied that salinity shifts could significantly alter the toxic effects of blooms. We repeated the original experiments separately in three different laboratories and found no evidence for increased brevetoxin production in response to low-salinity stress in any of the eight K. brevis strains we tested, including three used in the original study. Thus, we find no support for an osmoregulatory function of brevetoxins. The original publication also stated that there was no known cellular function for brevetoxins. However, there is increasing evidence that brevetoxins promote survival of the dinoflagellates by deterring grazing by zooplankton. Whether they have other as-yet-unidentified cellular functions is currently unknown.

Subcellular localization of dinoflagellate polyketide synthases and fatty acid synthase activity

Dinoflagellates are prolific producers of polyketide secondary metabolites. Dinoflagellate polyketide synthases (PKSs) have sequence similarity to Type I PKSs, megasynthases that encode all catalytic domains on a single polypeptide. However, in dinoflagellate PKSs identified to date, each catalytic domain resides on a separate transcript, suggesting multiprotein complexes similar to Type II PKSs. Here, we provide evidence through coimmunoprecipitation that single‐domain ketosynthase and ketoreductase proteins interact, suggesting a predicted multiprotein complex. In Karenia brevis (C.C. Davis) Gert Hansen & Ø. Moestrup, previously observed chloroplast localization of PKSs suggested that brevetoxin biosynthesis may take place in the chloroplast. Here, we report that PKSs are present in both cytosol and chloroplast. Furthermore, brevetoxin is not present in isolated chloroplasts, raising the question of what chloroplast‐localized PKS enzymes might be doing. Antibodies to K. brevis PKSs recognize cytosolic and chloroplast proteins in Ostreopsis cf. ovata Fukuyo, and Coolia monotis Meunier, which produce different suites of polyketide toxins, suggesting that these PKSs may share common pathways. Since PKSs are closely related to fatty acid synthases (FAS), we sought to determine if fatty acid biosynthesis colocalizes with either chloroplast or cytosolic PKSs. [3H]acetate labeling showed fatty acids are synthesized in the cytosol, with little incorporation in chloroplasts, consistent with a Type I FAS system. However, although 29 sequences in a K. brevis expressed sequence tag database have similarity (BLASTx e‐value <10−10) to PKSs, no transcripts for either Type I (cytosolic) or Type II (chloroplast) FAS are present. Further characterization of the FAS complexes may help to elucidate the functions of the PKS enzymes identified in dinoflagellates.

Global analysis of mRNA half-lives and de novo transcription in a dinoflagellate, Karenia brevis.

Dinoflagellates possess many physiological processes that appear to be under post-transcriptional control. However, the extent to which their genes are regulated post-transcriptionally remains unresolved. To gain insight into the roles of differential mRNA stability and de novo transcription in dinoflagellates, we biosynthetically labeled RNA with 4-thiouracil to isolate newly transcribed and pre-existing RNA pools in Karenia brevis. These isolated fractions were then used for analysis of global mRNA stability and de novo transcription by hybridization to a K. brevis microarray. Global K. brevis mRNA half-lives were calculated from the ratio of newly transcribed to pre-existing RNA for 7086 array features using the online software HALO (Half-life Organizer). Overall, mRNA half-lives were substantially longer than reported in other organisms studied at the global level, ranging from 42 minutes to greater than 144 h, with a median of 33 hours. Consistent with well-documented trends observed in other organisms, housekeeping processes, including energy metabolism and transport, were significantly enriched in the most highly stable messages. Shorter-lived transcripts included a higher proportion of transcriptional regulation, stress response, and other response/regulatory processes. One such family of proteins involved in post-transcriptional regulation in chloroplasts and mitochondria, the pentatricopeptide repeat (PPR) proteins, had dramatically shorter half-lives when compared to the arrayed transcriptome. As transcript abundances for PPR proteins were previously observed to rapidly increase in response to nutrient addition, we queried the newly synthesized RNA pools at 1 and 4 h following nitrate addition to N-depleted cultures. Transcriptome-wide there was little evidence of increases in the rate of de novo transcription during the first 4 h, relative to that in N-depleted cells, and no evidence for increased PPR protein transcription. These results lend support to the growing consensus of post-transcriptional control of gene expression in dinoflagellates.

Transcriptome remodeling associated with chronological aging in the dinoflagellate, Karenia brevis

The toxic dinoflagellate, Karenia brevis, forms dense blooms in the Gulf of Mexico that persist for many months in coastal waters, where they can cause extensive marine animal mortalities and human health impacts. The mechanisms that enable cell survival in high density, low growth blooms, and the mechanisms leading to often rapid bloom demise are not well understood. To gain an understanding of processes that underlie chronological aging in this dinoflagellate, a microarray study was carried out to identify changes in the global transcriptome that accompany the entry and maintenance of stationary phase up to the onset of cell death. The transcriptome of K. brevis was assayed using a custom 10,263 feature oligonucleotide microarray from mid-logarithmic growth to the onset of culture demise. A total of 2958 (29%) features were differentially expressed, with the mid-stationary phase timepoint demonstrating peak changes in expression. Gene ontology enrichment analyses identified a significant shift in transcripts involved in energy acquisition, ribosome biogenesis, gene expression, stress adaptation, calcium signaling, and putative brevetoxin biosynthesis. The extensive remodeling of the transcriptome observed in the transition into a quiescent non-dividing phase appears to be indicative of a global shift in the metabolic and signaling requirements and provides the basis from which to understand the process of chronological aging in a dinoflagellate.