Estimating RNA dynamics using one time point for one sample in a single-pulse metabolic labeling experiment.

BACKGROUND: Over the past decade, experimental procedures such as metabolic labeling for determining RNA turnover rates at the transcriptome-wide scale have been widely adopted and are now turning to single cell measurements. Several computational methods to estimate RNA synthesis, processing and degradation rates from such experiments have been suggested, but they all require several RNA sequencing samples. Here we present a method that can estimate those three rates from a single sample. METHODS: Our method relies on the analytical solution to the Zeisel model of RNA dynamics. It was validated on metabolic labeling experiments performed on mouse embryonic stem cells. Resulting degradation rates were compared both to previously published rates on the same system and to a state-of-the-art method applied to the same data. RESULTS: Our method is computationally efficient and outputs rates that correlate well with previously published data sets. Using it on a single sample, we were able to reproduce the observation that dynamic biological processes tend to involve genes with higher metabolic rates, while stable processes involve genes with lower rates. This supports the hypothesis that cells control not only the mRNA steady-state abundance, but also its responsiveness, i.e., how fast steady state is reached. Moreover, degradation rates obtained with our method compare favourably with the other tested method. CONCLUSIONS: In addition to saving experimental work and computational time, estimating rates for a single sample has several advantages. It does not require an error-prone normalization across samples and enables the use of replicates to estimate uncertainty and assess sample quality. Finally the method and theoretical results described here are general enough to be useful in other contexts such as nucleotide conversion methods and single cell metabolic labeling experiments.

Changes in resource partitioning between and within organs support growth adjustment to neighbor proximity in Brassicaceae seedlings.

In shade-intolerant plants, the perception of proximate neighbors rapidly induces architectural changes resulting in elongated stems and reduced leaf size. Sensing and signaling steps triggering this modified growth program have been identified. However, the underlying changes in resource allocation that fuel stem growth remain poorly understood. Through 14CO2 pulse labeling of Brassica rapa seedlings, we show that perception of the neighbor detection signal, low ratio of red to far-red light (R:FR), leads to increased carbon allocation from the major site of photosynthesis (cotyledons) to the elongating hypocotyl. While carbon fixation and metabolite levels remain similar in low R:FR, partitioning to all downstream carbon pools within the hypocotyl is increased. Genetic analyses using Arabidopsis thaliana mutants indicate that low-R:FR-induced hypocotyl elongation requires sucrose transport from the cotyledons and is regulated by a PIF7-dependent metabolic response. Moreover, our data suggest that starch metabolism in the hypocotyl has a growth-regulatory function. The results reveal a key mechanism by which metabolic adjustments can support rapid growth adaptation to a changing environment.

Seventy-five genetic loci influencing the human red blood cell.

Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth retardation and impaired physical capacity. To understand further the genetic factors influencing red blood cells, we carried out a genome-wide association study of haemoglobin concentration and related parameters in up to 135,367 individuals. Here we identify 75 independent genetic loci associated with one or more red blood cell phenotypes at P < 10(-8), which together explain 4-9% of the phenotypic variance per trait. Using expression quantitative trait loci and bioinformatic strategies, we identify 121 candidate genes enriched in functions relevant to red blood cell biology. The candidate genes are expressed preferentially in red blood cell precursors, and 43 have haematopoietic phenotypes in Mus musculus or Drosophila melanogaster. Through open-chromatin and coding-variant analyses we identify potential causal genetic variants at 41 loci. Our findings provide extensive new insights into genetic mechanisms and biological pathways controlling red blood cell formation and function.

A higher mutational burden in females supports a "female protective model" in neurodevelopmental disorders.

Increased male prevalence has been repeatedly reported in several neurodevelopmental disorders (NDs), leading to the concept of a "female protective model." We investigated the molecular basis of this sex-based difference in liability and demonstrated an excess of deleterious autosomal copy-number variants (CNVs) in females compared to males (odds ratio [OR] = 1.46, p = 8 × 10(-10)) in a cohort of 15,585 probands ascertained for NDs. In an independent autism spectrum disorder (ASD) cohort of 762 families, we found a 3-fold increase in deleterious autosomal CNVs (p = 7 × 10(-4)) and an excess of private deleterious single-nucleotide variants (SNVs) in female compared to male probands (OR = 1.34, p = 0.03). We also showed that the deleteriousness of autosomal SNVs was significantly higher in female probands (p = 0.0006). A similar bias was observed in parents of probands ascertained for NDs. Deleterious CNVs (>400 kb) were maternally inherited more often (up to 64%, p = 10(-15)) than small CNVs < 400 kb (OR = 1.45, p = 0.0003). In the ASD cohort, increased maternal transmission was also observed for deleterious CNVs and SNVs. Although ASD females showed higher mutational burden and lower cognition, the excess mutational burden remained, even after adjustment for those cognitive differences. These results strongly suggest that females have an increased etiological burden unlinked to rare deleterious variants on the X chromosome. Carefully phenotyped and genotyped cohorts will be required for identifying the symptoms, which show gender-specific liability to mutational burden.

