GEPIS--quantitative gene expression profiling in normal and cancer tissues

MOTIVATION: Expression profiling in diverse tissues is fundamental to understanding gene function as well as therapeutic target identification. The vast collection of expressed sequence tags (ESTs) and the associated tissue source information provides an attractive opportunity for studying gene expression. RESULTS: To facilitate EST-based expression analysis, we developed GEPIS (gene expression profiling in silico), a tool that integrates EST and tissue source information to compute gene expression patterns in a large panel of normal and tumor samples. We found EST-based expression patterns to be consistent with published papers as well as our own experimental results. We also built a GEPIS Regional Atlas that depicts expression characteristics of all genes in a selected genomic region. This program can be adapted for large-scale screening for genes with desirable expression patterns, as illustrated by our large-scale mining for tissue- and tumor-specific genes. AVAILABILITY: The email server version of the GEPIS application is freely available at http://share.gene.com/share/gepis. An interactive version of GEPIS will soon be freely available at http://www.cgl.ucsf.edu/Research/genentech/gepis/. The source code, modules, data and gene lists can be downloaded at http://share.gene.com/share/gepis.     

Steady-state expression of self-regulated genes

MOTIVATION: Regulatory gene networks contain generic modules such as feedback loops that are essential for the regulation of many biological functions. The study of the stochastic mechanisms of gene regulation is instrumental for the understanding of how cells maintain their expression at levels commensurate with their biological role, as well as to engineer gene expression switches of appropriate behavior. The lack of precise knowledge on the steady-state distribution of gene expression requires the use of Gillespie algorithms and Monte-Carlo approximations. Methodology: In this study, we provide new exact formulas and efficient numerical algorithms for computing/modeling the steady-state of a class of self-regulated genes, and we use it to model/compute the stochastic expression of a gene of interest in an engineered network introduced in mammalian cells. The behavior of the genetic network is then analyzed experimentally in living cells. RESULTS: Stochastic models often reveal counter-intuitive experimental behaviors, and we find that this genetic architecture displays a unimodal behavior in mammalian cells, which was unexpected given its known bimodal response in unicellular organisms. We provide a molecular rationale for this behavior, and we implement it in the mathematical picture to explain the experimental results obtained from this network.     

Classification using functional data analysis for temporal gene expression data

MOTIVATION: Temporal gene expression profiles provide an important characterization of gene function, as biological systems are predominantly developmental and dynamic. We propose a method of classifying collections of temporal gene expression curves in which individual expression profiles are modeled as independent realizations of a stochastic process. The method uses a recently developed functional logistic regression tool based on functional principal components, aimed at classifying gene expression curves into known gene groups. The number of eigenfunctions in the classifier can be chosen by leave-one-out cross-validation with the aim of minimizing the classification error. RESULTS: We demonstrate that this methodology provides low-error-rate classification for both yeast cell-cycle gene expression profiles and Dictyostelium cell-type specific gene expression patterns. It also works well in simulations. We compare our functional principal components approach with a B-spline implementation of functional discriminant analysis for the yeast cell-cycle data and simulations. This indicates comparative advantages of our approach which uses fewer eigenfunctions/base functions. The proposed methodology is promising for the analysis of temporal gene expression data and beyond. AVAILABILITY: MATLAB programs are available upon request.     

An efficient gene transfer method mediated by ultrasound and microbubbles into the kidney

