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

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

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 = 24h), and PGE2 (t = 24h). 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 (ntotal = 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

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

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 >1000 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

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

mTOR activation is known to generate polycystic kidneys, which show both increased proliferation and loss of oriented cell division (OCD). Here, Bonucci et al. show that loss of OCD is linked to S6K1 activation through its direct target Afadin and is separable from hyperproliferation

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

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 < similar to 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