The Genomic Threading Database: a comprehensive resource for structural annotations of the genomes from key organisms

Currently, the Genomic Threading Database (GTD) contains structural assignments for the proteins encoded within the genomes of nine eukaryotes and 101 prokaryotes. Structural annotations are carried out using a modified version of GenTHREADER, a reliable fold recognition method. The Gen THREADER annotation jobs are distributed across multiple clusters of processors using grid technology and the predictions are deposited in a relational database accessible via a web interface at http://bioinf.cs.ucl.ac.uk/GTD. Using this system, up to 84% of proteins encoded within a genome can be confidently assigned to known folds with 72% of the residues aligned. On average in the GTD, 64% of proteins encoded within a genome are confidently assigned to known folds and 58% of the residues are aligned to structures.

The Genomic Threading Database

The Genomic Threading Database currently contains structural annotations for the genomes of over 100 recently sequenced organisms. Annotations are carried out by using our modified GenTHREADER software and through implementing grid technology.

Improvement of the GenTHREADER method for genomic fold recognition

Motivation: In order to enhance genome annotation, the fully automatic fold recognition method GenTHREADER has been improved and benchmarked. The previous version of GenTHREADER consisted of a simple neural network which was trained to combine sequence alignment score, length information and energy potentials derived from threading into a single score representing the relationship between two proteins, as designated by CATH. The improved version incorporates PSI-BLAST searches, which have been jumpstarted with structural alignment profiles from FSSP, and now also makes use of PSIPRED predicted secondary structure and bi-directional scoring in order to calculate the final alignment score. Pairwise potentials and solvation potentials are calculated from the given sequence alignment which are then used as inputs to a multi-layer, feed-forward neural network, along with the alignment score, alignment length and sequence length. The neural network has also been expanded to accommodate the secondary structure element alignment (SSEA) score as an extra input and it is now trained to learn the FSSP Z-score as a measurement of similarity between two proteins. Results: The improvements made to GenTHREADER increase the number of remote homologues that can be detected with a low error rate, implying higher reliability of score, whilst also increasing the quality of the models produced. We find that up to five times as many true positives can be detected with low error rate per query. Total MaxSub score is doubled at low false positive rates using the improved method.

Traditional Chinese medicine research in the post-genomic era: good practice, priorities, challenges and opportunities

Background and aims: GP-TCM is the 1st EU-funded Coordination Action consortium dedicated to traditional Chinese medicine (TCM) research. This paper aims to summarise the objectives, structure and activities of the consortium and introduces the position of the consortium regarding good practice, priorities, challenges and opportunities in TCM research. Serving as the introductory paper for the GPTCM Journal of Ethnopharmacology special issue, this paper describes the roadmap of this special issue and reports how the main outputs of the ten GP-TCM work packages are integrated, and have led to consortium-wide conclusions. Materials and methods: Literature studies, opinion polls and discussions among consortium members and stakeholders. Results: By January 2012, through 3 years of team building, the GP-TCM consortium had grown into a large collaborative network involving ∼200 scientists from 24 countries and 107 institutions. Consortium members had worked closely to address good practice issues related to various aspects of Chinese herbal medicine (CHM) and acupuncture research, the focus of this Journal of Ethnopharmacology special issue, leading to state-of-the-art reports, guidelines and consensus on the application of omics technologies in TCM research. In addition, through an online survey open to GP-TCM members and non-members, we polled opinions on grand priorities, challenges and opportunities in TCM research. Based on the poll, although consortium members and non-members had diverse opinions on the major challenges in the field, both groups agreed that high-quality efficacy/effectiveness and mechanistic studies are grand priorities and that the TCM legacy in general and its management of chronic diseases in particular represent grand opportunities. Consortium members cast their votes of confidence in omics and systems biology approaches to TCM research and believed that quality and pharmacovigilance of TCM products are not only grand priorities, but also grand challenges. Non-members, however, gave priority to integrative medicine, concerned on the impact of regulation of TCM practitioners and emphasised intersectoral collaborations in funding TCM research, especially clinical trials. Conclusions: The GP-TCM consortium made great efforts to address some fundamental issues in TCM research, including developing guidelines, as well as identifying priorities, challenges and opportunities. These consortium guidelines and consensus will need dissemination, validation and further development through continued interregional, interdisciplinary and intersectoral collaborations. To promote this, a new consortium, known as the GP-TCM Research Association, is being established to succeed the 3-year fixed term FP7 GP-TCM consortium and will be officially launched at the Final GP-TCM Congress in Leiden, the Netherlands, in April 2012.

