Creating a health informatics data resource for hearing health research.

BACKGROUND: The National Institute of Health and Social Care Research (NIHR) Health Informatics Collaborative (HIC) for Hearing Health has been established in the UK to curate routinely collected hearing health data to address research questions. This study defines priority research areas, outlines its aims, governance structure and demonstrates how hearing health data have been integrated into a common data model using pure tone audiometry (PTA) as a case study. METHODS: After identifying key research aims in hearing health, the governance structure for the NIHR HIC for Hearing Health is described. The Observational Medical Outcomes Partnership (OMOP) was chosen as our common data model to provide a case study example. RESULTS: The NIHR HIC Hearing Health theme have developed a data architecture outlying the flow of data from all of the various siloed electronic patient record systems to allow the effective linkage of data from electronic patient record systems to research systems. Using PTAs as an example, OMOPification of hearing health data successfully collated a rich breadth of datapoints across multiple centres. CONCLUSION: This study identified priority research areas where routinely collected hearing health data could be useful. It demonstrates integration and standardisation of such data into a common data model from multiple centres. By describing the process of data sharing across the HIC, we hope to invite more centres to contribute and utilise data to address research questions in hearing health. This national initiative has the power to transform UK hearing research and hearing care using routinely collected clinical data.

An improved analysis methodology for translational profiling by microarray.

Translational regulation plays a central role in the global gene expression of a cell, and detection of such regulation has allowed deciphering of critical biological mechanisms. Genome-wide studies of the regulation of translation (translatome) performed on microarrays represent a substantial proportion of studies, alongside with recent advances in deep-sequencing methods. However, there has been a lack of development in specific processing methodologies that deal with the distinct nature of translatome array data. In this study, we confirm that polysome profiling yields skewed data and thus violates the conventional transcriptome analysis assumptions. Using a comprehensive simulation of translatome array data varying the percentage and symmetry of deregulation, we show that conventional analysis methods (Quantile and LOESS normalizations) and statistical tests failed, respectively, to correctly normalize the data and to identify correctly deregulated genes (DEGs). We thus propose a novel analysis methodology available as a CRAN package; Internal Control Analysis of Translatome (INCATome) based on a normalization tied to a group of invariant controls. We confirm that INCATome outperforms the other normalization methods and allows a stringent identification of DEGs. More importantly, INCATome implementation on a biological translatome data set (cells silenced for splicing factor PSF) resulted in the best normalization performance and an improved validation concordance for identification of true positive DEGs. Finally, we provide evidence that INCATome is able to infer novel biological pathways with superior discovery potential, thus confirming the benefits for researchers of implementing INCATome for future translatome studies as well as for existing data sets to generate novel avenues for research.

In situ growth in early lung adenocarcinoma may represent precursor growth or invasive clone outgrowth-a clinically relevant distinction.

The switch from in situ to invasive tumor growth represents a crucial stage in the evolution of lung adenocarcinoma. However, the biological understanding of this shift is limited, and 'Noguchi Type C' tumors, being early lung adenocarcinomas with mixed in situ and invasive growth, represent those that are highly valuable in advancing our understanding of this process. All Noguchi Type C adenocarcinomas (n = 110) from the LATTICE-A cohort were reviewed and two patterns of in situ tumor growth were identified: those deemed likely to represent a true shift from precursor in situ to invasive disease ('Noguchi C1') and those in which the lepidic component appeared to represent outgrowth of the invasive tumor along existing airspaces ('Noguchi C2'). Overall Ki67 fraction was greater in C2 tumors and only C1 tumors showed significant increasing Ki67 from in situ to invasive disease. P53 positivity was acquired from in situ to invasive disease in C1 tumors but both components were positive in C2 tumors. Likewise, vimentin expression was increased from in situ to invasive tumor in C1 tumors only. Targeted next generation sequencing of 18 C1 tumors identified four mutations private to the invasive regions, including two in TP53, while 6 C2 tumors showed no private mutations. In the full LATTICe-A cohort, Ki67 fraction classified as either less than or greater than 10% within the in situ component of lung adenocarcinoma was identified as a strong predictor of patient outcome. This supports the proposition that tumors of all stages that have 'high grade' in situ components represent those with aggressive lepidic growth of the invasive clone. Overall these data support that the combined growth of Noguchi C tumors can represent two differing biological states and that 'Noguchi C1' tumors represent the genuine biological shift from in situ to invasive disease.

