Formation of new genes explains lower intron density in mammalian Rhodopsin G protein-coupled receptors

Mammalian G protein-coupled receptor (GPCR) genes are characterised by a large proportion of intronless genes or a lower density of introns when compared with GPCRs of invertebrates. It is unclear which mechanisms have influenced intron density in this protein family, which is one of the largest in the mammalian genomes. We used a combination of Hidden Markov Models (HMM) and BLAST searches to establish the comprehensive repertoire of Rhodopsin GPCRs from seven species and performed overall alignments and phylogenetic analysis using the maximum parsimony method for over 1400 receptors in 12 subgroups. We identified 14 different Ancestral Receptor Groups (ARGs) that have members in both vertebrate and invertebrate species. We found that there exists a remarkable difference in the intron density among ancestral and new Rhodopsin GPCRs. The intron density among ARGs members was more than 3.5-fold higher than that within non-ARG members and more than 2-fold higher when considering only the 7TM region. This suggests that the new GPCR genes have been predominantly formed intronless while the ancestral receptors likely accumulated introns during their evolution. Many of the intron positions found in mammalian ARG receptor sequences were found to be present in orthologue invertebrate receptors suggesting that these intron positions are ancient. This analysis also revealed that one intron position is much more frequent than any other position and it is common for a number of phylogenetically different Rhodopsin GPCR groups. This intron position lies within a functionally important, conserved, DRY motif which may form a proto-splice site that could contribute to positional intron insertion. Moreover, we have found that other receptor motifs, similar to DRY, also contain introns between the second and third nucleotide of the arginine codon which also forms a proto-splice site. Our analysis presents compelling evidence that there was not a major loss of introns in mammalian GPCRs and formation of new GPCRs among mammals explains why these have fewer introns compared to invertebrate GPCRs. We also discuss and speculate about the possible role of different RNA- and DNA-based mechanisms of intron insertion and loss.     

The melanocortin system in Fugu: determination of POMC/AGRP/MCR gene repertoire and synteny, as well as pharmacology and anatomical distribution of the MCRs

The G-protein-coupled melanocortin receptors (MCRs) play an important role in a variety of essential functions such as the regulation of pigmentation, energy homeostasis, and steroid production. We performed a comprehensive characterization of the MC system in Fugu (Takifugu rubripes). We show that Fugu has an AGRP gene with high degree of conservation in the C-terminal region in addition to a POMC gene lacking gamma-MSH. The Fugu genome contains single copies of four MCRs, whereas the MC3R is missing. The MC2R and MC5R are found in tandem and remarkably contain one and two introns, respectively. We suggest that these introns were inserted through a reverse splicing mechanism into the DRY motif that is widely conserved through GPCRs. We were able to assemble large blocks around the MCRs in Fugu, showing remarkable synteny with human chromosomes 16 and 18. Detailed pharmacological characterization showed that ACTH had surprisingly high affinity for the Fugu MC1R and MC4R, whereas alpha-MSH had lower affinity. We also showed that the MC2R gene in Fugu codes for an ACTH receptor, which did not respond to alpha-MSH. All the Fugu receptors were able to couple functionally to cAMP production in line with the mammalian orthologs. The anatomical characterization shows that the MC2R is expressed in the brain in addition to the head-kidney, whereas the MC4R and MC5R are found in both brain regions and peripheral tissues. This is the first comprehensive genomic and functional characterization of a GPCR family within the Fugu genome. The study shows that some parts of the MC system are highly conserved through vertebrate evolution, such as regions in POMC coding for ACTH, alpha-MSH, and beta-MSH, the C-terminal region of AGRP, key binding units within the MC1R, MC2R, MC4R, and MC5R, synteny blocks around the MCRs, pharmacological properties of the MC2R, whereas other parts in the system are either missing, such as the MC3R and gamma-MSH, or different as compared to mammals, such as the affinity of ACTH and MSH peptides to MC1R and MC4R and the anatomical expression pattern of the MCRs.     

Major gender difference in association of FTO gene variant among severely obese children with obesity and obesity related phenotypes

Recent studies have shown that SNPs in the FTO gene predispose to childhood and adult obesity. In this study, we examined the association between variants in FTO and KIAA1005, a gene that maps closely to FTO, and obesity, as well as obesity related traits among 450 well characterised severely obese children and 512 normal weight controls. FTO showed significant association with several obesity related traits while SNPs in KIAA1005 did not. When stratified by gender, the FTO variant rs9939609 showed association with obesity and BMI among girls (P = 0.006 and 0.004, respectively) but not among boys. Gender differences were also found in the associations of the FTO rs9939609 with obesity related traits such as insulin sensitivity and plasma glucose. This study suggests that FTO may have an important role for gender specific development of severe obesity and insulin resistance in children.     

A Genome-Wide Analysis of Populations from European Russia Reveals a New Pole of Genetic Diversity in Northern Europe