TLR3-Mediated CD8+ Dendritic Cell Activation Is Coupled with Establishment of a Cell-Intrinsic Antiviral State.

Because of their unique capacity to cross-present Ags to CD8(+) T cells, mouse lymphoid tissue-resident CD8(+) dendritic cells (DCs) and their migratory counterparts are critical for priming antiviral T cell responses. High expression of the dsRNA sensor TLR3 is a distinctive feature of these cross-presenting DC subsets. TLR3 engagement in CD8(+) DCs promotes cross-presentation and the acquisition of effector functions required for driving antiviral T cell responses. In this study, we performed a comprehensive analysis of the TLR3-induced antiviral program and cell-autonomous immunity in CD8(+) DC lines and primary CD8(+) DCs. We found that TLR3-ligand polyinosinic-polycytidylic acid and human rhinovirus infection induced a potent antiviral protection against Sendai and vesicular stomatitis virus in a TLR3 and type I IFN receptor-dependent manner. Polyinosinic-polycytidylic acid-induced antiviral genes were identified by mass spectrometry-based proteomics and transcriptomics in the CD8(+) DC line. Nanostring nCounter experiments confirmed that these antiviral genes were induced by TLR3 engagement in primary CD8(+) DCs, and indicated that many are secondary TLR3-response genes requiring autocrine IFN-ß stimulation. TLR3-activation thus establishes a type I IFN-dependent antiviral program in a DC subtype playing crucial roles in priming adaptive antiviral immune responses. This mechanism is likely to shield the priming of antiviral responses against inhibition or abrogation by the viral infection. It could be particularly relevant for viruses detected mainly by TLR3, which may not trigger type I IFN production by DCs that lack TLR3, such as plasmacytoid DCs or CD8(-) DCs.

Light intensity modulates the regulatory network of the shade avoidance response in Arabidopsis.

Plants such as Arabidopsis thaliana respond to foliar shade and neighbors who may become competitors for light resources by elongation growth to secure access to unfiltered sunlight. Challenges faced during this shade avoidance response (SAR) are different under a light-absorbing canopy and during neighbor detection where light remains abundant. In both situations, elongation growth depends on auxin and transcription factors of the phytochrome interacting factor (PIF) class. Using a computational modeling approach to study the SAR regulatory network, we identify and experimentally validate a previously unidentified role for long hypocotyl in far red 1, a negative regulator of the PIFs. Moreover, we find that during neighbor detection, growth is promoted primarily by the production of auxin. In contrast, in true shade, the system operates with less auxin but with an increased sensitivity to the hormonal signal. Our data suggest that this latter signal is less robust, which may reflect a cost-to-robustness tradeoff, a system trait long recognized by engineers and forming the basis of information theory.

Phytochrome Kinase Substrate 4 is phosphorylated by the phototropin 1 photoreceptor.

Phototropism allows plants to redirect their growth towards the light to optimize photosynthesis under reduced light conditions. Phototropin 1 (phot1) is the primary low blue light-sensing receptor triggering phototropism in Arabidopsis. Light-induced autophosphorylation of phot1, an AGC-class protein kinase, constitutes an essential step for phototropism. However, apart from the receptor itself, substrates of phot1 kinase activity are less clearly established. Phototropism is also influenced by the cryptochromes and phytochromes photoreceptors that do not provide directional information but influence the process through incompletely characterized mechanisms. Here, we show that Phytochrome Kinase Substrate 4 (PKS4), a known element of phot1 signalling, is a substrate of phot1 kinase activity in vitro that is phosphorylated in a phot1-dependent manner in vivo. PKS4 phosphorylation is transient and regulated by a type 2-protein phosphatase. Moreover, phytochromes repress the accumulation of the light-induced phosphorylated form of PKS4 showing a convergence of photoreceptor activity on this signalling element. Our physiological analyses suggest that PKS4 phosphorylation is not essential for phototropism but is part of a negative feedback mechanism.

Pom1 gradient buffering through intermolecular auto-phosphorylation.

Concentration gradients provide spatial information for tissue patterning and cell organization, and their robustness under natural fluctuations is an evolutionary advantage. In rod-shaped Schizosaccharomyces pombe cells, the DYRK-family kinase Pom1 gradients control cell division timing and placement. Upon dephosphorylation by a Tea4-phosphatase complex, Pom1 associates with the plasma membrane at cell poles, where it diffuses and detaches upon auto-phosphorylation. Here, we demonstrate that Pom1 auto-phosphorylates intermolecularly, both in vitro and in vivo, which confers robustness to the gradient. Quantitative imaging reveals this robustness through two system's properties: The Pom1 gradient amplitude is inversely correlated with its decay length and is buffered against fluctuations in Tea4 levels. A theoretical model of Pom1 gradient formation through intermolecular auto-phosphorylation predicts both properties qualitatively and quantitatively. This provides a telling example where gradient robustness through super-linear decay, a principle hypothesized a decade ago, is achieved through autocatalysis. Concentration-dependent autocatalysis may be a widely used simple feedback to buffer biological activities.