Evaluating bias-reducing protocols for RNA sequencing library preparation.

BACKGROUND: Next-generation sequencing does not yield fully unbiased estimates for read abundance, which may impact on the conclusions that can be drawn from sequencing data. The ligation step in RNA sequencing library generation is a known source of bias, motivating developments in enzyme technology and library construction protocols. We present the first comparison of the standard duplex adaptor protocol supplied by Life Technologies for use on the Ion Torrent PGM with an alternate single adaptor approach involving CircLigase (CircLig protocol).A correlation between over-representation in sequenced libraries and degree of secondary structure has been reported previously, therefore we also investigated whether bias could be reduced by ligation with an enzyme that functions at a temperature not permissive for such structure. RESULTS: A pool of small RNA fragments of known composition was converted into a sequencing library using one of three protocols and sequenced on an Ion Torrent PGM. The CircLig protocol resulted in less over-representation of specific sequences than the standard protocol. Over-represented sequences are more likely to be predicted to have secondary structure and to co-fold with adaptor sequences. However, use of the thermostable ligase Methanobacterium thermoautotrophicum RNA ligase K97A (Mth K97A) was not sufficient to reduce bias. CONCLUSIONS: The single adaptor CircLigase-based approach significantly reduces, but does not eliminate, bias in Ion Torrent data. Ligases that function at temperatures to remove the possible influence of secondary structure on library generation may be of value, although Mth K97A is not effective in this case.

Exploring listening-related fatigue in children with and without hearing loss using self-report and parent-proxy measures.

Children with hearing loss appear to experience greater fatigue than children with normal hearing (CNH). Listening-related fatigue is often associated with an increase in effortful listening or difficulty in listening situations. This has been observed in children with bilateral hearing loss (CBHL) and, more recently, in children with unilateral hearing loss (CUHL). Available tools for measuring fatigue in children include general fatigue questionnaires such as the child self-report and parent-proxy versions of the PedsQLTM-Multidimensional Fatigue Scale (MFS) and the PROMIS Fatigue Scale. Recently, the Vanderbilt Fatigue Scale (VFS-C: child self-report; VFS-P: parent-proxy report) was introduced with a specific focus on listening-related fatigue. The aims of this study were to compare fatigue levels experienced by CNH, CUHL and CBHL using both generic and listening-specific fatigue measures and compare outcomes from the child self-report and parent-proxy reports. Eighty children aged 6-16 years (32 CNH, 19 CUHL, 29 CBHL), and ninety-nine parents/guardians (39 parents to CNH, 23 parents to CUHL, 37 parents to CBHL), completed the above fatigue questionnaires online. Kruskal-Wallis H tests were performed to compare fatigue levels between the CNH, CUHL and CBHL. To determine the agreement between parent-proxy and child self-report measures, Bland-Altman 95% limits of agreement were performed. All child self-report fatigue measures indicated that CBHL experience greater fatigue than CNH. Only the listening-specific tool (VFS-C) was sufficiently able to show greater fatigue in CUHL than in CNH. Similarly, all parent-proxy measures of fatigue indicated that CBHL experience significantly greater fatigue than CNH. The VFS-P and the PROMIS Fatigue Parent-Proxy also showed greater fatigue in CUHL than in CNH. Agreement between the parent-proxy and child self-report measures were found within the PedsQL-MFS and the PROMIS Fatigue Scale. Our results suggest that CBHL experience greater levels of daily-life fatigue compared to CNH. CUHL also appear to experience more fatigue than CNH, and listening-specific measures of fatigue may be better able to detect this effect. Further research is needed to understand the bases of fatigue in these populations and to clarify whether fatigue experienced by CBHL and CUHL is comparable in nature and degree.

PiRaNhA: a server for the computational prediction of RNA-binding residues in protein sequences.

