seekCRIT: Detecting and characterizing differentially expressed circular RNAs using high-throughput sequencing data.

Over the past two decades, researchers have discovered a special form of alternative splicing that produces a circular form of RNA. Although these circular RNAs (circRNAs) have garnered considerable attention in the scientific community for their biogenesis and functions, the focus of current studies has been on the tissue-specific circRNAs that exist only in one tissue but not in other tissues or on the disease-specific circRNAs that exist in certain disease conditions, such as cancer, but not under normal conditions. This approach was conducted in the relative absence of methods that analyze a group of common circRNAs that exist in both conditions, but are more abundant in one condition relative to another (differentially expressed). Studies of differentially expressed circRNAs (DECs) between two conditions would serve as a significant first step in filling this void. Here, we introduce a novel computational tool, seekCRIT (seek for differentially expressed CircRNAs In Transcriptome), that identifies the DECs between two conditions from high-throughput sequencing data. Using rat retina RNA-seq data from ischemic and normal conditions, we show that over 74% of identifiable circRNAs are expressed in both conditions and over 40 circRNAs are differentially expressed between two conditions. We also obtain a high qPCR validation rate of 90% for DECs with a FDR of < 5%. Our results demonstrate that seekCRIT is a novel and efficient approach to detect DECs using rRNA depleted RNA-seq data. seekCRIT is freely downloadable at https://github.com/UofLBioinformatics/seekCRIT. The source code is licensed under the MIT License. seekCRIT is developed and tested on Linux CentOS-7.

Choice of library size normalization and statistical methods for differential gene expression analysis in balanced two-group comparisons for RNA-seq studies.

BACKGROUND: High-throughput RNA sequencing (RNA-seq) has evolved as an important analytical tool in molecular biology. Although the utility and importance of this technique have grown, uncertainties regarding the proper analysis of RNA-seq data remain. Of primary concern, there is no consensus regarding which normalization and statistical methods are the most appropriate for analyzing this data. The lack of standardized analytical methods leads to uncertainties in data interpretation and study reproducibility, especially with studies reporting high false discovery rates. In this study, we compared a recently developed normalization method, UQ-pgQ2, with three of the most frequently used alternatives including RLE (relative log estimate), TMM (Trimmed-mean M values) and UQ (upper quartile normalization) in the analysis of RNA-seq data. We evaluated the performance of these methods for gene-level differential expression analysis by considering the factors, including: 1) normalization combined with the choice of a Wald test from DESeq2 and an exact test/QL (Quasi-likelihood) F-Test from edgeR; 2) sample sizes in two balanced two-group comparisons; and 3) sequencing read depths. RESULTS: Using the MAQC RNA-seq datasets with small sample replicates, we found that UQ-pgQ2 normalization combined with an exact test can achieve better performance in term of power and specificity in differential gene expression analysis. However, using an intra-group analysis of false positives from real and simulated data, we found that a Wald test performs better than an exact test when the number of sample replicates is large and that a QL F-test performs the best given sample sizes of 5, 10 and 15 for any normalization. The RLE, TMM and UQ methods performed similarly given a desired sample size. CONCLUSION: We found the UQ-pgQ2 method combined with an exact test/QL F-test is the best choice in order to control false positives when the sample size is small. When the sample size is large, UQ-pgQ2 with a QL F-test is a better choice for the type I error control in an intra-group analysis. We observed read depths have a minimal impact for differential gene expression analysis based on the simulated data.

Sample size calculations for the differential expression analysis of RNA-seq data using a negative binomial regression model.

High throughput RNA sequencing (RNA-seq) technology is increasingly used in disease-related biomarker studies. A negative binomial distribution has become the popular choice for modeling read counts of genes in RNA-seq data due to over-dispersed read counts. In this study, we propose two explicit sample size calculation methods for RNA-seq data using a negative binomial regression model. To derive these new sample size formulas, the common dispersion parameter and the size factor as an offset via a natural logarithm link function are incorporated. A two-sided Wald test statistic derived from the coefficient parameter is used for testing a single gene at a nominal significance level 0.05 and multiple genes at a false discovery rate 0.05. The variance for the Wald test is computed from the variance-covariance matrix with the parameters estimated from the maximum likelihood estimates under the unrestricted and constrained scenarios. The performance and a side-by-side comparison of our new formulas with three existing methods with a Wald test, a likelihood ratio test or an exact test are evaluated via simulation studies. Since other methods are much computationally extensive, we recommend our M1 method for quick and direct estimation of sample sizes in an experimental design. Finally, we illustrate sample sizes estimation using an existing breast cancer RNA-seq data.

