Quantification of G protein Gaalphas subunit splice variants in different human tissues and cells using pyrosequencing

The role of genetic and epigenetic GNAS alterations in the development of early-onset obesity

BACKGROUND

GNAS (guanine nucleotide-binding protein, alpha stimulating) is an imprinted gene that encodes Gsα, the α subunit of the heterotrimeric stimulatory G protein. This subunit mediates the signalling of a diverse array of G protein-coupled receptors (GPCRs), including the melanocortin 4 receptor (MC4R) that serves a pivotal role in regulating food intake, energy homoeostasis, and body weight. Genetic or epigenetic alterations in GNAS are known to cause pseudohypoparathyroidism in its different subtypes and have been recently associated with isolated, early-onset, severe obesity. Given the diverse biological functions that Gsα serves, multiple molecular mechanisms involving various GPCRs, such as MC4R, β2- and β3-adrenoceptors, and corticotropin-releasing hormone receptor, have been implicated in the pathophysiology of severe, early-onset obesity that results from genetic or epigenetic GNAS changes.

SCOPE OF REVIEW

This review examines the structure and function of GNAS and provides an overview of the disorders that are caused by defects in this gene and may feature early-onset obesity. Moreover, it elucidates the potential molecular mechanisms underlying Gsα deficiency-induced early-onset obesity, highlighting some of their implications for the diagnosis, management, and treatment of this complex condition.

MAJOR CONCLUSIONS

Gsα deficiency is an underappreciated cause of early-onset, severe obesity. Therefore, screening children with unexplained, severe obesity for GNAS defects is recommended, to enhance the molecular diagnosis and management of this condition.

Mapping genetic changes in the cAMP-signaling cascade in human atria

AIM

To obtain a quantitative expression profile of the main genes involved in the cAMP-signaling cascade in human control atria and in different cardiac pathologies.

METHODS AND RESULTS

Expression of 48 target genes playing a relevant role in the cAMP-signaling cascade was assessed by RT-qPCR. 113 samples were obtained from right atrial appendages (RAA) of patients in sinus rhythm (SR) with or without atrium dilation, paroxysmal atrial fibrillation (AF), persistent AF or heart failure (HF); and left atrial appendages (LAA) from patients in SR or with AF. Our results show that right and left atrial appendages in donor hearts or from SR patients have similar expression values except for AC7 and PDE2A. Despite the enormous chamber-dependent variability in the gene-expression changes between pathologies, several distinguishable patterns could be identified. PDE8A, PI3Kγ and EPAC2 were upregulated in AF. Different phosphodiesterase (PDE) families showed specific pathology-dependent changes.

CONCLUSION

By comparing mRNA-expression patterns of the cAMP-signaling cascade related genes in right and left atrial appendages of human hearts and across different pathologies, we show that 1) gene expression is not significantly affected by cardioplegic solution content, 2) it is appropriate to use SR atrial samples as controls, and 3) many genes in the cAMP-signaling cascade are affected in AF and HF but only few of them appear to be chamber (right or left) specific.

TOPIC

Genetic changes in human diseased atria.

TRANSLATIONAL PERSPECTIVE

The cyclic AMP signaling pathway is important for atrial function. However, expression patterns of the genes involved in the atria of healthy and diseased hearts are still unclear. We give here a general overview of how different pathologies affect the expression of key genes in the cAMP signaling pathway in human right and left atria appendages. Our study may help identifying new genes of interest as potential therapeutic targets or clinical biomarkers for these pathologies and could serve as a guide in future gene therapy studies.

The GNAS SNP c.393C>T (rs7121) as a marker for disease progression and survival in cancer

G-protein receptor signaling plays a key role in multiple signal transduction pathways. Aberrant activity of the stimulatory G_sα subunit has been frequently associated with cancer. GNAS sequence alterations and conformational changes of G_sα can both enhance or diminish its function and change downstream effects of G-protein receptor signaling. In this review and meta-analysis, we focus on the synonymous SNP rs7121 (FokI, c.393C>T), which is associated with either tumor progression or prolonged survival in cancer patients (overall hazard ratio = 2.256; p < 0.001). We finally point out the relevance of GNAS rs7121 as a promising biomarker and a prediction tool for therapy response and the need of further experiments to implement it into routine clinical diagnostics.

A novel aspect of GNAS imprinting: higher maternal expression of Gαs in human lymphoblasts, peripheral blood mononuclear cells, mammary adipose tissue, and heart

The human GNAS gene is imprinted in a tissue-specific manner, being expressed primarily from the maternal allele in pituitary, thyroid, renal proximal tubules, and gonads, but is supposed to be biallelically expressed with an equal allelic expression in most other tissues. We analysed allelic expression of Gαs using Pyrosequencing. By genotyping the GNAS T393C polymorphism we quantified mRNA transcripts in lymphoblasts (Ly, n=11), peripheral blood mononuclear cells (PBMC, n=18), mammary adipose tissue (MAT, n=23) and heart tissue (HT, n=44). Allelic expression analysis revealed an unequal allelic expression (ratio maternal/total×100±SEM: 55.7±1% (95% CI 53.4-58.1%) in Ly, 56.1±0.8 (95% CI 54.5-57.7%) in PBMC, 54.5±0.8% (95% CI 53-56.1%) in MAT and 54.1±0.6% (95% CI 53-55.3%) in HT). Maternal ratio differed significantly from the mean (p<0.0001). This phenomenon may be a general feature existing in all tissues.