BACKGROUND: Safety issues are of paramount importance in clinical human gene therapy. From this point of view, it would be better to develop a novel non-viral efficient gene transfer method. Recently, it was reported that ultrasound exposure could induce cell membrane permeabilization and enhance gene expression. METHODS: In this study, we examined the potential of ultrasound for gene transfer into the kidney. First, we transfected rat left kidney with luciferase plasmid mixed with microbubbles, Optison, to optimize the conditions (duration of ultrasound and concentration of Optison). Then, 4, 7, 14 and 21 days after gene transfer, luciferase activity was measured. Next, localization of gene expression was assessed by measuring luciferase activity and green fluorescent protein (GFP) expression. Expression of GFP plasmid was examined under a fluorescence microscope at 4 and 14 days after gene transfer. Finally, to examine the side effects of this gene transfer method, biochemical assays for aspartate aminotransferase (AST), alanine aminotransferase (ALT), blood urea nitrogen (BUN) and creatinine (Cre) were performed. RESULTS: Optison and/or ultrasound significantly enhanced the efficiency of gene transfer and expression in the kidney. Especially, 70-80% of total glomeruli could be transfected. Also, a significant dose-dependent effect of Optison was observed as assessed by luciferase assay (Optison 25%: 12.5 x 10(5) relative light units (RLU)/g tissue; 50%: 31.3 x 10(5) RLU/g tissue; 100%: 57.9 x 10(5) RLU/g tissue). GFP expression could be observed in glomeruli, tubules and interstitial area. Results of blood tests did not change significantly after gene transfer. CONCLUSIONS: Overall, an ultrasound-mediated gene transfer method with Optison enhanced the efficiency of gene transfer and expression in the rat kidney. This novel non-viral method may be useful for gene therapy for renal disease.     

Polynomial model approach for resynchronization analysis of cell-cycle gene expression data

MOTIVATION: Identification of genes expressed in a cell-cycle-specific periodical manner is of great interest to understand cyclic systems which play a critical role in many biological processes. However, identification of cell-cycle regulated genes by raw microarray gene expression data directly is complicated by the factor of synchronization loss, thus remains a challenging problem. Decomposing the expression measurements and extracting synchronized expression will allow to better represent the single-cell behavior and improve the accuracy in identifying periodically expressed genes. RESULTS: In this paper, we propose a resynchronization-based algorithm for identifying cell-cycle-related genes. We introduce a synchronization loss model by modeling the gene expression measurements as a superposition of different cell populations growing at different rates. The underlying expression profile is then reconstructed through resynchronization and is further fitted to the measurements in order to identify periodically expressed genes. Results from both simulations and real microarray data show that the proposed scheme is promising for identifying cyclic genes and revealing underlying gene expression profiles. AVAILABILITY: Contact the authors. SUPPLEMENTARY INFORMATION: Supplementary data are available at: http://dsplab.eng.umd.edu/~genomics/syn/     

Missing value estimation for DNA microarray gene expression data: local least squares imputation

MOTIVATION: Gene expression data often contain missing expression values. Effective missing value estimation methods are needed since many algorithms for gene expression data analysis require a complete matrix of gene array values. In this paper, imputation methods based on the least squares formulation are proposed to estimate missing values in the gene expression data, which exploit local similarity structures in the data as well as least squares optimization process. RESULTS: The proposed local least squares imputation method (LLSimpute) represents a target gene that has missing values as a linear combination of similar genes. The similar genes are chosen by k-nearest neighbors or k coherent genes that have large absolute values of Pearson correlation coefficients. Non-parametric missing values estimation method of LLSimpute are designed by introducing an automatic k-value estimator. In our experiments, the proposed LLSimpute method shows competitive results when compared with other imputation methods for missing value estimation on various datasets and percentages of missing values in the data. AVAILABILITY: The software is available at http://www.cs.umn.edu/~hskim/tools.html CONTACT: hpark@cs.umn.edu     

Development of an inducible gene expression system for Lactobacillus sakei

AIM: To develop an inducible gene expression system for Lactobacillus sakei, based on the regulatory system of sakacin A production. METHODS AND RESULTS: A Lactobacillus/Escherichia coli shuttle vector; pKRV3, was constructed including the signal transducing system genes of the bacteriocin sakacin A. The gusA gene fused to PsapA promoter, cloned in this vector allowed for inducible beta-glucuronidase expression in L. sakei and L. plantarum following the addition of the sakacin A inducing peptide. PsapA appeared to be a strong and tightly controlled promoter when compared with known promoters. CONCLUSION: The pKRV3 system can be used as an inducible gene expression system in lactobacilli. SIGNIFICANCE AND IMPACT OF THE STUDY: A novel, inducible gene expression system has been developed for lactic acid bacteria relevant in food fermentations.     

Web Tools for Rice Transcriptome Analyses