Urinary levels of the acrolein conjugates of carnosine are associated with inhaled toxicants.

OBJECTIVE: The inhalation of air-borne toxicants is associated with adverse health outcomes which can be somewhat mitigated by enhancing endogenous anti-oxidant capacity. Carnosine is a naturally occurring dipeptide (ß-alanine-L-histidine), present in high abundance in skeletal and cardiac muscle. This multi-functional dipeptide has anti-oxidant properties, can buffer intracellular pH, chelate metals, and sequester aldehydes such as acrolein. Due to these chemical properties, carnosine may be protective against inhaled pollutants which can contain metals and aldehydes and can stimulate the generation of electrophiles in exposed tissues. Thus, assessment of carnosine levels, or levels of its acrolein conjugates (carnosine-propanal and carnosine-propanol) may inform on level of exposure and risk assessment. METHODS: We used established mass spectroscopy methods to measure levels of urinary carnosine (n = 605) and its conjugates with acrolein (n = 561) in a subset of participants in the Louisville Healthy Heart Study (mean age = 51 ± 10; 52% male). We then determined associations between these measures and air pollution exposure and smoking behavior using statistical modeling approaches. RESULTS: We found that higher levels of non-conjugated carnosine, carnosine-propanal, and carnosine-propanol were significantly associated with males (p < 0.02) and those of Caucasian ethnicity (p < 0.02). Levels of carnosine-propanol were significantly higher in never-smokers (p = 0.001) but lower in current smokers (p = 0.037). This conjugate also demonstrated a negative association with mean-daily particulate air pollution (PM2.5) levels (p = 0.01). CONCLUSIONS: These findings suggest that urinary levels of carnosine-propanol may inform as to risk from inhaled pollutants.

Regulatory networks in retinal ischemia-reperfusion injury.

BACKGROUND: Retinal function is ordered by interactions between transcriptional and posttranscriptional regulators at the molecular level. These regulators include transcription factors (TFs) and posttranscriptional factors such as microRNAs (miRs). Some studies propose that miRs predominantly target the TFs rather than other types of protein coding genes and such studies suggest a possible interconnection of these two regulators in co-regulatory networks. RESULTS: Our lab has generated mRNA and miRNA microarray expression data to investigate time-dependent changes in gene expression, following induction of ischemia-reperfusion (IR) injury in the rat retina. Data from different reperfusion time points following retinal IR-injury were analyzed. Paired expression data for miRNA-target gene (TG), TF-TG, miRNA-TF were used to identify regulatory loop motifs whose expressions were altered by the IR injury paradigm. These loops were subsequently integrated into larger regulatory networks and biological functions were assayed. Systematic analyses of the networks have provided new insights into retinal gene regulation in the early and late periods of IR. We found both overlapping and unique patterns of molecular expression at the two time points. These patterns can be defined by their characteristic molecular motifs as well as their associated biological processes. We highlighted the regulatory elements of miRs and TFs associated with biological processes in the early and late phases of ischemia-reperfusion injury. CONCLUSIONS: The etiology of retinal ischemia-reperfusion injury is orchestrated by complex and still not well understood gene networks. This work represents the first large network analysis to integrate miRNA and mRNA expression profiles in context of retinal ischemia. It is likely that an appreciation of such regulatory networks will have prognostic potential. In addition, the computational framework described in this study can be used to construct miRNA-TF interactive systems networks for various diseases/disorders of the retina and other tissues.

Ischemia-reperfusion injury of the retina is linked to necroptosis via the ERK1/2-RIP3 pathway.