Skeletal progenitors and the GNAS gene: fibrous dysplasia of bone read through stem cells

Activating mutations of the GNAS gene, which causes fibrous dysplasia of bone (FD), lead to remarkable changes in the properties of skeletal progenitors, and it is these changes that mediate the pathological effect of this gene on bone. Mutated skeletal stem cells lose the ability to differentiate into adipocytes, and to maintain in situ, and transfer heterotopically, the hematopoietic microenvironment, leading to abnormal bone marrow histology in FD. They overexpress molecular effectors of osteoclastogenesis, thus promoting inappropriate bone resorption leading to fragility of FD bone. They express the phosphate-regulating hormone FGF-23 at normal levels, whose excess in the serum of FD patients correlates with the mass of osteogenic cells within FD lesions, leading to osteomalacia and deformity of the FD bone, and revealing that bone is an endocrine organ regulating renal handling of phosphate. Mechanisms of allelic selection and stem cell selection occur in mutated skeletal stem cells and contribute to the inherent diversity and evolution over time in FD. The definition of the etiological role of GNAS mutations marks the watershed between many decades of descriptive observation and the definition of cellular and molecular mechanisms that would explain and hopefully allow for a cure for the disease. Placing stem cells at center stage has permitted substantial advances in one decade, and promises more for the one to come.

Comparison of methods for quantification of subtle splice variants

Alternative splicing is capable of generating multiple mRNA variants from a single gene and is hence a key mediator of molecular diversity generated at the transcript level. Consequently, delivering quantitative information on the fractions of splice variants is essential for the understanding of their biological roles. Here we compare techniques for subtle splice variant quantification that are able to resolve length differences as small as one nucleotide: PAGE with ethidium-bromide densitometry, pyrosequencing, and CE-LIF. We give comprehensive descriptions of assay designs and calibration procedures and present an evaluation of these methods in terms of accuracy, reproducibility and applicability. We also examined template concentrations and reverse transcription-coupled PCR conditions as potential cause of biased results as they were observed for extreme low template concentrations and/or PCR amplicons with size differences of 195 nt. As proof of concept, we determine the splice ratios of variants differing by 3 and 12 nt in five human tissues. We demonstrate that CE-LIF is the most precise and also the most labor- and time-efficient method.

High resolution analysis of the human transcriptome: detection of extensive alternative splicing independent of transcriptional activity

BACKGROUND

Commercially available microarrays have been used in many settings to generate expression profiles for a variety of applications, including target selection for disease detection, classification, profiling for pharmacogenomic response to therapeutics, and potential disease staging. However, many commercially available microarray platforms fail to capture transcript diversity produced by alternative splicing, a major mechanism for driving proteomic diversity through transcript heterogeneity.

RESULTS

The human Genome-Wide SpliceArray(TM) (GWSA), a novel microarray platform, utilizes an existing probe design concept to monitor such transcript diversity on a genome scale. The human GWSA allows the detection of alternatively spliced events within the human genome through the use of exon body and exon junction probes to provide a direct measure of each transcript, through simple calculations derived from expression data. This report focuses on the performance and validation of the array when measured against standards recently published by the Microarray Quality Control (MAQC) Project. The array was shown to be highly quantitative, and displayed greater than 85% correlation with the HG-U133 Plus 2.0 array at the gene level while providing more extensive coverage of each gene. Almost 60% of splice events among genes demonstrating differential expression of greater than 3 fold also contained extensive splicing alterations. Importantly, almost 10% of splice events within the gene set displaying constant overall expression values had evidence of transcript diversity. Two examples illustrate the types of events identified: LIM domain 7 showed no differential expression at the gene level, but demonstrated deregulation of an exon skip event, while erythrocyte membrane protein band 4.1 -like 3 was differentially expressed and also displayed deregulation of a skipped exon isoform.

CONCLUSION

Significant changes were detected independent of transcriptional activity, indicating that the controls for transcript generation and transcription are distinct, and require novel tools in order to detect changes in specific transcript quantity. Our results demonstrate that the SpliceArray(TM) design will provide researchers with a robust platform to detect and quantify specific changes not only in overall gene expression, but also at the individual transcript level.

Analysis of splicing patterns by pyrosequencing

