Modeling dependent gene expression

In this paper we propose a Bayesian approach for inference about dependence of high throughput gene expression. Our goals are to use prior knowledge about pathways to anchor inference about dependence among genes; to account for this dependence while making inferences about differences in mean expression across phenotypes; and to explore differences in the dependence itself across phenotypes. Useful features of the proposed approach are a model-based parsimonious representation of expression as an ordinal outcome, a novel and flexible representation of prior information on the nature of dependencies, and the use of a coherent probability model over both the structure and strength of the dependencies of interest. We evaluate our approach through simulations and in the analysis of data on expression of genes in the Complement and Coagulation Cascade pathway in ovarian cancer.Comment: Published in at http://dx.doi.org/10.1214/11-AOAS525 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Generalized gene co-expression analysis via subspace clustering using low-rank representation

BACKGROUND: Gene Co-expression Network Analysis (GCNA) helps identify gene modules with potential biological functions and has become a popular method in bioinformatics and biomedical research. However, most current GCNA algorithms use correlation to build gene co-expression networks and identify modules with highly correlated genes. There is a need to look beyond correlation and identify gene modules using other similarity measures for finding novel biologically meaningful modules. RESULTS: We propose a new generalized gene co-expression analysis algorithm via subspace clustering that can identify biologically meaningful gene co-expression modules with genes that are not all highly correlated. We use low-rank representation to construct gene co-expression networks and local maximal quasi-clique merger to identify gene co-expression modules. We applied our method on three large microarray datasets and a single-cell RNA sequencing dataset. We demonstrate that our method can identify gene modules with different biological functions than current GCNA methods and find gene modules with prognostic values. CONCLUSIONS: The presented method takes advantage of subspace clustering to generate gene co-expression networks rather than using correlation as the similarity measure between genes. Our generalized GCNA method can provide new insights from gene expression datasets and serve as a complement to current GCNA algorithms

Stochastic Gene Expression in Single Gene Oscillator Variants

It is infeasible to understand all dynamics in cell, but we can aim to understand the impact of design choices under our control. Here we consider a single gene oscillator as a case study to understand the influence of DNA copy number and repressor choice on the resulting dynamics. We first switch the repressor in the oscillator from the originally published lacI to treRL, a chimeric repressor with a lacI DNA binding domain that is inducible by trehalose. This slightly modified system produces faster and more regular oscillations than the original lacI oscillator. We then compare the treRL oscillator at three different DNA copy numbers. The period and amplitude of oscillations increases as the copy number is decreased. We cannot explain the change in period with differential equation models without changing delays or degradation rates. The correlation and phase coherence between daughter cells after cell division also tend to fall off faster for the lower copy oscillator variants. These results suggest that lower copy number variants of our single gene oscillator produce more synchronized oscillations

A linear mixed model approach to gene expression-tumor aneuploidy association studies.

Aneuploidy, defined as abnormal chromosome number or somatic DNA copy number, is a characteristic of many aggressive tumors and is thought to drive tumorigenesis. Gene expression-aneuploidy association studies have previously been conducted to explore cellular mechanisms associated with aneuploidy. However, in an observational setting, gene expression is influenced by many factors that can act as confounders between gene expression and aneuploidy, leading to spurious correlations between the two variables. These factors include known confounders such as sample purity or batch effect, as well as gene co-regulation which induces correlations between the expression of causal genes and non-causal genes. We use a linear mixed-effects model (LMM) to account for confounding effects of tumor purity and gene co-regulation on gene expression-aneuploidy associations. When applied to patient tumor data across diverse tumor types, we observe that the LMM both accounts for the impact of purity on aneuploidy measurements and identifies a new association between histone gene expression and aneuploidy

Expression of Green Fluorescence Protein (GFP) in Zebrafish Muscle through Injection: A Gene Therapy Model

Expression of the target gene is important for gene therapy. Presently, localized transgenesis is used for gene therapy which can be achieved by a target gene expression. Here, we have reported the plasmid mediated gene therapy to zebrafish model. For this purpose, we have chosen green fluorescent protein (GFP) as a target gene because the expression can be detected easily. GFP was inserted in a plasmid vector, pQE30 to develop the vector pQE30GFP. The plasmid pQE30GFP was constructed form plasmid, pQE30 and pEGFPC2. pQE30GFP injected directly in one group of fish into the muscle where luciferase expression was noted. In another group, after injection electroporation was performed where we have also noted luciferase expression; but, electroporation cause muscle injury to the zebrafish. In our case, the expression was very strong at the site of injection in first group in compare to electroporation group and in both the cases expression was stable more than two weeks

