A blockchain-based framework to support pharmacogenetic data sharing.

The successful implementation of pharmacogenetics (PGx) into clinical practice requires patient genomic data to be shared between stakeholders in multiple settings. This creates a number of barriers to widespread adoption of PGx, including privacy concerns related to the storage and movement of identifiable genomic data. Informatic solutions that support secure and equitable data access for genomic data are therefore important to PGx. Here we propose a methodology that uses smart contracts implemented on a blockchain-based framework, PGxChain, to address this issue. The design requirements for PGxChain were identified through a systematic literature review, identifying technical challenges and barriers impeding the clinical implementation of pharmacogenomics. These requirements included security and privacy, accessibility, interoperability, traceability and legal compliance. A proof-of-concept implementation based on Ethereum was then developed that met the design requirements. PGxChain's performance was examined using Hyperledger Caliper for latency, throughput, and transaction success rate. The findings clearly indicate that blockchain technology offers considerable potential to advance pharmacogenetic data sharing, particularly with regard to PGx data security and privacy, large-scale accessibility of PGx data, PGx data interoperability between multiple health care providers and compliance with data-sharing laws and regulations.

Long-term follow-up of successful hepatitis C virus therapy: waning immune responses and disappearance of liver disease are consistent with cure.

BACKGROUND: A sustained viral response (SVR) after interferon-based therapy of chronic hepatitis C virus (HCV) infection is regarded to represent a cure. Previous studies have used different markers to clarify whether an SVR truly represents a cure, but no study has combined a clinical work-up with highly sensitive HCV RNA detection, and the determination of immune responses. AIM: To determine clinical, histological, virological and immunological markers 5-20 years after SVR. METHODS: In 54 patients, liver biochemistry, histology and elastography were evaluated. Liver biopsies, plasma and peripheral blood mononuclear cells (PBMCs) were tested for minute amounts of HCV RNA. HCV-specific T-cell responses were monitored by ELISpot and pentamer staining, and humoral responses by measuring HCV nonstructural (NS)3-specific antibodies and virus neutralisation. RESULTS: Liver disease regressed significantly in all patients, and 51 were HCV RNA-negative in all tissues tested. There was an inverse association between liver disease, HCV-specific T-cell responses and HCV antibody levels with time from SVR, supporting that the virus had been cleared. The three patients, who all lacked signs of liver disease, had HCV RNA in PBMCs 5-9 years after SVR. All three had HCV-specific T cells and NS3 antibodies, but no cross-neutralising antibodies. CONCLUSIONS: Our combined data confirm that a SVR corresponds to a long-term clinical cure. The waning immune responses support the disappearance of the antigenic stimulus. Transient HCV RNA traces may be detected in some patients up to 9 years after SVR, but no marker associates this with an increased risk for liver disease.

Vulnerability versus resilience to prenatal stress in male and female rats; implications from gene expression profiles in the hippocampus and frontal cortex.

Adverse life events during pregnancy may impact upon the developing fetus, predisposing prenatally stressed offspring to the development of psychopathology. In the present study, we examined the effects of prenatal restraint stress (PS) on anxiety- and depression-related behavior in both male and female adult Sprague-Dawley rats. In addition, gene expression profiles within the hippocampus and frontal cortex (FC) were examined in order to gain more insight into the molecular mechanisms that mediate the behavioral effects of PS exposure. PS significantly increased anxiety-related behavior in male, but not female offspring. Likewise, depression-related behavior was increased in male PS rats only. Further, male PS offspring showed increased basal plasma corticosterone levels in adulthood, whereas both PS males and females had lower stress-induced corticosterone levels when compared to controls. Microarray-based profiling of the hippocampus and FC showed distinct sex-dependent changes in gene expression after PS. Biological processes and/or signal transduction cascades affected by PS included glutamatergic and GABAergic neurotransmission, mitogen-activated protein kinase (MAPK) signaling, neurotrophic factor signaling, phosphodiesterase (PDE)/ cyclic nucleotide signaling, glycogen synthase kinase 3 (GSK3) signaling, and insulin signaling. Further, the data indicated that epigenetic regulation is affected differentially in male and female PS offspring. These sex-specific alterations may, at least in part, explain the behavioral differences observed between both sexes, i.e. relative vulnerability versus resilience to PS in male versus female rats, respectively. These data reveal novel potential targets for antidepressant and mood stabilizing drug treatments including PDE inhibitors and histone deacetylase (HDAC) inhibitors.

