Nutrient-regulated control of lysosome function by signaling lipid conversion.

Lysosomes serve dual antagonistic functions in cells by mediating anabolic growth signaling and the catabolic turnover of macromolecules. How these janus-faced activities are regulated in response to cellular nutrient status is poorly understood. We show here that lysosome morphology and function are reversibly controlled by a nutrient-regulated signaling lipid switch that triggers the conversion between peripheral motile mTOR complex 1 (mTORC1) signaling-active and static mTORC1-inactive degradative lysosomes clustered at the cell center. Starvation-triggered relocalization of phosphatidylinositol 4-phosphate (PI(4)P)-metabolizing enzymes reshapes the lysosomal surface proteome to facilitate lysosomal proteolysis and to repress mTORC1 signaling. Concomitantly, lysosomal phosphatidylinositol 3-phosphate (PI(3)P), which marks motile signaling-active lysosomes in the cell periphery, is erased. Interference with this PI(3)P/PI(4)P lipid switch module impairs the adaptive response of cells to altering nutrient supply. Our data unravel a key function for lysosomal phosphoinositide metabolism in rewiring organellar membrane dynamics in response to cellular nutrient status.

mTOR and S6K1 drive polycystic kidney by the control of Afadin-dependent oriented cell division.

mTOR activation is essential and sufficient to cause polycystic kidneys in Tuberous Sclerosis Complex (TSC) and other genetic disorders. In disease models, a sharp increase of proliferation and cyst formation correlates with a dramatic loss of oriented cell division (OCD). We find that OCD distortion is intrinsically due to S6 kinase 1 (S6K1) activation. The concomitant loss of S6K1 in Tsc1-mutant mice restores OCD but does not decrease hyperproliferation, leading to non-cystic harmonious hyper growth of kidneys. Mass spectrometry-based phosphoproteomics for S6K1 substrates revealed Afadin, a known component of cell-cell junctions required to couple intercellular adhesions and cortical cues to spindle orientation. Afadin is directly phosphorylated by S6K1 and abnormally decorates the apical surface of Tsc1-mutant cells with E-cadherin and α-catenin. Our data reveal that S6K1 hyperactivity alters centrosome positioning in mitotic cells, affecting oriented cell division and promoting kidney cysts in conditions of mTOR hyperactivity.

ZRF1 is a novel S6 kinase substrate that drives the senescence programme.

The inactivation of S6 kinases mimics several aspects of caloric restriction, including small body size, increased insulin sensitivity and longevity. However, the impact of S6 kinase activity on cellular senescence remains to be established. Here, we show that the constitutive activation of mammalian target of rapamycin complex 1 (mTORC1) by tuberous sclerosis complex (TSC) mutations induces a premature senescence programme in fibroblasts that relies on S6 kinases. To determine novel molecular targets linking S6 kinase activation to the control of senescence, we set up a chemical genetic screen, leading to the identification of the nuclear epigenetic factor ZRF1 (also known as DNAJC2, MIDA1, Mpp11). S6 kinases phosphorylate ZRF1 on Ser47 in cultured cells and in mammalian tissues in vivo Knock-down of ZRF1 or expression of a phosphorylation mutant is sufficient to blunt the S6 kinase-dependent senescence programme. This is traced by a sharp alteration in p16 levels, the cell cycle inhibitor and a master regulator of senescence. Our findings reveal a mechanism by which nutrient sensing pathways impact on cell senescence through the activation of mTORC1-S6 kinases and the phosphorylation of ZRF1.

Regulation of YAP by mTOR and autophagy reveals a therapeutic target of tuberous sclerosis complex.

Genetic studies have shown that the tuberous sclerosis complex (TSC) 1-TSC2-mammalian target of Rapamycin (mTOR) and the Hippo-Yes-associated protein 1 (YAP) pathways are master regulators of organ size, which are often involved in tumorigenesis. The crosstalk between these signal transduction pathways in coordinating environmental cues, such as nutritional status and mechanical constraints, is crucial for tissue growth. Whether and how mTOR regulates YAP remains elusive. Here we describe a novel mouse model of TSC which develops renal mesenchymal lesions recapitulating human perivascular epithelioid cell tumors (PEComas) from patients with TSC. We identify that YAP is up-regulated by mTOR in mouse and human PEComas. YAP inhibition blunts abnormal proliferation and induces apoptosis of TSC1-TSC2-deficient cells, both in culture and in mosaic Tsc1 mutant mice. We further delineate that YAP accumulation in TSC1/TSC2-deficient cells is due to impaired degradation of the protein by the autophagosome/lysosome system. Thus, the regulation of YAP by mTOR and autophagy is a novel mechanism of growth control, matching YAP activity with nutrient availability under growth-permissive conditions. YAP may serve as a potential therapeutic target for TSC and other diseases with dysregulated mTOR activity.

