Predicting interpretability of metabolome models based on behavior, putative identity, and biological relevance of explanatory signals

Enot, D. P., Beckmann, M., Overy, D., Draper, J. (2006). Predicting interpretability of metabolome models based on behavior, putative identity, and biological relevance of explanatory signals. Proceedings of the National Academy of Sciences of the USA, 103(40), 14865-14870. Sponsorship: BBSRC RAE2008

Autoimmunity in Chagas disease cardiopathy: biological relevance of a cardiac myosin-specific epitope crossreactive to an immunodominant Trypanosoma cruzi antigen.

Heart tissue destruction in chronic Chagas disease cardiopathy (CCC) may be caused by autoimmune recognition of heart tissue by a mononuclear cell infiltrate decades after Trypanosoma cruzi infection. Indirect evidence suggests that there is antigenic crossreactivity between T. cruzi and heart tissue. As there is evidence for immune recognition of cardiac myosin in CCC, we searched for a putative myosin-crossreactive T. cruzi antigen. T. cruzi lysate immunoblots were probed with anti-cardiac myosin heavy chain IgG antibodies (AMA) affinity-purified from CCC or asymptomatic Chagas disease patient-seropositive sera. A 140/116-kDa doublet was predominantly recognized by AMA from CCC sera. Further, recombinant T. cruzi protein B13--whose native protein is also a 140- and 116-kDa double band--was identified by crossreactive AMA. Among 28 sera tested in a dot-blot assay, AMA from 100% of CCC sera but only 14% of the asymptomatic Chagas disease sera recognized B13 protein (P = 2.3 x 10(-6)). Sequence homology to B13 protein was found at positions 8-13 and 1442-1447 of human cardiac myosin heavy chain. Competitive ELISA assays that used the correspondent myosin synthetic peptides to inhibit serum antibody binding to B13 protein identified the heart-specific AAALDK (1442-1447) sequence of human cardiac myosin heavy chain and the homologous AAAGDK B13 sequence as the respective crossreactive epitopes. The recognition of a heart-specific T. cruzi crossreactive epitope, in strong association with the presence of chronic heart lesions, suggests the involvement of crossreactivity between cardiac myosin and B13 in the pathogenesis of CCC.

The synthesis and application of novel profluorescent nitroxides as probes for polymer degradation

