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.     

Variation in the promoter of the human hormone sensitive lipase gene shows gender specific effects on insulin and lipid levels: results from the Ely study.

We previously identified a hormone sensitive lipase (HSL) promoter variant, -60C>G, which in vitro exhibits 40% reduced promoter activity. In this study we examined the effect of the -60C>G on glycemic and lipid measures in the population based Ely study of metabolic function and insulin resistance in 218 middle-aged men and 276 middle-aged women. Adipose tissue HSL is the rate-limiting step in triglyceride lipolysis, generating free fatty acids for energy utilization. HSL is also expressed in pancreatic beta-cells where its activity therefore may affect insulin secretion. In the women, carriers of the HSL -60G allele had significantly lower fasting insulin levels (P=0.0005) and a lower total area under the curve for insulin during the oral glucose tolerance test (P=0.005). There was no demonstrable association in men with these measures of insulin sensitivity but carriers of the -60G allele had significantly lower fasting non-esterified fatty acid (NEFA) levels (P=0.025) and higher low density lipoprotein cholesterol levels (P=0.02) than men who were non-carriers. This study provides additional evidence for a role for HSL in the development of insulin resistance, from which carriers of the -60G allele, associated here with markers of insulin sensitivity in women, and with lower NEFA levels in men, might be protected.     

DNA ligase 1 deficient plants display severe growth defects and delayed repair of both DNA single and double strand breaks

Background  

DNA ligase enzymes catalyse the joining of adjacent polynucleotides and as such play important roles in DNA replication and repair pathways. Eukaryotes possess multiple DNA ligases with distinct roles in DNA metabolism, with clear differences in the functions of DNA ligase orthologues between animals, yeast and plants. DNA ligase 1, present in all eukaryotes, plays critical roles in both DNA repair and replication and is indispensable for cell viability.     

Urotensin-II System in Genetic Control of Blood Pressure and Renal Function

Urotensin-II controls ion/water homeostasis in fish and vascular tone in rodents. We hypothesised that common genetic variants in urotensin-II pathway genes are associated with human blood pressure or renal function. We performed family-based analysis of association between blood pressure, glomerular filtration and genes of the urotensin-II pathway (urotensin-II, urotensin-II related peptide, urotensin-II receptor) saturated with 28 tagging single nucleotide polymorphisms in 2024 individuals from 520 families; followed by an independent replication in 420 families and 7545 unrelated subjects. The expression studies of the urotensin-II pathway were carried out in 97 human kidneys. Phylogenetic evolutionary analysis was conducted in 17 vertebrate species. One single nucleotide polymorphism (rs531485 in urotensin-II gene) was associated with adjusted estimated glomerular filtration rate in the discovery cohort (p = 0.0005). It showed no association with estimated glomerular filtration rate in the combined replication resource of 8724 subjects from 6 populations. Expression of urotensin-II and its receptor showed strong linear correlation (r = 0.86, p<0.0001). There was no difference in renal expression of urotensin-II system between hypertensive and normotensive subjects. Evolutionary analysis revealed accumulation of mutations in urotensin-II since the divergence of primates and weaker conservation of urotensin-II receptor in primates than in lower vertebrates. Our data suggest that urotensin-II system genes are unlikely to play a major role in genetic control of human blood pressure or renal function. The signatures of evolutionary forces acting on urotensin-II system indicate that it may have evolved towards loss of function since the divergence of primates.     

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.      

Report on injuries sustained by patients treated at the Birmingham General Hospital following the recent bomb explosions.

As a result of recent bomb explosions a total of 82 patients were treated at the Birmingham General Hospital, 61 with minor injuries. Bomb injuries may be divided into three main groups due to the blast effect (such as blast lung and ruptured tympanic membranes), the flash (such as burns to the exposed part of the body), and shrapnel (which may cause a wide variety of injuries). The amount of warning of such explosions is usually minimal, and so the prepared accident schemes of most hospitals are inappropriate. If the disaster occurs outside normal working hours much responsibility initially falls on the resident staff. This report gives some idea of the type of injuries they are likely to see.