Search

Introduction: Shrunken pore syndrome (SPS), originally defined by cystatin C−based estimated glomerular filtration rate (eGFRcystatin C) being less than 60% of creatinine-based estimated glomerular filtration rate (eGFRcreatinine) in the absence of extrarenal influences on the plasma levels of cystatin C or creatinine, is associated with a high increase in mortality, even in the absence of reduced glomerular filtration rate (GFR). The objective of the present study was to determine whether the proteome of patients with SPS shows differences from that of patients with normal or reduced measured GFR (mGFR) without SPS. Methods: Four patient cohorts were included: 1 cohort with normal mGFR without SPS, 1 with normal mGFR with SPS, 1 with reduced mGFR without SPS, and 1 with reduced mGFR with SPS. The plasma levels of 177 selected proteins were analyzed. Results: Differences in the levels of 30 proteins were specific for SPS; 31 differences were specific for patients with both SPS and reduced mGFR; and 27 were specific for reduced mGFR. Eighteen of the differences specific for SPS concerned proteins described as promoting, or being associated with, atherosclerosis. Twelve of the differences specific for patients with both SPS and reduced mGFR and 10 of the differences specific for reduced mGFR also concerned proteins described as promoting, or being associated with, atherosclerosis. Almost all (82 of 88) of the concentration differences represented increased levels. For SPS, but not for reduced mGFR, a correlation between protein size and increase in level was observed, with smaller proteins being associated with higher levels. Conclusion: The high mortality in shrunken pore syndrome might be caused by the accumulation of atherosclerosis-promoting proteins in this condition. Keywords: atherosclerosis, creatinine, cystatin C, GFR, kidney, mortality

Summary A previous study revealed a strong association between the DMRT3:Ser301STOP mutation in horses and alternate gaits as well as performance in harness racing. Several follow-up studies have confirmed a high frequency of the mutation in gaited horse breeds and an effect on gait quality. The aim of this study was to determine when and where the mutation arose, to identify additional potential causal mutations and to determine the coalescence time for contemporary haplotypes carrying the stop mutation. We utilized sequences from 89 horses representing 26 breeds to identify 102 SNPs encompassing the DMRT3 gene that are in strong linkage disequilibrium with the stop mutation. These 102 SNPs were genotyped in an additional 382 horses representing 72 breeds, and we identified 14 unique haplotypes. The results provided conclusive evidence that DMRT3:Ser301STOP is causal, as no other sequence polymorphisms showed an equally strong association to locomotion traits. The low sequence diversity among mutant chromosomes demonstrated that they must have diverged from a common ancestral sequence within the last 10 000 years. Thus, the mutation occurred either just before domestication or more likely some time after domestication and then spread across the world as a result of selection on locomotion traits.

Linked loci and Galapagos finch size Observations of parallel evolution in the finches of the Galapagos, including body and beak size, contributed to Darwin's theories. Lamichhaney et al. carried out whole-genome sequencing of 60 Darwin's finches. These included small, medium, and large ground finches as well as small, medium, and large tree finches. A genomic region containing the HMGA2 gene correlated strongly with beak size across different species. This locus appears to have played a role in beak diversification throughout the radiation of Darwin's finches. Science , this issue p. 470

Single nucleotide polymorphisms (SNPs) within a genetic region including the first two introns of the gene encoding FTO have consistently been shown to be the strongest genetic factors influencing body mass index (BMI). However, this same also contains several regulatory DNA elements that affect the expression of IRX3 and IRX5, which respectively, are located approximately 500 kb and 1.2 Mbp downstream from the BMI-associated FTO locus. Through these affected regulatory elements, genetic variation at the FTO locus influences adipocyte development leading to decreased thermogenesis and increased lipid storage. These findings provide a genomic model for the functional implications of genetic variations at this locus, and also demonstrate the importance of accounting for chromatin-chromatin interactions when constructing hypotheses for the mechanisms of trait and disease-associated common genetic variants. Several consortia have generated genome-wide datasets describing different aspects of chromatin biology which can be utilized to predict functionality and propose biologically relevant descriptions of specific DNA regions. Here, we review some of the publically available data resources on genome function and organization that can be used to gain an overview of genetic regions of interest and to generate testable hypotheses for future studies. We use the BMI- and obesity-associated FTO locus as a subject as it poses an illustrative example on the value of these resources. We find that public databases strongly support long-range interactions between regulatory elements in the FTO locus with the IRXB cluster genes IRX3 and IRX5. Chromatin configuration capture data also support interactions across a large region stretching across from the RPGRIP1L gene, FTO and the IRXB gene cluster.

Darwin's finches, inhabiting the Galápagos archipelago and Cocos Island, constitute an iconic model for studies of speciation and adaptive evolution. Here we report the results of whole-genome re-sequencing of 120 individuals representing all of the Darwin's finch species and two close relatives. Phylogenetic analysis reveals important discrepancies with the phenotype-based taxonomy. We find extensive evidence for interspecific gene flow throughout the radiation. Hybridization has given rise to species of mixed ancestry. A 240 kilobase haplotype encompassing the ALX1 gene that encodes a transcription factor affecting craniofacial development is strongly associated with beak shape diversity across Darwin's finch species as well as within the medium ground finch (Geospiza fortis), a species that has undergone rapid evolution of beak shape in response to environmental changes. The ALX1 haplotype has contributed to diversification of beak shapes among the Darwin's finches and, thereby, to an expanded utilization of food resources.

Recent advances in comparative genomics have revealed that divergence between populations is not necessarily uniform across all parts of the genome. There are examples of regions with divergent haplotypes that are substantially more different from each other that the genomic average.Typically, these regions are of interest, as their persistence over long periods of time may reflect balancing selection. However, they are hard to detect unless the divergent sub-populations are known prior to analysis.Here, we introduce HaploDistScan, an R-package implementing model-free detection of deep-divergence genomic regions based on the distribution of pair-wise haplotype distances, and show that it can detect such regions without use of a priori information about population sub-division. We apply the method to real-world data sets, from ruff and Darwin’s finches, and show that we are able to recover known instances of balancing selection – originally identified in studies reliant on detailed phenotyping – using only genotype data. Furthermore, in addition to replicating previously known divergent haplotypes as a proof-of-concept, we identify novel regions of interest in the Darwin’s finch genome and propose a plausible, data-driven evolutionary history for each novel locus individually.In conclusion, HaploDistScan requires neither phenotypic nor demographic input data, thus filling a gap in the existing set of methods for genome scanning, and provides a useful tool for identification of regions under balancing selection or similar evolutionary processes.

