tigaR: integrative significance analysis of temporal differential gene expression induced by genomic abnormalities

Background

To determine which changes in the host cell genome are crucial for cervical carcinogenesis, a longitudinal in vitro model system of HPV-transformed keratinocytes was profiled in a genome-wide manner. Four cell lines affected with either HPV16 or HPV18 were assayed at 8 sequential time points for gene expression (mRNA) and gene copy number (DNA) using high-resolution microarrays. Available methods for temporal differential expression analysis are not designed for integrative genomic studies.

Results

Here, we present a method that allows for the identification of differential gene expression associated with DNA copy number changes over time. The temporal variation in gene expression is described by a generalized linear mixed model employing low-rank thin-plate splines. Model parameters are estimated with an empirical Bayes procedure, which exploits integrated nested Laplace approximation for fast computation. Iteratively, posteriors of hyperparameters and model parameters are estimated. The empirical Bayes procedure shrinks multiple dispersion-related parameters. Shrinkage leads to more stable estimates of the model parameters, better control of false positives and improvement of reproducibility. In addition, to make estimates of the DNA copy number more stable, model parameters are also estimated in a multivariate way using triplets of features, imposing a spatial prior for the copy number effect.

Conclusion

With the proposed method for analysis of time-course multilevel molecular data, more profound insight may be gained through the identification of temporal differential expression induced by DNA copy number abnormalities. In particular, in the analysis of an integrative oncogenomics study with a time-course set-up our method finds genes previously reported to be involved in cervical carcinogenesis. Furthermore, the proposed method yields improvements in sensitivity, specificity and reproducibility compared to existing methods. Finally, the proposed method is able to handle count (RNAseq) data from time course experiments as is shown on a real data set.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-327) contains supplementary material, which is available to authorized users.     

Nucleolar Integrity Is Required for the Maintenance of Long-Term Synaptic Plasticity

Long-term memory (LTM) formation requires new protein synthesis and new gene expression. Based on our work in Aplysia, we hypothesized that the rRNA genes, stimulation-dependent targets of the enzyme Poly(ADP-ribose) polymerase-1 (PARP-1), are primary effectors of the activity-dependent changes in synaptic function that maintain synaptic plasticity and memory. Using electrophysiology, immunohistochemistry, pharmacology and molecular biology techniques, we show here, for the first time, that the maintenance of forskolin-induced late-phase long-term potentiation (L-LTP) in mouse hippocampal slices requires nucleolar integrity and the expression of new rRNAs. The activity-dependent upregulation of rRNA, as well as L-LTP expression, are poly(ADP-ribosyl)ation (PAR) dependent and accompanied by an increase in nuclear PARP-1 and Poly(ADP) ribose molecules (pADPr) after forskolin stimulation. The upregulation of PARP-1 and pADPr is regulated by Protein kinase A (PKA) and extracellular signal-regulated kinase (ERK)—two kinases strongly associated with long-term plasticity and learning and memory. Selective inhibition of RNA Polymerase I (Pol I), responsible for the synthesis of precursor rRNA, results in the segmentation of nucleoli, the exclusion of PARP-1 from functional nucleolar compartments and disrupted L-LTP maintenance. Taken as a whole, these results suggest that new rRNAs (28S, 18S, and 5.8S ribosomal components)—hence, new ribosomes and nucleoli integrity—are required for the maintenance of long-term synaptic plasticity. This provides a mechanistic link between stimulation-dependent gene expression and the new protein synthesis known to be required for memory consolidation.     

Cognitive Behavioral Therapy and Fluoxetine as Adjuncts to Group Behavioral Therapy for Binge Eating Disorder