The PiRaNhA web server is a publicly available online resource that automatically predicts the location of RNA-binding residues (RBRs) in protein sequences. The goal of functional annotation of sequences in the field of RNA binding is to provide predictions of high accuracy that require only small numbers of targeted mutations for verification. The PiRaNhA server uses a support vector machine (SVM), with position-specific scoring matrices, residue interface propensity, predicted residue accessibility and residue hydrophobicity as features. The server allows the submission of up to 10 protein sequences, and the predictions for each sequence are provided on a web page and via email. The prediction results are provided in sequence format with predicted RBRs highlighted, in text format with the SVM threshold score indicated and as a graph which enables users to quickly identify those residues above any specific SVM threshold. The graph effectively enables the increase or decrease of the false positive rate. When tested on a non-redundant data set of 42 protein sequences not used in training, the PiRaNhA server achieved an accuracy of 85%, specificity of 90% and a Matthews correlation coefficient of 0.41 and outperformed other publicly available servers. The PiRaNhA prediction server is freely available at http://www.bioinformatics.sussex.ac.uk/PIRANHA.

Mitochondrial ROS signalling requires uninterrupted electron flow and is lost during ageing in flies.

Mitochondrial reactive oxygen species (mtROS) are cellular messengers essential for cellular homeostasis. In response to stress, reverse electron transport (RET) through respiratory complex I generates high levels of mtROS. Suppression of ROS production via RET (ROS-RET) reduces survival under stress, while activation of ROS-RET extends lifespan in basal conditions. Here, we demonstrate that ROS-RET signalling requires increased electron entry and uninterrupted electron flow through the electron transport chain (ETC). We find that in old fruit flies, ROS-RET is abolished when electron flux is decreased and that their mitochondria produce consistently high levels of mtROS. Finally, we demonstrate that in young flies, limiting electron exit, but not entry, from the ETC phenocopies mtROS generation observed in old individuals. Our results elucidate the mechanism by which ROS signalling is lost during ageing.

The pathogenesis of mesothelioma is driven by a dysregulated translatome.

Malignant mesothelioma (MpM) is an aggressive, invariably fatal tumour that is causally linked with asbestos exposure. The disease primarily results from loss of tumour suppressor gene function and there are no 'druggable' driver oncogenes associated with MpM. To identify opportunities for management of this disease we have carried out polysome profiling to define the MpM translatome. We show that in MpM there is a selective increase in the translation of mRNAs encoding proteins required for ribosome assembly and mitochondrial biogenesis. This results in an enhanced rate of mRNA translation, abnormal mitochondrial morphology and oxygen consumption, and a reprogramming of metabolic outputs. These alterations delimit the cellular capacity for protein biosynthesis, accelerate growth and drive disease progression. Importantly, we show that inhibition of mRNA translation, particularly through combined pharmacological targeting of mTORC1 and 2, reverses these changes and inhibits malignant cell growth in vitro and in ex-vivo tumour tissue from patients with end-stage disease. Critically, we show that these pharmacological interventions prolong survival in animal models of asbestos-induced mesothelioma, providing the basis for a targeted, viable therapeutic option for patients with this incurable disease.

Mitochondrial complex I derived ROS regulate stress adaptation in Drosophila melanogaster.

Reactive Oxygen Species (ROS) are essential cellular messengers required for cellular homeostasis and regulate the lifespan of several animal species. The main site of ROS production is the mitochondrion, and within it, respiratory complex I (CI) is the main ROS generator. ROS produced by CI trigger several physiological responses that are essential for the survival of neurons, cardiomyocytes and macrophages. Here, we show that CI produces ROS when electrons flow in either the forward (Forward Electron Transport, FET) or reverse direction (Reverse Electron Transport, RET). We demonstrate that ROS production via RET (ROS-RET) is activated under thermal stress conditions and that interruption of ROS-RET production, through ectopic expression of the alternative oxidase AOX, attenuates the activation of pro-survival pathways in response to stress. Accordingly, we find that both suppressing ROS-RET signalling or decreasing levels of mitochondrial H2O2 by overexpressing mitochondrial catalase (mtCAT), reduces survival dramatically in flies under stress. Our results uncover a specific ROS signalling pathway where hydrogen peroxide (H2O2) generated by CI via RET is required to activate adaptive mechanisms, maximising survival under stress conditions.

Nolz1 expression is required in dopaminergic axon guidance and striatal innervation.