The dramatic increase in membrane proteome complexity is arguably one of the most pivotal evolutionary events that underpins the origin of multicellular animals. However, the origin of a significant number of membrane families involved in metazoan development has not been clarified. In this study, we have manually curated the membrane proteomes of 22 metazoan and 2 unicellular holozoan species. We identify 123,014 membrane proteins in these 24 eukaryotic species and classify 86% of the dataset. We determine 604 functional clusters that are present from the last holozoan common ancestor (LHCA) through many metazoan species. Intriguingly, we show that more than 70% of the metazoan membrane protein families have a premetazoan origin. The data show that enzymes are more highly represented in the LHCA and expand less than threefold throughout metazoan species; in contrast to receptors that are relatively few in the LHCA but expand nearly eight fold within metazoans. Expansions related to cell adhesion, communication, immune defence, and developmental processes are shown in conjunction with emerging biological systems, such as neuronal development, cytoskeleton organization, and the adaptive immune response. This study defines the possible LHCA membrane proteome and describes the fundamental functional clusters that underlie metazoan diversity and innovation.

The type and the amount of dietary fat have a significant influence on the metabolic pathways involved in the development of obesity, metabolic syndrome, diabetes type 2 and cardiovascular diseases. However, it is unknown to what extent this modulation is achieved through DNA methylation. We assessed the effects of cholesterol intake, the proportion of energy intake derived from fat, the ratio of polyunsaturated fatty acids (PUFA) to saturated fatty acids (SFA), the ratio of monounsaturated fatty acids (MUFA) to SFA, and the ratio of MUFA+PUFA to SFA on genome-wide DNA methylation patterns in normal-weight and obese children. We determined the genome-wide methylation profile in the blood of 69 Greek preadolescents (∼10 years old) as well as their dietary intake for two consecutive weekdays and one weekend day. The methylation levels of one CpG island shore and four sites were significantly correlated with total fat intake. The methylation levels of 2 islands, 11 island shores and 16 sites were significantly correlated with PUFA/SFA; of 9 islands, 26 island shores and 158 sites with MUFA/SFA; and of 10 islands, 40 island shores and 130 sites with (MUFA+PUFA)/SFA. We found significant gene enrichment in 34 pathways for PUFA/SFA, including the leptin pathway, and a significant enrichment in 5 pathways for (MUFA+PUFA)/SFA. Our results suggest that specific changes in DNA methylation may have an important role in the mechanisms involved in the physiological responses to different types of dietary fat.

Variations in the FTO gene and near the TMEM18 gene are risk factors for common form of obesity, but have also been linked with type 2 diabetes (T2D). Our aim was to investigate the contribution of these variants to risk of T2D in a population in Latvia. Four single nucleotide polymorphisms (SNP) in the first and fourth intronic regions of FTO and one close to TMEM18 were genotyped in 987 patients with T2D and 1080 controls selected from the Latvian Genome Data Base (LGDB). We confirmed association of SNPs in the first intron (rs11642015, rs62048402 and rs9939609) of FTO and rs7561317 representing the TMEM18 locus with T2D. Association between SNP in FTO and T2D remained significant after correction for body mass index (BMI). The rs57103849 located in the fourth intron of FTO and rs7561317 in TMEM18 showed BMI independent association with younger age at diagnosis of T2D. Our results add to the evidence that BMI related variants in and near FTO and TMEM18 may increase the risk for T2D not only through secondary effects of obesity. The influence of variants in the fourth intron of the FTO gene on development of T2D may be mediated by mechanisms other than those manifested by SNPs in the first intron of the same gene.

Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation. DOI: http://dx.doi.org/10.7554/eLife.12081.001, eLife digest The Atlantic herring is one of the most common fish in the world and has been a crucial food resource in northern Europe. One school of herring may comprise billions of fish, but previous studies had only revealed very few genetic differences in herring from different geographic regions. This was unexpected since Atlantic herring is one of the few marine species that can reproduce throughout the brackish Baltic Sea, which can be about a tenth as salty as the Atlantic Ocean. This unexpected finding could be explained in at least two different ways. Firstly, perhaps Atlantic herring are flexible enough to adapt to very different environments (i.e. high or low salinity) without much genetic change. Secondly, the previous studies only looked at a handful of sites in the Atlantic herring’s genome and so it is possible that genetic differences at other genes control this fish’s adaptation instead. Now, Martinez Barrio, Lamichhaney, Fan, Rafati et al. have sequenced entire genomes from groups of Atlantic herring and revealed hundreds of sites that are associated with adaptation to the Baltic Sea. The analysis also identified a number of genes that control when these fish reproduce by comparing herring that spawn in the autumn with those that spawn in spring. This is important because natural populations must carefully time when they reproduce to maximize the survival of their young. These new findings provide compelling evidence that changes in protein-coding genes and stretches of DNA that regulate the expression of other genes both contribute to adaptation in herrings. The analysis also clearly shows that variants of genes that contribute to adaptation were likely to evolve over time by accumulating multiple sequence changes affecting the same gene. Furthermore, these gene variants essentially form a rich “tool-box” that underlies the Atlantic herring’s adaptation to its environment, and different subpopulations of herring were found to have their own optimal sets of gene variants. For instance, autumn-spawning herring and spring-spawning herring from the Baltic Sea both have gene variants that favor adaptation to low salinity. However, autumn-spawning Baltic herring also share gene variants that favor spawning in the autumn with autumn-spawning herring from the North Sea, but not with spring-spawning Baltic herring. The next step will be to study how the 500 or so genes identified affect adaptation at the molecular level. This will likely involve experiments with other model fish such as zebrafish and sticklebacks. Finally, these new findings can be directly applied to monitor stocks of herring to make herring fisheries more sustainable. DOI: http://dx.doi.org/10.7554/eLife.12081.002

Genome-wide association studies have identified a number of genes associated with human body weight. While some of these genes are large fields within obesity research, such as MC4R, POMC, FTO and BDNF, the majority do not have a clearly defined functional role explaining why they may affect body weight. Here, we searched biological databases and discovered 33 additional genes associated with human obesity (CADM2, GIPR, GPCR5B, LRP1B, NEGR1, NRXN3, SH2B1, FANCL, GNPDA2, HMGCR, MAP2K5, NUDT3, PRKD1, QPCTL, TNNI3K, MTCH2, DNAJC27, SLC39A8, MTIF3, RPL27A, SEC16B, ETV5, HMGA1, TFAP2B, TUB, ZNF608, FAIM2, KCTD15, LINGO2, POC5, PTBP2, TMEM18, TMEM160). We find that the majority have orthologues in distant species, such as D. melanogaster and C. elegans, suggesting that they are important for the biology of most bilateral species. Intriguingly, signalling cascade genes and transcription factors are enriched among these obesity genes, and several of the genes show properties that could be useful for potential drug discovery. In this review, we demonstrate how information from several distant model species, interactomics and signalling pathway analysis represents an important way to better understand the functional diversity of the surprisingly high number of molecules that seem to be important for human obesity.