PURPOSE: Ischemia-reperfusion (IR) injury is involved in the pathology of many retinal disorders since it contributes to the death of retinal neurons and the subsequent decline in vision. We determined the molecular patterns of some of the principal molecules involved in necroptosis and investigated whether IR retinal injury is associated with the extracellular signal-regulated kinase-1/2- receptor-interacting protein kinase 3 (ERK1/2-RIP3) pathway. METHODS: The cellular and subcellular localization of molecules involved in the cell death pathway, including RAGE, ERK1/2, FLIP, and RIP3, was determined with immunohistochemistry of cryosections of IR-injured retina from 2-month-old Long Evans rats. The total and phosphorylated protein levels were analyzed with quantitative western blots. ERK1/2 activity was inhibited by intravitreal injection of U0126, a highly selective inhibitor of mitogen-activated protein kinase 1/2 (MEK1/2). RESULTS: The IR-injured rat retinas expressed two RAGE isoforms with different intracellular localizations at early time points after injury. At that time point, frame inhibition of ERK activation decreased RIP3 accumulation and cell death. FLIP was detected in the IR-injured rat retinas at early time points after ischemia reperfusion. CONCLUSIONS: We report that the necroptotic cell death mechanism is executed by an ERK1/2-RIP3 pathway in the retinal ganglion cells in early stages after IR injury. Inhibition of ERK1/2 activity increased retinal ganglion cell (RGC) survival indicating that targeting of this pathway within the initial 12 h after IR injury can be used to inhibit the necroptosis pathway. We also provide evidence that a novel RAGE isoform is expressed in the early stages in rat retinal RGCs.

Transcriptomic analysis of the zebrafish inner ear points to growth hormone mediated regeneration following acoustic trauma.

BACKGROUND: Unlike mammals, teleost fishes are capable of regenerating sensory inner ear hair cells that have been lost following acoustic or ototoxic trauma. Previous work indicated that immediately following sound exposure, zebrafish saccules exhibit significant hair cell loss that recovers to pre-treatment levels within 14 days. Following acoustic trauma in the zebrafish inner ear, we used microarray analysis to identify genes involved in inner ear repair following acoustic exposure. Additionally, we investigated the effect of growth hormone (GH) on cell proliferation in control zebrafish utricles and saccules, since GH was significantly up-regulated following acoustic trauma. RESULTS: Microarray analysis, validated with the aid of quantitative real-time PCR, revealed several genes that were highly regulated during the process of regeneration in the zebrafish inner ear. Genes that had fold changes of ≥ 1.4 and P -values ≤ 0.05 were considered significantly regulated and were used for subsequent analysis. Categories of biological function that were significantly regulated included cancer, cellular growth and proliferation, and inflammation. Of particular significance, a greater than 64-fold increase in growth hormone (gh1) transcripts occurred, peaking at 2 days post-sound exposure (dpse) and decreasing to approximately 5.5-fold by 4 dpse. Pathway Analysis software was used to reveal networks of regulated genes and showed how GH affected these networks. Subsequent experiments showed that intraperitoneal injection of salmon growth hormone significantly increased cell proliferation in the zebrafish inner ear. Many other gene transcripts were also differentially regulated, including heavy and light chain myosin transcripts, both of which were down-regulated following sound exposure, and major histocompatability class I and II genes, several of which were significantly regulated on 2 dpse. CONCLUSIONS: Transcripts for GH, MHC Class I and II genes, and heavy- and light-chain myosins, as well as many others genes, were differentially regulated in the zebrafish inner ear following overexposure to sound. GH injection increased cell proliferation in the inner ear of non-sound-exposed zebrafish, suggesting that GH could play an important role in sensory hair cell regeneration in the teleost ear.

Characterization of the transcripts and protein isoforms for cytoplasmic polyadenylation element binding protein-3 (CPEB3) in the mouse retina.

BACKGROUND: Cytoplasmic polyadenylation element binding proteins (CPEBs) regulate translation by binding to regulatory motifs of defined mRNA targets. This translational mechanism has been shown to play a critical role in oocyte maturation, early development, and memory formation in the hippocampus. Little is known about the presence or functions of CPEBs in the retina. The purpose of the current study is to investigate the alternative splicing isoforms of a particular CPEB, CPEB3, based on current databases, and to characterize the expression of CPEB3 in the retina. RESULTS: In this study, we have characterized CPEB3, whose putative role is to regulate the translation of GluR2 mRNA. We identify the presence of multiple alternative splicing isoforms of CPEB3 transcripts and proteins in the current databases. We report the presence of eight alternative splicing patterns of CPEB3, including a novel one, in the mouse retina. All but one of the patterns appear to be ubiquitous in 13 types of tissue examined. The relative abundance of the patterns in the retina is demonstrated. Experimentally, we show that CPEB3 expression is increased in a time-dependent manner during the course of postnatal development, and CPEB3 is localized mostly in the inner retina, including retinal ganglion cells. CONCLUSION: The level of CPEB3 was up-regulated in the retina during development. The presence of multiple CPEB3 isoforms indicates remarkable complexity in the regulation and function of CPEB3.