Comparative genomic indexing reveals the phylogenomics of Escherichia coli pathogens

The Escherichia coli O26 serogroup includes important food-borne pathogens associated with human and animal diarrheal disease. Current typing methods have revealed great genetic heterogeneity within the O26 group; the data are often inconsistent and focus only on verotoxin (VT)-positive O26 isolates. To improve current understanding of diversity within this serogroup, the genomic relatedness of VT-positive and -negative O26 strains was assessed by comparative genomic indexing. Our results clearly demonstrate that irrespective of virulence characteristics and pathotype designation, the O26 strains show greater genomic similarity to each other than to any other strain included in this study. Our data suggest that enteropathogenic and VT-expressing E. coli O26 strains represent the same clonal lineage and that W-expressing E. coli O26 strains have gained additional virulence characteristics. Using this approach, we established the core genes which are central to the E. coli species and identified regions of variation from the E. coli K-12 chromosomal backbone.

Mapping the platelet profile for functional genomic studies and demonstration of the effect size of the GP6 locus

BACKGROUND: Evidence suggests the wide variation in platelet response within the population is genetically controlled. Unraveling the complex relationship between sequence variation and platelet phenotype requires accurate and reproducible measurement of platelet response. OBJECTIVE: To develop a methodology suitable for measuring signaling pathway-specific platelet phenotype, to use this to measure platelet response in a large cohort, and to demonstrate the effect size of sequence variation in a relevant model gene. METHODS: Three established platelet assays were evaluated: mobilization of [Ca(2+)](i), aggregometry and flow cytometry, each in response to adenosine 5'-diphosphate (ADP) or the glycoprotein (GP) VI-specific crosslinked collagen-related peptide (CRP). Flow cytometric measurement of fibrinogen binding and P-selectin expression in response to a single, intermediate dose of each agonist gave the best combination of reproducibility and inter-individual variability and was used to measure the platelet response in 506 healthy volunteers. Pathway specificity was ensured by blocking the main subsidiary signaling pathways. RESULTS: Individuals were identified who were hypo- or hyper-responders for both pathways, or who had differential responses to the two agonists, or between outcomes. 89 individuals, retested three months later using the same methodology, showed high concordance between the two visits in all four assays (r(2) = 0.872, 0.868, 0.766 and 0.549); all subjects retaining their phenotype at recall. The effect of sequence variation at the GP6 locus accounted for approximately 35% of the variation in the CRP-XL response. CONCLUSION: Genotyping-phenotype association studies in a well-characterized, large cohort provides a powerful strategy to measure the effect of sequence variation in genes regulating the platelet response.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Genomic signals of migration and continuity in Britain before the Anglo-Saxons

The purported migrations that have formed the peoples of Britain have been the focus of generations of scholarly controversy. However, this has not benefited from direct analyses of ancient genomes. Here we report nine ancient genomes (~1 x) of individuals from northern Britain: seven from a Roman era York cemetery, bookended by earlier Iron-Age and later Anglo-Saxon burials. Six of the Roman genomes show affinity with modern British Celtic populations, particularly Welsh, but significantly diverge from populations from Yorkshire and other eastern English samples. They also show similarity with the earlier Iron-Age genome, suggesting population continuity, but differ from the later Anglo-Saxon genome. This pattern concords with profound impact of migrations in the Anglo-Saxon period. Strikingly, one Roman skeleton shows a clear signal of exogenous origin, with affinities pointing towards the Middle East, confirming the cosmopolitan character of the Empire, even at its northernmost fringes.

Genomic evidence for the Pleistocene and recent population history of Native Americans

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

Early membrane estrogenic effects required for full expression of slower genomic actions in a nerve cell line.

Interpretations of steroid hormone actions as slow, nuclear, transcriptional events have frequently been seen as competing against inferences of rapid membrane actions. We have discovered conditions where membrane-limited effects potentiate later transcriptional actions in a nerve cell line. Making use of a two-pulse hormonal schedule in a transfection system, early and brief administration of conjugated, membrane-limited estradiol was necessary but not sufficient for full transcriptional potency of the second estrogen pulse. Efficacy of the first pulse depended on intact signal transduction pathways. Surprisingly, the actions of both pulses were blocked by a classical estrogen receptor (ER) antagonist. Thus, two different modes of steroid hormone action can synergize.

Integration of steroid hormone initiated membrane action to genomic function in the brain