Midbrain dopaminergic (DA) axons make long longitudinal projections towards the striatum. Despite the importance of DA striatal innervation, processes involved in establishment of DA axonal connectivity remain largely unknown. Here we demonstrate a striatal-specific requirement of transcriptional regulator Nolz1 in establishing DA circuitry formation. DA projections are misguided and fail to innervate the striatum in both constitutive and striatal-specific Nolz1 mutant embryos. The lack of striatal Nolz1 expression results in nigral to pallidal lineage conversion of striatal projection neuron subtypes. This lineage switch alters the composition of secreted factors influencing DA axonal tract formation and renders the striatum non-permissive for dopaminergic and other forebrain tracts. Furthermore, transcriptomic analysis of wild-type and Nolz1-/- mutant striatal tissue led to the identification of several secreted factors that underlie the observed guidance defects and proteins that promote DA axonal outgrowth. Together, our data demonstrate the involvement of the striatum in orchestrating dopaminergic circuitry formation.

CD40L/IL-4-stimulated CLL demonstrates variation in translational regulation of DNA damage response genes including ATM.

CD40L/interleukin-4 (IL-4) stimulation occurs in vivo in the tumor microenvironment and induces global translation to varying degrees in individuals with chronic lymphocytic leukemia (CLL) in vitro. However, the implications of CD40L/IL-4 for the translation of specific genes is not known. To determine the most highly translationally regulated genes in response to CD40L/IL-4, we carried out ribosome profiling, a next-generation sequencing method. Significant differences in the translational efficiency of DNA damage response genes, specifically ataxia-telangiectasia-mutated kinase (ATM) and the MRE11/RAD50/NBN (MRN) complex, were observed between patients, suggesting different patterns of translational regulation. We confirmed associations between CD40L/IL-4 response and baseline ATM levels, induction of ATM, and phosphorylation of the ATM targets, p53 and H2AX. X-irradiation was used to demonstrate that CD40L/IL-4 stimulation tended to improve DNA damage repair. Baseline ATM levels, independent of the presence of 11q deletion, correlated with overall survival (OS). Overall, we suggest that there are individual differences in translation of specific genes, including ATM, in response to CD40L/IL-4 and that these interpatient differences might be clinically important.

RTN3 Is a Novel Cold-Induced Protein and Mediates Neuroprotective Effects of RBM3.

Cooling and hypothermia are profoundly neuroprotective, mediated, at least in part, by the cold shock protein, RBM3. However, the neuroprotective effector proteins induced by RBM3 and the mechanisms by which mRNAs encoding cold shock proteins escape cooling-induced translational repression are unknown. Here, we show that cooling induces reprogramming of the translatome, including the upregulation of a new cold shock protein, RTN3, a reticulon protein implicated in synapse formation. We report that this has two mechanistic components. Thus, RTN3 both evades cooling-induced translational elongation repression and is also bound by RBM3, which drives the increased expression of RTN3. In mice, knockdown of RTN3 expression eliminated cooling-induced neuroprotection. However, lentivirally mediated RTN3 overexpression prevented synaptic loss and cognitive deficits in a mouse model of neurodegeneration, downstream and independently of RBM3. We conclude that RTN3 expression is a mediator of RBM3-induced neuroprotection, controlled by novel mechanisms of escape from translational inhibition on cooling.

A ligand-centric analysis of the diversity and evolution of protein-ligand relationships in E.coli.

As enzymes evolve and diverge from common ancestor sequences, they often keep their overall reaction chemistry but specialize in the binding of different cognate ligands. This study borrows methods for the computational assessment of 2D similarity of small molecules from the field of chemoinformatics, to examine the extent of structure conservation of cognate ligands binding to similar proteins. Proteins from 87 structural superfamilies from Escherichia coli form the core dataset, which is extended using homologues with functional assignments from any organism. We find that correlation of the substrate similarity with protein similarity (measured by either sequence-based or structure-based scores) can only be clearly established for very similar proteins. At low sequence identities, the superfamily to which a protein belongs can give helpful clues to its function, and more importantly, the confidence attached to such clues is superfamily-dependent. Our data indicate that only a few superfamilies show great substrate diversity, and that most exhibit conservation of at least part of the structural scaffold of the substrate.