Objective: Although binge eating disorder is a common and distressing concomitant of obesity, it has not yet been established whether affected individuals presenting to behavioral weight control programs should receive specialized treatments to supplement standard treatment. This study was designed to examine the added benefit of two adjunctive interventions, individual cognitive behavioral therapy (CBT) and fluoxetine, offered in the context of group behavioral weight control treatment. Research Methods and Procedures: One hundred sixteen overweight/obese women and men with binge eating disorder were all assigned to receive a 16‐session group behavioral weight control treatment over 20 weeks. Simultaneously, subjects were randomly assigned to receive CBT + fluoxetine, CBT + placebo, fluoxetine, or placebo in a two‐by‐two factorial design. Outcome measures, assessed at the end of the 16‐session acute treatment phase, included binge frequency, weight, and measures of eating‐related and general psychopathology. Results: Overall, subjects showed substantial improvement in binge eating and both general and eating‐related psychopathology, but little weight loss. Subjects who received individual CBT improved more in binge frequency than did those not receiving CBT (p < 0.001), and binge abstinence was significantly more common in subjects receiving CBT vs. those who did not (62% vs. 33%, p < 0.001). Fluoxetine treatment was associated with greater reduction in depressive symptoms (p < 0.05). The 54 subjects who achieved binge abstinence improved more on all measures than the 62 subjects who did not. In particular, these subjects lost, on average, 6.2 kg compared with a gain of 0.7 kg among non‐abstainers. Discussion: Adjunctive individual CBT results in significant additional binge reduction in obese binge eaters receiving standard behavioral weight control treatment.     

Conduction of pEC22, a plasmid coding for MR.EcoT22I, mediated by a resident Tn3-like transposon, Tn5396.

pEC22 is a small plasmid that encodes the restriction-modification system MR.EcoT22I. Restriction and functional analysis of the plasmid identified the positions of genes encoding that system. The plasmid is able to be conducted by conjugal plasmids, a process mediated by a transposon contained within pEC22. This cryptic transposon, called Tn5396, was isolated from pEC22 and partially sequenced. The sequence of Tn5396 is for the most part typical of transposons of the Tn3 family and is most similar to that of Tn1000. The transposon differs from closely related transposons in that it lacks well-conserved sequences in the inverted-repeat region and has an unusually long terminal inverted repeat. Consideration of regions of internal sequence similarity in this and other transposons in the Tn3 family supports a theory of the mechanism by which the ends of Tn3-like transposons may maintain substantial identity between their inverted repeats over the course of evolutionary time.     

Isolation and complementation of mutants of Anabaena sp. strain PCC 7120 unable to grow aerobically on dinitrogen.

Mutants of Anabaena sp. strain PCC 7120 unable to grow aerobically on dinitrogen were isolated by mutagenesis with UV irradiation, followed by a period of incubation in yellow light and then by penicillin enrichment. A cosmid vector, pRL25C, containing replicons functional in Escherichia coli and in Anabaena species was constructed. DNA from wild-type Anabaena sp. strain PCC 7120 was partially digested with Sau3AI, and size-fractionated fragments about 40 kilobases (kb) in length were ligated into the phosphatase-treated unique BamHI site of pRL25C. A library of 1,054 cosmid clones was generated in E. coli DH1 bearing helper plasmid pDS4101. A derivative of conjugative plasmid RP-4 was transferred to this library by conjugation, and the library was replicated to lawns of mutant Anabaena strains with defects in the polysaccharide layer of the envelopes of the heterocysts. Mutant EF116 was complemented by five cosmids, three of which were subjected to detailed restriction mapping; a 2.8-kb fragment of DNA derived from one of the cosmids was found to complement EF116. Mutant EF113 was complemented by a single cosmid, which was also restriction mapped, and was shown to be complemented by a 4.8-kb fragment of DNA derived from this cosmid.     

Ribosomal RNA secondary structure: compensatory mutations and implications for phylogenetic analysis

Using sequence data from the 28S ribosomal RNA (rRNA) genes of selected vertebrates, we investigated the effects that constraints imposed by secondary structure have on the phylogenetic analysis of rRNA sequence data. Our analysis indicates that characters from both base-pairing regions (stems) and non-base-pairing regions (loops) contain phylogenetic information, as judged by the level of support of the phylogenetic results compared with a well-established tree based on both morphological and molecular data. The best results (the greatest level of support of well-accepted nodes) were obtained when the complete data set was used. However, some previously supported nodes were resolved using either the stem or loop bases alone. Stem bases sustain a greater number of compensatory mutations than would be expected at random, but the number is < 40% of that expected under a hypothesis of perfect compensation to maintain secondary structure. Therefore, we suggest that in phylogenetic analyses, the weighting of stem characters be reduced by no more than 20%, relative to that of loop characters. In contrast to previous suggestions, we do not recommend weighting of stem positions by one-half, compared with that of loop positions, because this overcompensates for the constraints that selection imposes on the secondary structure of rRNA.