Genetic variations at the human growth hormone receptor (GHR) gene locus are associated with idiopathic short stature.

GH plays an essential role in the growing child by binding to the growth hormone receptor (GHR) on target cells and regulating multiple growth promoting and metabolic effects. Mutations in the GHR gene coding regions result in GH insensitivity (dwarfism) due to a dysfunctional receptor protein. However, children with idiopathic short stature (ISS) show growth impairment without GH or GHR defects. We hypothesized that decreased expression of the GHR gene may be involved. To test this, we investigated whether common genetic variants (microsatellites, SNPs) in regulatory regions of the GHR gene region were associated with the ISS phenotype. Genotyping of a GT-repeat microsatellite in the GHR 5'UTR in a Montreal ISS cohort (n = 37 ISS, n = 105 controls) revealed that the incidence of the long/short (L/S) genotype was 3.3× higher in ISS children than controls (P = 0.04, OR = 3.85). In an Italian replication cohort (n = 143 ISS, n = 282 controls), the medium/short (M/S) genotype was 1.9× more frequent in the male ISS than controls (P = 0.017, OR = 2.26). In both ISS cohorts, logistic regression analysis of 27 SNPs showed an association of ISS with rs4292454, while haplotype analysis revealed specific risk haplotypes in the 3' haploblocks. In contrast, there were no differences in GT genotype frequencies in a cohort of short stature (SS) adults versus controls (CARTaGENE: n = 168 SS, n = 207 controls) and the risk haplotype in the SS cohort was located in the most 5' haploblock. These data suggest that the variants identified are potentially genetic markers specifically associated with the ISS phenotype.

Fine-scale variation and genetic determinants of alternative splicing across individuals.

Recently, thanks to the increasing throughput of new technologies, we have begun to explore the full extent of alternative pre-mRNA splicing (AS) in the human transcriptome. This is unveiling a vast layer of complexity in isoform-level expression differences between individuals. We used previously published splicing sensitive microarray data from lymphoblastoid cell lines to conduct an in-depth analysis on splicing efficiency of known and predicted exons. By combining publicly available AS annotation with a novel algorithm designed to search for AS, we show that many real AS events can be detected within the usually unexploited, speculative majority of the array and at significance levels much below standard multiple-testing thresholds, demonstrating that the extent of cis-regulated differential splicing between individuals is potentially far greater than previously reported. Specifically, many genes show subtle but significant genetically controlled differences in splice-site usage. PCR validation shows that 42 out of 58 (72%) candidate gene regions undergo detectable AS, amounting to the largest scale validation of isoform eQTLs to date. Targeted sequencing revealed a likely causative SNP in most validated cases. In all 17 incidences where a SNP affected a splice-site region, in silico splice-site strength modeling correctly predicted the direction of the micro-array and PCR results. In 13 other cases, we identified likely causative SNPs disrupting predicted splicing enhancers. Using Fst and REHH analysis, we uncovered significant evidence that 2 putative causative SNPs have undergone recent positive selection. We verified the effect of five SNPs using in vivo minigene assays. This study shows that splicing differences between individuals, including quantitative differences in isoform ratios, are frequent in human populations and that causative SNPs can be identified using in silico predictions. Several cases affected disease-relevant genes and it is likely some of these differences are involved in phenotypic diversity and susceptibility to complex diseases.

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.

Alternative splicing is frequent during early embryonic development in mouse.

BACKGROUND: Alternative splicing is known to increase the complexity of mammalian transcriptomes since nearly all mammalian genes express multiple pre-mRNA isoforms. However, our knowledge of the extent and function of alternative splicing in early embryonic development is based mainly on a few isolated examples. High throughput technologies now allow us to study genome-wide alternative splicing during mouse development. RESULTS: A genome-wide analysis of alternative isoform expression in embryonic day 8.5, 9.5 and 11.5 mouse embryos and placenta was carried out using a splicing-sensitive exon microarray. We show that alternative splicing and isoform expression is frequent across developmental stages and tissues, and is comparable in frequency to the variation in whole-transcript expression. The genes that are alternatively spliced across our samples are disproportionately involved in important developmental processes. Finally, we find that a number of RNA binding proteins, including putative splicing factors, are differentially expressed and spliced across our samples suggesting that such proteins may be involved in regulating tissue and temporal variation in isoform expression. Using an example of a well characterized splicing factor, Fox2, we demonstrate that changes in Fox2 expression levels can be used to predict changes in inclusion levels of alternative exons that are flanked by Fox2 binding sites. CONCLUSIONS: We propose that alternative splicing is an important developmental regulatory mechanism. We further propose that gene expression should routinely be monitored at both the whole transcript and the isoform level in developmental studies.

Sex-specific effects of a microsatellite polymorphism on human growth hormone receptor gene expression.

Growth hormone (GH) binds to its specific receptor (GHR) at the surface of target cells activating multiple signaling pathways implicated in growth and metabolism. Dysregulation of GHRs leads to pathophysiological states that most commonly affect stature. We previously showed the association of a polymorphic (n = 15-37) GT microsatellite in the human GHR gene promoter with short stature in a sex-specific manner. In the present study we evaluated the functional relevance of this polymorphism in regulating GHR expression. Using luciferase reporter assays, we found that the GT repeat had a significant cis regulatory effect in response to HIF1α and a potential repressor role following C/EBPß stimulation. Using a digital PCR application to measure allelic imbalance (AI), we showed a high prevalence of AI (∼76%) at the GHR locus in lymphoblastoid cell lines (LCLs), with a significantly higher degree of imbalance in LCLs derived from males. Examination of expression of GHR as well as other members of the GH-IGF1 axis in the LCLs revealed significant associations of GHR, IGF1 and BCL2 expression with GT genotype in a sex-specific manner. Our results suggest that this GT microsatellite exerts both cis and trans effects in a sex-specific context, revealing a new mechanism by which GHR gene expression is regulated.

