Individual genetic variation affects gene expression and cell phenotype by acting within complex molecular circuits, but this relationship is still largely unknown. Here, we combine genomic and meso-scale profiling with novel computational methods to detect genetic variants that affect the responsiveness of gene expression to stimulus (responsiveness QTLs) and position them in circuit diagrams. We apply this approach to study individual variation in transcriptional responsiveness to three different pathogen components in the model response of primary bone marrow dendritic cells (DCs) from recombinant inbred mice strains. We show that reQTLs are common both in cis (affecting a single target gene) and in trans (pleiotropically affecting co-regulated gene modules) and are specific to some stimuli but not others. Leveraging the stimulus-specific activity of reQTLs and the differential responsiveness of their associated targets, we show how to position reQTLs within the context of known pathways in this regulatory circuit. For example, we find that a pleiotropic trans-acting genetic factor in chr1:129-165Mb affects the responsiveness of 35 anti-viral genes only during an anti-viral like stimulus. Using RNAi we uncover RGS16 the likely causal gene in this interval, and an activator of the antiviral response. Our approach charts an experimental and analytic path to decipher the mechanisms underlying genetic variation in other complex circuits in primary mammalian cells.
Deciphering molecular circuits from genetic variation underlying transcriptional responsiveness to stimuli.
Age, Specimen part
View SamplesHundreds of immune cell types work in coordination to maintain tissue homeostasis. Upon infection, dramatic changes occur with the localization, migration and proliferation of the immune cells to first alert the body of the danger, confine it to limit spreading, and finally extinguish the threat and bring the tissue back to homeostasis. Since current technologies can follow the dynamics of only a limited number of cell types, we have yet to grasp the full complexity of global in vivo cell dynamics in normal developmental processes and disease. Here we devise a computational method, digital cell quantification (DCQ), which combines genomewide gene expression data with an immune cell compendium to infer in vivo dynamical changes in the quantities of 213 immune cell subpopulations. DCQ was applied to study global immune cell dynamics in mice lungs at ten time points during a 7-day time course of flu infection. We find dramatic changes in quantities of 70 immune cell types, including various innate, adaptive and progenitor immune cells. We focus on the previously unreported dynamics of four immune dendritic cell subtypes, and suggest a specific role for CD103+CD11b- cDCs in early stages of disease and CD8+ pDC in late stages of flu infection. Overall design: To study pathogenesis of Influenza infection, C57BL/6 mice (5 weeks) were infected intranasally with 4x103 PFU of influenza PR8 virus. We measured using RNA-Seq global gene expression in lung tissue at ten time points during a 7-day time course of infection, two infected individuals in each time point and four un-infected individuals as control. The lung organ was removed and transferred immediately into RNA Latter solution (Invitrogen).
Digital cell quantification identifies global immune cell dynamics during influenza infection.
Age, Specimen part, Cell line, Subject, Time
View SamplesHundreds of immune cell types work in coordination to maintain tissue homeostasis. Upon infection, dramatic changes occur with the localization, migration and proliferation of the immune cells to first alert the body of the danger, confine it to limit spreading, and finally extinguish the threat and bring the tissue back to homeostasis. Since current technologies can follow the dynamics of only a limited number of cell types, we have yet to grasp the full complexity of global in vivo cell dynamics in normal developmental processes and disease. Here we devise a computational method, digital cell quantification (DCQ), which combines genomewide gene expression data with an immune cell compendium to infer in vivo dynamical changes in the quantities of 213 immune cell subpopulations. DCQ was applied to study global immune cell dynamics in mice lungs at ten time points during a 7-day time course of flu infection. We find dramatic changes in quantities of 70 immune cell types, including various innate, adaptive and progenitor immune cells. We focus on the previously unreported dynamics of four immune dendritic cell subtypes, and suggest a specific role for CD103+CD11b- cDCs in early stages of disease and CD8+ pDC in late stages of flu infection. Overall design: To better understand the physiological role of these differential dynamic changes in the DCs, we measured the genome-wide RNA expression of all four DC subpopulations from lung of influenza infected mice at four time points following infections (two mice per time-point). For sorting dendritic cells from lungs, the lungs from infected and control uninfected C57BL/6J mice were immersed in cold PBS, cut into small pieces in 5 ml DMEM media containing 10% Bovine Fetal Serum, the cell suspensions were grinded using 1ml syringe cup on a 70 µm cell strainers (BD Falcon). The cells were washed with ice cold PBS. Remaining red blood cells were lysed using ammonium chloride solution (Sigma). Cells were harvested, immersed 1ml FACS buffer [PBS+2% FBS, 1mM EDTA], Fc receptors were blocked with anti-mouse CD16/CD32, washed with FACS buffer and divided into two tubes for sorting cDC and pDC cells.
Digital cell quantification identifies global immune cell dynamics during influenza infection.
Age, Specimen part, Cell line, Subject, Time
View SamplesThe host immune response against an infection requires the coordinated action of many diverse cell subsets that dynamically adapt to the pathogen threat. Here we combined WGCNA and DCQ to analyse time-resolved mouse splenic transcriptomes in acute and chronic LCMV infections. This approach allowed to better characterize the dynamic cell events occurring in complex tissues such as the induction of the adaptive T cell response which requires the coordination of monocytes/macrophages and CD8+ T cells. Overall design: mRNA profiles of CD8 T cells and macrophages (in duplicate days 0 and 7 post-infection) from C57BL/6 mice infected with 2x10E2 pfu of LCMV strain Docile, generated by deep sequencing.
Linking Cell Dynamics With Gene Coexpression Networks to Characterize Key Events in Chronic Virus Infections.
No sample metadata fields
View SamplesSingle-cell expression profiling is a rich resource of cellular heterogeneity. While profiling every sample under study is advantageous, such workflow is time consuming and costly. We devised CPM - a deconvolution algorithm in which cellular heterogeneity is inferred from bulk expression data based on pre-existing collection of single-cell RNA-seq profiles. We applied CPM to investigate individual variation in heterogeneity of murine lung cells during in vivo influenza virus infection, revealing that the relations between cell quantities and clinical outcomes varies in a gradual manner along the cellular activation process. Validation experiments confirmed these gradual changes along the cellular activation trajectory. Additional analysis suggests that clinical outcomes relate to the rate of cell activation at the early stages of this process. These findings demonstrate the utility of CPM as a mapping deconvolution tool at single-cell resolution, and highlight the importance of such fine cell landscape for understanding diversity of clinical outcomes. Overall design: Lungs gene expression of Collaborative Cross mice taken 48h after the infection with either the influenza virus or PBS.
Cell composition analysis of bulk genomics using single-cell data.
Specimen part, Subject, Time
View SamplesA central challenge in pharmaceutical research is to investigate genetic variation in response to drugs. The Collaborative Cross (CC) mouse reference population is a promising model for pharmacogenomic studies because of its large amount of genetic variation, genetic reproducibility, and dense recombination sites. While the CC lines are phenotypically diverse, their genetic diversity in drug disposition processes, such as detoxification reactions, is still largely uncharacterized. Here we systematically measured RNA-sequencing expression profiles from livers of 29 CC lines under baseline conditions. We then leveraged a reference collection of metabolic biotransformation pathways to map potential relations between drugs and their underlying expression quantitative trait loci (eQTLs). By applying this approach on proximal eQTLs, including eQTLs acting on the overall expression of genes and on the expression of particular transcript isoforms, we were able to construct the organization of hepatic eQTL-drug connectivity across the CC population. The analysis revealed a substantial impact of genetic variation acting on drug biotransformation, allowed mapping of potential joint genetic effects in the context of individual drugs, and demonstrated crosstalk between drug metabolism and lipid metabolism. Our findings provide a resource for investigating drug disposition in the CC strains, and offer a new paradigm for integrating biotransformation reactions to corresponding variations in DNA sequences. Overall design: This dataset includes RNA-Seq data of mRNA that were extracted from the liver of 55 male mice. The 55 mice belong to 29 different collaborative cross strains. The number of individual mice per strains is 3 for 3 strains, 2 for 16 strains, and 1 for 8 strains. All the mice are naïve without any special treatment.
Dissecting the Effect of Genetic Variation on the Hepatic Expression of Drug Disposition Genes across the Collaborative Cross Mouse Strains.
Specimen part, Cell line, Subject
View SamplesIschemic tolerance can be induced by numerous preconditioning stimuli, including various Toll-like receptor (TLR) ligands. We have shown previously that systemic administration of the TLR4 ligand, lipopolysaccharide (LPS) or the TLR9 ligand, unmethylated CpG ODNs prior to transient brain ischemia in mice confers substantial protection against ischemic damage. To elucidate the molecular mechanisms of preconditioning, we compared brain and blood genomic profiles in response to preconditioning with these TLR ligands and to preconditioning via exposure to brief ischemia.
Multiple preconditioning paradigms converge on interferon regulatory factor-dependent signaling to promote tolerance to ischemic brain injury.
Specimen part, Treatment
View SamplesWe used microarrays to detail the global program of gene expression during early hESC differentiation to mesendoderm using FBS, with and without RUNX1 depletion.
Transient RUNX1 Expression during Early Mesendodermal Differentiation of hESCs Promotes Epithelial to Mesenchymal Transition through TGFB2 Signaling.
Specimen part, Cell line
View SamplesThe signaling molecule retinoic acid (RA) regulates rod and cone photoreceptor fate, differentiation, and survival. The purpose of this study was to identify eye-specific genes controlled by RA during photoreceptor differentiation in the zebrafish.
Retinoic Acid Signaling Regulates Differential Expression of the Tandemly-Duplicated Long Wavelength-Sensitive Cone Opsin Genes in Zebrafish.
Specimen part
View SamplesThis SuperSeries is composed of the SubSeries listed below.
Anti-diabetic rosiglitazone remodels the adipocyte transcriptome by redistributing transcription to PPARγ-driven enhancers.
Cell line, Treatment, Time
View Samples