The complement of enzymatic sets in different species.

We present here a comprehensive analysis of the complement of enzymes in a large variety of species. As enzymes are a relatively conserved group there are several classification systems available that are common to all species and link a protein sequence to an enzymatic function. Enzymes are therefore an ideal functional group to study the relationship between sequence expansion, functional divergence and phenotypic changes. By using information retrieved from the well annotated SWISS-PROT database together with sequence information from a variety of fully sequenced genomes and information from the EC functional scheme we have aimed here to estimate the fraction of enzymes in genomes, to determine the extent of their functional redundancy in different domains of life and to identify functional innovations and lineage specific expansions in the metazoa lineage. We found that prokaryote and eukaryote species differ both in the fraction of enzymes in their genomes and in the pattern of expansion of their enzymatic sets. We observe an increase in functional redundancy accompanying an increase in species complexity. A quantitative assessment was performed in order to determine the degree of functional redundancy in different species. Finally, we report a massive expansion in the number of mammalian enzymes involved in signalling and degradation.

Temporal Regulation of Distinct Internal Ribosome Entry Sites of the Dicistroviridae Cricket Paralysis Virus.

Internal ribosome entry is a key mechanism for viral protein synthesis in a subset of RNA viruses. Cricket paralysis virus (CrPV), a member of Dicistroviridae, has a positive-sense single strand RNA genome that contains two internal ribosome entry sites (IRES), a 5'untranslated region (5'UTR) and intergenic region (IGR) IRES, that direct translation of open reading frames (ORF) encoding the viral non-structural and structural proteins, respectively. The regulation of and the significance of the CrPV IRESs during infection are not fully understood. In this study, using a series of biochemical assays including radioactive-pulse labelling, reporter RNA assays and ribosome profiling, we demonstrate that while 5'UTR IRES translational activity is constant throughout infection, IGR IRES translation is delayed and then stimulated two to three hours post infection. The delay in IGR IRES translation is not affected by inhibiting global translation prematurely via treatment with Pateamine A. Using a CrPV replicon that uncouples viral translation and replication, we show that the increase in IGR IRES translation is dependent on expression of non-structural proteins and is greatly stimulated when replication is active. Temporal regulation by distinct IRESs within the CrPV genome is an effective viral strategy to ensure optimal timing and expression of viral proteins to facilitate infection.

SCOPEC: a database of protein catalytic domains.

MOTIVATION: Domains are the units of protein structure, function and evolution. It is therefore essential to utilize knowledge of domains when studying the evolution of function, or when assigning function to genome sequence data. For this purpose, we have developed a database of catalytic domains, SCOPEC, by combining structural domain information from SCOP, full-length sequence information from Swiss-Prot, and verified functional information from the Enzyme Classification (EC) database. Two major problems need to be overcome to create a database of domain-function relationships; (1) for sequences, EC numbers are typically assigned to whole sequences rather than the functional unit, and (2) The Protein Data Bank (PDB) structures elucidated from a larger multi-domain protein will often have EC annotation although the relevant catalytic domain may lie elsewhere. RESULTS: SCOPEC entries have high quality enzyme assignments; having passed both computational and manual checks. SCOPEC currently contains entries for 75% of all EC annotations in the PDB. Overall, EC number is fairly well conserved within a superfamily, even when the proteins are distantly related. Initial analysis is encouraging; suggesting that there is a 50:50 chance of conserved function in distant homologues first detected by a third iteration PSI-BLAST search. Therefore, we envisage that a knowledge-based approach to function assignment using the domain-EC relationships in SCOPEC will gain a marked improvement over this base line. AVAILABILITY: The SCOPEC database is a valuable resource in the analysis and prediction of protein structure and function. It can be obtained or queried at our website http://www.enzome.com

TRANS PROF DB: A new resource for sharing translational profiles.