Pleiotropic effects of negative energy balance in the postpartum dairy cow on splenic gene expression: repercussions for innate and adaptive immunity.

Increased energy demands to support lactation, coupled with lowered feed intake capacity results in negative energy balance (NEB) and is typically characterized by extensive mobilization of body energy reserves in the early postpartum dairy cow. The catabolism of stored lipid leads to an increase in the systemic concentrations of nonesterified fatty acids (NEFA) and beta-hydroxy butyrate (BHB). Oxidation of NEFA in the liver result in the increased production of reactive oxygen species and the onset of oxidative stress and can lead to disruption of normal metabolism and physiology. The immune system is depressed in the peripartum period and early lactation and dairy cows are therefore more vulnerable to bacterial infections causing mastitis and or endometritis at this time. A bovine Affymetrix oligonucleotide array was used to determine global gene expression in the spleen of dairy cows in the early postpartum period. Spleen tissue was removed post mortem from five severe NEB (SNEB) and five medium NEB (MNEB) cows 15 days postpartum. SNEB increased systemic concentrations of NEFA and BHB, and white blood cell and lymphocyte numbers were decreased in SNEB animals. A total of 545 genes were altered by SNEB. Network analysis using Ingenuity Pathway Analysis revealed that SNEB was associated with NRF2-mediated oxidative stress, mitochondrial dysfunction, endoplasmic reticulum stress, natural killer cell signaling, p53 signaling, downregulation of IL-15, BCL-2, and IFN-gamma; upregulation of BAX and CHOP and increased apoptosis with a potential negative impact on innate and adaptive immunity.

Genomics approaches to study the biology underlying resistance to trypanosomiasis--some unexpected lessons.

An international multidisciplinary consortium is conducting a programme of research on the host response to trypanosome infection. This builds upon quantitative trait loci (QTL) mapping which identified genome regions influencing susceptibility to pathology following T. congolense infection in both cattle and mice. The approach uses expression analysis to examine the response of both susceptible and resistant strains and a series of novel informatics tools to identify pathways which are activated as a result of challenge, and those which are differentially used by resistant and susceptible strains. Of particular interest are those pathways which simultaneously satisfy both criteria, i.e. are significantly differentially activated and contain genes within QTL regions. However, it is important to stress that it is not required that the genes within the QTL region are differentially expressed themselves.

A systematic, data-driven approach to the combined analysis of microarray and QTL data.