Transmembrane proteins are gatekeepers of the cells; controlling the transport of substrates as well as communicating signals among cells and between the organelles and cytosol. Solute carriers (SLC) and G protein-coupled receptors (GPCR) are the largest family of membrane transporters and membrane receptors respectively. The overall aim of this thesis was to provide a basic understanding of some of the novel SLCs and GPCRs with emphasis on expression, transport property, evolution and probable function. The first part of the thesis directs towards the study of some novel solute carriers. In an initial study, we provided an overall picture of the sequence relationship and tissue expression of 14 diverse atypical SLCs confirming some of their evolutionary conservation and highly specific expression pattern. The focus then was on the SLC17 family (mainly vesicular proteins) and a novel member named Slc17a9. This study revealed that SLC17 family could be divided into four main phylogenetic clades which were all present before the divergence of the insect lineage with Slc17a9 having the most restricted evolutionary history. Detailed expression study of Slc17a9 in the mouse brain suggests that it is also expressed in some regions important for purinergic neurotransmission. Further, we deorphanised an aminoacid transporter Slc38a7 which was expressed in a majority of neurons in the CNS and showed that it preferably mediate transport of L–glutamine and L–histidine. The second part of the thesis focuses on the study of two GPCRs belonging to the Rhodopsin superfamily, Gpr162 and Gpr153. A phylogenetic analysis revealed that both Gpr153 and Gpr162 originated from a common ancestor before the radiation of the mammalian lineage. Expression study revealed that Gpr162 had a predominant expression in the CNS and relatively lower expression in the other tissue tested whereas Gpr153 had a more widespread and similar expression pattern in both CNS and peripheral tissues. The functional studies of the two GPCRs were done using the antisense oligodeoxynucleotide knockdown rat model. These studies provided evidence linking the orphan Gpr162 gene with the regulation of food intake– related behaviour whereas Gpr153 gene caused only a slight reduction in food intake.

The discovery and exploitation of new drug targets is a key focus for both the pharmaceutical industry and academic biomedical research. To provide an insight into trends in the exploitation of new drug targets, we have analysed the drugs that were approved by the US Food and Drug Administration during the past three decades and examined the interactions of these drugs with therapeutic targets that are encoded by the human genome, using the DrugBank database and extensive manual curation. We have identified 435 effect-mediating drug targets in the human genome, which are modulated by 989 unique drugs, through 2,242 drug-target interactions. We also analyse trends in the introduction of drugs that modulate previously unexploited targets, and discuss the network pharmacology of the drugs in our data set.

We recently used genome sequencing to study the evolutionary history of the Darwin's finches. A prominent feature of our data was that different polymorphic sites in the genome tended to indicate different genetic relationships among these closely related species. Such patterns are expected in recently diverged genomes as a result of incomplete lineage sorting. However, we uncovered conclusive evidence that these patterns have also been influenced by interspecies hybridisation, a process that has likely played an important role in the radiation of Darwin's finches. A major discovery was that segregation of two haplotypes at the ALX1 locus underlies variation in beak shape among the Darwin's finches, and that differences between the two haplotypes in a 240 kb region in blunt and pointed beaked birds involve both coding and regulatory changes. As we review herein, the evolution of such adaptive haplotypes comprising multiple causal changes appears to be an important mechanism contributing to the evolution of biodiversity.

Background The mechanisms by which genetic variants, such as single nucleotide polymorphisms (SNPs), identified in genome-wide association studies act to influence body mass remain unknown for most of these SNPs, which continue to puzzle the scientific community. Recent evidence points to the epigenetic and chromatin states of the genome as having important roles. Methods We genotyped 355 healthy young individuals for 52 known obesity-associated SNPs and obtained DNA methylation levels in their blood using the Illumina 450 K BeadChip. Associations between alleles and methylation at proximal cytosine residues were tested using a linear model adjusted for age, sex, weight category, and a proxy for blood cell type counts. For replication in other tissues, we used two open-access datasets (skin fibroblasts, n = 62; four brain regions, n = 121–133) and an additional dataset in subcutaneous and visceral fat (n = 149). Results We found that alleles at 28 of these obesity-associated SNPs associate with methylation levels at 107 proximal CpG sites. Out of 107 CpG sites, 38 are located in gene promoters, including genes strongly implicated in obesity (MIR148A, BDNF, PTPMT1, NR1H3, MGAT1, SCGB3A1, HOXC12, PMAIP1, PSIP1, RPS10-NUDT3, RPS10, SKOR1, MAP2K5, SIX5, AGRN, IMMP1L, ELP4, ITIH4, SEMA3G, POMC, ADCY3, SSPN, LGR4, TUFM, MIR4721, SULT1A1, SULT1A2, APOBR, CLN3, SPNS1, SH2B1, ATXN2L, and IL27). Interestingly, the associated SNPs are in known eQTLs for some of these genes. We also found that the 107 CpGs are enriched in enhancers in peripheral blood mononuclear cells. Finally, our results indicate that some of these associations are not blood-specific as we successfully replicated four associations in skin fibroblasts. Conclusions Our results strongly suggest that many obesity-associated SNPs are associated with proximal gene regulation, which was reflected by association of obesity risk allele genotypes with differential DNA methylation. This study highlights the importance of DNA methylation and other chromatin marks as a way to understand the molecular basis of genetic variants associated with human diseases and traits. Electronic supplementary material The online version of this article (doi:10.1186/s13073-015-0225-4) contains supplementary material, which is available to authorized users.

SummaryGenome-wide association studies (GWAS) have revealed association of a locus approximately 25b downstream of the TMEM18 gene with body mass and obesity. We utilized targeted re-sequencing of the body mass associated locus in proximity of TMEM18 in a case-control population of severely obese children and adolescents from the Stockholm area. We expanded our study to include the TMEM18 gene itself, with the aim of identifying body mass associated genetic variants. Sequencing was performed on the SOLiD platform, on long-range PCR fragments generated through targeted amplification of the regions of interest. Candidate single nucleotide polymorphisms (SNPs) were validated by TaqMan genotyping. We were able to observe 131 SNPs across the re-sequenced regions. Chi squared tests comparing the allele frequencies between cases and controls revealed 57 SNPs as candidates for association with obesity. Validation and replication genotyping revealed robust associations for SNPs within the haplotype block region located downstream from the TMEM18 gene. This study provides a high resolution map of the genetic variation pattern in the TMEM18 gene, as well as the associated haplotype block, and further strengthens the association of variants within the proximal haplotype block with obesity and body mass.