A combined approach with gene-wise normalization improves the analysis of RNA-seq data in human breast cancer subtypes.

Breast cancer (BC) is increasing in incidence and resistance to treatment worldwide. The challenges in limited therapeutic options and poor survival outcomes in BC subtypes persist because of its molecular heterogeneity and resistance to standard endocrine therapy. Recently, high throughput RNA sequencing (RNA-seq) has been used to identify biomarkers of disease progression and signaling pathways that could be amenable to specific therapies according to the BC subtype. However, there is no single generally accepted pipeline for the analysis of RNA-seq data in biomarker discovery due, in part, to the needs of simultaneously satisfying constraints of sensitivity and specificity. We proposed a combined approach using gene-wise normalization, UQ-pgQ2, followed by a Wald test from DESeq2. Our approach improved the analysis based on within-group comparisons in terms of the specificity when applied to publicly available RNA-seq BC datasets. In terms of identifying differentially expressed genes (DEGs), we combined an optimized log2 fold change cutoff with a nominal false discovery rate of 0.05 to further minimize false positives. Using this method in the analysis of two GEO BC datasets, we identified 797 DEGs uniquely expressed in triple negative BC (TNBC) and significantly associated with T cell and immune-related signaling, contributing to the immunotherapeutic efficacy in TNBC patients. In contrast, we identified 1403 DEGs uniquely expressed in estrogen positive and HER2 negative BC (ER+HER2-BC) and significantly associated with eicosanoid, notching and FAK signaling while a common set of genes was associated with cellular growth and proliferation. Thus, our approach to control for false positives identified two distinct gene expression profiles associated with these two subtypes of BC which are distinguishable by their molecular and functional attributes.

Analysis of probe level patterns in Affymetrix microarray data.

BACKGROUND: Microarrays have been used extensively to analyze the expression profiles for thousands of genes in parallel. Most of the widely used methods for analyzing Affymetrix Genechip microarray data, including RMA, GCRMA and Model Based Expression Index (MBEI), summarize probe signal intensity data to generate a single measure of expression for each transcript on the array. In contrast, other methods are applied directly to probe intensities, negating the need for a summarization step. RESULTS: In this study, we used the Affymetrix rat genome Genechip to explore variability in probe response patterns within transcripts. We considered a number of possible sources of variability in probe sets including probe location within the transcript, middle base pair of the probe sequence, probe overlap, sequence homology and affinity. Although affinity, middle base pair and probe location effects may be seen at the gross array level, these factors only account for a small proportion of the variation observed at the gene level. A BLAST search and the presence of probe by treatment interactions for selected differentially expressed genes showed high sequence homology for many probes to non-target genes. CONCLUSION: We suggest that examination and modeling of probe level intensities can be used to guide researchers in refining their conclusions regarding differentially expressed genes. We discuss implications for probe sequence selection for confirmatory analysis using real time PCR.

CaMKIIalphaB mediates a survival response in retinal ganglion cells subjected to a glutamate stimulus.

