A population genomics overview of the neuronal nitric oxide synthase (nNOS) gene and its relationship to migraine susceptibility

The ubiquitous chemical messenger molecule nitric oxide (NO) has been implicated in a diverse range of biological activities including neurotransmission, smooth muscle motility and mediation of nociception. Endogenous synthesis of NO by the neuronal isoform of the nitric oxide synthase gene family has an essential role within the central and peripheral nervous systems in addition to the autonomic innervation of cerebral blood vessels. To investigate the potential role of NO and more specifically the neuronal nitric oxide synthase (nNOS) gene in migraine susceptibility, we investigated two microsatellite repeat variants residing within the 5′ and 3′ regions of the nNOS gene. Population genomic evaluation of the two nNOS repeat variants indicated significant linkage disequilibrium between the two loci. Z-DNA conformational sequence structures within the 5′ region of the nNOS gene have the potential to enhance or repress gene promoter activity. We suggest that genetic analysis of this 5′ repeat variant is the more functional variant expressing gene wide information that could affect endogenous NO synthesis and potentially result in diseased states. However, no association with migraine (with or without aura) was seen in our extensive case-control cohort (n = 579 affected with matched controls), when both the 5′ and 3′ genetic variants were investigated.

Comparative kinomics of human and chimpanzee reveal unique kinship and functional diversity generated by new domain combinations

Background: Phosphorylation by protein kinases is a common event in many cellular processes. Further, many kinases perform specialized roles and are regulated by non-kinase domains tethered to kinase domain. Perturbation in the regulation of kinases leads to malignancy. We have identified and analysed putative protein kinases encoded in the genome of chimpanzee which is a close evolutionary relative of human. Result: The shared core biology between chimpanzee and human is characterized by many orthologous protein kinases which are involved in conserved pathways. Domain architectures specific to chimp/human kinases have been observed. Chimp kinases with unique domain architectures are characterized by deletion of one or more non-kinase domains in the human kinases. Interestingly, counterparts of some of the multi-domain human kinases in chimp are characterized by identical domain architectures but with kinase-like non-kinase domain. Remarkably, out of 587 chimpanzee kinases no human orthologue with greater than 95% sequence identity could be identified for 160 kinases. Variations in chimpanzee kinases compared to human kinases are brought about also by differences in functions of domains tethered to the catalytic kinase domain. For example, the heterodimer forming PB1 domain related to the fold of ubiquitin/Ras-binding domain is seen uniquely tethered to PKC-like chimpanzee kinase. Conclusion: Though the chimpanzee and human are evolutionary very close, there are chimpanzee kinases with no close counterpart in the human suggesting differences in their functions. This analysis provides a direction for experimental analysis of human and chimpanzee protein kinases in order to enhance our understanding on their specific biological roles.

Identification of putative regulatory motifs in the upstream regions of co-expressed functional groups of genes in Plasmodium falciparum

Background: Regulation of gene expression in Plasmodium falciparum (Pf) remains poorly understood. While over half the genes are estimated to be regulated at the transcriptional level, few regulatory motifs and transcription regulators have been found. Results: The study seeks to identify putative regulatory motifs in the upstream regions of 13 functional groups of genes expressed in the intraerythrocytic developmental cycle of Pf. Three motif-discovery programs were used for the purpose, and motifs were searched for only on the gene coding strand. Four motifs – the 'G-rich', the 'C-rich', the 'TGTG' and the 'CACA' motifs – were identified, and zero to all four of these occur in the 13 sets of upstream regions. The 'CACA motif' was absent in functional groups expressed during the ring to early trophozoite transition. For functional groups expressed in each transition, the motifs tended to be similar. Upstream motifs in some functional groups showed 'positional conservation' by occurring at similar positions relative to the translational start site (TLS); this increases their significance as regulatory motifs. In the ribonucleotide synthesis, mitochondrial, proteasome and organellar translation machinery genes, G-rich, C-rich, CACA and TGTG motifs, respectively, occur with striking positional conservation. In the organellar translation machinery group, G-rich motifs occur close to the TLS. The same motifs were sometimes identified for multiple functional groups; differences in location and abundance of the motifs appear to ensure different modes of action. Conclusion: The identification of positionally conserved over-represented upstream motifs throws light on putative regulatory elements for transcription in Pf.

Comparative and functional genome analysis of fungi for development of the protein production host Trichoderma reesei