Wild and domesticated Atlantic salmon males display large variation for sea age at sexual maturation, which varies between 1-5 years. Previous studies have uncovered a genetic predisposition for age at maturity with moderate heritability, thus suggesting a polygenic or complex nature of this trait. The aim of this study was to identify associated genetic loci, genes and ultimately specific sequence variants conferring sea age at maturity in salmon. We performed a GWAS using a pool sequencing approach (20 individuals per river and trait) of salmon returning to rivers as sexually mature either after one sea winter (2009) or three sea winters (2011) in six rivers in Norway. The study revealed one major selective sweep, which covered 76 significant SNP in which 74 were found in a 370 kb region of chromosome 25. Genotyping other smolt year classes of wild salmon and domesticated salmon confirmed this finding. Genotyping domesticated fish narrowed the haplotype region to four SNPs covering 2386 bp, containing thevgll3gene, including two missense mutations explaining 33-36% phenotypic variation. This study demonstrates a single locus playing a highly significant role in governing sea age at maturation in this species. The SNPs identified may be both used as markers to guide breeding for late maturity in salmon aquaculture and in monitoring programs of wild salmon. Interestingly, a SNP in proximity of the VGLL3 gene in human (Homo sapiens), has previously been linked to age at puberty suggesting a conserved mechanism for timing of puberty in vertebrates.

Background Kiwi, comprising five species from the genus Apteryx, are endangered, ground-dwelling bird species endemic to New Zealand. They are the smallest and only nocturnal representatives of the ratites. The timing of kiwi adaptation to a nocturnal niche and the genomic innovations, which shaped sensory systems and morphology to allow this adaptation, are not yet fully understood. Results We sequenced and assembled the brown kiwi genome to 150-fold coverage and annotated the genome using kiwi transcript data and non-redundant protein information from multiple bird species. We identified evolutionary sequence changes that underlie adaptation to nocturnality and estimated the onset time of these adaptations. Several opsin genes involved in color vision are inactivated in the kiwi. We date this inactivation to the Oligocene epoch, likely after the arrival of the ancestor of modern kiwi in New Zealand. Genome comparisons between kiwi and representatives of ratites, Galloanserae, and Neoaves, including nocturnal and song birds, show diversification of kiwi’s odorant receptors repertoire, which may reflect an increased reliance on olfaction rather than sight during foraging. Further, there is an enrichment of genes influencing mitochondrial function and energy expenditure among genes that are rapidly evolving specifically on the kiwi branch, which may also be linked to its nocturnal lifestyle. Conclusions The genomic changes in kiwi vision and olfaction are consistent with changes that are hypothesized to occur during adaptation to nocturnal lifestyle in mammals. The kiwi genome provides a valuable genomic resource for future genome-wide comparative analyses to other extinct and extant diurnal ratites. Electronic supplementary material The online version of this article (doi:10.1186/s13059-015-0711-4) contains supplementary material, which is available to authorized users.

