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 (C(t)), 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

Transcriptomic analysis of polyketide synthases in a highly ciguatoxic dinoflagellate, Gambierdiscus polynesiensis and low toxicity Gambierdiscus pacificus, from French Polynesia.

Marine dinoflagellates produce a diversity of polyketide toxins that are accumulated in marine food webs and are responsible for a variety of seafood poisonings. Reef-associated dinoflagellates of the genus Gambierdiscus produce toxins responsible for ciguatera poisoning (CP), which causes over 50,000 cases of illness annually worldwide. The biosynthetic machinery for dinoflagellate polyketides remains poorly understood. Recent transcriptomic and genomic sequencing projects have revealed the presence of Type I modular polyketide synthases in dinoflagellates, as well as a plethora of single domain transcripts with Type I sequence homology. The current transcriptome analysis compares polyketide synthase (PKS) gene transcripts expressed in two species of Gambierdiscus from French Polynesia: a highly toxic ciguatoxin producer, G. polynesiensis, versus a non-ciguatoxic species G. pacificus, each assembled from approximately 180 million Illumina 125 nt reads using Trinity, and compares their PKS content with previously published data from other Gambierdiscus species and more distantly related dinoflagellates. Both modular and single-domain PKS transcripts were present. Single domain β-ketoacyl synthase (KS) transcripts were highly amplified in both species (98 in G. polynesiensis, 99 in G. pacificus), with smaller numbers of standalone acyl transferase (AT), ketoacyl reductase (KR), dehydratase (DH), enoyl reductase (ER), and thioesterase (TE) domains. G. polynesiensis expressed both a larger number of multidomain PKSs, and larger numbers of modules per transcript, than the non-ciguatoxic G. pacificus. The largest PKS transcript in G. polynesiensis encoded a 10,516 aa, 7 module protein, predicted to synthesize part of the polyether backbone. Transcripts and gene models representing portions of this PKS are present in other species, suggesting that its function may be performed in those species by multiple interacting proteins. This study contributes to the building consensus that dinoflagellates utilize a combination of Type I modular and single domain PKS proteins, in an as yet undefined manner, to synthesize polyketides

RNA-Seq analysis of seasonal and individual variation in blood transcriptomes of healthy managed bottlenose dolphins

All globin depleted samples and associated physical or hematological parameters measured. (XLS 43 kb

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 ≤ 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| ≥ 1.5 and p ≤ 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°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

Modeling population effects of the Deepwater Horizon oil spill on a long-lived species

This research was enabled partly by a grant from The Gulf of Mexico Research Initiative (GOMRI).The 2010 Deepwater Horizon (DWH) oil spill exposed common bottlenose dolphins (Tursiops truncatus) in Barataria Bay, Louisiana to heavy oiling that caused increased mortality and chronic disease and impaired reproduction in surviving dolphins. We conducted photographic surveys and veterinary assessments in the decade following the spill. We assigned a prognostic score (good, fair, guarded, poor, or grave) for each dolphin to provide a single integrated indicator of overall health, and we examined temporal trends in prognostic scores. We used expert elicitation to quantify the implications of trends for the proportion of the dolphins that would recover within their lifetime. We integrated expert elicitation, along with other new information, in a population dynamics model to predict the effects of observed health trends on demography. We compared the resulting population trajectory with that predicted under baseline (no spill) conditions. Disease conditions persisted and have recently worsened in dolphins that were presumably exposed to DWH oil: 78% of those assessed in 2018 had a guarded, poor, or grave prognosis. Dolphins born after the spill were in better health. We estimated that the population declined by 45% (95% CI 14–74) relative to baseline and will take 35 years (95% CI 18–67) to recover to 95% of baseline numbers. The sum of annual differences between baseline and injured population sizes (i.e., the lost cetacean years) was 30,993 (95% CI 6607–94,148). The population is currently at a minimum point in its recovery trajectory and is vulnerable to emerging threats, including planned ecosystem restoration efforts that are likely to be detrimental to the dolphins’ survival. Our modeling framework demonstrates an approach for integrating different sources and types of data, highlights the utility of expert elicitation for indeterminable input parameters, and emphasizes the importance of considering and monitoring long-term health of long-lived species subject to environmental disasters. Article impact statement: Oil spills can have long-term consequences for the health of long-lived species; thus, effective restoration and monitoring are needed.Publisher PDFPeer reviewe

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

<p>Abstract</p> <p>Background</p> <p>Ciguatoxins (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.</p> <p>Results</p> <p>A 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 (Spearman's Rho, p < 0.0001) with microarray results.</p> <p>Conclusions</p> <p>Many 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.</p

Transcriptomic response of the red tide dinoflagellate, Karenia brevis, to nitrogen and phosphorus depletion and addition

<p>Abstract</p> <p>Background</p> <p>The role of coastal nutrient sources in the persistence of <it>Karenia brevis </it>red tides in coastal waters of Florida is a contentious issue that warrants investigation into the regulation of nutrient responses in this dinoflagellate. In other phytoplankton studied, nutrient status is reflected by the expression levels of N- and P-responsive gene transcripts. In dinoflagellates, however, many processes are regulated post-transcriptionally. All nuclear encoded gene transcripts studied to date possess a 5' <it>trans</it>-spliced leader (SL) sequence suggestive, based on the trypanosome model, of post-transcriptional regulation. The current study therefore sought to determine if the transcriptome of <it>K. brevis </it>is responsive to nitrogen and phosphorus and is informative of nutrient status.</p> <p>Results</p> <p>Microarray analysis of N-depleted <it>K. brevis </it>cultures revealed an increase in the expression of transcripts involved in N-assimilation (nitrate and ammonium transporters, glutamine synthetases) relative to nutrient replete cells. In contrast, a transcriptional signal of P-starvation was not apparent despite evidence of P-starvation based on their rapid growth response to P-addition. To study transcriptome responses to nutrient addition, the limiting nutrient was added to depleted cells and changes in global gene expression were assessed over the first 48 hours following nutrient addition. Both N- and P-addition resulted in significant changes in approximately 4% of genes on the microarray, using a significance cutoff of 1.7-fold and p ≤ 10^-4. By far, the earliest responding genes were dominated in both nutrient treatments by pentatricopeptide repeat (PPR) proteins, which increased in expression up to 3-fold by 1 h following nutrient addition. PPR proteins are nuclear encoded proteins involved in chloroplast and mitochondria RNA processing. Correspondingly, other functions enriched in response to both nutrients were photosystem and ribosomal genes.</p> <p>Conclusions</p> <p>Microarray analysis provided transcriptomic evidence for N- but not P-limitation in <it>K. brevis</it>. Transcriptomic responses to the addition of either N or P suggest a concerted program leading to the reactivation of chloroplast functions. Even the earliest responding PPR protein transcripts possess a 5' SL sequence that suggests post-transcriptional control. Given the current state of knowledge of dinoflagellate gene regulation, it is currently unclear how these rapid changes in such transcript levels are achieved.</p

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

Array features with significantly different abundance during the fourth hour post-nitrate addition, relative to N-depleted cultures.

<p>Contigs in bold are also significantly changing during the first hour post-nitrate addition.</p