Filamentous fungi of the subphylum Pezizomycotina are well known as protein and secondary metabolite producers. Various industries take advantage of these capabilities. However, the molecular biology of yeasts, i.e. Saccharomycotina and especially that of Saccharomyces cerevisiae, the baker's yeast, is much better known. In an effort to explain fungal phenotypes through their genotypes we have compared protein coding gene contents of Pezizomycotina and Saccharomycotina. Only biomass degradation and secondary metabolism related protein families seem to have expanded recently in Pezizomycotina. Of the protein families clearly diverged between Pezizomycotina and Saccharomycotina, those related to mitochondrial functions emerge as the most prominent. However, the primary metabolism as described in S. cerevisiae is largely conserved in all fungi. Apart from the known secondary metabolism, Pezizomycotina have pathways that could link secondary metabolism to primary metabolism and a wealth of undescribed enzymes. Previous studies of individual Pezizomycotina genomes have shown that regardless of the difference in production efficiency and diversity of secreted proteins, the content of the known secretion machinery genes in Pezizomycotina and Saccharomycotina appears very similar. Genome wide analysis of gene products is therefore needed to better understand the efficient secretion of Pezizomycotina. We have developed methods applicable to transcriptome analysis of non-sequenced organisms. TRAC (Transcriptional profiling with the aid of affinity capture) has been previously developed at VTT for fast, focused transcription analysis. We introduce a version of TRAC that allows more powerful signal amplification and multiplexing. We also present computational optimisations of transcriptome analysis of non-sequenced organism and TRAC analysis in general. Trichoderma reesei is one of the most commonly used Pezizomycotina in the protein production industry. In order to understand its secretion system better and find clues for improvement of its industrial performance, we have analysed its transcriptomic response to protein secretion stress conditions. In comparison to S. cerevisiae, the response of T. reesei appears different, but still impacts on the same cellular functions. We also discovered in T. reesei interesting similarities to mammalian protein secretion stress response. Together these findings highlight targets for more detailed studies.

Structure-Based Phylogeny as a Diagnostic for Functional Characterization of Proteins with a Cupin Fold

Background: The members of cupin superfamily exhibit large variations in their sequences, functions, organization of domains, quaternary associations and the nature of bound metal ion, despite having a conserved beta-barrel structural scaffold. Here, an attempt has been made to understand structure-function relationships among the members of this diverse superfamily and identify the principles governing functional diversity. The cupin superfamily also contains proteins for which the structures are available through world-wide structural genomics initiatives but characterized as ``hypothetical''. We have explored the feasibility of obtaining clues to functions of such proteins by means of comparative analysis with cupins of known structure and function. Methodology/Principal Findings: A 3-D structure-based phylogenetic approach was undertaken. Interestingly, a dendrogram generated solely on the basis of structural dissimilarity measure at the level of domain folds was found to cluster functionally similar members. This clustering also reflects an independent evolution of the two domains in bicupins. Close examination of structural superposition of members across various functional clusters reveals structural variations in regions that not only form the active site pocket but are also involved in interaction with another domain in the same polypeptide or in the oligomer. Conclusions/Significance: Structure-based phylogeny of cupins can influence identification of functions of proteins of yet unknown function with cupin fold. This approach can be extended to other proteins with a common fold that show high evolutionary divergence. This approach is expected to have an influence on the function annotation in structural genomics initiatives.

Functional diversity of human protein kinase splice variants marks significant expansion of human kinome

Background: Protein kinases are involved in diverse spectrum of cellular processes. Availability of draft version of the human genomic data in the year 2001 enabled recognition of repertoire of protein kinases. However, over the years the human genomic data is being refined and the current release of human genomic data has helped us to recognize a larger repertoire of over 900 human protein kinases represented mainly by splice variants. Results: Many of these identified protein kinases are alternatively spliced products. Interestingly, some of the human kinase splice variants appear to be significantly diverged in terms of their functional properties as represented by incorporation or absence of one or more domains. Many sets of protein kinase splice variants have substantially different domain organization and in a few sets of splice variants kinase domains belong to different subfamilies of kinases suggesting potential participation in different signal transduction pathways. Conclusions: Addition or deletion of a domain between splice variants of multi-domain kinases appears to be a means of generating differences in the functional features of otherwise similar kinases. It is intriguing that marked sequence diversity within the catalytic regions of some of the splice variant kinases result in kinases belonging to different subfamilies. These human kinase splice variants with different functions might contribute to diversity of eukaryotic cellular signaling.

Peptidomics and genomics analysis of novel antimicrobial peptides from the frog, Rana nigrovittata