Context DNA methylation has been proposed to play a critical role in many cellular and biological processes. Objective To examine the influence of Roux-en-Y gastric bypass (RYGB) surgery on genome-wide promoter-specific DNA methylation in obese patients. Promoters are involved in the initiation and regulation of gene transcription. Methods Promoter-specific DNA methylation in whole blood was measured in 11 obese patients (presurgery BMI >35 kg/m2, 4 females), both before and 6 months after RYGB surgery, as well as once only in a control group of 16 normal-weight men. In addition, body weight and fasting plasma glucose were measured after an overnight fast. Results The mean genome-wide distance between promoter-specific DNA methylation of obese patients at six months after RYGB surgery and controls was shorter, as compared to that at baseline (p

Heterotrimeric G proteins perform a crucial role as molecular switches controlling various cellular responses mediated by G protein-coupled receptor (GPCR) signaling pathway. Recent data have shown that the vertebrate-like G protein families are found across metazoans and their closest unicellular relatives. However, an overall evolutionary hierarchy of vertebrate-like G proteins, including gene family annotations and in particular mapping individual gene gain/loss events across diverse holozoan lineages is still incomplete. Here, with more expanded invertebrate taxon sampling, we have reconstructed phylogenetic trees for each of the G protein classes/families and provide a robust classification and hierarchy of vertebrate-like heterotrimeric G proteins. Our results further extend the evidence that the common ancestor (CA) of holozoans had at least five ancestral Gα genes corresponding to all major vertebrate Gα classes and contain a total of eight genes including two Gβ and one Gγ. Our results also indicate that the GNAI/O-like gene likely duplicated in the last CA of metazoans to give rise to GNAI- and GNAO-like genes, which are conserved across invertebrates. Moreover, homologs of GNB1-4 paralogon- and GNB5 family-like genes are found in most metazoans and that the unicellular holozoans encode two ancestral Gβ genes. Similarly, most bilaterian invertebrates encode two Gγ genes which include a representative of the GNG gene cluster and a putative homolog of GNG13. Interestingly, our results also revealed key evolutionary events such as the Drosophila melanogaster eye specific Gβ subunit that is found conserved in most arthropods and several previously unidentified species specific expansions within Gαi/o, Gαs, Gαq, Gα12/13 classes and the GNB1-4 paralogon. Also, we provide an overall proposed evolutionary scenario on the expansions of all G protein families in vertebrate tetraploidizations. Our robust classification/hierarchy is essential to further understand the differential roles of GPCR/G protein mediated intracellular signaling system across various metazoan lineages.

Genome-wide association studies (GWAS) have identified common genetic factors influencing body mass as well as body adiposity. The functional implications of these loci are currently under investigation. Intense scrutiny of the body mass-associated FTO locus revealed age-specific effects, or a weakened effect in elderly populations. In this study, we aimed to determine the effects of single nucleotide polymorphisms (SNPs) representing 35 GWAS-identified body mass- and adiposity-associated genetic loci. In our analysis, 949 participants of the Prospective Investigation of the Vasculature in Uppsala Seniors cohort were included. All participants were born between 1920 and 1924. Data were available for 474 male and 475 female participants at age 70 and 380 male and 390 female participants at age 75. Genetic associations with BMI and change in BMI from age 70 to 75 were analyzed. In our analysis, rs10968576, an intronic SNP within the LINGO2 (LERN3, LRRN6C) gene, was associated with body mass in a cross section of elderly Swedes at age 70. This is the first study to replicate the association of a LINGO2-related genetic variant with body mass in an independent cohort of elderly citizens.

The combination of the obesity epidemic and an aging population presents growing challenges for the healthcare system. Obesity and aging are major risk factors for a diverse number of diseases and it is of importance to understand their interaction and the underlying molecular mechanisms. Herein the authors examined the methylation levels of 27578 CpG sites in 46 samples from adult peripheral blood. The effect of obesity and aging was ascertained with general linear models. More than one hundred probes were correlated to aging, nine of which belonged to the KEGG group map04080. Additionally, 10 CpG sites had diverse methylation profiles in obese and lean individuals, one of which was the telomerase catalytic subunit (TERT). In eight of ten cases the methylation change was reverted between obese and lean individuals. One region proved to be differentially methylated with obesity (LINC00304) independent of age. This study provides evidence that obesity influences age driven epigenetic changes, which provides a molecular link between aging and obesity. This link and the identified markers may prove to be valuable biomarkers for the understanding of the molecular basis of aging, obesity and associated diseases.

Spinal muscular atrophy (SMA) is a monogenic disorder that is subdivided into four different types and caused by survival motor neuron gene 1 (SMN1) deletion. Discordant cases of SMA suggest that there exist additional severity modifying factors, apart from the SMN2 gene copy number. Here we performed the first genome-wide methylation profiling of SMA patients and healthy individuals to study the association of DNA methylation status with the severity of the SMA phenotype. We identified strong significant differences in methylation level between SMA patients and healthy controls in CpG sites close to the genes CHML, ARHGAP22, CYTSB, CDK2AP1 and SLC23A2. Interestingly, the CHML and ARHGAP22 genes are associated with the activity of Rab and Rho GTPases, which are important regulators of vesicle formation, actin dynamics, axonogenesis, processes that could be critical for SMA development. We suggest that epigenetic modifications may influence the severity of SMA and that these novel genetic positions could prove to be valuable biomarkers for the understanding of SMA pathogenesis.

The melanocortin 2 receptor (MC2R) accessory proteins, MRAP, along with its homolog, MRAP2, are two among a growing number of G protein-coupled receptor accessory proteins that have been identified in recent years. These proteins interact directly with MC2R and are essential for trafficking of this receptor from the endoplasmic reticulum to the cell surface, where it mediates the effects of ACTH. lthough earlier studies have identified MRAP and MRAP2 subtypes in distant species, an overall evolutionary analysis of these families is still missing. Here, we performed a comprehensive evolutionary analysis of the MRAP and MRAP2 homologs based on whole genome sequences. We systematically mined and analyzed the genomes of metazoans to identify these genes. Overall, we identified 70 sequences of MRAP and MRAP2 from 44 species belonging to several vertebrate lineages, including at least 40 new sequences previously not reported in the literature. Herein, we provide evidence that MRAP2 is likely to be the ancestor of the MRAP family because MRAP2-like protein, but not MRAP, was identified in Petromyzon marinus (sea lamprey), which belong to an ancient basal vertebrate lineage. Later in vertebrate evolution, MRAP2 duplicated and gave rise to MRAP in an event before the emergence of actinopterygii (ray-finned fishes). However, we observed losses of MRAP in sarcopterygii (lobe-finned fish), amphibians and reptiles while both subtypes are present in chicken and most mammals studied. Synteny analysis showed a conserved synteny within same lineages and an inversion of gene order between lineages. An evolutionary rate shift analysis indicated that these genes were under high purifying selection. Overall, this study provides a comprehensive analysis of the evolution and gene repertoire of MRAP and MRAP2.

Saccoglossus kowalevskii (the acorn worm) is a hemichordate belonging to the superphylum of deuterostome bilateral animals. Hemichordates are sister group to echinoderms, and closely related to chordates. S. kowalevskii has chordate like morphological traits and serves as an important model organism, helping developmental biologists to understand the evolution of the central nervous system (CNS). Despite being such an important model organism, the signalling system repertoire of the largest family of integral transmembrane receptor proteins, G protein-coupled receptors (GPCRs) is largely unknown in S. kowalevskii. Here, we identified 260 unique GPCRs and classified as many as 257 of them into five main mammalian GPCR families; Glutamate (23), Rhodopsin (212), Adhesion (18), Frizzled (3) and Secretin (1). Despite having a diffuse nervous system, the acorn worm contains well conserved orthologues for human Adhesion and Glutamate family members, with a similar N-terminal domain architecture. This is particularly true for genes involved in CNS development and regulation in vertebrates. The average sequence identity between the GPCR orthologues in human and S. kowalevskii is around 47%, and this is same as observed in couple of the closest vertebrate relatives, Ciona intestinalis (41%) and Branchiostoma floridae (~47%). The Rhodopsin family has fewer members than vertebrates and lacks clear homologues for 6 of the 13 subgroups, including olfactory, chemokine, prostaglandin, purine, melanocyte concentrating hormone receptors and MAS-related receptors. However, the peptide and somatostatin binding receptors have expanded locally in the acorn worm. Overall, this study is the first large scale analysis of a major signalling gene superfamily in the hemichordate lineage. The establishment of orthologue relationships with genes involved in neurotransmission and development of the CNS in vertebrates provides a foundation for understanding the evolution of signal transduction and allows for further investigation of the hemichordate neurobiology.

The Interferon induced transmembrane proteins (IFITM) are a family of transmembrane proteins that is known to inhibit cell invasion of viruses such as HIV-1 and influenza. We show that the IFITM genes are a subfamily in a larger family of transmembrane (TM) proteins that we call Dispanins, which refers to a common 2TM structure. We mined the Dispanins in 36 eukaryotic species, covering all major eukaryotic groups, and investigated their evolutionary history using Bayesian and maximum likelihood approaches to infer a phylogenetic tree. We identified ten human genes that together with the known IFITM genes form the Dispanin family. We show that the Dispanins first emerged in eukaryotes in a common ancestor of choanoflagellates and metazoa, and that the family later expanded in vertebrates where it forms four subfamilies (A-D). Interestingly, we also find that the family is found in several different phyla of bacteria and propose that it was horizontally transferred to eukaryotes from bacteria in the common ancestor of choanoflagellates and metazoa. The bacterial and eukaryotic sequences have a considerably conserved protein structure. In conclusion, we introduce a novel family, the Dispanins, together with a nomenclature based on the evolutionary origin.

Background The Fat mass and obesity gene (FTO) has been identified through genome wide association studies as an important genetic factor contributing to a higher body mass index (BMI). However, the molecular context in which this effect is mediated has yet to be determined. We investigated the potential molecular network for FTO by analyzing co-expression and protein-protein interaction databases, Coxpresdb and IntAct, as well as the functional coupling predicting multi-source database, FunCoup. Hypothalamic expression of FTO-linked genes defined with this bioinformatics approach was subsequently studied using quantitative real time-PCR in mouse feeding models known to affect FTO expression. Results We identified several candidate genes for functional coupling to FTO through database studies and selected nine for further study in animal models. We observed hypothalamic expression of Profilin 2 (Pfn2), cAMP-dependent protein kinase catalytic subunit beta (Prkacb), Brain derived neurotrophic factor (Bdnf), neurotrophic tyrosine kinase, receptor, type 2 (Ntrk2), Signal transducer and activator of transcription 3 (Stat3), and Btbd12 to be co-regulated in concert with Fto. Pfn2 and Prkacb have previously not been linked to feeding regulation. Conclusions Gene expression studies validate several candidates generated through database studies of possible FTO-interactors. We speculate about a wider functional role for FTO in the context of current and recent findings, such as in extracellular ligand-induced neuronal plasticity via NTRK2/BDNF, possibly via interaction with the transcription factor CCAAT/enhancer binding protein β (C/EBPβ).

Solute carriers (SLCs), the second largest super-family of membrane proteins in the human genome, transport amino acids, sugars, fatty acids, inorganic ions, essential metals and drugs over membranes. To date no study has provided a comprehensive analysis of SLC localization along the entire GI tract. The aim of the present study was to provide a comprehensive, segment-specific description of the localization of SLC genes along the rat GI tract by employing bioinformatics and molecular biology methods. The Unigene database was screened for rat SLC entries in the intestinal tissue. Using qPCR we measured expression of the annotated genes in the GI tract divided into the following segments: the esophagus, the corpus and the antrum of the stomach, the proximal and distal parts of the duodenum, ileum, jejunum and colon, and the cecum. Our Unigene-derived gene pool was expanded with data from in-house tissue panels and a literature search. We found 44 out of 78 (56%) of gut SLC transcripts to be expressed in all GI tract segments, whereas the majority of remaining SLCs were detected in more than five segments. SLCs are predominantly expressed in gut regions with absorptive functions although expression was also found in segments unrelated to absorption. The proximal jejunum had the highest number of differentially expressed SLCs. In conclusion, SLCs are a crucial molecular component of the GI tract, with many of them expressed along the entire GI tract. This work presents the first overall road map of localization of transporter genes in the GI tract.

Background The drug/metabolite transporter superfamily comprises a diversity of protein domain families with multiple functions including transport of nucleotide sugars. Drug/metabolite transporter domains are contained in both solute carrier families 30, 35 and 39 proteins as well as in acyl-malonyl condensing enzyme proteins. In this paper, we present an evolutionary analysis of nucleotide sugar transporters in relation to the entire superfamily of drug/metabolite transporters that considers crucial intra-protein duplication events that have shaped the transporters. We use a method that combines the strengths of hidden Markov models and maximum likelihood to find relationships between drug/metabolite transporter families, and branches within families. Results We present evidence that the triose-phosphate transporters, domain unknown function 914, uracil-diphosphate glucose-N-acetylglucosamine, and nucleotide sugar transporter families have evolved from a domain duplication event before the radiation of Viridiplantae in the EamA family (previously called domain unknown function 6). We identify previously unknown branches in the solute carrier 30, 35 and 39 protein families that emerged simultaneously as key physiological developments after the radiation of Viridiplantae, including the "35C/E" branch of EamA, which formed in the lineage of T. adhaerens (Animalia). We identify a second cluster of DMTs, called the domain unknown function 1632 cluster, which has non-cytosolic N- and C-termini, and thus appears to have been formed from a different domain duplication event. We identify a previously uncharacterized motif, G-X(6)-G, which is overrepresented in the fifth transmembrane helix of C-terminal domains. We present evidence that the family called fatty acid elongases are homologous to transporters, not enzymes as had previously been thought. Conclusions The nucleotide sugar transporters families were formed through differentiation of the gene cluster EamA (domain unknown function 6) before Viridiplantae, showing for the first time the significance of EamA.

Several families of G protein-coupled receptors (GPCRs) show no significant sequence similarities to each other, and it has been debated which of them share a common origin. We developed and performed integrated and independent HHsearch, Needleman--Wunsch-based and motif analyses on more than 6,600 unique GPCRs from 12 species. Moreover, we mined the evolutionary important Trichoplax adhaerens, Nematostella vectensis, Thalassiosira pseudonana, and Strongylocentrotus purpuratus genomes, revealing remarkably rich vertebrate-like GPCR repertoires already in the early Metazoan species. We found strong evidence that the Adhesion and Frizzled families are children to the cyclic AMP (cAMP) family with HHsearch homology probabilities of 99.8% and 99.4%, respectively, also supported by the Needleman--Wunsch analysis and several motifs. We also found that the large Rhodopsin family is likely a child of the cAMP family with an HHsearch homology probability of 99.4% and conserved motifs. Therefore, we suggest that the Adhesion and Frizzled families originated from the cAMP family in an event close to that which gave rise to the Rhodopsin family. We also found convincing evidence that the Rhodopsin family is parent to the important sensory families; Taste 2 and Vomeronasal type 1 as well as the Nematode chemoreceptor families. The insect odorant, gustatory, and Trehalose receptors, frequently referred to as GPCRs, form a separate cluster without relationship to the other families, and we propose, based on these and others' results, that these families are ligand-gated ion channels rather than GPCRs. Overall, we suggest common descent of at least 97% of the GPCRs sequences found in humans.

Obesity is a complex and a highly individualized disease and the molecular mechanisms behind this disorder need to be better elucidated. Identification of genes and genetic variants that are involved provide opportunities to establish a genetic understanding of the disease. These findings may also provide more rational approaches to therapy, either by identifying underlying causes or point out the need for different treatments. In addition, the timing and severity of obesity may provide insights into the aetiology of obesity and also identify age-specific determinants of weight gain. Recently, genome-wide association studies have led to a rapid progress in our understanding of the genetic basis of various diseases and candidate genes for obesity have been identified. The overall aim of this thesis was to investigate the genetic impact on severity of childhood obesity and the associations between obesity and genetic variants in the fat mass and obesity associated gene, FTO, and MGAT1, the gene encoding mannosyl (α-1,3-)-glycoprotein β-1,2-N-acetyl-glucosaminyltransferase. We show that the impact of parental body mass index (BMI) on the severity of obesity in children is strengthened as the child grows older, whereas the age at obesity onset is of limited importance. By association studies, we show that single nucleotide polymorphisms downstream MGAT1 influence susceptibility to obesity. Moreover, these variants affect the levels of unsaturated fatty acids and desaturase indices, variables previously shown to correlate with obesity. Furthermore, one variant in the first intronic region of FTO is associated with obesity among children but not with BMI or other measures of adiposity at older ages. However, this variant shows a weight-dependent association with cognitive function among elderly men. By direct sequencing, we identified novel variants in FTO, affecting glucose homeostasis in a BMI-independent manner. Furthermore, we found gender specific effects for FTO, both regarding obesity susceptibility and related phenotypes.

Background TMEM18 is a hypothalamic gene that has recently been linked to obesity and BMI in genome wide association studies. However, the functional properties of TMEM18 are obscure. Methods The evolutionary history of TMEM18 was inferred using phylogenetic and bioinformatic methods. The gene's expression profile was investigated with real-time PCR in a panel of rat and mouse tissues and with immunohistochemistry in the mouse brain. Also, gene expression changes were analyzed in three feeding-related mouse models: food deprivation, reward and diet-induced increase in body weight. Finally, we genotyped 502 severely obese and 527 healthy Swedish children for two SNPs near TMEM18 (rs6548238 and rs756131). Results TMEM18 was found to be remarkably conserved and present in species that diverged from the human lineage over 1500 million years ago. The TMEM18 gene was widely expressed and detected in the majority of cells in all major brain regions, but was more abundant in neurons than other cell types. We found no significant changes in the hypothalamic and brainstem expression in the feeding-related mouse models. There was a strong association for two SNPs (rs6548238 and rs756131) of the TMEM18 locus with an increased risk for obesity (p = 0.001 and p = 0.002). Conclusion We conclude that TMEM18 is involved in both adult and childhood obesity. It is one of the most conserved human obesity genes and it is found in the majority of all brain sites, including the hypothalamus and the brain stem, but it is not regulated in these regions in classical energy homeostatic models.

Background Membrane proteins form key nodes in mediating the cell's interaction with the surroundings, which is one of the main reasons why the majority of drug targets are membrane proteins. Results Here we mined the human proteome and identified the membrane proteome subset using three prediction tools for alpha-helices: Phobius, TMHMM, and SOSUI. This dataset was reduced to a non-redundant set by aligning it to the human genome and then clustered with our own interactive implementation of the ISODATA algorithm. The genes were classified and each protein group was manually curated, virtually evaluating each sequence of the clusters, applying systematic comparisons with a range of databases and other resources. We identified 6,718 human membrane proteins and classified the majority of them into 234 families of which 151 belong to the three major functional groups: receptors (63 groups, 1,352 members), transporters (89 groups, 817 members) or enzymes (7 groups, 533 members). Also, 74 miscellaneous groups with 697 members were determined. Interestingly, we find that 41% of the membrane proteins are singlets with no apparent affiliation or identity to any human protein family. Our results identify major differences between the human membrane proteome and the ones in unicellular organisms and we also show a strong bias towards certain membrane topologies for different functional classes: 77% of all transporters have more than six helices while 60% of proteins with an enzymatic function and 88% receptors, that are not GPCRs, have only one single membrane spanning α-helix. Further, we have identified and characterized new gene families and novel members of existing families. Conclusion Here we present the most detailed roadmap of gene numbers and families to our knowledge, which is an important step towards an overall classification of the entire human proteome. We estimate that 27% of the total human proteome are alpha-helical transmembrane proteins and provide an extended classification together with in-depth investigations of the membrane proteome's functional, structural, and evolutionary features.

We studied the genomic positions of 38,129 putative ncRNAs from the RIKEN dataset in relation to protein-coding genes. We found that the dataset has 41% sense, 6% antisense, 24% intronic and 29% intergenic transcripts. Interestingly, 17,678 (47%) of the FANTOM3 transcripts were found to potentially be internally primed from longer transcripts. The highest fraction of these transcripts was found among the intronic transcripts and as many as 77% or 6929 intronic transcripts were both internally primed and unspliced. We defined a filtered subset of 8535 transcripts that did not overlap with protein-coding genes, did not contain ORFs longer than 100 residues and were not internally primed. This dataset contains 53% of the FANTOM3 transcripts associated to known ncRNA in RNAdb and expands previous similar efforts with 6523 novel transcripts. This bioinformatic filtering of the FANTOM3 non-coding dataset has generated a lead dataset of transcripts without signs of being artefacts, providing a suitable dataset for investigation with hybridization-based techniques.

Coat color in Holstein dairy cattle is primarily controlled by the melanocortin 1 receptor (MC1R) gene, a central determinant of black (eumelanin) vs. red/brown pheomelanin synthesis across animal species. The major MC1R alleles in Holsteins are Dominant Black (MC1RD) and Recessive Red (MC1Re). A novel form of dominant red coat color was first observed in an animal born in 1980. The mutation underlying this phenotype was named Dominant Red and is epistatic to the constitutively activated MC1RD. Here we show that a missense mutation in the coatomer protein complex, subunit alpha (COPA), a gene with previously no known role in pigmentation synthesis, is completely associated with Dominant Red in Holstein dairy cattle. The mutation results in an arginine to cysteine substitution at an amino acid residue completely conserved across eukaryotes. Despite this high level of conservation we show that both heterozygotes and homozygotes are healthy and viable. Analysis of hair pigment composition shows that the Dominant Red phenotype is similar to the MC1R Recessive Red phenotype, although less effective at reducing eumelanin synthesis. RNA-seq data similarly show that Dominant Red animals achieve predominantly pheomelanin synthesis by downregulating genes normally required for eumelanin synthesis. COPA is a component of the coat protein I seven subunit complex that is involved with retrograde and cis-Golgi intracellular coated vesicle transport of both protein and RNA cargo. This suggests that Dominant Red may be caused by aberrant MC1R protein or mRNA trafficking within the highly compartmentalized melanocyte, mimicking the effect of the Recessive Red loss of function MC1R allele. Published version

Nematode chemosensory GPCRs in Caenorhabditis elegans (NemChRs) are classified into 19 gene families, and are initially thought to have split from the ancestral Rhodopsin family of GPCRs. However, earlier studies have shown that among all 19 NemChR gene families, only the srw family has a clear sequence relationship to the ancestral Rhodopsin GPCR family. Yet, the phylogenetic relationships between the srw family of NemChRs and the Rhodopsin subfamilies are not fully understood. Also, a widespread search was not previously performed to check for the presence of putative srw family-like sequences or the other 18 NemChR families in several new protostome species outside the nematode lineage. In this study, we have investigated for the presence of 19 NemChR families across 26 eukaryotic species, covering basal eukaryotic branches and provide the first evidence that the srw family of NemChRs is indeed present across several phyla of protostomes. We could identify 29 putative orthologs of the srw family in insects (15 genes), molluscs (11 genes) and Schistosoma mansoni (3 genes). Furthermore, using HMM-HMM profile based comparisons and phylogenetic analysis we show that among all Rhodopsin subfamilies, the peptide and SOG (somatostatin/opioid/galanin) subfamilies are phylogenetically the closest relatives to the srw family of NemChRs. Taken together, we demonstrate that the srw family split from the large Rhodopsin family, possibly from the peptide and/or SOG subfamilies, well before the split of the nematode lineage, somewhere close to the divergence of the common ancestor of protostomes. Our analysis also suggests that the srsx family of NemChRs shares a clear sequence homology with the Rhodopsin subfamilies, as well as with few of the vertebrate olfactory receptors. Overall, this study provides further insights into the evolutionary events that shaped the GPCR chemosensory system in protostome species.

Humans carrying the prevalent rs9939609 A allele of the fat mass and obesity-associated (FTO) gene are more susceptible to developing obesity than noncarries. Recently, polymorphisms in the FTO gene of elderly subjects have also been linked to a reduced volume in the frontal lobe as well as increased risk for incident Alzheimer disease. However, so far there is no evidence directly linking the FTO gene to functional cognitive processes. Here we examined whether the FTO rs9939609 A allele is associated with verbal fluency performance in 355 elderly men at the age of 82 years who have no clinically apparent cognitive impairment. Retrieval of verbal memory is a good surrogate measure reflecting frontal lobe functioning. Here we found that obese and overweight but not normal weight FTO A allele carriers showed a lower performance on verbal fluency than non-carriers (homozygous for rs9939609 T allele). This effect was not observed for a measure of general cognitive performance (i.e., Mini-Mental State Examination score), thereby indicating that the FTO gene primarily affects frontal lobe-dependent cognitive processes in elderly men.

Variants of the FTO gene show strong association with obesity, but the mechanisms behind this association remain unclear. We determined the genome wide DNA methylation profile in blood from 47 female preadolescents. We identified sites associated with the genes KARS, TERF2IP, DEXI, MSI1, STON1 and BCAS3 that had a significant differential methylation level in the carriers of the FTO risk allele (rs9939609). In addition, we identified 20 differentially methylated sites associated with obesity. Our findings suggest that the effect of the FTO obesity risk allele may be mediated through epigenetic changes. Further, these sites might prove to be valuable biomarkers for the understanding of obesity and its comorbidites.

Background Membrane proteins are key components in a large spectrum of diverse functions and thus account for the major proportion of the drug-targeted portion of the genome. From a structural perspective, the α-helical transmembrane proteins can be categorized into major groups based on the number of transmembrane helices and these groups are often associated with specific functions. When compared to the well-characterized seven-transmembrane containing proteins (7TM), other TM groups are less explored and in particular the 4TM group. In this study, we identify the complete 4TM complement from the latest release of the human genome and assess the 4TM structure group as a whole. We functionally characterize this dataset and evaluate the resulting groups and ubiquitous functions, and furthermore describe disease and drug target involvement. Results We classified 373 proteins, which represents ~7 % of the human membrane proteome, and includes 69 more proteins than our previous estimate. We have characterized the 4TM dataset based on functional, structural, and/or evolutionary similarities. Proteins that are involved in transport activity constitute 37 % of the dataset, 23 % are receptor-related, and 13 % have enzymatic functions. Intriguingly, proteins involved in transport are more than double the 15 % of transporters in the entire human membrane proteome, which might suggest that the 4TM topological architecture is more favored for transporting molecules over other functions. Moreover, we found an interesting exception to the ubiquitous intracellular N- and C-termini localization that is found throughout the entire membrane proteome and 4TM dataset in the neurotransmitter gated ion channel families. Overall, we estimate that 58 % of the dataset has a known association to disease conditions with 19 % of the genes possibly involved in different types of cancer. Conclusions We provide here the most robust and updated classification of the 4TM complement of the human genome as a platform to further understand the characteristics of 4TM functions and to explore pharmacological opportunities. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2592-7) contains supplementary material, which is available to authorized users.