Interaction between genetic and epigenetic variation defines gene expression patterns at the asthma-associated locus 17q12-q21 in lymphoblastoid cell lines.

Phenotypic variation results from variation in gene expression, which is modulated by genetic and/or epigenetic factors. To understand the molecular basis of human disease, interaction between genetic and epigenetic factors needs to be taken into account. The asthma-associated region 17q12-q21 harbors three genes, the zona pellucida binding protein 2 (ZPBP2), gasdermin B (GSDMB) and ORM1-like 3 (ORMDL3), that show allele-specific differences in expression levels in lymphoblastoid cell lines (LCLs) and CD4+ T cells. Here, we report a molecular dissection of allele-specific transcriptional regulation of the genes within the chromosomal region 17q12-q21 combining in vitro transfection, formaldehyde-assisted isolation of regulatory elements, chromatin immunoprecipitation and DNA methylation assays in LCLs. We found that a single nucleotide polymorphism rs4795397 influences the activity of ZPBP2 promoter in vitro in an allele-dependent fashion, and also leads to nucleosome repositioning on the asthma-associated allele. However, variable methylation of exon 1 of ZPBP2 masks the strong genetic effect on ZPBP2 promoter activity in LCLs. In contrast, the ORMDL3 promoter is fully unmethylated, which allows detection of genetic effects on its transcription. We conclude that the cis-regulatory effects on 17q12-q21 gene expression result from interaction between several regulatory polymorphisms and epigenetic factors within the cis-regulatory haplotype region.

S6 kinase 1 is required for rapamycin-sensitive liver proliferation after mouse hepatectomy.

Rapamycin is an antibiotic inhibiting eukaryotic cell growth and proliferation by acting on target of rapamycin (TOR) kinase. Mammalian TOR (mTOR) is thought to work through 2 independent complexes to regulate cell size and cell replication, and these 2 complexes show differential sensitivity to rapamycin. Here we combine functional genetics and pharmacological treatments to analyze rapamycin-sensitive mTOR substrates that are involved in cell proliferation and tissue regeneration after partial hepatectomy in mice. After hepatectomy, hepatocytes proliferated rapidly, correlating with increased S6 kinase phosphorylation, while treatment with rapamycin derivatives impaired regeneration and blocked S6 kinase activation. In addition, genetic deletion of S6 kinase 1 (S6K1) caused a delay in S phase entry in hepatocytes after hepatectomy. The proliferative defect of S6K1-deficient hepatocytes was cell autonomous, as it was also observed in primary cultures and hepatic overexpression of S6K1-rescued proliferation. We found that S6K1 controlled steady-state levels of cyclin D1 (Ccnd1) mRNA in liver, and cyclin D1 expression was required to promote hepatocyte cell cycle. Notably, in vivo overexpression of cyclin D1 was sufficient to restore the proliferative capacity of S6K-null livers. The identification of an S6K1-dependent mechanism participating in cell proliferation in vivo may be relevant for cancer cells displaying high mTOR complex 1 activity and cyclin D1 accumulation.

Genome-wide assessment of imprinted expression in human cells.

BACKGROUND: Parent-of-origin-dependent expression of alleles, imprinting, has been suggested to impact a substantial proportion of mammalian genes. Its discovery requires allele-specific detection of expressed transcripts, but in some cases detected allelic expression bias has been interpreted as imprinting without demonstrating compatible transmission patterns and excluding heritable variation. Therefore, we utilized a genome-wide tool exploiting high density genotyping arrays in parallel measurements of genotypes in RNA and DNA to determine allelic expression across the transcriptome in lymphoblastoid cell lines (LCLs) and skin fibroblasts derived from families. RESULTS: We were able to validate 43% of imprinted genes with previous demonstration of compatible transmission patterns in LCLs and fibroblasts. In contrast, we only validated 8% of genes suggested to be imprinted in the literature, but without clear evidence of parent-of-origin-determined expression. We also detected five novel imprinted genes and delineated regions of imprinted expression surrounding annotated imprinted genes. More subtle parent-of-origin-dependent expression, or partial imprinting, could be verified in four genes. Despite higher prevalence of monoallelic expression, immortalized LCLs showed consistent imprinting in fewer loci than primary cells. Random monoallelic expression has previously been observed in LCLs and we show that random monoallelic expression in LCLs can be partly explained by aberrant methylation in the genome. CONCLUSIONS: Our results indicate that widespread parent-of-origin-dependent expression observed recently in rodents is unlikely to be captured by assessment of human cells derived from adult tissues where genome-wide assessment of both primary and immortalized cells yields few new imprinted loci.