Gene expression databases provide profiling data for the expression of thousands of genes to researchers worldwide. Oligonucleotide microarray technology is a useful tool that has been employed to produce gene expression profiles in most species. In rice, there are five genome-wide DNA microarray platforms: NSF 45K, BGI/Yale 60K, Affymetrix, Agilent Rice 44K, and NimbleGen 390K. Presently, more than 1,700 hybridizations of microarray gene expression data are available from public microarray depositing databases such as NCBI gene expression omnibus and Arrayexpress at EBI. More processing or reformatting of public gene expression data is required for further applications or analyses. Web-based databases for expression meta-analyses are useful for guiding researchers in designing relevant research schemes. In this review, we summarize various databases for expression meta-analyses of rice genes and web tools for further applications, such as the development of co-expression network or functional gene network.     

Development of a noninvasive reporter system for gene expression in Porphyromonas gingivalis

Several reports have supported the association of Porphyromonas gingivalis with periodontal disease. Genetic studies are vital for understanding the relative importance of virulence factors in this organism. Thus, gene reporters may prove useful for the study of gene expression in this organism. We have investigated the use of the green fluorescent protein (GFP), bacterial luciferase, and bifunctional xylosidase/arabinosidase enzyme (XA) as reporters of gene expression in P. gingivalis. Fusion cassettes containing the promoterless tetracycline resistant gene [tetA(A)Q2] and the promoterless gfp, luxAB, or xa gene were placed under the control of the rgpA promoter in P. gingivalis W83 using recombinational allelic exchange. The rgpA gene encodes for an arginine-specific protease in P. gingivalis. No GFP activity was detected in P. gingivalis isogenic mutants carrying the rgpA::gfp-tetA(Q)2 fusion construct. Luciferase activity in P. gingivalis mutants carrying the rgpA::luxAB-tetA(Q)2 fusion was only detected in the presence of exogenous FMNH(2). xa gene expression in P. gingivalis with the rgpA::xa-tetA(Q)2 fusion construct was detected in crude extracts using rho-nitrophenol derivatives as substrate and on agar plates with methylumbelliferyl derivatives under long-wave ultraviolet light. This indicates that both luxAB and xa genes can be used as reporters of gene expression in P. gingivalis. However, only the xa gene can be used as a noninvasive reporter gene.     

A phosphate-regulated promoter for fine-tuned and reversible overexpression in Ostreococcus: application to circadian clock functional analysis

BACKGROUND: The green picoalga Ostreococcus tauri (Prasinophyceae), which has been described as the smallest free-living eukaryotic organism, has minimal cellular ultra-structure and a very small genome. In recent years, O. tauri has emerged as a novel model organism for systems biology approaches that combine functional genomics and mathematical modeling, with a strong emphasis on light regulated processes and circadian clock. These approaches were made possible through the implementation of a minimal molecular toolbox for gene functional analysis including overexpression and knockdown strategies. We have previously shown that the promoter of the High Affinity Phosphate Transporter (HAPT) gene drives the expression of a luciferase reporter at high and constitutive levels under constant light. METHODOLOGY/PRINCIPAL FINDINGS: Here we report, using a luciferase reporter construct, that the HAPT promoter can be finely and reversibly tuned by modulating the level and nature of phosphate in culture medium. This HAPT regulation was additionally used to analyze the circadian clock gene Time of Cab expression 1 (TOC1). The phenotype of a TOC1ox/CCA1:Luc line was reverted from arrhythmic to rhythmic simply by adding phosphate to the culture medium. Furthermore, since the time of phosphate injection had no effect on the phase of CCA1:Luc expression, this study suggests further that TOC1 is a central clock gene in Ostreococcus. CONCLUSIONS/PERSPECTIVES: We have developed a phosphate-regulated expression system that allows fine gene function analysis in Ostreococcus. Recently, there has been a growing interest in microalgae as cell factories. This non-toxic phosphate-regulated system may prove useful in tuning protein expression levels quantitatively and temporally for biotechnological applications.     

EXP-PAC: providing comparative analysis and storage of next generation gene expression data