Nuclear phytochrome A signaling promotes phototropism in Arabidopsis.

Phototropin photoreceptors (phot1 and phot2 in Arabidopsis thaliana) enable responses to directional light cues (e.g., positive phototropism in the hypocotyl). In Arabidopsis, phot1 is essential for phototropism in response to low light, a response that is also modulated by phytochrome A (phyA), representing a classical example of photoreceptor coaction. The molecular mechanisms underlying promotion of phototropism by phyA remain unclear. Most phyA responses require nuclear accumulation of the photoreceptor, but interestingly, it has been proposed that cytosolic phyA promotes phototropism. By comparing the kinetics of phototropism in seedlings with different subcellular localizations of phyA, we show that nuclear phyA accelerates the phototropic response, whereas in the fhy1 fhl mutant, in which phyA remains in the cytosol, phototropic bending is slower than in the wild type. Consistent with this data, we find that transcription factors needed for full phyA responses are needed for normal phototropism. Moreover, we show that phyA is the primary photoreceptor promoting the expression of phototropism regulators in low light (e.g., PHYTOCHROME KINASE SUBSTRATE1 [PKS1] and ROOT PHOTO TROPISM2 [RPT2]). Although phyA remains cytosolic in fhy1 fhl, induction of PKS1 and RPT2 expression still occurs in fhy1 fhl, indicating that a low level of nuclear phyA signaling is still present in fhy1 fhl.

Mechanical forces drive ordered patterning of hair cells in the mammalian inner ear.

Periodic organization of cells is required for the function of many organs and tissues. The development of such periodic patterns is typically associated with mechanisms based on intercellular signaling such as lateral inhibition and Turing patterning. Here we show that the transition from disordered to ordered checkerboard-like pattern of hair cells and supporting cells in the mammalian hearing organ, the organ of Corti, is likely based on mechanical forces rather than signaling events. Using time-lapse imaging of mouse cochlear explants, we show that hair cells rearrange gradually into a checkerboard-like pattern through a tissue-wide shear motion that coordinates intercalation and delamination events. Using mechanical models of the tissue, we show that global shear and local repulsion forces on hair cells are sufficient to drive the transition from disordered to ordered cellular pattern. Our findings suggest that mechanical forces drive ordered hair cell patterning in a process strikingly analogous to the process of shear-induced crystallization in polymer and granular physics.

High capacity in G protein-coupled receptor signaling.

G protein-coupled receptors (GPCRs) constitute a large family of receptors that activate intracellular signaling pathways upon detecting specific extracellular ligands. While many aspects of GPCR signaling have been uncovered through decades of studies, some fundamental properties, like its channel capacity-a measure of how much information a given transmission system can reliably transduce-are still debated. Previous studies concluded that GPCRs in individual cells could transmit around one bit of information about the concentration of the ligands, allowing only for a reliable on or off response. Using muscarinic receptor-induced calcium response measured in individual cells upon repeated stimulation, we show that GPCR signaling systems possess a significantly higher capacity. We estimate the channel capacity of this system to be above two, implying that at least four concentration levels of the agonist can be distinguished reliably. These findings shed light on the basic principles of GPCR signaling.

Cell-Cell Contact Area Affects Notch Signaling and Notch-Dependent Patterning.

During development, cells undergo dramatic changes in their morphology. By affecting contact geometry, these morphological changes could influence cellular communication. However, it has remained unclear whether and how signaling depends on contact geometry. This question is particularly relevant for Notch signaling, which coordinates neighboring cell fates through direct cell-cell signaling. Using micropatterning with a receptor trans-endocytosis assay, we show that signaling between pairs of cells correlates with their contact area. This relationship extends across contact diameters ranging from micrometers to tens of micrometers. Mathematical modeling predicts that dependence of signaling on contact area can bias cellular differentiation in Notch-mediated lateral inhibition processes, such that smaller cells are more likely to differentiate into signal-producing cells. Consistent with this prediction, analysis of developing chick inner ear revealed that ligand-producing hair cell precursors have smaller apical footprints than non-hair cells. Together, these results highlight the influence of cell morphology on fate determination processes.

Replication and Characterization of Association between ABO SNPs and Red Blood Cell Traits by Meta-Analysis in Europeans.