Global analysis of the impact of environmental perturbation on cis-regulation of gene expression.

Genetic variants altering cis-regulation of normal gene expression (cis-eQTLs) have been extensively mapped in human cells and tissues, but the extent by which controlled, environmental perturbation influences cis-eQTLs is unclear. We carried out large-scale induction experiments using primary human bone cells derived from unrelated donors of Swedish origin treated with 18 different stimuli (7 treatments and 2 controls, each assessed at 2 time points). The treatments with the largest impact on the transcriptome, verified on two independent expression arrays, included BMP-2 (t = 2h), dexamethasone (DEX) (t = 24 h), and PGE2 (t = 24 h). Using these treatments and control, we performed expression profiling for 18,144 RefSeq transcripts on biological replicates of the complete study cohort of 113 individuals (n(total) = 782) and combined it with genome-wide SNP-genotyping data in order to map treatment-specific cis-eQTLs (defined as SNPs located within the gene ± 250 kb). We found that 93% of cis-eQTLs at 1% FDR were observed in at least one additional treatment, and in fact, on average, only 1.4% of the cis-eQTLs were considered as treatment-specific at high confidence. The relative invariability of cis-regulation following perturbation was reiterated independently by genome-wide allelic expression tests where only a small proportion of variance could be attributed to treatment. Treatment-specific cis-regulatory effects were, however, 2- to 6-fold more abundant among differently expressed genes upon treatment. We further followed-up and validated the DEX-specific cis-regulation of the MYO6 and TNC loci and found top cis-regulatory variants located 180 kb and 250 kb upstream of the transcription start sites, respectively. Our results suggest that, as opposed to tissue-specificity of cis-eQTLs, the interactions between cellular environment and cis-variants are relatively rare (∼1.5%), but that detection of such specific interactions can be achieved by a combination of functional genomic approaches as described here.

Global patterns of cis variation in human cells revealed by high-density allelic expression analysis.

Cis-acting variants altering gene expression are a source of phenotypic differences. The cis-acting components of expression variation can be identified through the mapping of differences in allelic expression (AE), which is the measure of relative expression between two allelic transcripts. We generated a map of AE associated SNPs using quantitative measurements of AE on Illumina Human1M BeadChips. In 53 lymphoblastoid cell lines derived from donors of European descent, we identified common cis variants affecting 30% (2935/9751) of the measured RefSeq transcripts at 0.001 permutation significance. The pervasive influence of cis-regulatory variants, which explain 50% of population variation in AE, extend to full-length transcripts and their isoforms as well as to unannotated transcripts. These strong effects facilitate fine mapping of cis-regulatory SNPs, as demonstrated by dissection of heritable control of transcripts in the systemic lupus erythematosus-associated C8orf13-BLK region in chromosome 8. The dense collection of associations will facilitate large-scale isolation of cis-regulatory SNPs.

Allele-specific chromatin remodeling in the ZPBP2/GSDMB/ORMDL3 locus associated with the risk of asthma and autoimmune disease.

Common SNPs in the chromosome 17q12-q21 region alter the risk for asthma, type 1 diabetes, primary biliary cirrhosis, and Crohn disease. Previous reports by us and others have linked the disease-associated genetic variants with changes in expression of GSDMB and ORMDL3 transcripts in human lymphoblastoid cell lines (LCLs). The variants also alter regulation of other transcripts, and this domain-wide cis-regulatory effect suggests a mechanism involving long-range chromatin interactions. Here, we further dissect the disease-linked haplotype and identify putative causal DNA variants via a combination of genetic and functional analyses. First, high-throughput resequencing of the region and genotyping of potential candidate variants were performed. Next, additional mapping of allelic expression differences in Yoruba HapMap LCLs allowed us to fine-map the basis of the cis-regulatory differences to a handful of candidate functional variants. Functional assays identified allele-specific differences in nucleosome distribution, an allele-specific association with the insulator protein CTCF, as well as a weak promoter activity for rs12936231. Overall, this study shows a common disease allele linked to changes in CTCF binding and nucleosome occupancy leading to altered domain-wide cis-regulation. Finally, a strong association between asthma and cis-regulatory haplotypes was observed in three independent family-based cohorts (p = 1.78 x 10(-8)). This study demonstrates the requirement of multiple parallel allele-specific tools for the investigation of noncoding disease variants and functional fine-mapping of human disease-associated haplotypes.

Tissue effect on genetic control of transcript isoform variation.