This PhD project has expanded the knowledge in the area of profluorescent nitroxides with regard to the synthesis and characterisations of novel profluorescent nitroxide probes as well as physical characterisation of the probe molecules in various polymer/physical environments. The synthesis of the first example of an azaphenalene-based fused aromatic nitroxide TMAO, [1,1,3,3-tetramethyl-2,3-dihydro-2-azaphenalen-2-yloxyl, was described. This novel nitroxide possesses some of the structural rigidity of the isoindoline class of nitroxides, as well as some properties akin to TEMPO nitroxides. Additionally, the integral aromatic ring imparts fluorescence that is switched on by radical scavenging reactions of the nitroxide, which makes it a sensitive probe for polymer degradation. In addition to the parent TMAO, 5 other azaphenalene derivatives were successfully synthesised. This new class of nitroxide was expected to have interesting redox properties when the structure was investigated by high-level ab initio molecular orbitals theory. This was expected to have implications with biological relevance as the calculated redox potentials for the azaphenalene ring class would make them potent antioxidant compounds. The redox potentials of 25 cyclic nitroxides from four different structural classes (pyrroline, piperidine, isoindoline and azaphenalene) were determined by cyclic voltammetry in acetonitrile. It was shown that potentials related to the one electron processes of the nitroxide were influenced by the type of ring system, ring substituents or groups surrounding the moiety. Favourable comparisons were found between theoretical and experimental potentials for pyrroline, piperidine and isoindoline ring classes. Substitution of these ring classes, were correctly calculated to have a small yet predictable effect on the potentials. The redox potentials of the azaphenalene ring class were underestimated by the calculations in all cases by at least a factor of two. This is believed to be due to another process influencing the redox potentials of the azaphenalene ring class which is not taken into account by the theoretical model. It was also possible to demonstrate the use of both azaphenalene and isoindoline nitroxides as additives for monitoring radical mediated damage that occurs in polypropylene as well as in more commercially relevant polyester resins. Polymer sample doped with nitroxide were exposed to both thermo-and photo-oxidative conditions with all nitroxides showing a protective effect. It was found that isoindoline nitroxides were able to indicate radical formation in polypropylene aged at elevated temperatures via fluorescence build-up. The azaphenalene nitroxide TMAO showed no such build-up of fluorescence. This was believed to be due to the more labile bond between the nitroxide and macromolecule and the protection may occur through a classical Denisov cycle, as is expected for commercially available HAS units. Finally, A new profluorescent dinitroxide, BTMIOA (9,10-bis(1,1,3,3- tetramethylisoindolin-2-yloxyl-5-yl)anthracene), was synthesised and shown to be a powerful probe for detecting changes during the initial stages of thermo-oxidative degradation of polypropylene. This probe, which contains a 9,10-diphenylanthracene core linked to two nitroxides, possesses strongly suppressed fluorescence due to quenching by the two nitroxide groups. This molecule also showed the greatest protective effect on thermo-oxidativly aged polypropylene. Most importantly, BTMIOA was found to be a valuable tool for imaging and mapping free-radical generation in polypropylene using fluorescence microscopy.

Exploring visual features through Gabor representations for facial expression detection

Gabor representations have been widely used in facial analysis (face recognition, face detection and facial expression detection) due to their biological relevance and computational properties. Two popular Gabor representations used in literature are: 1) Log-Gabor and 2) Gabor energy filters. Even though these representations are somewhat similar, they also have distinct differences as the Log-Gabor filters mimic the simple cells in the visual cortex while the Gabor energy filters emulate the complex cells, which causes subtle differences in the responses. In this paper, we analyze the difference between these two Gabor representations and quantify these differences on the task of facial action unit (AU) detection. In our experiments conducted on the Cohn-Kanade dataset, we report an average area underneath the ROC curve (A`) of 92.60% across 17 AUs for the Gabor energy filters, while the Log-Gabor representation achieved an average A` of 96.11%. This result suggests that small spatial differences that the Log-Gabor filters pick up on are more useful for AU detection than the differences in contours and edges that the Gabor energy filters extract.

Clinical dyslipidaemia is associated with changes in the lipid composition and inflammatory properties of apolipoprotein-B-containing lipoproteins from women with type 2 diabetes

The aim of this study was to use lipidomics to determine if the lipid composition of apolipoprotein-B-containing lipoproteins is modified by dyslipidaemia in type 2 diabetes and if any of the identified changes potentially have biological relevance in the pathophysiology of type 2 diabetes. VLDL and LDL from normolipidaemic and dyslipidaemic type 2 diabetic women and controls were isolated and quantified with HPLC and mass spectrometry. A detailed molecular characterisation of VLDL triacylglycerols (TAG) was also performed using the novel ozone-induced dissociation method, which allowed us to distinguish vaccenic acid (C18:1 n-7) from oleic acid (C18:1 n-9) in specific TAG species. Lipid class composition was very similar in VLDL and LDL from normolipidaemic type 2 diabetic and control participants. By contrast, dyslipidaemia was associated with significant changes in both lipid classes (e.g. increased diacylglycerols) and lipid species (e.g. increased C16:1 and C20:3 in phosphatidylcholine and cholesteryl ester and increased C16:0 [palmitic acid] and vaccenic acid in TAG). Levels of palmitic acid in VLDL and LDL TAG correlated with insulin resistance, and VLDL TAG enriched in palmitic acid promoted increased secretion of proinflammatory mediators from human smooth muscle cells. We showed that dyslipidaemia is associated with major changes in both lipid class and lipid species composition in VLDL and LDL from women with type 2 diabetes. In addition, we identified specific molecular lipid species that both correlate with clinical variables and are proinflammatory. Our study thus shows the potential of advanced lipidomic methods to further understand the pathophysiology of type 2 diabetes.

Phenotypical changes of osteocytes in osteoarthritis

This thesis focuses on the role of a bone cell, osteocytes in the progression of osteoarthritis. The biological relevance of this study has opened doors to new possibilities of understanding the pathogenesis of osteoarthritis.

Genetic variants underlying risk of endometriosis: Insights from meta-analysis of eight genome-wide association and replication datasets

BACKGROUND Endometriosis is a heritable common gynaecological condition influenced by multiple genetic and environmental factors. Genome-wide association studies (GWASs) have proved successful in identifying common genetic variants of moderate effects for various complex diseases. To date, eight GWAS and replication studies from multiple populations have been published on endometriosis. In this review, we investigate the consistency and heterogeneity of the results across all the studies and their implications for an improved understanding of the aetiology of the condition. METHODS Meta-analyses were conducted on four GWASs and four replication studies including a total of 11 506 cases and 32 678 controls, and on the subset of studies that investigated associations for revised American Fertility Society (rAFS) Stage III/IV including 2859 cases. The datasets included 9039 cases and 27 343 controls of European (Australia, Belgium, Italy, UK, USA) and 2467 cases and 5335 controls of Japanese ancestry. Fixed and Han and Elkin random-effects models, and heterogeneity statistics (Cochran's Q test), were used to investigate the evidence of the nine reported genome-wide significant loci across datasets and populations. RESULTS Meta-analysis showed that seven out of nine loci had consistent directions of effect across studies and populations, and six out of nine remained genome-wide significant (P < 5 × 10(-8)), including rs12700667 on 7p15.2 (P = 1.6 × 10(-9)), rs7521902 near WNT4 (P = 1.8 × 10(-15)), rs10859871 near VEZT (P = 4.7 × 10(-15)), rs1537377 near CDKN2B-AS1 (P = 1.5 × 10(-8)), rs7739264 near ID4 (P = 6.2 × 10(-10)) and rs13394619 in GREB1 (P = 4.5 × 10(-8)). In addition to the six loci, two showed borderline genome-wide significant associations with Stage III/IV endometriosis, including rs1250248 in FN1 (P = 8 × 10(-8)) and rs4141819 on 2p14 (P = 9.2 × 10(-8)). Two independent inter-genic loci, rs4141819 and rs6734792 on chromosome 2, showed significant evidence of heterogeneity across datasets (P < 0.005). Eight of the nine loci had stronger effect sizes among Stage III/IV cases, implying that they are likely to be implicated in the development of moderate to severe, or ovarian, disease. While three out of nine loci were inter-genic, the remaining were in or near genes with known functions of biological relevance to endometriosis, varying from roles in developmental pathways to cellular growth/carcinogenesis. CONCLUSIONS Our meta-analysis shows remarkable consistency in endometriosis GWAS results across studies, with little evidence of population-based heterogeneity. They also show that the phenotypic classifications used in GWAS to date have been limited. Stronger associations with Stage III/IV disease observed for most loci emphasize the importance for future studies to include detailed sub-phenotype information. Functional studies in relevant tissues are needed to understand the effect of the variants on downstream biological pathways.

Molecular cloning, expression and immunological characterisation of Pas n 1, the major allergen of Bahia grass Paspalum notatum pollen

Bahia grass, Paspalum notatum, is a clinically important subtropical grass with a prolonged pollination season from spring to autumn. We aimed to clone and characterise the major Bahia grass pollen allergen, Pas n 1. Grass pollen-allergic patients presenting to a tertiary hospital allergy clinic were tested for IgE reactivity with Bahia grass pollen extract by skin prick testing, ImmunoCAP, ELISA and immunoblotting. Using primers deduced from the N-terminal peptide sequence of a group 1 allergen of Bahia grass pollen extract separated by two-dimensional gel electrophoresis, the complete Pas n 1 cDNA was obtained by rapid amplification of cDNA ends and cloned. Biological relevance of recombinant Pas n 1 expressed in Escherichia coli was assessed by serum IgE reactivity and basophil activation. Twenty-nine of 34 (85%) consecutive patients presenting with grass pollen allergy were skin prick test positive to Bahia grass pollen. The Pas n 1 cDNA has sequence homology with the β-expansin 1 glycoprotein family and is more closely related to the maize pollen group 1 allergen (85% identity) than to ryegrass Lol p 1 or Timothy grass Phl p 1 (64 and 66% identity, respectively). rPas n 1 reacted with serum IgE in 47 of 55 (85%) Bahia grass pollen-allergic patients, activated basophils and inhibited serum IgE reactivity with the 29 kDa band of Bahia grass pollen extract. In conclusion the cDNA for the major group 1 allergen of the subtropical Bahia grass pollen, Pas n 1, was identified and cloned. rPas n 1 is immunologically active and is a valuable reagent for diagnosis and specific immunotherapy of grass pollen allergy.

Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs

An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.

The Pseudomonas syringae-derived HrpA pilins - molecular characterization and biotechnological application of the transcripts

Many Gram-negative bacteria pathogenic to plants and animals possess type III secretion systems that are used to cause disease. Effector proteins are injected into host cells using the type III secretion machineries. Despite vigorous studies, the nature of the secretion signal for type III secreted proteins still remains elusive. Both mRNA and proteinaceous signals have been proposed. Findings on coupling of translation to secretion by the type III secretion systems are also still contradictory. This study dealt with the secretion signal of HrpA from Pseudomonas syringae pathovar tomato. HrpA is the major component of the type III secretion system-associated Hrp pilus and a substrate for the type III secretion systems. The secretion signal was shown to reside in the first 15 codons or amino acids, a location typical for type III secretion signals. Translation of HrpA in the absence of a functional type III secretion system was established, but it does not exclude the possibility of coupling of translation to secretion when the secretion apparatus is present. The hrpA transcripts from various unrelated plant pathogenic bacteria were shown to be extremely stable. The biological relevance of this observation is unknown, but possible explanations include the high prevalence of HrpA protein, an mRNA secretion signal or timing of secretion. The hrpA mRNAs are stable over a wide range of temperatures, in the absence of translating ribosomes and even in the heterologous host Escherichia coli. The untranslated regions (UTRs) of hrpA transcripts from at least 20 pathovars of Pseudomonas syringae are highly homologous, whilst their coding regions exhibit low similarity. The stable nature of hrpA messenger RNAs is likely to be due to the folding of their 5 and 3 UTRs. In silico the UTRs seem to form stem-loop structures, the hairpin structures in the 3 UTRs being rich in guanidine and cytosine residues. The stable nature of the hrpA transcript redirected the studies to the stabilization of heterologous transcripts and to the use of stable messenger RNAs in recombinant protein production. Fragments of the hrpA transcript can be used to confer stability on heterologous transcripts from several sources of bacterial and eukaryotic origin, and to elevate the levels of production of the corresponding recombinant proteins several folds. hrpA transcript stabilizing elements can be used for improving the yields of recombinant proteins even in Escherichia coli, one of the most commonly used industrial protein production hosts.

Prediction of protein-protein interactions between Helicobacter pylori and a human host

A lack of information on protein-protein interactions at the host-pathogen interface is impeding the understanding of the pathogenesis process. A recently developed, homology search-based method to predict protein-protein interactions is applied to the gastric pathogen, Helicobacter pylori to predict the interactions between proteins of H. pylori and human proteins in vitro. Many of the predicted interactions could potentially occur between the pathogen and its human host during pathogenesis as we focused mainly on the H. pylori proteins that have a transmembrane region or are encoded in the pathogenic island and those which are known to be secreted into the human host. By applying the homology search approach to protein-protein interaction databases DIP and iPfam, we could predict in vitro interactions for a total of 623 H. pylori proteins with 6559 human proteins. The predicted interactions include 549 hypothetical proteins of as yet unknown function encoded in the H. pylori genome and 13 experimentally verified secreted proteins. We have recognized 833 interactions involving the extracellular domains of transmembrane proteins of H. pylori. Structural analysis of some of the examples reveals that the interaction predicted by us is consistent with the structural compatibility of binding partners. Examples of interactions with discernible biological relevance are discussed.

Effect of deletion of Ghrelin-O-Acyltransferase on the pulsatile release of growth hormone in mice

Ghrelin, a gut hormone originating from the post-translational cleavage of preproghrelin, is the endogenous ligand of growth hormone secretagogue receptor 1a (GHS-R1a). Within the growth hormone (GH) axis, the biological activity of ghrelin requires octanoylation by ghrelin-O-acyltransferase (GOAT), conferring selective binding to the GHS-R1a receptor via acylated ghrelin. Complete loss of preproghrelin-derived signalling (through deletion of the Ghrl gene) contributes to a decline in peak GH release; however, the selective contribution of endogenous acyl-ghrelin to pulsatile GH release remains to be established. We assessed the pulsatile release of GH in ad lib. fed male germline goat−/− mice, extending measures to include mRNA for key hypothalamic regulators of GH release, and peripheral factors that are modulated relative to GH release. The amount of GH released was reduced in young goat−/− mice compared to age-matched wild-type mice, whereas pulse frequency and irregularity increased. Altered GH release did not coincide with alterations in hypothalamic Ghrh, Srif, Npy or Ghsr mRNA expression, or pituitary GH content, suggesting that loss of Goat does not compromise canonical mechanisms that contribute to pituitary GH production and release. Although loss of Goat resulted in an irregular pattern of GH release (characterised by an increase in the number of GH pulses observed during extended secretory events), this did not contribute to a change in the expression of sexually dimorphic GH-dependent liver genes. Of interest, circulating levels of insulin-like growth factor (IGF)-1 were elevated in goat−/− mice. This rise in circulating levels of IGF-1 was correlated with an increase in GH pulse frequency, suggesting that sustained or increased IGF-1 release in goat−/− mice may occur in response to altered GH release patterning. Our observations demonstrate that germline loss of Goat alters GH release and patterning. Although the biological relevance of altered GH secretory patterning remains unclear, we propose that this may contribute to sustained IGF-1 release and growth in goat−/− mice.

Molecular features of colorectal cancer predisposition and progression

Both inherited genetic variations and somatically acquired mutations drive cancer development. The aim of this thesis was to gain insight into the molecular mechanisms underlying colorectal cancer (CRC) predisposition and tumor progression. Whereas one-third of CRC may develop in the context of hereditary predisposition, the known highly penetrant syndromes only explain a small fraction of all cases. Genome-wide association studies have shown that ten common single nucleotide polymorphisms (SNPs) modestly predispose to CRC. Our population-based sample series of around thousand CRC cases and healthy controls was genotyped for these SNPs. Tumors of heterozygous patients were analyzed for allelic imbalance, in an attempt to reveal the role of these SNPs in somatic tumor progression. The risk allele of rs6983267 at 8q24 was favored in the tumors significantly more often than the neutral allele, indicating that this germline variant is somatically selected for. No imbalance targeting the risk allele was observed in the remaining loci, suggesting that most of the low-penetrance CRC SNPs mainly play a role in the early stages of the neoplastic process. The ten SNPs were further analyzed in 788 CRC cases, 97 of which had a family history of CRC, to evaluate their combined contribution. A significant association appeared between the overall number of risk alleles and familial CRC and these ten SNPs seem to explain around 9% of the familial clustering of CRC. Finding more CRC susceptibility alleles may facilitate individualized risk prediction and cancer prevention in the future. Microsatellite instability (MSI), resulting from defective mismatch repair function, is a hallmark of Lynch syndrome and observed in a subset of all CRCs. Our aim was to identify microsatellite frameshift mutations that inactivate tumor suppressor genes in MSI CRCs. By sequencing microsatellite repeats of underexpressed genes we found six novel MSI target genes that were frequently mutated in 100 MSI CRCs: 51% in GLYR1, 47% in ABCC5, 43% in WDTC1, 33% in ROCK1, 30% in OR51E2, and 28% in TCEB3. Immunohistochemical staining of GLYR1 revealed defective protein expression in homozygously mutated tumors, providing further support for the loss of function hypothesis. Another mutation screening effort sought to identify MSI target genes with putative oncogenic functions. Microsatellites were similarly sequenced in genes that were overexpressed and, upon mutation, predicted to avoid nonsense-mediated mRNA decay. The mitotic checkpoint kinase TTK harbored protein-elongating mutations in 59% of MSI CRCs and the mutant protein was detected in heterozygous MSI CRC cells. No checkpoint dysregulation or defective protein localization was observable however, and the biological relevance of this mutation may hence be related to other mechanisms. In conclusion, these two large-scale and unbiased efforts identified frequently mutated genes that are likely to contribute to the development of this cancer type and may be utilized in developing diagnostic and therapeutic applications.