Background G protein-coupled receptors (GPCRs) play a central role in eukaryotic signal transduction. However, the GPCR component of this signalling system, at the early origins of metazoans is not fully understood. Here we aim to identify and classify GPCRs in Amphimedon queenslandica (sponge), a member of an earliest diverging metazoan lineage (Porifera). Furthermore, phylogenetic comparisons of sponge GPCRs with eumetazoan and bilaterian GPCRs will be essential to our understanding of the GPCR system at the roots of metazoan evolution. Results We present a curated list of 220 GPCRs in the sponge genome after excluding incomplete sequences and false positives from our initial dataset of 282 predicted GPCR sequences obtained using Pfam search. Phylogenetic analysis reveals that the sponge genome contains members belonging to four of the five major GRAFS families including Glutamate (33), Rhodopsin (126), Adhesion (40) and Frizzled (3). Interestingly, the sponge Rhodopsin family sequences lack orthologous relationships with those found in eumetazoan and bilaterian lineages, since they clustered separately to form sponge specific groups in the phylogenetic analysis. This suggests that sponge Rhodopsins diverged considerably from that found in other basal metazoans. A few sponge Adhesions clustered basal to Adhesion subfamilies commonly found in most vertebrates, suggesting some Adhesion subfamilies may have diverged prior to the emergence of Bilateria. Furthermore, at least eight of the sponge Adhesion members have a hormone binding motif (HRM domain) in their N-termini, although hormones have yet to be identified in sponges. We also phylogenetically clarified that sponge has homologs of metabotropic glutamate (mGluRs) and GABA receptors. Conclusion Our phylogenetic comparisons of sponge GPCRs with other metazoan genomes suggest that sponge contains a significantly diversified set of GPCRs. This is evident at the family/subfamily level comparisons for most GPCR families, in particular for the Rhodopsin family of GPCRs. In summary, this study provides a framework to perform future experimental and comparative studies to further verify and understand the roles of GPCRs that predates the divergence of bilaterian and eumetazoan lineages. Electronic supplementary material The online version of this article (doi:10.1186/s12862-014-0270-4) contains supplementary material, which is available to authorized users.