Red blood cell (RBC) traits are routinely measured in clinical practice as important markers of health. Deviations from the physiological ranges are usually a sign of disease, although variation between healthy individuals also occurs, at least partly due to genetic factors. Recent large scale genetic studies identified loci associated with one or more of these traits; further characterization of known loci and identification of new loci is necessary to better understand their role in health and disease and to identify potential molecular mechanisms. We performed meta-analysis of Metabochip association results for six RBC traits-hemoglobin concentration (Hb), hematocrit (Hct), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), mean corpuscular volume (MCV) and red blood cell count (RCC)-in 11 093 Europeans from seven studies of the UCL-LSHTM-Edinburgh-Bristol (UCLEB) Consortium. We identified 394 non-overlapping SNPs in five loci at genome-wide significance: 6p22.1-6p21.33 (with HFE among others), 6q23.2 (with HBS1L among others), 6q23.3 (contains no genes), 9q34.3 (only ABO gene) and 22q13.1 (with TMPRSS6 among others), replicating previous findings of association with RBC traits at these loci and extending them by imputation to 1000 Genomes. We further characterized associations between ABO SNPs and three traits: hemoglobin, hematocrit and red blood cell count, replicating them in an independent cohort. Conditional analyses indicated the independent association of each of these traits with ABO SNPs and a role for blood group O in mediating the association. The 15 most significant RBC-associated ABO SNPs were also associated with five cardiometabolic traits, with discordance in the direction of effect between groups of traits, suggesting that ABO may act through more than one mechanism to influence cardiometabolic risk.

Distinct levels in Pom1 gradients limit Cdr2 activity and localization to time and position division.

Where and when cells divide are fundamental questions. In rod-shaped fission yeast cells, the DYRK-family kinase Pom1 is organized in concentration gradients from cell poles and controls cell division timing and positioning. Pom1 gradients restrict to mid-cell the SAD-like kinase Cdr2, which recruits Mid1/Anillin for medial division. Pom1 also delays mitotic commitment through Cdr2, which inhibits Wee1. Here, we describe quantitatively the distributions of cortical Pom1 and Cdr2. These reveal low profile overlap contrasting with previous whole-cell measurements and Cdr2 levels increase with cell elongation, raising the possibility that Pom1 regulates mitotic commitment by controlling Cdr2 medial levels. However, we show that distinct thresholds of Pom1 activity define the timing and positioning of division. Three conditions-a separation-of-function Pom1 allele, partial downregulation of Pom1 activity, and haploinsufficiency in diploid cells-yield cells that divide early, similar to pom1 deletion, but medially, like wild-type cells. In these cells, Cdr2 is localized correctly at mid-cell. Further, Cdr2 overexpression promotes precocious mitosis only in absence of Pom1. Thus, Pom1 inhibits Cdr2 for mitotic commitment independently of regulating its localization or cortical levels. Indeed, we show Pom1 restricts Cdr2 activity through phosphorylation of a C-terminal self-inhibitory tail. In summary, our results demonstrate that distinct levels in Pom1 gradients delineate a medial Cdr2 domain, for cell division placement, and control its activity, for mitotic commitment.

Mapping genetic variants associated with beta-adrenergic responses in inbred mice.

ß-blockers and ß-agonists are primarily used to treat cardiovascular diseases. Inter-individual variability in response to both drug classes is well recognized, yet the identity and relative contribution of the genetic players involved are poorly understood. This work is the first genome-wide association study (GWAS) addressing the values and susceptibility of cardiovascular-related traits to a selective ß(1)-blocker, Atenolol (ate), and a ß-agonist, Isoproterenol (iso). The phenotypic dataset consisted of 27 highly heritable traits, each measured across 22 inbred mouse strains and four pharmacological conditions. The genotypic panel comprised 79922 informative SNPs of the mouse HapMap resource. Associations were mapped by Efficient Mixed Model Association (EMMA), a method that corrects for the population structure and genetic relatedness of the various strains. A total of 205 separate genome-wide scans were analyzed. The most significant hits include three candidate loci related to cardiac and body weight, three loci for electrocardiographic (ECG) values, two loci for the susceptibility of atrial weight index to iso, four loci for the susceptibility of systolic blood pressure (SBP) to perturbations of the ß-adrenergic system, and one locus for the responsiveness of QTc (p<10(-8)). An additional 60 loci were suggestive for one or the other of the 27 traits, while 46 others were suggestive for one or the other drug effects (p<10(-6)). Most hits tagged unexpected regions, yet at least two loci for the susceptibility of SBP to ß-adrenergic drugs pointed at members of the hypothalamic-pituitary-thyroid axis. Loci for cardiac-related traits were preferentially enriched in genes expressed in the heart, while 23% of the testable loci were replicated with datasets of the Mouse Phenome Database (MPD). Altogether these data and validation tests indicate that the mapped loci are relevant to the traits and responses studied.