Several different mRNAs can be produced from a given pre-mRNA by regulated alternative splicing, or as the result of deregulations that may lead to pathological states. Analysing splicing patterns is therefore of importance to describe and understand developmental programs, cellular responses to internal or external cues, or human diseases. We describe here a method, Pyrosequencing Analysis of Splicing Patterns (PASP), that combines RT–PCR and pyrosequencing of PCR products. We demonstrated that: (i) Ratios of two pure RNAs mixed in various proportions were accurately measured by PASP; (ii) PASP can be adapted to virtually any splicing event, including mutually exclusive exons, complex patterns of exon skipping or inclusion, and alternative 3′ terminal exons; (iii) In extracts from different organs, the proportions of RNA isoforms measured by PASP reflected those measured by other methods. The PASP method is therefore reliable for analysing splicing patterns. All steps are done in 96-wells microplates, without gel electrophoresis, opening the way to high-throughput comparisons of RNA from several sources.

Violating the splicing rules: TG dinucleotides function as alternative 3' splice sites in U2-dependent introns

TG dinucleotides functioning as alternative 3' splice sites were identified and experimentally verified in 36 human genes.

Background

Despite some degeneracy of sequence signals that govern splicing of eukaryotic pre-mRNAs, it is an accepted rule that U2-dependent introns exhibit the 3' terminal dinucleotide AG. Intrigued by anecdotal evidence for functional non-AG 3' splice sites, we carried out a human genome-wide screen.

Results

We identified TG dinucleotides functioning as alternative 3' splice sites in 36 human genes. The TG-derived splice variants were experimentally validated with a success rate of 92%. Interestingly, ratios of alternative splice variants are tissue-specific for several introns. TG splice sites and their flanking intron sequences are substantially conserved between orthologous vertebrate genes, even between human and frog, indicating functional relevance. Remarkably, TG splice sites are exclusively found as alternative 3' splice sites, never as the sole 3' splice site for an intron, and we observed a distance constraint for TG-AG splice site tandems.

Conclusion

Since TGs splice sites are exclusively found as alternative 3' splice sites, the U2 spliceosome apparently accomplishes perfect specificity for 3' AGs at an early splicing step, but may choose 3' TGs during later steps. Given the tiny fraction of TG 3' splice sites compared to the vast amount of non-viable TGs, cis-acting sequence signals must significantly contribute to splice site definition. Thus, we consider TG-AG 3' splice site tandems as promising subjects for studies on the mechanisms of 3' splice site selection.

High qualitative and quantitative conservation of alternative splicing in Caenorhabditis elegans and Caenorhabditis briggsae

Alternative splicing (AS) is an important contributor to proteome diversity and is regarded as an explanatory factor for the relatively low number of human genes compared with less complex animals. To assess the evolutionary conservation of AS and its developmental regulation, we have investigated the qualitative and quantitative expression of 21 orthologous alternative splice events through the development of 2 nematode species separated by 85-110 Myr of evolutionary time. We demonstrate that most of these alternative splice events present in Caenorhabditis elegans are conserved in Caenorhabditis briggsae. Moreover, we find that relative isoform expression levels vary significantly during development for 78% of the AS events and that this quantitative variation is highly conserved between the 2 species. Our results suggest that AS is generally tightly regulated through development and that the regulatory mechanisms controlling AS are to a large extent conserved during the evolution of Caenorhabditis. This strong conservation indicates that both major and minor splice forms have important functional roles and that the relative quantities in which they are expressed are crucial. Our results therefore suggest that the quantitative regulation of isoform expression levels is an intrinsic part of most AS events. Moreover, our results indicate that AS contributes little to transcript variation in Caenorhabditis genes and that gene duplication may be the major evolutionary mechanism for the origin of novel transcripts in these 2 species.

Methods and platforms for the quantification of splice variants' expression

The relatively limited number of human protein encoding genes highlights the importance of the diversity generated at the level of the mRNA transcripts. As alternative RNA splicing plays a key role in mediating this diversity, it becomes critical to develop the tools and platforms that will deliver quantitative information on the specific expression levels associated with splice isoforms. This chapter describes the constraints generated by this global transcriptome analysis and the state-of-the-art techniques and products available to the scientific community.

Pyrosequencing protocol requiring a unique biotinylated primer

BACKGROUND

DNA sequencing has markedly changed the nature of biomedical research. Large-scale sequencing projects have generated several millions of potential polymorphisms widespread in the human genome requiring validation and incorporation into screening panels. As a consequence, high-throughput analysis of these variants in different populations of interest is now the cornerstone of structural genomics. Pyrosequencing is a versatile technique allowing an easy 96-well typing format. However, every polymorphism requires a specific labeled primer to generate a single-stranded DNA fragment containing the region of interest.

METHODS

We describe how with an adjusted primer stoichiometry we can standardize the labeling of every amplicon with a single biotinylated universal primer (BM13S).

RESULTS

We circumvent the need for specific biotinylated primers for each single-nucleotide polymorphism (SNP) under study. As an example, we assessed this novel protocol by genotyping three SNPs mapping calpain-10, caveolin-1 and CYP19A1.

CONCLUSION

The present approach represents an alternative to standard pyrosequencing protocols, since it requires a single biotinylated primer that is suitable for each SNP under study.