The translational efficiency of individual mRNAs can be measured on a genome-wide scale using translational profiling techniques. Data from such experiments are an enormously important resource in the quest to understand the impact of cellular state on gene expression. To improve our understanding of these data, we have created TRANS PROF DB, a manually curated resource containing the translational status of human mRNAs under defined conditions. Results are provided at the level of an annotated conclusion for each gene, e.g. "Translation up-regulated", and also, where available, in a rawer form such as the normalized analyzed microarray output. TRANS PROF DB aims to provide a central resource for the sharing of translational profiles to facilitate reuse of published data and to enable meta-analyses across data sets. As the database expands, it will provide an easily searchable archive of publicly available translational profiling data sets. We encourage all researchers to deposit their translational profiling data into TRANS PROF DB to enable us to create a truly comprehensive resource. TRANS PROF DB is available without restriction at mrctools.mrctox.le.ac.uk/TRANS_PROF_DB.

Glycogen synthase kinase-3 inhibition enhances translation of pluripotency-associated transcription factors to contribute to maintenance of mouse embryonic stem cell self-renewal.

Maintenance of embryonic stem cell (ESC) self-renewal and pluripotency are controlled by extrinsic factors, molecular signaling pathways and transcriptional regulators. While many of the key players have been studied in depth, how the molecular signals interact with transcription factors of the pluripotency network to regulate their action remains less well understood. Inhibition of glycogen synthase kinase 3 (Gsk-3) has been implicated in the maintenance of mouse ESC pluripotency, although there is contradictory data on its role, with enhancement of cell survival and metabolism, stabilisation of c-Myc and activation of Wnt signalling proposed as potential mechanisms. We have discovered that suppression of Gsk-3 activity leads to enhanced protein levels of key transcriptional regulators of the pluripotency network, notably Nanog, Tbx3 and c-Myc. Protein stability was unchanged following Gsk-3 inhibition, although interestingly, Nanog and Tbx3 proteins were found to have half-lives of 1-3 h, while that of Oct4 protein was longer, at 6 h. We demonstrate that the effects on protein levels seen following inhibition of Gsk-3 are due to both enhanced de novo synthesis of Nanog protein and increases in the proportion of Nanog and Tbx3 RNAs bound to polysomes, findings consistent with Gsk-3 regulating translation of these factors. These effects were not due to changes in regulators of general translation initiation machinery nor mediated via the 5' or 3' UTR sequences of Nanog alone. The data we present provide both new conceptual insight into the mechanisms regulated by Gsk-3 that may contribute to ESC self-renewal and, importantly, establish control of protein translation as an additional mechanism involved in modulation of ESC pluripotency.

RNA-binding residues in sequence space: conservation and interaction patterns.

RNA-binding proteins (RBPs) perform fundamental and diverse functions within the cell. Approximately 15% of proteins sequences are annotated as RNA-binding, but with a significant number of proteins without functional annotation, many RBPs are yet to be identified. A percentage of uncharacterised proteins can be annotated by transferring functional information from proteins sharing significant sequence homology. However, genomes contain a significant number of orphan open reading frames (ORFs) that do not share significant sequence similarity to other ORFs, but correspond to functional proteins. Hence methods for protein function annotation that go beyond sequence homology are essential. One method of annotation is the identification of ligands that bind to proteins, through the characterisation of binding site residues. In the current work RNA-binding residues (RBRs) are characterised in terms of their evolutionary conservation and the patterns they form in sequence space. The potential for such characteristics to be used to identify RBPs from sequence is then evaluated. In the current work the conservation of residues in 261 RBPs is compared for (a) RBRs vs. non-RBRs surface residues, and for (b) specific and non-specific RBRs. The analysis shows that RBRs are more conserved than other surface residues, and RBRs hydrogen-bonded to the RNA backbone are more conserved than those making hydrogen bonds to RNA bases. This observed conservation of RBRs was then used to inform the construction of RBR sequence patterns from known protein-RNA structures. A series of RBR patterns were generated for a case study protein aspartyl-tRNA synthetase bound to tRNA; and used to differentiate between RNA-binding and non-RNA-binding protein sequences. Six sequence patterns performed with high precision values of >80% and recall values 7 times that of an homology search. When the method was expanded to the complete dataset of 261 proteins, many patterns were of poor predictive value, as they had not been manipulated on a family-specific basis. However, two patterns with precision values > or = 85% were used to make function predictions for a set of hypothetical proteins. This revealed a number of potential RBPs that require experimental verification.