Adenovirus-delivered antisense RNA and shRNA exhibit different silencing efficiencies for the endogenous transforming growth factor-beta (TGF-beta) type II receptor.

Gene silencing is an essential tool in gene discovery and gene therapy. Traditionally, viral delivery of antisense RNA and, more recently, small interfering RNA (siRNA) molecules in the form of small hairpin RNAs (shRNA) has been used as a strategy to achieve gene silencing. Nevertheless, the enduring challenge is to identify molecules that specifically and optimally silence a given target gene. In this study, we tested a set of adenovirus-delivered antisense RNA fragments and adenovirus-delivered shRNA molecules for their ability to target human transforming growth factor-beta type II receptor (TGFbetaRII). We used a dicistronic reporter, consisting of the coding sequences for TGFbetaRII and green fluorescent protein (GFP) to screen for optimal silencing agents targeting TGFbetaRII. Our results show, for both antisense RNA and shRNA molecules, that their effectiveness in the GFP screen correlated directly with their ability to reduce exogenously expressed TGFbetaRII. Unexpectedly, the antisense RNAs were unable to silence endogenous TGFbetaRII. In contrast, the shRNAs were able to silence endogenous TGFbetaRII. The shRNA that demonstrated the most pronounced effect on the dicistronic TGFbetaRII/GFP reporter reduced endogenous TGFbetaRII protein expression by 70% in A549 cells and reduced TGFbeta signaling by >80% in HeLa cells.

Exon-level transcriptome profiling in murine breast cancer reveals splicing changes specific to tumors with different metastatic abilities.

BACKGROUND: Breast cancer is the second most frequent type of cancer affecting women. We are increasingly aware that changes in mRNA splicing are associated with various characteristics of cancer. The most deadly aspect of cancer is metastasis, the process by which cancer spreads from the primary tumor to distant organs. However, little is known specifically about the involvement of alternative splicing in the formation of macroscopic metastases. Our study investigates transcript isoform changes that characterize tumors of different abilities to form growing metastases. METHODS AND FINDINGS: To identify alternative splicing events (ASEs) that are associated with the fully metastatic phenotype in breast cancer, we used Affymetrix Exon Microarrays to profile mRNA isoform variations genome-wide in weakly metastatic (168FARN and 4T07) and highly metastatic (4T1) mammary carcinomas. Statistical analysis identified significant expression changes in 7606 out of 155,994 (4%) exons and in 1725 out of 189,460 (1%) intronic regions, which affect 2623 out of 16,654 (16%) genes. These changes correspond to putative alternative isoforms-several of which are novel-that are differentially expressed between tumors of varying metastatic phenotypes. Gene pathway analysis showed that 1224 of genes expressing alternative isoforms were involved in cell growth, cell interactions, cell proliferation, cell migration and cell death and have been previously linked to cancers and genetic disorders. We chose ten predicted splice variants for RT-PCR validation, eight of which were successfully confirmed (MED24, MFI2, SRRT, CD44, CLK1 and HNRNPH1). These include three novel intron retentions in CD44, a gene in which isoform variations have been previously associated with the metastasis of several cancers. CONCLUSION: Our findings reveal that various genes are differently spliced and/or expressed in association with the metastatic phenotype of tumor cells. Identification of metastasis-specific isoforms may contribute to the development of improved breast cancer stage identification and targeted therapies.

Genome-wide analysis of transcript isoform variation in humans.

We have performed a genome-wide analysis of common genetic variation controlling differential expression of transcript isoforms in the CEU HapMap population using a comprehensive exon tiling microarray covering 17,897 genes. We detected 324 genes with significant associations between flanking SNPs and transcript levels. Of these, 39% reflected changes in whole gene expression and 55% reflected transcript isoform changes such as splicing variants (exon skipping, alternative splice site use, intron retention), differential 5' UTR (initiation of transcription) use, and differential 3' UTR (alternative polyadenylation) use. These results demonstrate that the regulatory effects of genetic variation in a normal human population are far more complex than previously observed. This extra layer of molecular diversity may account for natural phenotypic variation and disease susceptibility.

Heritability of alternative splicing in the human genome.

Alternative pre-mRNA splicing increases proteomic diversity and provides a potential mechanism underlying both phenotypic diversity and susceptibility to genetic disorders in human populations. To investigate the variation in splicing among humans on a genome-wide scale, we use a comprehensive exon-targeted microarray to examine alternative splicing in lymphoblastoid cell lines (LCLs) derived from the CEPH HapMap population. We show the identification of transcripts containing sequence verified exon skipping, intron retention, and cryptic splice site usage that are specific between individuals. A number of novel alternative splicing events with no previous annotations in either the RefSeq and EST databases were identified, indicating that we are able to discover de novo splicing events. Using family-based linkage analysis, we demonstrate Mendelian inheritance and segregation of specific splice isoforms with regulatory haplotypes for three genes: OAS1, CAST, and CRTAP. Allelic association was further used to identify individual SNPs or regulatory haplotype blocks linked to the alternative splicing event, taking advantage of the high-resolution genotype information from the CEPH HapMap population. In one candidate, we identified a regulatory polymorphism that disrupts a 5' splice site of an exon in the CAST gene, resulting in its exclusion in the mutant allele. This report illustrates that our approach can detect both annotated and novel alternatively spliced variants, and that such variation among individuals is heritable and genetically controlled.