Several studies examined the fine-scale structure of human genetic variation in Europe. However, the European sets analyzed represent mainly northern, western, central, and southern Europe. Here, we report an analysis of approximately 166,000 single nucleotide polymorphisms in populations from eastern (northeastern) Europe: four Russian populations from European Russia, and three populations from the northernmost Finno-Ugric ethnicities (Veps and two contrast groups of Komi people). These were compared with several reference European samples, including Finns, Estonians, Latvians, Poles, Czechs, Germans, and Italians. The results obtained demonstrated genetic heterogeneity of populations living in the region studied. Russians from the central part of European Russia (Tver, Murom, and Kursk) exhibited similarities with populations from central–eastern Europe, and were distant from Russian sample from the northern Russia (Mezen district, Archangelsk region). Komi samples, especially Izhemski Komi, were significantly different from all other populations studied. These can be considered as a second pole of genetic diversity in northern Europe (in addition to the pole, occupied by Finns), as they had a distinct ancestry component. Russians from Mezen and the Finnic-speaking Veps were positioned between the two poles, but differed from each other in the proportions of Komi and Finnic ancestries. In general, our data provides a more complete genetic map of Europe accounting for the diversity in its most eastern (northeastern) populations.     

Evidence for constitutive dimerization of niacin receptor subtypes

The recently deorphanized niacin receptor subtypes NIACR1 (GPR109A) and NIACR2 (GPR109B) play an essential role in the regulation of metabolic processes and immune reactions. Both receptors belong to the G-protein-coupled receptor (GPCR) family, whose members have traditionally been treated as monomeric entities, but now appear to exist and function as both homodimers and heterodimers. In this study, a close physical interaction is shown between the highly homologous niacin receptor subtypes, NIACR1 and NIACR2, using bioluminescence resonance energy transfer (BRET²) in living cells. The extent of homo- and hetero-dimerization of the niacin receptors did not vary after activation of the receptors with selective agonists, indicating that the dimerization state of NIACR1 and NIACR2 is not regulated by ligand binding. Moreover, detection of niacin receptor dimers in both plasma membrane- and endoplasmic reticulum-enriched fractions suggests that they are formed early in the biosynthetic pathway. Taken together, these results demonstrate that niacin receptor dimerization is a constitutive process occurring early during biosynthesis.     

A long-range pseudoknot in Qbeta RNA is essential for replication

A puzzling aspect of replication of bacteriophage Qbeta RNA has always been that replicase binds at an internal segment, the M-site, some 1450 nt away from the 3' end. Here, we report on the existence of a long-range pseudoknot, base-pairing eight nt in the loop of the 3' terminal hairpin to a single-stranded interdomain sequence located about 1200 nt upstream, close to the internal replicase binding site. Introduction of a single mismatch into this pseudoknot is sufficient to abolish replication, but the inhibition is fully reversed by a second-site substitution that restores the pairing. The pseudoknot is part of an elaborate structure that seems to hold the 3' end in a fixed position vis a vis the replicase binding site. Our results imply that the shape of the RNA confers the functonality. We discuss the possible relevance of our findings for replication of other viral RNAs.     

Identification and analysis of functionally important amino acids in human purinergic 12 receptor using a Saccharomyces cerevisiae expression system

The purinergic 12 receptor (P2Y12) is a major drug target for anticoagulant therapies, but little is known about the regions involved in ligand binding and activation of this receptor. We generated four randomized P2Y12 libraries and investigated their ligand binding characteristics. P2Y12 was expressed in a Saccharomyces cerevisiae model system. Four libraries were generated with randomized amino acids at positions 181, 256, 265 and 280. Mutant variants were screened for functional activity in yeast using the natural P2Y12 ligand ADP. Activation results were investigated using quantitative structure-activity relationship (QSAR) models and ligand-receptor docking. We screened four positions in P2Y12 for functional activity by substitution with amino acids with diverse physiochemical properties. This analysis revealed that positions E181, R256 and R265 alter the functional activity of P2Y12 in a specific manner. QSAR models for E181 and R256 mutant libraries strongly supported the experimental data. All substitutions of amino acid K280 were completely inactive, highlighting the crucial role of this residue in P2Y12 function. Ligand-receptor docking revealed that K280 is likely to be a key element in the ligand-binding pocket of P2Y12. The results of this study demonstrate that positions 181, 256, 265 and 280 of P2Y12 are important for the functional integrity of the receptor. Moreover, K280 appears to be a crucial feature of the P2Y12 ligand-binding pocket. These results are important for rational design of novel antiplatelet agents.     

Expression of human melanocortin 4 receptor in <Emphasis Type=Italic>Saccharomyces cerevisiae</Emphasis>

The melanocortin 4 receptor (MC4R) is involved in the regulation of energy homeostasis and is known as one of the major hypothalamic regulators of food intake. Several studies have shown that replacement of aspartic acid at position 126 of the MC4R abolishes the ligand binding. We used the modified yeast Saccharomyces cerevisiae strain MMY28 to functionally express the MC4R and characterise the importance of this amino acid for ligand based activation of the receptor. The efficiency of the functional expression system was estimated by activation with αMSH, ACTH and THIQ and compared with cAMP response in mammalian cells. We generated the library of MC4R mutants randomised at the amino acid position 126. Recombinant MC4R clones were screened for the αMSH induced activity in yeast. From 9 different amino acids obtained only the natural aspartic acid displayed the ligand dependent activity of MC4R. The MC4R variants with glutamic acid and leucine at position 126, however, displayed higher background activity than other amino acid substitutions. The results suggest that the yeast expression system is suitable for screening of the MC4R receptor ligands and that the substitution of aspartic acid at position 126 of MC4R by different amino acids functionally inactivates the receptor.