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.     

To cell cycle, swing the APC/C

For successful mitosis, Cyclin B1 and Securin must be degraded efficiently before anaphase. Destruction of these mitotic regulators by the 26S proteasome is the result of their poly-ubiquitination by a multi-subunit E3 ligase: the Anaphase-Promoting Complex or Cyclosome (APC/C). Clearly, the APC/C is not just important for mitosis. Destruction of APC/C substrates such as Cdc20, Plk1, Aurora A and Skp2 directs events in G1. Strikingly, the APC/C needs to stay active even in quiescent cells to keep them out of the cell cycle and forms an intriguing link with pRb. An inactive APC/C stabilizes Geminin, Cyclin A and Cyclin B1, thereby securing completion of DNA synthesis and progression through G2-phase. In prometaphase the APC/C becomes active again, but is controlled by the spindle assembly checkpoint. Here we discuss how the APC/C is either held in check or released. We argue that shedding more light on the APC/C is also important to understand cancer and could help the design of treatment.     

Carbon for soils,not soils for carbon

The role of soil organic carbon (SOC) sequestration as a ‘win-win’ solution to both climate change and food insecurity receives an increasing promotion. The opportunity may be too good to be missed! Yet the tremendous complexity of the two issues at stake calls for a detailed and nuanced examination of any potential solution, no matter how appealing. Here, we critically re-examine the benefits of global SOC sequestration strategies on both climate change mitigation and food production. While estimated contributions of SOC sequestration to climate change vary, almost none take SOC saturation into account. Here, we show that including saturation in estimations decreases any potential contribution of SOC sequestration to climate change mitigation by 53%–81% towards 2100. In addition, reviewing more than 21 meta-analyses, we found that observed yield effects of increasing SOC are inconsistent, ranging from negative to neutral to positive. We find that the promise of a win-win outcome is confirmed only when specific land management practices are applied under specific conditions. Therefore, we argue that the existing knowledge base does not justify the current trend to set global agendas focusing first and foremost on SOC sequestration. Away from climate-smart soils, we need a shift towards soil-smart agriculture, adaptative and adapted to each local context, and where multiple soil functions are quantified concurrently. Only such comprehensive assessments will allow synergies for land sustainability to be maximised and agronomic requirements for food security to be fulfilled. This implies moving away from global targets for SOC in agricultural soils. SOC sequestration may occur along this pathway and contribute to climate change mitigation and should be regarded as a co-benefit.     

FOSSILS AND A LARGE MOLECULAR PHYLOGENY SHOW THAT THE EVOLUTION OF SPECIES RICHNESS,GENERIC DIVERSITY,AND TURNOVER RATES ARE DISCONNECTED

The magnitude and extent of global change during the Cenozoic is remarkable, yet the impacts of these global changes on the biodiversity and evolutionary dynamics of species diversification remain poorly understood. To investigate this question, we combine paleontological and neontological data for the angiosperm order Fagales, an ecologically important clade of about 1370 species of trees with an exceptional fossil record. We show differences in patterns of accumulation of generic diversity, species richness, and turnover rates for Fagales. Generic diversity evolved rapidly since the Late Cretaceous and peaked during the Eocene or Oligocene. Turnover rates were high during periods of extreme global climate change, but relatively low when the climate remained stable. Species richness accumulated gradually throughout the Cenozoic, possibly at an accelerated pace after the Middle Miocene. Species diversification occurred in new environments: Quercoids radiating in Oligocene subtropical seasonally arid habitats, Casuarinaceae in Australian pyrophytic biomes, and Betula in Late Neogene holarctic habitats. These radiations were counterbalanced by regional extinctions in Late Neogene mesic warm‐temperate forests. Thus, the overall diversification at species level is linked to regional radiations of clades with appropriate ecologies exploiting newly available habitats.     

Prognostic significance of tumor-infiltrating T-lymphocytes in primary and metastatic lesions of advanced stage ovarian cancer

Purpose  Ovarian cancer patients with intra-tumoral CD3⁺ T-lymphocytes in primary tumor tissue have a better prognosis. This study aims to analyze the presence and relative influence of three important T-lymphocyte subsets, tumor-infiltrating CD8⁺ cytotoxic T-lymphocytes (CTL), CD45R0⁺ memory T-lymphocytes, and FoxP3⁺ regulatory T-lymphocytes (Treg), in primary tumor tissue and omental metastases of patients with ovarian cancer. Experimental design  The number of CD8⁺, CD45R0⁺, and FoxP3⁺ T-lymphocytes was determined by immunohistochemistry on a tissue micro array containing ovarian tumor tissue and/or omental metastases obtained at primary debulking surgery from 306 FIGO stage I–IV ovarian cancer patients. Immunohistochemistry data were correlated to clinicopathological parameters and survival data. Results  High number of CD8⁺ CTL and a high CD8⁺/FoxP3⁺ ratio in ovarian-derived tumor tissue were associated with increased disease-specific survival and proved to be independent prognostic factors in multivariate analyses. In advanced stage patients, the presence of CD8⁺ CTL, CD45R0⁺ memory T-lymphocytes, FoxP3⁺ Treg or a high CD8⁺/FoxP3⁺ ratio in ovarian-derived tumor tissue was associated with an increased disease specific survival in univariate analysis, as was the presence of CD45R0⁺ memory T-lymphocytes and FoxP3⁺ Treg in omental metastases. Furthermore, in advanced stage patients CD8⁺ cytotoxic and FoxP3⁺ regulatory T-lymphocytes infiltrating ovarian-derived tumor tissue were independent predictors of increased prognosis. Conclusions  T-lymphocytes infiltrating primary and metastatic ovarian cancer sites are associated with improved prognosis. These associations are especially distinct in advanced stage patients, underlining the potential for immunotherapy as a broadly applicable therapeutic strategy. Ninke Leffers and Marloes J. M. Gooden contributed equally.     

