REMOVING TECHNICAL VARIABILITY IN RNA-SEQ DATA USING CONDITIONAL QUANTILE NORMALIZATION

The ability to measure gene expression on a genome-wide scale is one of the most promising accomplishments in molecular biology. Microarrays, the technology that first permitted this, were riddled with problems due to unwanted sources of variability. Many of these problems are now mitigated, after a decade’s worth of statistical methodology development. The recently developed RNA sequencing (RNA-seq) technology has generated much excitement in part due to claims of reduced variability in comparison to microarrays. However, we show RNA-seq data demonstrates unwanted and obscuring variability similar to what was first observed in microarrays. In particular, we find GC-content has a strong sample specific effect on gene expression measurements that, if left uncorrected, leads to false positives in downstream results. We also report on commonly observed data distortions that demonstrate the need for data normalization. Here we describe statistical methodology that improves precision by 42% without loss of accuracy. Our resulting conditional quantile normalization (CQN) algorithm combines robust generalized regression to remove systematic bias introduced by deterministic features such as GC-content, and quantile normalization to correct for global distortions.

Establishment of murine cell lines by constitutive and conditional immortalization

Mouse cell lines were immortalized by introduction of specific immortalizing genes. Embryonic and adult animals and an embryonal stem cell line were used as a source of primary cells. The immortalizing genes were either introduced by DNA transfection or by ecotropic retrovirus transduction. Fibroblasts were obtained by expression of SV40 virus large T antigen (TAg). The properties of the resulting fibroblast cell lines were reproducible, independent of the donor mouse strains employed and the cells showed no transformed properties in vitro and did not form tumors in vivo. Endothelial cell lines were generated by Polyoma virus middle T antigen expression in primary embryonal cells. These cell lines consistently expressed relevant endothelial cell surface markers. Since the expression of the immortalizing genes was expected to strongly influence the cellular characteristics fibroblastoid cells were reversibly immortalized by using a vector that allows conditional expression of the TAg. Under inducing conditions, these cells exhibited properties that were highly similar to the properties of constitutively immortalized cells. In the absence of TAg expression, cell proliferation stops. Cell growth is resumed when TAg expression is restored. Gene expression profiling indicates that TAg influences the expression levels of more than 1000 genes that are involved in diverse cellular processes. The data show that conditionally immortalized cell lines have several advantageous properties over constitutively immortalized cells.

Transgenic system for conditional induction and rescue of chronic myocardial hibernation provides insights into genomic programs of hibernation

A key energy-saving adaptation to chronic hypoxia that enables cardiomyocytes to withstand severe ischemic insults is hibernation, i.e., a reversible arrest of contractile function. Whereas hibernating cardiomyocytes represent the critical reserve of dysfunctional cells that can be potentially rescued, a lack of a suitable animal model has hampered insights on this medically important condition. We developed a transgenic mouse system for conditional induction of long-term hibernation and a system to rescue hibernating cardiomyocytes at will. Via myocardium-specific induction (and, in turn, deinduction) of a VEGF-sequestering soluble receptor, we show that VEGF is indispensable for adjusting the coronary vasculature to match increased oxygen consumption and exploit this finding to generate a hypoperfused heart. Importantly, ensuing ischemia is tunable to a level at which large cohorts of cardiomyocytes are driven to enter a hibernation mode, without cardiac cell death. Relieving the VEGF blockade even months later resulted in rapid revascularization and full recovery of contractile function. Furthermore, we show that left ventricular remodeling associated with hibernation is also fully reversible. The unique opportunity to uncouple hibernation from other ischemic heart phenotypes (e.g., infarction) was used to determine the genetic program of hibernation; uncovering hypoxia-inducible factor target genes associated with metabolic adjustments and induced expression of several cardioprotective genes. Autophagy, specifically self-digestion of mitochondria, was identified as a key prosurvival mechanism in hibernating cardiomyocytes. This system may lend itself for examining the potential utility of treatments to rescue dysfunctional cardiomyocytes and reverse maladaptive remodeling.

Multivariate modelling of responses to conditional items: New possibilities for latent class analysis

Questionnaire data may contain missing values because certain questions do not apply to all respondents. For instance, questions addressing particular attributes of a symptom, such as frequency, triggers or seasonality, are only applicable to those who have experienced the symptom, while for those who have not, responses to these items will be missing. This missing information does not fall into the category 'missing by design', rather the features of interest do not exist and cannot be measured regardless of survey design. Analysis of responses to such conditional items is therefore typically restricted to the subpopulation in which they apply. This article is concerned with joint multivariate modelling of responses to both unconditional and conditional items without restricting the analysis to this subpopulation. Such an approach is of interest when the distributions of both types of responses are thought to be determined by common parameters affecting the whole population. By integrating the conditional item structure into the model, inference can be based both on unconditional data from the entire population and on conditional data from subjects for whom they exist. This approach opens new possibilities for multivariate analysis of such data. We apply this approach to latent class modelling and provide an example using data on respiratory symptoms (wheeze and cough) in children. Conditional data structures such as that considered here are common in medical research settings and, although our focus is on latent class models, the approach can be applied to other multivariate models.

Do Taxes Matter in the CAPM?

The traditional literature on the CAPM assumes that investor's tax payments simply vanish from the model. This assumption is not at all consistent with the actual behavior of the Treasury. The theory of general equilibrium states that an interest rate rf = 0 will not affect prices if taxes are introduced. We show that this result can be extended to the CAPM if the tax payments are redistributed among investors.

In vitro differentiation of near-unlimited numbers of functional mouse basophils using conditional Hoxb8

Background: Basophils constitute a rare leukocyte population known for their effector functions in inflammation and allergy, as well as more recently described immunoregulatory roles. Besides their low frequency, functional analysis of basophils is hindered by a short life span, inefficient ex vivo differentiation protocols, and lack of suitable cell models. A method to produce large quantities of basophils in vitro would facilitate basophil research and constitute a sought-after tool for diagnostic and drug testing purposes. Methods: A method is described to massively expand bone marrow–derived basophils in vitro. Myeloid progenitors are conditionally immortalized using Hoxb8 in the presence of interleukin-3 (IL-3) and outgrowing cell lines selected for their potential to differentiate into basophils upon shutdown of Hoxb8 expression. Results: IL-3-dependent, conditional Hoxb8-immortalized progenitor cell lines can be expanded and maintained in culture for prolonged periods. Upon shutdown of Hoxb8 expression, near-unlimited numbers of mature functional basophils can be differentiated in vitro within six days. The cells are end-differentiated and short-lived and express basophil-specific surface markers and proteases. Upon IgE- as well as C5a-mediated activation, differentiated basophils release granule enzymes and histamine and secrete Th2-type cytokines (IL-4, IL-13) and leukotriene C4. IL-3-deprivation induces apoptosis correlating with upregulation of the BH3-only proteins BCL-2-interacting mediator of cell death (BIM) and p53 upregulated modulator of apoptosis (PUMA) and downregulation of proviral integration site for Moloney murine leukemia virus 1 kinase (PIM-1). Conclusion: A novel method is presented to generate quantitative amounts of mouse basophils in vitro, which moreover allows genetic manipulation of conditionally immortalized progenitors. This approach may represent a useful alternative method to isolating primary basophils.