PURPOSE: During N-methyl-d-aspartate-induced cell death in the neural retina, levels of the nuclear isoform of CaMKIIalpha, CaMKIIalphaB, previously reported to be detected only in the midbrain and diencephalon, become elevated. The purpose of this study was to investigate whether CaMKIIalphaB is present specifically in retinal ganglion cells (RGCs) and to determine whether it can be implicated in the cell death or cell survival of signal transduction pathways. METHODS: Pan-purified RGCs were obtained from the retinas of postnatal day (P)6 to P8 Sprague-Dawley rats. The expression level of CaMKIIalphaB was investigated in RGCs with the aid of RT-PCR and immunostaining under normal and glutamate-stressed conditions. siRNA targeted to CaMKIIalphaB was used to knock down the level of endogenous mRNA in RGCs, and cell viability was tested. The putative role of CaMKIIalphaB in the downstream expression of survival genes such as BDNF was evaluated in CaMKIIalphaB knocked-down RGCs with the aid of RT-PCR, real-time PCR, and immunofluorescence microscopy. RESULTS: Basal levels of CaMKIIalphaB were expressed in RGCs. Expression levels became increased in response to glutamate treatment and were translocated to the nuclei after a glutamate stimulus. In pan-purified RGCs with knocked down levels of CaMKIIalphaB, a glutamate stimulus led to an increase in cell death. When CaMKIIalphaB was knocked down in RGCs, a corresponding decrease occurred in the level of BDNF expression. CONCLUSIONS: These data indicate that the presence of basal levels of CaMKIIalphaB in RGCs may afford them some ongoing protection from a stressful environment. In response to the glutamate stimulus, the expression of survival genes such as BDNF may be enhanced through elevation of this particular isoform of the CaMKIIalpha gene.

Microarray reveals complement components are regulated in the serum-deprived rat retinal ganglion cell line.

PURPOSE: Glaucoma is a progressive eye disease that leads to blindness due to loss of retinal ganglion cells (RGCs). There are difficulties in using primary cultures of purified RGC to study this pathophysiology. RGC-5, a transformed not RGC line, expresses several markers characteristic of the RGCs. The aim of this study was to generate a genome-wide gene expression of RGC-5 following serum deprivation and to identify candidate genes that may be involved in the signal transduction pathways. METHODS: Apoptosis in the transformed rat RGC-5 was induced by serum deprivation for 0, 8, 24, 48, and 96 h. Briefly, 400 ng of RNA from each sample was reverse transcribed and labeled with Cy3 dye. Fragmented fluorescent cRNA was mixed with hybridization buffer and incubated at 60 degrees C for 16 h. Labeled cRNA was hybridized to Rat Genome Oligonucleotide Arrays. These arrays contain 22,775 transcripts with one oligonucleotide per transcript (60-mer). Gene expression from scanned images was quantified and analyzed using ArrayVision software. Reproducibility among triplicate arrays was determined by ANOVA statistical analysis. Significant differences in gene expression between apoptotic and nonapoptotic cells were determined based on p-values. RESULTS: Of the 22,775 transcripts present on the arrays (Agilent rat genome, 60-mer), 713 (8 h), 1,967 (24 h), 1,011 (48 h), and 1,161 (96 h) were differentially expressed relative to the 0 h time point (p-values <0.05). Twenty-three transcripts were common to 8, 24, 48, and 96 h and 130 transcripts were common to the 24, 48, and 96 h time points. The two most highly upregulated genes were Fdft1 and Lgals3 (8 h), C3 and Fcgrt (24 h), C and Lcn2 (48 h), and Mgp and C3 (96 h). A subset of the differentially expressed genes identified in microarray data (Ftl1, C3, C1s, Neu1, Polr2g, Acadm, Nupr1, Gch, Dia1, DNase1, Tgfb2, and Cyr61) were validated using quantitative real time polymerase chain reaction (QRT-PCR). Here we show that complement factor H (CFH), the major inhibitor of the alternative complement pathway is downregulated in serum-deprived RGC-5. CFH protein was detected within RGC-5 cells as well as the rat retina with the aid of immunocytochemistry and confocal microscopy. CONCLUSIONS: This study was undertaken to generate a genome-wide gene expression profile of RGC-5 after serum deprivation, and to identify candidate and novel genes that may be involved in the signal transduction pathways leading to apoptosis. RGC-5 serum deprivation revealed up-and downregulation in gene expression profiles. The data gathered from this study was the first report that the genes identified in microarray data and validated by real-time RT-PCR may play an important role in RGC-5 cell death. Among the validated genes, C3 and C1s showed significant upregulation of the complement component pathway. The results further indicate that components of the complement pathway are present in neurons of the rat retina. The data indicated that complement factors are likely involved in the pathway leading to ganglion cell death in the serum-deprivation paradigm, which may be similar to the mechanism of cell death in glaucoma.

A causal mediation model of ischemia reperfusion injury in the retina.