Microarrays and more recently RNA sequencing has led to an increase in available gene expression data. How to manage and store this data is becoming a key issue. In response we have developed EXP-PAC, a web based software package for storage, management and analysis of gene expression and sequence data. Unique to this package is SQL based querying of gene expression data sets, distributed normalization of raw gene expression data and analysis of gene expression data across experiments and species. This package has been populated with lactation data in the international milk genomic consortium web portal (http://milkgenomics.org/). Source code is also available which can be hosted on a Windows, Linux or Mac APACHE server connected to a private or public network (http://mamsap.it.deakin.edu.au/~pcc/Release/EXP_PAC.html).     

A database for management of gene expression data in situ

MOTIVATION: To create a spatiotemporal atlas of Drosophila segmentation gene expression at cellular resolution. RESULTS: The expression of segmentation genes plays a crucial role in the establishment of the Drosophila body plan. Using the IBM DB2 Relational Database Management System we have designed and implemented the FlyEx database. FlyEx contains 2832 images of 14 segmentation gene expression patterns obtained from 954 embryos and 2,073,662 quantitative data records. The averaged data is available for most of segmentation genes at eight time points. FlyEx supports operations on images of gene expression patterns. The database can be used to examine the quality of data, analyze the dynamics of formation of segmentation gene expression domains, as well as estimate the variability of gene expression patterns. We also provide the capability to download data of interest. AVAILABILITY: http://urchin.spbcas.ru/flyex, http://flyex.ams.sunysb.edu/flyex     

Changes in the Pattern of Twisted Gastrulation Gene Expression Among Drosophila Species

A long-standing hypothesis posits that morphological changes may be more likely to result from changes in regulation of gene expression than from changes in the protein coding sequences of genes. We have compared the expression pattern of the twisted gastrulation (tsg) gene among five Drosophila species: D. melanogaster, D. simulans, D. subobscura, D. mojavensis, and D. virilis. The tsg gene encodes a secreted protein that is required for the specification of dorsal midline fates in the Drosophila early embryo. TSG is unlike other secreted growth and differentiation factors in Drosophila in that its expression pattern can be experimentally varied and still result in normal development. Because of this, its regulatory region may be freer to diverge than that of other developmental genes whose misexpression may lead to lethal defects. Thus, the tsg gene may be a good indicator of the frequency and nature of evolutionary changes affecting patterns of gene expression. Over ∼60 million years (Myr), the tsg gene has retained a dorsal-on/ventral-off pattern and a middorsal region of expression; but there have been marked changes in the middorsal domain of expression as well as the appearance/loss of other domains of expression along the anterior/posterior axis. Changes between closely related species (∼2–5 Myr since divergence) that are not reflected among more distantly related species suggest frequent changes in gene expression over evolutionary time. These changes in gene expression may serve as the raw material for eventual evolutionary changes in morphology. Received: 24 March 1997 / Accepted: 20 June 1997     

Stimulatory effect of polyethylene glycol (PEG) on gene expression in mouse liver following hydrodynamics-based transfection

BACKGROUND: Rapid intravenous injection of a large volume of plasmid DNA (pDNA), i.e. a transfection procedure based on hydrodynamics, is known to be an efficient and liver-specific method of in vivo gene delivery. However, the gene expression is transient. METHODS: We investigated the effect of addition of polyethylene glycol (PEG) to a solution of naked pDNA (luciferase) on the expression of the gene in mouse liver following transfection by the hydrodynamics-based technique. In addition, the mechanism leading to the enhancement of the gene expression was studied. RESULTS: The addition of 1% (w/v) PEG2000 to the pDNA solution enhanced the resulting gene expression in the liver. Increasing the PEG2000 concentration to more than 1 and up to 10% (w/v) rather diminished the gene expression level. By contrast, increasing the molecular weight of PEG to over 2000 up to 10 000 did not affect the level of gene expression. Histopathological and serum-chemistry examinations indicated that hydrostatic or osmotic pressure increased tissue and hepatocellular damage in a PEG-concentration-dependent manner, and resulted in a decrease in gene expression. Quantitative evaluation showed that the enhanced gene expression resulted from stabilization of the pDNA introduced into the hepatocytes and an enhancement of the transport of intact pDNA to the nucleus. CONCLUSIONS: For most gene therapy applications and gene function studies, sustained expression of the introduced gene(s) is necessary. This simple method to achieve enhanced gene expression in liver may have a great potential for a wide variety of laboratory studies in molecular and cellular biology as well as possibly for future clinical applications in humans.     

Highly efficient gene transfer with degradable poly(ester amine) based on poly(ethylene glycol) diacrylate and polyethylenimine in vitro and in vivo

BACKGROUND: Polyethylenimine (PEI) is toxic although it is one of the most successful and widely used gene delivery polymers with the aid of the proton sponge effect. Therefore, development of new novel gene delivery carriers having high efficiency with less toxicity is necessary. METHODS: In this study, a degradable poly(ester amine) carrier based on poly(ethylene glycol) diacrylate (PEGDA) and low molecular weight linear PEI was prepared. Furthermore, we compared the gene expression of the polymer/DNA complexes using two delivery methods: intravenous administration as an invasive method and aerosol as a non-invasive method. RESULTS: The synthesized polymer had a relatively small molecular weight (MW = 7980) with 25 h half-life in vitro. The polymer/DNA complexes were formed at an N/P ratio of 9. The particle sizes and zeta-potentials of the complexes were dependent on N/P ratio. Compared to PEI 25K, the newly synthesized polymer exhibited high transfection efficiency with low toxicity. Poly(ester amine)-mediated gene expression in the lung and liver was higher than that of the conventional PEI carrier. Interestingly, non-invasive aerosol delivery induced higher gene expression in all organs compared to intravenous method in an in vivo mice study. Such an expressed gene via a single aerosol administration in the lung and liver remained unchanged for 7 days. CONCLUSIONS: Our study demonstrates that poly(ester amine) may be applied as an useful gene carrier.     

Chitosan-g-PEG/DNA complexes deliver gene to the rat liver via intrabiliary and intraportal infusions

BACKGROUND: Chitosan has been shown to be a non-toxic and efficient vector for in vitro gene transfection and in vivo gene delivery through pulmonary and oral administrations. Recently, we have shown that chitosan/DNA nanoparticles could mediate high levels of gene expression following intrabiliary infusion 1. In this study, we have examined the possibility of using polyethylene glycol (PEG)-grafted chitosan/DNA complexes to deliver genes to the liver through bile duct and portal vein infusions. METHODS: PEG (Mw: 5 kDa) was grafted onto chitosan (Mw: 47 kDa, deacetylation degree: 94%) with grafting degrees of 3.6% and 9.6% (molar percentage of chitosan monosaccharide units grafted with PEG). The stability of chitosan-g-PEG/DNA complexes was studied by measuring the change in particle size and by agarose gel electrophoresis against bile or serum challenge. The influence of PEG grafting on gene transfection efficiency was evaluated in HepG2 cells using luciferase reporter gene. Chitosan and chitosan-g-PEG/DNA complexes were delivered to the liver through bile duct and portal vein infusions with a syringe pump. Gene expression in the liver and the distribution of gene expression in other organs were evaluated. The acute liver toxicity of chitosan and chitosan-g-PEG/DNA complexes was examined by measuring serum alanine aminotranferase (ALT) and aspartate aminotransferase (AST) activities as a function of time. RESULTS: Both chitosan and chitosan-g-PEG displayed comparable gene transfection efficiency in HepG2 cells. After challenge with serum and bile, chitosan-g-PEG/DNA complexes, especially those prepared with chitosan-g-PEG (GD = 9.6%), did not form large aggregates like chitosan/DNA complexes but remained stable for up to 30 min. In addition, chitosan-g-PEG prevented the degradation of DNA in the presence of serum and bile. On day 3 after bile duct infusion, chitosan-g-PEG (GD = 9.6%)/DNA complexes mediated three times higher gene expression in the liver than chitosan/DNA complexes and yielded background levels of gene expression in other organs. On day 1 following portal vein infusion, gene expression level induced by chitosan/DNA complexes was hardly detectable but chitosan-g-PEG (GD = 9.6%) mediated significant transgene expression. Interestingly, transgene expression by chitosan-g-PEG/DNA complexes in other organs after portal vein infusion increased with increasing grafting degree of PEG. The ALT and AST assays indicated that grafting of PEG to chitosan reduced the acute liver toxicity towards the complexes. CONCLUSION: This study demonstrated the potential of chitosan-g-PEG as a safe and more stable gene carrier to the liver.