RTTN Mutations Link Primary Cilia Function to Organization of the Human Cerebral Cortex

Polymicrogyria is a malformation of the developing cerebral cortex caused by abnormal organization and characterized by many small gyri and fusion of the outer molecular layer. We have identified autosomal-recessive mutations in RTTN, encoding Rotatin, in individuals with bilateral diffuse polymicrogyria from two separate families. Rotatin determines early embryonic axial rotation, as well as anteroposterior and dorsoventral patterning in the mouse. Human Rotatin has recently been identified as a centrosome-associated protein. The Drosophila melanogaster homolog of Rotatin, Ana3, is needed for structural integrity of centrioles and basal bodies and maintenance of sensory neurons. We show that Rotatin colocalizes with the basal bodies at the primary cilium. Cultured fibroblasts from affected individuals have structural abnormalities of the cilia and exhibit downregulation of BMP4, WNT5A, and WNT2B, which are key regulators of cortical patterning and are expressed at the cortical hem, the cortex-organizing center that gives rise to Cajal-Retzius (CR) neurons. Interestingly, we have shown that in mouse embryos, Rotatin colocalizes with CR neurons at the subpial marginal zone. Knockdown experiments in human fibroblasts and neural stem cells confirm a role for RTTN in cilia structure and function. RTTN mutations therefore link aberrant ciliary function to abnormal development and organization of the cortex in human individuals.     

Interaction between Sleep and Metabolism in Drosophila with Altered Octopamine Signaling

Sleep length and metabolic dysfunction are correlated, but the causal relationship between these processes is unclear. Octopamine promotes wakefulness in the fly by acting through the insulin-producing cells (IPCs) in the fly brain. To determine if insulin signaling mediates the effects of octopamine on sleep:wake behavior, we assayed flies in which insulin signaling activity was genetically altered. We found that increasing insulin signaling does not promote wake, nor does insulin appear to mediate the wake-promoting effects of octopamine. Octopamine also affects metabolism in invertebrate species, including, as we show here, Drosophila melanogaster. Triglycerides are decreased in mutants with compromised octopamine signaling and elevated in flies with increased activity of octopaminergic neurons. Interestingly, this effect is mediated at least partially by insulin, suggesting that effects of octopamine on metabolism are independent of its effects on sleep. We further investigated the relative contribution of metabolic and sleep phenotypes to the starvation response of flies with altered octopamine signaling. Hyperactivity (indicative of foraging) induced by starvation was elevated in octopamine receptor mutants, despite their high propensity for sleep, indicating that their metabolic state dictates their behavioral response under these conditions. Moreover, flies with increased octopamine signaling do not suppress sleep in response to starvation, even though they are normally hyper-aroused, most likely because of their high triglyceride levels. Together, these data suggest that observed correlations between sleep and metabolic phenotypes can result from shared molecular pathways rather than causality, and environmental conditions can lead to the dominance of one phenotype over the other.     

Diversification in evolutionary arenas—Assessment and synthesis

Understanding how and why rates of evolutionary diversification vary is a key issue in evolutionary biology, ecology, and biogeography. Evolutionary rates are the net result of interacting processes summarized under concepts such as adaptive radiation and evolutionary stasis. Here, we review the central concepts in the evolutionary diversification literature and synthesize these into a simple, general framework for studying rates of diversification and quantifying their underlying dynamics, which can be applied across clades and regions, and across spatial and temporal scales. Our framework describes the diversification rate (d) as a function of the abiotic environment (a), the biotic environment (b), and clade‐specific phenotypes or traits (c); thus, d ~ a,b,c. We refer to the four components (a–d) and their interactions collectively as the “Evolutionary Arena.” We outline analytical approaches to this framework and present a case study on conifers, for which we parameterize the general model. We also discuss three conceptual examples: the Lupinus radiation in the Andes in the context of emerging ecological opportunity and fluctuating connectivity due to climatic oscillations; oceanic island radiations in the context of island formation and erosion; and biotically driven radiations of the Mediterranean orchid genus Ophrys. The results of the conifer case study are consistent with the long‐standing scenario that low competition and high rates of niche evolution promote diversification. The conceptual examples illustrate how using the synthetic Evolutionary Arena framework helps to identify and structure future directions for research on evolutionary radiations. In this way, the Evolutionary Arena framework promotes a more general understanding of variation in evolutionary rates by making quantitative results comparable between case studies, thereby allowing new syntheses of evolutionary and ecological processes to emerge.