Much attention has been paid on amphibian peptides for their wide-ranging pharmacological properties, clinical potential, and gene-encoded origin. More than 300 antimicrobial peptides (AMPs) from amphibians have been studied. Peptidomics and genomics analysis combined with functional test including microorganism killing, histamine-releasing, and mast cell degranulation was used to investigate antimicrobial peptide diversity. Thirty-four novel AMPs from skin secretions of Rana nigrovittata were identified in current work, and they belong to 9 families, including 6 novel families. Other three families are classified into rugosin, gaegurin, and temporin family of amphibian AMP, respectively. These AMPs share highly conserved preproregions including signal peptides and spacer acidic peptides, while greatly diversified on mature peptides structures. In this work, peptidomics combined with genomics analysis was confirmed to be an effective way to identify amphibian AMPs, especially novel families. Some AMPs reported here will provide leading molecules for designing novel antimicrobial agents. (C) 2009 Elsevier Inc. All rights reserved

Comparative genomics of Bifidobacteria

Chapter 2 of this thesis describes the sequence analysis of 14 bifidobacterial genomes from various species of the genus Bifidobacterium, and the determination of their open pan-genome trend. This analysis first determined the total number of genes to be considered as the reservoir of functions available to representatives of this genus. Many identified genes are still uncharacterized, but may be involved in the adaptation to the gut environment. This comparative genomic analysis also determined a pool of ortholog functions used to infer their phylogenetic relationship, thereby providing a more robust approach compared to that based solely on the16S rRNA-encoding gene. The genome sequence of an isolate from the insect hindgut Bifidobacterium asteroides PRL2011 was fully characterized in Chapter 3, surprisingly revealing a putative respiratory metabolism, which was also found to be present in other insect isolates, suggesting that respiration was an ancient feature of this genus, but also an adaptative trait to different atmosferic oxygen levels. Chapter 4 of this thesis outlines a comparative study which focused on the analysis of representatives of the Bifidobacterium breve species, revealing that the genetic variability among members of this species principally consists of genes with a role in adaptation to host environment and gut colonization. Finally, Chapter 5 describes the analysis of the genome sequence of Bifidobacterium animalis subsp. lactis BLC1, a probiotic bacterium widely used in food industries as an ingredient of functional foods, providing information that will allow future investigations of this species.

Functional annotation signatures of disease susceptibility loci improve SNP association analysis.

BACKGROUND: Genetic association studies are conducted to discover genetic loci that contribute to an inherited trait, identify the variants behind these associations and ascertain their functional role in determining the phenotype. To date, functional annotations of the genetic variants have rarely played more than an indirect role in assessing evidence for association. Here, we demonstrate how these data can be systematically integrated into an association study's analysis plan. RESULTS: We developed a Bayesian statistical model for the prior probability of phenotype-genotype association that incorporates data from past association studies and publicly available functional annotation data regarding the susceptibility variants under study. The model takes the form of a binary regression of association status on a set of annotation variables whose coefficients were estimated through an analysis of associated SNPs in the GWAS Catalog (GC). The functional predictors examined included measures that have been demonstrated to correlate with the association status of SNPs in the GC and some whose utility in this regard is speculative: summaries of the UCSC Human Genome Browser ENCODE super-track data, dbSNP function class, sequence conservation summaries, proximity to genomic variants in the Database of Genomic Variants and known regulatory elements in the Open Regulatory Annotation database, PolyPhen-2 probabilities and RegulomeDB categories. Because we expected that only a fraction of the annotations would contribute to predicting association, we employed a penalized likelihood method to reduce the impact of non-informative predictors and evaluated the model's ability to predict GC SNPs not used to construct the model. We show that the functional data alone are predictive of a SNP's presence in the GC. Further, using data from a genome-wide study of ovarian cancer, we demonstrate that their use as prior data when testing for association is practical at the genome-wide scale and improves power to detect associations. CONCLUSIONS: We show how diverse functional annotations can be efficiently combined to create 'functional signatures' that predict the a priori odds of a variant's association to a trait and how these signatures can be integrated into a standard genome-wide-scale association analysis, resulting in improved power to detect truly associated variants.

AutoFACT: An Automatic Functional Annotation and Classification Tool

Affiliation: Centre Robert-Cedergren de l'Université de Montréal en bio-informatique et génomique & Département de biochimie, Université de Montréal

Automatic prediction of functional site regions in low-resolution protein structures

