An atlas of human kinase regulation

The coordinated regulation of protein kinases is a rapid mechanism that integrates diverse cues and swiftly determines appropriate cellular responses. However, our understanding of cellular decision‐making has been limited by the small number of simultaneously monitored phospho‐regulatory events. Here, we have estimated changes in activity in 215 human kinases in 399 conditions derived from a large compilation of phosphopeptide quantifications. This atlas identifies commonly regulated kinases as those that are central in the signaling network and defines the logic relationships between kinase pairs. Co‐regulation along the conditions predicts kinase–complex and kinase–substrate associations. Additionally, the kinase regulation profile acts as a molecular fingerprint to identify related and opposing signaling states. Using this atlas, we identified essential mediators of stem cell differentiation, modulators of Salmonella infection, and new targets of AKT1. This provides a global view of human phosphorylation‐based signaling and the necessary context to better understand kinase‐driven decision‐making.     

The shape of human gene family phylogenies

Background  

The shape of phylogenetic trees has been used to make inferences about the evolutionary process by comparing the shapes of actual phylogenies with those expected under simple models of the speciation process. Previous studies have focused on speciation events, but gene duplication is another lineage splitting event, analogous to speciation, and gene loss or deletion is analogous to extinction. Measures of the shape of gene family phylogenies can thus be used to investigate the processes of gene duplication and loss. We make the first systematic attempt to use tree shape to study gene duplication using human gene phylogenies.     

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.     

Non-native interactions play an effective role in protein folding dynamics

Systematic Monte Carlo simulations of simple lattice models show that the final stage of protein folding is an ordered process where native contacts get locked (i.e., the residues come into contact and remain in contact for the duration of the folding process) in a well‐defined order. The detailed study of the folding dynamics of protein‐like sequences designed as to exhibit different contact energy distributions, as well as different degrees of sequence optimization (i.e., participation of non‐native interactions in the folding process), reveals significant differences in the corresponding locking scenarios—the collection of native contacts and their average locking times, which are largely ascribable to the dynamics of non‐native contacts. Furthermore, strong evidence for a positive role played by non‐native contacts at an early folding stage was also found. Interestingly, for topologically simple target structures, a positive interplay between native and non‐native contacts is observed also toward the end of the folding process, suggesting that non‐native contacts may indeed affect the overall folding process. For target models exhibiting clear two‐state kinetics, the relation between the nucleation mechanism of folding and the locking scenario is investigated. Our results suggest that the stabilization of the folding transition state can be achieved through the establishment of a very small network of native contacts that are the first to lock during the folding process.     

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.     

Reproductive biology of Myrtillocactus geometrizans (Cactaceae) flowers with contrasting orientation

In columnar cacti, a higher production of reproductive structures on branches oriented towards the Equator has been explained by their higher interception of photosynthetic active radiation (PAR) as well as resource availability. The goal of this study was to evaluate the effect of orientation on diverse aspects of the reproductive biology of Myrtillocactus geometrizans. Phenology was studied in north- and south-facing branches. Floral cycle events, floral visitors, reproductive traits associated with sexual and attraction functions, and reproductive success were estimated from reproductive structures with contrasting orientation. Pollination experiments were conducted to evaluate the effect of orientation on mating system. Our results showed that south-facing branches had a longer duration of the mature fruit phenophase. Moreover, flower synchrony, production of reproductive structures, and floral traits associated with the male (number of anthers and pollen grains per floral bud), female (number and size of ovules and dimensions of both ovary and ovary cavity), and attraction (petal size) functions had higher values in south-facing flowers. The beginning and ending of the male function and the end of flower anthesis occurred earlier in south-facing flowers. Diversity of floral visitors was similar between orientations, except for beetles whose abundance was greater in flowers oriented towards the south. North- and south-facing flowers had a mixed mating system, with similar reproductive success. Our results showed strong differences in the reproductive biology of an intertropical columnar cactus, probably in response to the uneven PAR interception and resource availability in branches and flowers with contrasting orientation.     

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.