The goal of this study is to develop a model that explains the relationship between microRNAs, transcription factors, and their co-target genes. This relationship was previously reported in gene regulatory loops associated with 24 hour (24h) and 7 day (7d) time periods following ischemia-reperfusion injury in a rat's retina. Using a model system of retinal ischemia-reperfusion injury, we propose that microRNAs first influence transcription factors, which in turn act as mediators to influence transcription of genes via triadic regulatory loops. Analysis of the relative contributions of direct and indirect regulatory influences on genes revealed that a substantial fraction of the regulatory loops (69% for 24 hours and 77% for 7 days) could be explained by causal mediation. Over 40% of the mediated loops in both time points were regulated by transcription factors only, while about 20% of the loops were regulated entirely by microRNAs. The remaining fractions of the mediated regulatory loops were cooperatively mediated by both microRNAs and transcription factors. The results from these analyses were supported by the patterns of expression of the genes, transcription factors, and microRNAs involved in the mediated loops in both post-ischemic time points. Additionally, network motif detection for the mediated loops showed a handful of time specific motifs related to ischemia-reperfusion injury in a rat's retina. In summary, the effects of microRNAs on genes are mediated, in large part, via transcription factors.

A comparison of per sample global scaling and per gene normalization methods for differential expression analysis of RNA-seq data.

Normalization is an essential step with considerable impact on high-throughput RNA sequencing (RNA-seq) data analysis. Although there are numerous methods for read count normalization, it remains a challenge to choose an optimal method due to multiple factors contributing to read count variability that affects the overall sensitivity and specificity. In order to properly determine the most appropriate normalization methods, it is critical to compare the performance and shortcomings of a representative set of normalization routines based on different dataset characteristics. Therefore, we set out to evaluate the performance of the commonly used methods (DESeq, TMM-edgeR, FPKM-CuffDiff, TC, Med UQ and FQ) and two new methods we propose: Med-pgQ2 and UQ-pgQ2 (per-gene normalization after per-sample median or upper-quartile global scaling). Our per-gene normalization approach allows for comparisons between conditions based on similar count levels. Using the benchmark Microarray Quality Control Project (MAQC) and simulated datasets, we performed differential gene expression analysis to evaluate these methods. When evaluating MAQC2 with two replicates, we observed that Med-pgQ2 and UQ-pgQ2 achieved a slightly higher area under the Receiver Operating Characteristic Curve (AUC), a specificity rate > 85%, the detection power > 92% and an actual false discovery rate (FDR) under 0.06 given the nominal FDR (≤0.05). Although the top commonly used methods (DESeq and TMM-edgeR) yield a higher power (>93%) for MAQC2 data, they trade off with a reduced specificity (<70%) and a slightly higher actual FDR than our proposed methods. In addition, the results from an analysis based on the qualitative characteristics of sample distribution for MAQC2 and human breast cancer datasets show that only our gene-wise normalization methods corrected data skewed towards lower read counts. However, when we evaluated MAQC3 with less variation in five replicates, all methods performed similarly. Thus, our proposed Med-pgQ2 and UQ-pgQ2 methods perform slightly better for differential gene analysis of RNA-seq data skewed towards lowly expressed read counts with high variation by improving specificity while maintaining a good detection power with a control of the nominal FDR level.

The role of CaMKII in BDNF-mediated neuroprotection of retinal ganglion cells (RGC-5).

The purpose of the study is to determine if expression or secretion of brain-derived neurotrophic factor (BDNF) in retinal ganglion cells (RGC-5) is mediated by NFkappaB or Ca2+/calmodulin-dependent protein kinase II (CaMKII). RGC-5 cells were exposed to 1 mM glutamate for various periods of time, in the presence or absence of prospective regulatory molecules. BDNF mRNA and protein expression were assessed with the aid of real-time PCR and immunoblots, respectively, and BDNF secretion was determined by ELISA. The NFkappaB inhibitor (TLCK and PTD-p65), or a specific CaMKII inhibitor (m-AIP), was used to study association of NFkappaB or CaMKII with BDNF expression/secretion in RGC-5 cells. Glutamate stimulated a transient increase in BDNF mRNA and protein in RGC-5 cells, and also stimulated an early release of BDNF into the culture media. Neutralizing the BDNF or blocking the TrkB receptor enhanced the glutamate-induced cytotoxicity. NFkappaB nuclear translocation was revealed in response to glutamate treatment. Application of TLCK or PTD-p65 inhibited the glutamate-induced BDNF expression and secretion. Inhibition of CaMKII by m-AIP did not affect expression but significantly enhanced the release of BDNF from glutamate challenged cells. Our data suggest that glutamate treatment may stimulate expression of BDNF in RGC-5 cells through NFkappaB activation. A novel mechanism for neuroprotection is proposed for the CaMKII inhibitor, AIP, which appears to protect RGC-5 cells from cytotoxicity by enhancing the release of BDNF from glutamate challenged cells.