High-throughputtechnologies inevitably produce vast quantities of data. This presents challenges in terms of developing effective analysis methods, particularly where the analysis involves combining data derived from different experimental technologies. In this investigation, a systematic approach was applied to combine microarray gene expression data, quantitative trait loci (QTL) data and pathway analysis resources in order to identify functional candidate genes underlying tolerance to Trypanosoma congolense infection in cattle. We automated much of the analysis using Taverna workflows previously developed for the study of trypanotolerance in the mouse model. Pathways represented by genes within the QTL regions were identified, and this list was subsequently ranked according to which pathways were over-represented in the set of genes that were differentially expressed (over time or between tolerant N'dama and susceptible Boran breeds) at various timepoints after T. congolense infection. The genes within the QTLthat played a role in the highest ranked pathways were flagged as good targets for further investigation and experimental confirmation.

Gene expression profiling in a mouse model for African trypanosomiasis.

This study aimed to provide the foundation for an integrative approach to the identification of the mechanisms underlying the response to infection with Trypanosoma congolense, and to identify pathways that have previously been overlooked. We undertook a large-scale gene expression analysis study comparing susceptible A/J and more tolerant C57BL/6 mice. In an initial time course experiment, we monitored the development of parasitaemia and anaemia in every individual. Based on the kinetics of disease progression, we extracted total RNA from liver at days 0, 4, 7, 10 and 17 post infection and performed a microarray analysis. We identified 64 genes that were differentially expressed in the two strains in non-infected animals, of which nine genes remained largely unaffected by the disease. Gene expression profiling at stages of low, peak, clearance and recurrence of parasitaemia suggest that susceptibility is associated with high expression of genes coding for chemokines (e.g. Ccl24, Ccl27 and Cxcl13), complement components (C1q and C3) and interferon receptor alpha (Ifnar1). Additionally, susceptible A/J mice expressed higher levels of some potassium channel genes. In contrast, messenger RNA levels of a few immune response, metabolism and protease genes (e.g. Prss7 and Mmp13) were higher in the tolerant C57BL/6 strain as compared to A/J.

Protein classification using ontology classification.

MOTIVATION: The classification of proteins expressed by an organism is an important step in understanding the molecular biology of that organism. Traditionally, this classification has been performed by human experts. Human knowledge can recognise the functional properties that are sufficient to place an individual gene product into a particular protein family group. Automation of this task usually fails to meet the 'gold standard' of the human annotator because of the difficult recognition stage. The growing number of genomes, the rapid changes in knowledge and the central role of classification in the annotation process, however, motivates the need to automate this process. RESULTS: We capture human understanding of how to recognise members of the protein phosphatases family by domain architecture as an ontology. By describing protein instances in terms of the domains they contain, it is possible to use description logic reasoners and our ontology to assign those proteins to a protein family class. We have tested our system on classifying the protein phosphatases of the human and Aspergillus fumigatus genomes and found that our knowledge-based, automatic classification matches, and sometimes surpasses, that of the human annotators. We have made the classification process fast and reproducible and, where appropriate knowledge is available, the method can potentially be generalised for use with any protein family. AVAILABILITY: All components described in this paper are freely available. OWL ontology http://www.bioinf.man.ac.uk/phosphabase myGrid http://www.mygrid.org.uk Instance Store http://instancestore.man.ac.uk.

Identification of novel genes in intestinal tissue that are regulated after infection with an intestinal nematode parasite.

Infection of resistant or susceptible mice with Trichuris muris provokes mesenteric lymph node responses which are polarized towards Th2 or Th1, respectively. These responses are well documented in the literature. In contrast, little is known about the local responses occurring within the infected intestine. Through microarray analyses, we demonstrate that the gene expression profile of infected gut tissue differs according to whether the parasite is expelled or not. Genes differentially regulated postinfection in resistant BALB/c mice include several antimicrobial genes, in particular, intelectin (Itln). In contrast, analyses in AKR mice which ultimately progress to chronic infection provide evidence for a Th1-dominated mucosa with up-regulated expression of genes regulated by gamma interferon. Increases in the expression of genes associated with tryptophan metabolism were also apparent with the coinduction of tryptophanyl tRNA synthetase (Wars) and indoleamine-2,3-dioxygenase (Indo). With the emerging literature on the role of these gene products in the suppression of T-cell responses in vitro and in vivo, their up-regulated expression here may suggest a role for tryptophan metabolism in the parasite survival strategy.

Constructing ontology-driven protein family databases.

MOTIVATION: Protein family databases provide a central focus for scientific communities as well as providing useful resources to aide research. However, such resources require constant curation and often become outdated and discontinued. We have developed an ontology-driven system for capturing and managing protein family data that addresses the problems of maintenance and sustainability. RESULTS: Using protein phosphatases and ABC transporters as model protein families, we constructed two protein family database resources around a central DAML+OIL ontology. Each resource contains specialist information about each protein family, providing specialized domain-specific resources based on the same template structure. The formal structure, combined with the extraction of biological data using GO terms, allows for automated update strategies. Despite the functional differences between the two protein families, the ontology model was equally applicable to both, demonstrating the generic nature of the system. AVAILABILITY: The protein phosphatase resource, PhosphaBase, is freely available on the internet (http://www.bioinf.man.ac.uk/phosphabase). The DAML+OIL ontology for the protein phosphatases and the ABC transporters is available on request from the authors. CONTACT: kwolstencroft@cs.man.ac.uk.

PhosphaBase: an ontology-driven database resource for protein phosphatases.

PhosphaBase is an ontology-driven database resource containing information on the protein phosphatase family. It is the first public resource dedicated to protein phosphatases, which are enzymes that perform dephosphorylation reactions. In conjunction with the phosphorylation action of protein kinases, phosphatases are involved in important control and communication mechanisms in the cell. They have also been implicated in many human diseases, including diabetes and obesity, cancers, and neurodegenerative conditions. PhosphaBase aims to centralize the growing base of knowledge in the phosphatase research domain. The resource is built around a formal, domain-specific DAML+OIL ontology, and the data are collected from heterogeneous biological sources using Gene Ontology terms as a means of data extraction. The overall ontology-driven architecture provides a robust structure with distinct advantages for sustainability and provides the potential for the development of diagnostic tools, as well as a data repository.

Exploring Williams-Beuren syndrome using myGrid.

MOTIVATION: In silico experiments necessitate the virtual organization of people, data, tools and machines. The scientific process also necessitates an awareness of the experience base, both of personal data as well as the wider context of work. The management of all these data and the co-ordination of resources to manage such virtual organizations and the data surrounding them needs significant computational infra-structure support. RESULTS: In this paper, we show that (my)Grid, middleware for the Semantic Grid, enables biologists to perform and manage in silico experiments, then explore and exploit the results of their experiments. We demonstrate (my)Grid in the context of a series of bioinformatics experiments focused on a 1.5 Mb region on chromosome 7 which is deleted in Williams-Beuren syndrome (WBS). Due to the highly repetitive nature of sequence flanking/in the WBS critical region (WBSCR), sequencing of the region is incomplete leaving documented gaps in the released sequence. (my)Grid was used in a series of experiments to find newly sequenced human genomic DNA clones that extended into these 'gap' regions in order to produce a complete and accurate map of the WBSCR. Once placed in this region, these DNA sequences were analysed with a battery of prediction tools in order to locate putative genes and regulatory elements possibly implicated in the disorder. Finally, any genes discovered were submitted to a range of standard bioinformatics tools for their characterization. We report how (my)Grid has been used to create workflows for these in silico experiments, run those workflows regularly and notify the biologist when new DNA and genes are discovered. The (my)Grid services collect and co-ordinate data inputs and outputs for the experiment, as well as much provenance information about the performance of experiments on WBS. AVAILABILITY: The (my)Grid software is available via http://www.mygrid.org.uk

Investigating semantic similarity measures across the Gene Ontology: the relationship between sequence and annotation.

MOTIVATION: Many bioinformatics data resources not only hold data in the form of sequences, but also as annotation. In the majority of cases, annotation is written as scientific natural language: this is suitable for humans, but not particularly useful for machine processing. Ontologies offer a mechanism by which knowledge can be represented in a form capable of such processing. In this paper we investigate the use of ontological annotation to measure the similarities in knowledge content or 'semantic similarity' between entries in a data resource. These allow a bioinformatician to perform a similarity measure over annotation in an analogous manner to those performed over sequences. A measure of semantic similarity for the knowledge component of bioinformatics resources should afford a biologist a new tool in their repertoire of analyses. RESULTS: We present the results from experiments that investigate the validity of using semantic similarity by comparison with sequence similarity. We show a simple extension that enables a semantic search of the knowledge held within sequence databases. AVAILABILITY: Software available from http://www.russet.org.uk.

Semantic similarity measures as tools for exploring the gene ontology.

Many bioinformatics resources hold data in the form of sequences. Often this sequence data is associated with a large amount of annotation. In many cases this data has been hard to model, and has been represented as scientific natural language, which is not readily computationally amenable. The development of the Gene Ontology provides us with a more accessible representation of some of this data. However it is not clear how this data can best be searched, or queried. Recently we have adapted information content based measures for use with the Gene Ontology (GO). In this paper we present detailed investigation of the properties of these measures, and examine various properties of GO, which may have implications for its future design.

Mutagenic and chemical modification of the ABA-1 allergen of the nematode Ascaris: consequences for structure and lipid binding properties.

The polyprotein allergens/antigens of nematodes (NPAs) are the only lipid binding proteins known to be produced as polyproteins. Cleavage of the large polyprotein precursors at regularly spaced proteinase cleavage sites produces 10 or 11 individual protein units of approximately 15 kDa. The sequences of these units are highly diverse within and between species, but there are five absolutely or strongly conserved amino acid positions (Trp15, Gln20, Leu42, Cys64, and Cys120). We have tested the role of these signature amino acids by mutational or chemical alteration of the ABA-1 protein of Ascaris, and examined the resulting modified proteins for perturbations of their lipid binding activities and structural integrity. Substitution of Trp15 and Gln20 both affect the stability of the protein in terms of resistance to thermal or chemical denaturation, but the ligand binding function is unaffected. Mutation of Leu42, however, disrupts both the protein's structural stability and functional integrity, as does chemical disruption of the disulfide bridge formed between Cys64 and Cys120. We also find that the C-terminal, but not the N-terminal, half of the protein binds fatty acids, indicating that the binding site may be confined to this part of the protein. This also supports the idea that the NPA units are themselves derived from an ancient duplication event, and that they may comprise two functionally distinct domains.

Detection of herpes viruses in clinical samples using real-time PCR.

Herpes viruses including cytomegalovirus, varicella zoster and herpes simplex are an important cause of morbidity and mortality, especially in immunocompromised patients. Real-time PCR assays were developed with the aim of introducing a rapid and sensitive test to replace culture, and as a surveillance system for high-risk patients. The assays were optimised using cell culture derived material, and the sensitivity ascertained using cloned product before applying to extracted and non-extracted clinical samples. The sensitivity was between 1--100 virus copies with increased sensitivity to detect less than 10 copies possible when an initial round of amplification was carried out using external primers. Results were available within four hours of receipt compared with a median of 4.4 days for culture and immunofluorescence. Real-time PCR was found to be a sensitive and rapid method of detecting these viruses and will be a valuable tool for the surveillance of immunosuppressed patients.

A classification of tasks in bioinformatics.

MOTIVATION: This paper reports on a survey of bioinformatics tasks currently undertaken by working biologists. The aim was to find the range of tasks that need to be supported and the components needed to do this in a general query system. This enabled a set of evaluation criteria to be used to assess both the biology and mechanical nature of general query systems. RESULTS: A classification of the biological content of the tasks gathered offers a checklist for those tasks (and their specialisations) that should be offered in a general bioinformatics query system. This semantic analysis was contrasted with a syntactic analysis that revealed the small number of components required to describe all bioinformatics questions. Both the range of biological tasks and syntactic task components can be seen to provide a set of bioinformatics requirements for general query systems. These requirements were used to evaluate two bioinformatics query systems.

Rapid real time PCR to distinguish between high risk human papillomavirus types 16 and 18.

AIMS: To assess the validity and practicality of real time polymerase chain reaction (PCR) for human papillomavirus (HPV) testing in combination with liquid based cytology samples for cervical screening. METHODS: Real time PCR using consensus (GPS+/6+) and type specific primers was developed to detect genital HPV types. This provides rapid, efficient amplification followed by denaturation of the product and computer analysis of the kinetics data that are generated. Liquid based cytology samples were obtained from patients attending routine cervical screening clinics. DNA was extracted from the residual cellular suspension after cytology using spin columns. RESULTS: Real time PCR successfully distinguished between HPV-16 and HPV-18 on the basis of amplification with consensus primers followed by DNA melting temperature (Tm) analysis. Sensitivities of one to 10 copies of HPV-16 (mean Tm = 79.4 degrees C; 2 SD, 0.8) and four to 40 copies of HPV-18 (mean Tm = 80.4 degrees C; 2 SD, 0.4) were obtained. In a mixed population of SiHa and HeLa cells containing known copy numbers of HPV-16 and HPV-18 genomes, HPV-16 and HPV-18 products were clearly separated by Tm analysis in mixtures varying from equivalence to 111000. Together with detailed melt analysis, type specific primers from the same region of the L1 gene confirmed the differential ability of this system. The method was applied to 100 liquid based cytology samples where HPV status using conventional GP5+/6+ PCR was already known. There was 95% agreement between the methods, with 55 positives detected by conventional PCR and 59 with real time PCR. The method was then tested on 200 routine liquid based cytology samples. Approximately 10% were positive by real time PCR, most of which were classified as HPV-16 by detailed melt analysis. Thirteen (6.8%) HPV positives were identified in 189 samples showing no evidence of cervical cytological abnormality. CONCLUSIONS: Real time PCR is a rapid, efficient method for the detection of HPV with the separation of HPV-16 and HPV-18 on the basis of differential Tm. Preliminary results suggest it could prove

Apoptosis signals in atopy and asthma measured with cDNA arrays.

A variety of studies have stressed the importance of the control of inflammatory cell longevity and the balance of pro-survival and pro-apoptotic signalling. Recently, asthma was found to be associated with reduced apoptosis of inflammatory cells in lung tissue. The aim of the study was to investigate the systemic activation of apoptosis pathways using cDNA array technology in atopy and asthma. Eighteen atopic asthmatics (AA), eight atopic non-asthmatic (AN) and 14 healthy control subjects (C) were included in the study. Peripheral blood mononuclear cells were separated with gradient centrifugation, mRNA purified and the reverse-transcribed probes hybridized to cDNA arrays. The signals were compared by standardizing to the 100 most expressed genes and group differences assessed with the Mann-Whitney U-test. We found a concerted up-regulation of several pro-survival cytokines and growth factors in AN and AA. FAS and FASL were not differentially expressed, but FAST kinase was over-expressed in AN and AA. The tumour necrosis factor pathway was activated in AN and AA with increased cytokine and receptor levels and increased TRAF2, an intracellular signalling product. There were indications of a down-regulated p53 system. In contrast, the Bcl-2 family of genes showed a net pro-apoptotic profile in AN and AA. The group of caspases showed a constant gene expression pattern in all groups. In conclusion, significant differences in the expression of apoptosis-related genes were found in peripheral blood of atopic individuals with and without asthma. cDNA array technology proved to be useful and may be complementary to DNA-based studies in order to analyse interactive and multidimensional pathways as shown here for apoptosis.

Conceptual modelling of genomic information.

MOTIVATION: Genome sequencing projects are making available complete records of the genetic make-up of organisms. These core data sets are themselves complex, and present challenges to those who seek to store, analyse and present the information. However, in addition to the sequence data, high throughput experiments are making available distinctive new data sets on protein interactions, the phenotypic consequences of gene deletions, and on the transcriptome, proteome, and metabolome. The effective description and management of such data is of considerable importance to bioinformatics in the post-genomic era. The provision of clear and intuitive models of complex information is surprisingly challenging, and this paper presents conceptual models for a range of important emerging information resources in bioinformatics. It is hoped that these can be of benefit to bioinformaticians as they attempt to integrate genetic and phenotypic data with that from genomic sequences, in order to both assign gene functions and elucidate the different pathways of gene action and interaction. RESULTS: This paper presents a collection of conceptual (i.e. implementation-independent) data models for genomic data. These conceptual models are amenable to (more or less direct) implementation on different computing platforms.