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.     

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.      

Dynamic compression counteracts IL-1β induced inducible nitric oxide synthase and cyclo-oxygenase-2 expression in chondrocyte/agarose constructs

Background  

Nitric oxide and prostaglandin E₂ (PGE₂play pivotal roles in both the pathogenesis of osteoarthritis and catabolic processes in articular cartilage. These mediators are influenced by both IL-1β and mechanical loading, and involve alterations in the inducible nitric oxide synthase (iNOS) and cyclo-oxygenase (COX)-2 enzymes. To identify the specific interactions that are activated by both types of stimuli, we examined the effects of dynamic compression on levels of expression of iNOS and COX-2 and involvement of the p38 mitogen-activated protein kinase (MAPK) pathway.     

Focusing on the genetics of hearing: you ain't heard nothin' yet

The complexity of genetic pathways for hearing is beginning to be amenable to unraveling by systematic functional genomic analysis. Genome-wide mutagenesis studies in the mouse are beginning to shed further light on the structure and regulation of the machinery of hearing.     

Crystal structure of the inclusion complex of the antibacterial agent triclosan with cyclomaltoheptaose and NMR study of its molecular encapsulation in positively and negatively charged cyclomaltoheptaose derivatives

The inclusion complexes of triclosan with native cyclomaltoheptaose (beta-cyclodextrin, betaCD) as well as with negatively and positively charged derivatives are studied. The structure of the inclusion complex betaCD/triclosan in the crystalline state [P1, a=15.189(5), b=15.230(6), c=16.293(6), alpha=91.07(4), beta=91.05(3) gamma=100.71(3)] comprises two crystallographically independent host macrocycles A and B. The packing results in betaCD dimers that align head-to-head and form infinite channels along the c-axis. Only one guest molecule statistically disordered over two positions, (the dichlorophenyl ring in the cavities of either A or B) corresponds to each dimer (a 2:1 host/guest complex). The enclosed dichlorophenyl ring enters the dimer through the primary side, whereas the hydrophilic chlorophenol ring extends in the space between dimers. Water molecules in five positions are also enclosed in the intradimer region, arranged on a plane perpendicular to the sevenfold axis of betaCD. The NMR spectroscopic studies in aqueous solution show the presence of both 1:1 and 2:1 betaCD/triclosan complexes. In the first case, two different 1:1 complexes are simultaneously present, each with either ring entering the narrow primary side of one betaCD molecule. In the 2:1 complex both rings of triclosan are included in two independent betaCD hosts, a precursor to the supramolecular arrangement found in the crystalline form. In the case of the negatively charged sodium heptakis[6-deoxy-6-(3-thiopropionate)]-betaCD, the NMR studies at pH 7.9 show a complete inclusion of triclosan inside the host in two orientations, one for the non-ionized (phenol) and reverse for the ionized (phenolate) form. Finally, for the positively charged heptakis(6-aminoethylamino-6-deoxy)-betaCD, inclusion of triclosan is possible only when the pH is raised to 10 and it is concluded that both aromatic rings are alternatively inside the cavity. However in that case also, inclusion of the entire guest in the elongated cavity is suggested.     

Monitoring ungulates in steep non-forest habitat: a comparison of faecal pellet and helicopter counts

Faecal pellet counts have been widely used to monitor the abundances of introduced ungulates in New Zealand, but ground-based sampling cannot be conducted safely in the steep non-forest habitats that are common in New Zealand's Southern Alps. Helicopter counts may be an effective technique for monitoring ungulates in steep non-forest habitat. We evaluated the relationship between faecal pellet and helicopter counts of ungulates (primarily feral goat Capra hircus) at 12 non-forest sites in the Southern Alps. Within each site we counted the numbers of ungulates from a helicopter on three occasions and the number of intact faecal pellets along 30 transects. Mean observed densities of feral goats derived from helicopter counts ranged from 0.0 to 20.2 km^?2. There was a positive curvilinear (concave down) relationship between faecal pellet and helicopter counts. Compared with faecal pellet counts, helicopter counts were cheaper, could identify ungulate species and provided estimates of absolute density. Helicopter counts are a cost-effective method for monitoring ungulates in the steep non-forest habitats of New Zealand's Southern Alps.