Building a glaucoma interaction network using a text mining approach.

BACKGROUND: The volume of biomedical literature and its underlying knowledge base is rapidly expanding, making it beyond the ability of a single human being to read through all the literature. Several automated methods have been developed to help make sense of this dilemma. The present study reports on the results of a text mining approach to extract gene interactions from the data warehouse of published experimental results which are then used to benchmark an interaction network associated with glaucoma. To the best of our knowledge, there is, as yet, no glaucoma interaction network derived solely from text mining approaches. The presence of such a network could provide a useful summative knowledge base to complement other forms of clinical information related to this disease. RESULTS: A glaucoma corpus was constructed from PubMed Central and a text mining approach was applied to extract genes and their relations from this corpus. The extracted relations between genes were checked using reference interaction databases and classified generally as known or new relations. The extracted genes and relations were then used to construct a glaucoma interaction network. Analysis of the resulting network indicated that it bears the characteristics of a small world interaction network. Our analysis showed the presence of seven glaucoma linked genes that defined the network modularity. A web-based system for browsing and visualizing the extracted glaucoma related interaction networks is made available at http://neurogene.spd.louisville.edu/GlaucomaINViewer/Form1.aspx. CONCLUSIONS: This study has reported the first version of a glaucoma interaction network using a text mining approach. The power of such an approach is in its ability to cover a wide range of glaucoma related studies published over many years. Hence, a bigger picture of the disease can be established. To the best of our knowledge, this is the first glaucoma interaction network to summarize the known literature. The major findings were a set of relations that could not be found in existing interaction databases and that were found to be new, in addition to a smaller subnetwork consisting of interconnected clusters of seven glaucoma genes. Future improvements can be applied towards obtaining a better version of this network.

MPrime: efficient large scale multiple primer and oligonucleotide design for customized gene microarrays.

BACKGROUND: Enhancements in sequencing technology have recently yielded assemblies of large genomes including rat, mouse, human, fruit fly, and zebrafish. The availability of large-scale genomic and genic sequence data coupled with advances in microarray technology have made it possible to study the expression of large numbers of sequence products under several different conditions in days where traditional molecular biology techniques might have taken months, or even years. Therefore, to efficiently study a number of gene products associated with a disease, pathway, or other biological process, it is necessary to be able to design primer pairs or oligonucleotides en masse rather than using a time consuming and laborious gene-by-gene method. RESULTS: We have developed an integrated system, MPrime, in order to efficiently calculate primer pairs or specific oligonucleotides for multiple genic regions based on a keyword, gene name, accession number, or sequence fasta format within the rat, mouse, human, fruit fly, and zebrafish genomes. A set of products created for mouse housekeeping genes from MPrime-designed primer pairs has been validated using both PCR-amplification and DNA sequencing. CONCLUSION: These results indicate MPrime accurately incorporates standard PCR primer design characteristics to produce high scoring primer pairs for genes of interest. In addition, sequence similarity for a set of oligonucleotides constructed for the same set of genes indicates high specificity in oligo design.

Retinal ganglion cell death and neuroprotection: Involvement of the CaMKIIalpha gene.