Modelling the burden caused by gene expression: an in silico investigation into the interactions between synthetic gene circuits and their chassis cell

In this paper we motivate and develop a model of gene expression for the purpose of studying the interaction between synthetic gene circuits and the chassis cell within which they are in- serted. This model focuses on the translational aspect of gene expression as this is where the literature suggests the crucial interaction between gene expression and shared resources lies

Positive selection underlies Faster-Z evolution of gene expression in birds.

The elevated rate of evolution for genes on sex chromosomes compared to autosomes (Fast-X or Fast-Z evolution) can result either from positive selection in the heterogametic sex, or from non-adaptive consequences of reduced relative effective population size. Recent work in birds suggests that Fast-Z of coding sequence is primarily due to relaxed purifying selection resulting from reduced relative effective population size. However, gene sequence and gene expression are often subject to distinct evolutionary pressures, therefore we tested for Fast-Z in gene expression using next-generation RNA-sequencing data from multiple avian species. Similar to studies of Fast-Z in coding sequence, we recover clear signatures of Fast-Z in gene expression, however in contrast to coding sequence, our data indicate that Fast-Z in expression is due to positive selection acting primarily in females. In the soma, where gene expression is highly correlated between the sexes, we detected Fast-Z in both sexes, although at a higher rate in females, suggesting that many positively selected expression changes in females are also expressed in males. In the gonad, where inter-sexual correlations in expression are much lower, we detected Fast-Z for female gene expression, but crucially, not males. This suggests that a large amount of expression variation is sex-specific in its effects within the gonad. Taken together, our results indicate that Fast-Z evolution of gene expression is the product of positive selection acting on recessive beneficial alleles in the heterogametic sex. More broadly, our analysis suggests that the adaptive potential of Z chromosome gene expression may be much greater than that of gene sequence, results which have important implications for the role of sex chromosomes in speciation and sexual selection

Digital gene expression analysis of the zebra finch genome

Background: In order to understand patterns of adaptation and molecular evolution it is important to quantify both variation in gene expression and nucleotide sequence divergence. Gene expression profiling in non-model organisms has recently been facilitated by the advent of massively parallel sequencing technology. Here we investigate tissue specific gene expression patterns in the zebra finch (Taeniopygia guttata) with special emphasis on the genes of the major histocompatibility complex (MHC). Results: Almost 2 million 454-sequencing reads from cDNA of six different tissues were assembled and analysed. A total of 11,793 zebra finch transcripts were represented in this EST data, indicating a transcriptome coverage of about 65%. There was a positive correlation between the tissue specificity of gene expression and non-synonymous to synonymous nucleotide substitution ratio of genes, suggesting that genes with a specialised function are evolving at a higher rate (or with less constraint) than genes with a more general function. In line with this, there was also a negative correlation between overall expression levels and expression specificity of contigs. We found evidence for expression of 10 different genes related to the MHC. MHC genes showed relatively tissue specific expression levels and were in general primarily expressed in spleen. Several MHC genes, including MHC class I also showed expression in brain. Furthermore, for all genes with highest levels of expression in spleen there was an overrepresentation of several gene ontology terms related to immune function. Conclusions: Our study highlights the usefulness of next-generation sequence data for quantifying gene expression in the genome as a whole as well as in specific candidate genes. Overall, the data show predicted patterns of gene expression profiles and molecular evolution in the zebra finch genome. Expression of MHC genes in particular, corresponds well with expression patterns in other vertebrates

The role of spatial averaging in the precision of gene expression patterns

During embryonic development, differentiating cells respond via gene expression to positional cues from morphogen gradients. While gene expression is often highly erratic, embryonic development is precise. We show by theory and simulations that diffusion of the expressed protein can enhance the precision of its expression domain. While diffusion lessens the sharpness of the expression boundary, it also reduces super-Poissonian noise by washing out bursts of gene expression. Balancing these effects yields an optimal diffusion constant maximizing the precision of the expression domain.Comment: 5 pages, 4 EPS figure