World-wide structural genomics initiatives are rapidly accumulating structures for which limited functional information is available. Additionally, state-of-the art structural prediction programs are now capable of generating at least low resolution structural models of target proteins. Accurate detection and classification of functional sites within both solved and modelled protein structures therefore represents an important challenge. We present a fully automatic site detection method, FuncSite, that uses neural network classifiers to predict the location and type of functionally important sites in protein structures. The method is designed primarily to require only backbone residue positions without the need for specific side-chain atoms to be present. In order to highlight effective site detection in low resolution structural models FuncSite was used to screen model proteins generated using mGenTHREADER on a set of newly released structures. We found effective metal site detection even for moderate quality protein models illustrating the robustness of the method.

Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs

Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca2+ transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca2+ transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.

Contributions and perspectives of chicken genomics in Brazil: from biological model to export commodity

Chicken is one of the most important sources of animal protein for human consumption, and breeding programmes have been responsible for constant improvements in production efficiency and product quality. Furthermore, chicken has largely contributed to fundamental discoveries in biology for the last 100 years. In this article we review recent developments in poultry genomics and their contribution to adding functional information to the already existing structural genomics, including the availability of the complete genome sequence, a comprehensive collection of mRNA sequences ( ESTs), microarray platforms, and their use to complement QTL mapping strategies in the identification of genes that underlie complex traits. Efforts of the Brazilian Poultry Genomics Programme in this area resulted in generation of a resource population, which was used for identification of Quantitative Trait Loci ( QTL) regions, generation of ESTs and candidate gene studies that contributed to furthering our understanding of the complex biological processes involved in growth and muscular development in chicken.

Comparative genomics and phylogenetics of the vertebrate CYP3 family

CYP3A verstoffwechselt mehr als 50% aller gegenwärtig in der Therapie eingesetzten Wirkstoffe, die häufig an klinisch relevanten Arzneimitttel-Wechselwirkungen beteiligt sind. Das Verständnis über die Bedeutung und die Regulation von einzelnen CYP3A Genen in der Pharmakologie und Physiologie ist unvollständig. Wir untersuchten die Evolution des CYP3 Genlokus über einen Zeitraum von 450 Millionen Jahre mittels genomischer Sequenzen von 16 Tierarten. Neue CYP3 Unterfamilien (CYP3B, C und D) entstanden über eine beschleunigte Evolution aus CYP3A Vorstufen von Clupeocephala Spezies. Ausgeprägte funktionelle Unterschiede traten zwischen CYP3A in Säugern und Clupeocephala CYP3 auf. Alle amnioten CYP3A Gene entwickelten sich aus zwei CYP3A Urgenen. Aufgrund der Entstehung von Säugern mit Plazenta ging eines von ihnen verloren während das andere eine neue genomische Umgebung infolge einer Translokation erlangte. In Primaten unterzog sich CYP3A mit mehreren Genduplikationen, Deletionen, Pseudogenisierung und Genkonversionen einer raschen evolutionären Veränderung. Die Entwicklung von CYP3A in Schmalnasenaffen (Alte Welt Affen, große Menschenaffen und Menschen) unterschieden sich wesentlich von Neue Welt Primaten (z.B. gewöhnlichen Krallenaffen) und Feuchtnasenaffen (z.B. Galago). Stellvertretend für die CYP3A Protein-codierende Sequenz entdeckten wir zwei frühe Episoden von besonders starker positiver Selektion: (1) auf CYP3A7 in der frühen hominoiden Evolution, welche im fetalen Zeitraum von einer Einschränkung der hepatischen Expression begleitet war, und (2) auf humanes CYP3A4 im Anschluss an die Teilung der Abstammungslinie in Schimpansen und Mensch. In Übereinstimmung mit diesen Befunden beeinflussen drei von vier positiv ausgewählten Aminosäuren, die in früheren biochemischen CYP3A Studien untersucht wurden, die Aktivität und Regioselektivität. Es ist somit naheliegend, dass CYP3A7 und CYP3A4 katalytische Funktionen erworben haben können, die besonders wichtig waren für die Evolution von Hominoiden und Menschen. Die Charakterisierung von CYP3A Promotoren in Primaten zeigte eine Anreicherung von ER6 Elementen in CYP3A Promotoren von Primaten und einen Trend in Richtung Erhöhung der ER6 Enstehung entlang den Abstammungslinien, die zu humanen und Schimpansen CYP3A4 führten. Die steigende Anzahl an ER6 Elementen kann durch die ausgeprägte CYP3A4 Induzierbarkeit und Expressionsvariabilität im Menschen verursacht sein.