Current genome-wide association studies (GWAS) are moving towards the use of large cohorts of primary cell lines to study a disease of interest and to assign biological relevance to the genetic signals identified. Here, we use a panel of human osteoblasts (HObs) to carry out a transcriptomic survey, similar to recent studies in lymphoblastoid cell lines (LCLs). The distinct nature of HObs and LCLs is reflected by the preferential grouping of cell type-specific genes within biologically and functionally relevant pathways unique to each tissue type. We performed cis-association analysis with SNP genotypes to identify genetic variations of transcript isoforms, and our analysis indicates that differential expression of transcript isoforms in HObs is also partly controlled by cis-regulatory genetic variants. These isoforms are regulated by genetic variants in both a tissue-specific and tissue-independent fashion, and these associations have been confirmed by RT-PCR validation. Our study suggests that multiple transcript isoforms are often present in both tissues and that genetic control may affect the relative expression of one isoform to another, rather than having an all-or-none effect. Examination of the top SNPs from a GWAS of bone mineral density show overlap with probeset associations observed in this study. The top hit corresponding to the FAM118A gene was tested for association studies in two additional clinical studies, revealing a novel transcript isoform variant. Our approach to examining transcriptome variation in multiple tissue types is useful for detecting the proportion of genetic variation common to different cell types and for the identification of cell-specific isoform variants that may be functionally relevant, an important follow-up step for GWAS.

Population genomics in a disease targeted primary cell model.

The common genetic variants associated with complex traits typically lie in noncoding DNA and may alter gene regulation in a cell type-specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an expression quantitative trait loci (eQTL) study of primary human osteoblasts (HOb) derived from 95 unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top 2000 BMD-associated SNPs (P < approximately 10(-3)) were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment (threefold) of converging cis-eQTLs as compared to LCLs. The top 10 BMD loci with SNPs showing strong cis-effects on gene expression in HObs (P = 6 x 10(-10) - 7 x 10(-16)) were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All 10 variants were tested in the Swedish MrOS Cohort (n = 3014), providing evidence for two novel BMD loci (SRR and MSH3). These variants were then tested in the Rotterdam Study (n = 2090), yielding converging evidence for BMD association at the 17p13.3 SRR locus (P(combined) = 5.6 x 10(-5)). The cis-regulatory effect was further fine-mapped to the proximal promoter of the SRR gene (rs3744270, r(2) = 0.5, P = 2.6 x 10(-15)). Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.

Targeted screening of cis-regulatory variation in human haplotypes.

Regulatory cis-acting variants account for a large proportion of gene expression variability in populations. Cis-acting differences can be specifically measured by comparing relative levels of allelic transcripts within a sample. Allelic expression (AE) mapping for cis-regulatory variant discovery has been hindered by the requirements of having informative or heterozygous single nucleotide polymorphisms (SNPs) within genes in order to assign the allelic origin of each transcript. In this study we have developed an approach to systematically screen for heritable cis-variants in common human haplotypes across >1,000 genes. In order to achieve the highest level of information per haplotype studied, we carried out allelic expression measurements by using both intronic and exonic SNPs in primary transcripts. We used a novel RNA pooling strategy in immortalized lymphoblastoid cell lines (LCLs) and primary human osteoblast cell lines (HObs) to allow for high-throughput AE. Screening hits from RNA pools were further validated by performing allelic expression mapping in individual samples. Our results indicate that >10% of expressed genes in human LCLs show genotype-linked AE. In addition, we have validated cis-acting variants in over 20 genes linked with common disease susceptibility in recent genome-wide studies. More generally, our results indicate that RNA pooling coupled with AE read-out by second generation sequencing or by other methods provides a high-throughput tool for cataloging the impact of common noncoding variants in the human genome.

Tunable sustained release properties of "onion-like" phospholipids multilamellar vesicles.

"Onion-type" multilamellar micro-vesicles of phospholipids (spherulites) were doped with different amounts of a cationic cosurfactant ((-)N-dodecyl-N-methylephedrinium bromide) for the purpose of controlling the sustained release of anionic drugs. Three weak acid probes (methyl red, chlorophenol red, and ibuprofen) were encapsulated in the vesicles as drug models. The kinetics and rate of release were studied by absorption spectroscopy and HPLC. The effect of probe charge (pH above and below pKa of the probes), of cosurfactant concentration and of added salt was investigated. It was found that, above pKa (i.e., when the probes are anionic), the release can be almost totally inhibited by doping the vesicles with 2.4 wt% of cationic cosurfactant. The release properties can even be finely tuned by controlling the amounts of the cosurfactant. Salt and pH effects demonstrate the role of electrostatic interactions in sustaining the release.