The purpose of this study is to determine if calcium/calmodulin-dependent protein kinase-II (CaMKII) plays a role in neuronal cell death and if inhibition of this kinase affords some neuroprotection in the RGC-5 retinal ganglion cell line. The RGC-5 cells were treated with glutamate at various concentrations for increasing increments of time. Cytotoxicity was assayed by measuring the lactate dehydrogenase (LDH) leakage from non-viable cells and TUNEL assays. The involvement of caspase-3, Bcl-2 and caspase-8 in glutamate-induced cytotoxicity was determined by immunoblots and/or real time RT-PCR. In addition, the autocamtide-2-related inhibitory peptide (AIP), a specific inhibitor of CaMKII, was used to determine the involvement of CaMKII in glutamate-induced RGC-5 cell death. Application of increasing concentrations of glutamate to RGC-5 cells caused a dose-dependent increase in the level of cell death after 24 h. There was a glutamate-stimulated increase in the expression of caspase-8 and caspase-3 and a corresponding decrease in Bcl-2. The active fragment of caspase-3 increased in glutamate-treated cells. An early transient increase in the expression of CaMKIIalpha(B) gene and a corresponding CaMKIIalpha nuclear translocation was found in glutamate-treated cells. Treatment with AIP blocked the activation of caspase-3 and protected RGC from glutamate-mediated cell death but did not alter the glutamate-enhanced expression levels of caspase-8 or caspase-3. This report shows the likely involvement of a transcript of the CaMKIIalpha gene in the cytotoxicity response of RGC-5 cells similar to previous reports in the neural retina. AIP is shown to be a neuroprotectant for RGC-5 cells as was reported for the neural retina.

Human adult olfactory neuroepithelial derived progenitors retain telomerase activity and lack apoptotic activity.

Olfactory epithelium (OE) contains a population of progenitors responsible for its life-long regenerative capacity. Procedures for the isolation of these progenitors have been established [F.J. Roisen, K.M. Klueber, C.L. Lu, L.M. Hatcher, A. Dozier, C.B. Shields, Adult human olfactory stem cells, Brain Res., 890 (2001) 11-12.] and over 40 patient-specific cell lines from adult postmortem OE and endoscopic biopsy from patients undergoing nasal sinus surgery have been obtained. As these cells emerged in primary cultures, they formed neurospheres (NSFCs). The purpose of the present study was to further characterize these adult human olfactory-derived progenitors. Subcultures of the NSFCs have been passaged nearly 200 times, with a mitotic cycle of 18-20 h. Telomerase activity remains in stem cells; therefore, ELISA was employed to determine the telomerase activity of different lines and passages. Since progenitors undergo low levels of apoptosis, the levels of apoptosis were also examined in these populations. The levels of telomerase and apoptotic activity in 12 NSFC lines remained relatively constant irrespective of donor age, culture duration, or sex. To further study the apoptotic characteristics of the NSFCs, nine different caspases (cysteine proteases) known to be critical in apoptosis were evaluated using gene-microarrays comparing cells from a single line at passages 14, 88, and 183. No increases were found in caspase activity in all passages studied. ELISA confirmed the absence of caspase activity over the entire range of passages. This study further suggests that NSFCs can be obtained and used from patients, irrespective of age, sex, or time in culture without altered viability expanding the potential utility of these cells for autologous transplantation and possible diagnostic testing.

Enhanced expression of TNF-R1 protein in NMDA-mediated cell death in the retina.

Multiple apoptosis-related genes are expressed in the retina after exposure to N-methyl-D-aspartic acid (NMDA). For example, the mRNAs for TNF-R1, FasL, and Nur77 are up-regulated between 2.8 and 7-fold [Mol. Brain Res. 91 (2001) 34-42]. The purpose of the present study is to examine prospective changes in protein expression for these genes and to determine their cellular localizations subsequent to NMDA stimulation. Following anesthesia, a single intravitreal injection of 4 mM NMDA was administered into the right eye of anesthetized rats. The left eye was injected with phosphate-buffered saline. Retinae were harvested at 2 and 24 h postinjection. Western-blot and immunocytochemical techniques were used to detect changes in protein expression levels, and to localize their distributions within the retina. Analyses of Western blots demonstrated a significant increase in TNF-R1 (100 and 80%) compared to the sham-controls at 2 and 24 h postinjection with NMDA. Immunolabeling of TNF-R1 was observed in the inner nuclear layer (INL) at 2 h postinjection with NMDA. TNF-R1 was also clearly evident in cells within the INL and ganglion cell layers (GCL) at 24 h post-injection with NMDA. In contrast to these changes in TNF-R1 there were no significant changes in the levels or distributions of FasL or Nur77 in NMDA-stimulated animals at either 2 or 24 h when compared to the sham-controls. These results implicate the TNF-R1 signal transduction pathway in NMDA-induced cell death in the INL and GCL of the retina.