Estrogen is a ligand for the estrogen receptor (ER), which on binding 17beta-estradiol, functions as a ligand-activated transcription factor and regulates the transcription of target genes. This is the slow genomic mode of action. However, rapid non-genomic actions of estrogen also exist at the cell membrane. Using a novel two-pulse paradigm in which the first pulse rapidly initiates non-genomic actions using a membrane-limited estrogen conjugate (E-BSA), while the second pulse promotes genomic transcription from a consensus estrogen response element (ERE), we have demonstrated that rapid actions of estrogen potentiate the slower transcriptional response from an ERE-reporter in neuroblastoma cells. Since rapid actions of estrogen activate kinases, we used selective inhibitors in the two-pulse paradigm to determine the intracellular signaling cascades important in such potentiation. Inhibition of protein kinase A (PKA), PKC, mitogen activated protein kinase (MAPK) or phosphatidylinositol 3-OH kinase (PI-3K) in the first pulse decreases potentiation of transcription. Also, our data with both dominant negative and constitutive mutants of Galpha subunits show that Galpha(q) initiates the rapid signaling cascade at the membrane in SK-N-BE(2)C neuroblastoma cells. We discuss two models of multiple kinase activation at the membrane Pulses of estrogen induce lordosis behavior in female rats. Infusion of E-BSA into the ventromedial hypothalamus followed by 17beta-estradiol in the second pulse could induce lordosis behavior, demonstrating the applicability of this paradigm in vivo. A model where non-genomic actions of estrogen couple to genomic actions unites both aspects of hormone action.

Calcium Flux in Neuroblastoma Cells Is a Coupling Mechanism between Non-Genomic and Genomic Modes of Estrogens

Estrogens have been demonstrated to rapidly modulate calcium levels in a variety of cell types. However, the significance of estrogen-mediated calcium flux in neuronal cells is largely unknown. The relative importance of intra- and extracellular sources of calcium in estrogenic effects on neurons is also not well understood. Previously, we have demonstrated that membrane-limited estrogens, such as E-BSA given before an administration of a 2-hour pulse of 17beta-estradiol (E(2)), can potentiate the transcription mediated by E(2) from a consensus estrogen response element (ERE)-driven reporter gene. Inhibitors to signal transduction cascades given along with E-BSA or E(2) demonstrated that calcium flux is important for E-BSA-mediated potentiation of transcription in a transiently transfected neuroblastoma cell line. In this report, we have used inhibitors to different voltage-gated calcium channels (VGCCs) and to intracellular store receptors along with E-BSA in the first pulse or with E(2) in the second pulse to investigate the relative importance of these channels to estrogen-mediated transcription. Neither L- nor P-type VGCCs seem to play a role in estrogen action in these cells; while N-type VGCCs are important in both the non-genomic and genomic modes of estrogen action. Specific inhibitors also showed that the ryanodine receptor and the inositol trisphosphate receptor are important to E-BSA-mediated transcriptional potentiation. This report provides evidence that while intracellular stores of calcium are required to couple non-genomic actions of estrogen initiated at the membrane to transcription in the nucleus, extracellular sources of calcium are also important in both non-genomic and genomic actions of estrogens. Copyright (c) 2005 S. Karger AG, Basel.

Estrogens and thyroid hormone: Non-genomic effects are coupled to transcription.

Estrogens and thyroid hormones are regulators of important diverse physiological processes such as reproduction, thermogenesis, neural development, neural differentiation and cardiovascular functions. Both are ligands for receptors in the nuclear receptor superfamily, which act as ligand-dependent transcription factors, regulating transcription. However, estrogens and thyroid hormones also rapidly (within minutes or seconds) activate kinase cascades and calcium increases, presumably initiated at the cell membrane. We discuss the relevance of both modes of hormone action, including the membrane estrogen receptor, to physiology, with particular reference to lordosis behavior. We first showed that estrogen restricted to the membrane can, in fact, lead to subsequent increases in transcription from a consensus estrogen response element-based reporter in the neuroblastoma cell line, SK-N-BE(2)C. Using a novel hormonal paradigm, we also showed that the activation of protein kinase A, protein kinase C, mitogen activated protein kinase and increases in calcium were important in the ability of the membrane-limited estrogen to potentiate transcription. We discuss the source of calcium important in transcriptional potentiation. Since estrogens and thyroid hormones have common effects on neuroprotection, cognition and mood, we also hypothesized that crosstalk could occur between the rapid actions of thyroid hormones and the genomic actions of estrogens. In neural cells, we showed that triiodothyronine acting rapidly via MAPK can increase transcription by the nuclear estrogen receptor ERa from a consensus estrogen response element, possibly by the phosphorylation of the ERa. Novel mechanisms that link signals initiated by hormones from the membrane to the nucleus are physiologically relevant and can achieve neuroendocrine integration