Permutation Tests for Classification

We introduce and explore an approach to estimating statistical significance of classification accuracy, which is particularly useful in scientific applications of machine learning where high dimensionality of the data and the small number of training examples render most standard convergence bounds too loose to yield a meaningful guarantee of the generalization ability of the classifier. Instead, we estimate statistical significance of the observed classification accuracy, or the likelihood of observing such accuracy by chance due to spurious correlations of the high-dimensional data patterns with the class labels in the given training set. We adopt permutation testing, a non-parametric technique previously developed in classical statistics for hypothesis testing in the generative setting (i.e., comparing two probability distributions). We demonstrate the method on real examples from neuroimaging studies and DNA microarray analysis and suggest a theoretical analysis of the procedure that relates the asymptotic behavior of the test to the existing convergence bounds.

A tecnologia de microarranjos na identificação de genes de interesse na bovinocultura.

Como ferramenta de análise da expressão gênica, a tecnologia dos microarranjos de DNA permite a investigação de milhares de transcritos de um tecido, de maneira simultânea. Essa metodologia tem revolucionado diversas áreas do conhecimento, por meio do aumento substancial da capacidade analítica dos processos moleculares. Hoje, dentro da ciência animal, a disponibilidade desse método de investigação tem permitido aos pesquisadores identificar variações na expressão de determinados genes que possam ocorrer como respostas biológicas devido à condição experimental (idade, raça, estado fisiológico). Sob um aspecto prático, essa tecnologia tem seu papel na identificação de genes envolvidos na determinação de características fenotípicas de interesse econômico, os quais podem ser usados no desenvolvimento de métodos para a seleção de genótipos superiores dentro dos programas de Melhoramento Genético. Esse é um dos motivos de sua adoção em um dos projetos que compõe a carteira dos Macroprogramas da Embrapa, o MP1 Rede Genômica Animal. Assim, esta publicação tem como propósito posicionar o leitor a respeito da contribuição que a tecnologia de microarranjos de expressão tem trazido com respeito à identificação de genes e processos biológicos que atuam na manifestação de características de interesse econômico, em especial na bovinocultura. Além disso, pretende despertar o interesse da comunidade científica para uma área que ainda demanda muita pesquisa e à qual a Embrapa Informática Agropecuária tem dedicado um grande esforço: a análise dos dados gerados dos experimentos com microarranjos da Rede Genômica Animal.

Regulation of yeast oscillatory dynamics

Douglas B. Murray, Manfred Beckmann, and Hiroaki Kitano. (2007). Regulation of yeast oscillatory dynamics. Proceedings of the National Academy of Sciences of the USA, 104 (7), 2241-2246 Sponsorship: Solution-Oriented Research for Science and Technology Agency to the Systems Biology Institute /21st Century Center of Excellence Program and Special Coordination Program of the Ministry of Education, Sports, Culture, Science, and Technology to Keio University RAE2008

Porównawcza analiza udziału pokrewnych kinaz MAPK (MKK1, MKK2, MKK4 i MKK5) Arabidopsis thaliana w odpowiedzi na stres suszy

Wydział Biologii: Instytut Biologii Molekularnej i Biotechnologii

SIEGE: Smoking Induced Epithelial Gene Expression Database

The SIEGE (Smoking Induced Epithelial Gene Expression) database is a clinical resource for compiling and analyzing gene expression data from epithelial cells of the human intra-thoracic airway. This database supports a translational research study whose goal is to profile the changes in airway gene expression that are induced by cigarette smoke. RNA is isolated from airway epithelium obtained at bronchoscopy from current-, former- and never-smoker subjects, and hybridized to Affymetrix HG-U133A Genechips, which measure the level of expression of ~22 500 human transcripts. The microarray data generated along with relevant patient information is uploaded to SIEGE by study administrators using the database's web interface, found at http://pulm.bumc.bu.edu/siegeDB. PERL-coded scripts integrated with SIEGE perform various quality control functions including the processing, filtering and formatting of stored data. The R statistical package is used to import database expression values and execute a number of statistical analyses including t-tests, correlation coefficients and hierarchical clustering. Values from all statistical analyses can be queried through CGI-based tools and web forms found on the �Search� section of the database website. Query results are embedded with graphical capabilities as well as with links to other databases containing valuable gene resources, including Entrez Gene, GO, Biocarta, GeneCards, dbSNP and the NCBI Map Viewer.

Epidemiology and comparative analysis of Yersinia in Ireland

Yersiniosis is an acute or chronic enteric zoonosis caused by enteropathogenic Yersinia species. Although yersiniosis is predominantly associated with gastroenteric forms of infection, extraintestinal forms are often reported from the elderly or patients with predisposing factors. Yersiniosis is often reported in countries with cold and mild climates (Northern and Central Europe, New Zealand and North of Russian Federation). The Irish Health Protection Surveillance Centre (HPSC) currently records only 3-7 notified cases of yersiniosis per year. At the same time pathogenic Yersinia enterocolitica is recovered from pigs (main source of pathogenic Y. enterocolitica) at the levels similar to that observed in Yersinia endemic countries. Introduction of Yersinia selective culture procedures may increase Yersinia isolation rates. To establish whether the small number of notifications of human disease was an underestimate due to lack of specific selective culture for Yersinia we carried out a prospective culture study of faecal samples from outpatients with diarrhoea, with additional culture of appendix and throat swabs. Higher levels of anti-Yersinia seroprevalence than yersiniosis notification rates in endemic countries suggests that most yersiniosis cases are unrecognised by culture. Subsequently, in addition to a prospective culture study of clinical specimens, we carried out serological screening of Irish blood donors and environmental screening of human sewage. Pathogenic Yersinia strains were not isolated from 1,189 faeces samples, nor from 297 throat swabs, or 23 appendix swabs. This suggested that current low notification rates in Ireland are not due to the lack of specific Yersinia culture procedures. Molecular screening detected a wider variety of Y. enterocolitica-specific targets in pig slurry than in human sewage. A serological survey for antibodies against Yersinia YOP (Yersinia Outer Proteins) proteins in Irish blood donors found antibodies in 25%, with an age-related trend to increased seropositivity, compatible with the hypothesis that yersiniosis may have been more prevalent in Ireland in the recent past. Y. enterocolitica is a heterogeneous group of microorganisms that comprises strains with different degree of pathogenicity. Although non-pathogenic Y. enterocolitica lack conventional virulence factors, these strains can be isolated from patients with diarrhoea. Insecticidal Toxin Complex (ITC) and Cytolethal Distending Toxins can potentially contribute to the virulence of non-pathogenic Y. enterocolitica in the absence of other virulence factors. We compared distribution of ITC and CDT loci among pathogenic and non-pathogenic Y. enterocolitica. Additionally, to demonstrate potential pathogenicity of non-pathogenic Y. enterocolitica we compared their virulence towards Galleria mellonella larvae (a non-mammalian model of human bacterial infections) with the virulence of highly and mildly pathogenic Y. enterocolitica strains. Surprisingly, virulence of pathogenic and non-pathogenic Y. enterocolitica in Galleria mellonella larvae observed at 37°C did not correlate with their pathogenic potential towards humans. Comparative phylogenomic analysis detects predicted coding sequences (CDSs) that define host-pathogen interactions and hence providing insights into molecular evolution of bacterial virulence. Comparative phylogenomic analysis of microarray data generated in Y. enterocolitica strains isolated in the Great Britain from humans with diarrhoea and domestic animals revealed high genetic heterogeneity of these species. Because of the extensive human, animal and food exchanges between the UK and Ireland the objective of this study was to gain further insight into genetic heterogeneity and relationships among clinical and non-clinical Y. enterocolitica strains of various pathogenic potential isolated in Ireland and Great Britain. No evidence of direct transfer of strains between the two countries was found.

Integrated genetic analysis systems

The overall objective of this thesis is to integrate a number of micro/nanotechnologies into integrated cartridge type systems to implement such biochemical protocols. Instrumentation and systems were developed to interface such cartridge systems: (i) implementing microfluidic handling, (ii) executing thermal control during biochemical protocols and (iii) detection of biomolecules associated with inherited or infectious disease. This system implements biochemical protocols for DNA extraction, amplification and detection. A digital microfluidic chip (ElectroWetting on Dielectric) manipulated droplets of sample and reagent implementing sample preparation protocols. The cartridge system also integrated a planar magnetic microcoil device to generate local magnetic field gradients, manipulating magnetic beads. For hybridisation detection a fluorescence microarray, screening for mutations associated with CFTR gene is printed on a waveguide surface and integrated within the cartridge. A second cartridge system was developed to implement amplification and detection screening for DNA associated with disease-causing pathogens e.g. Escherichia coli. This system incorporates (i) elastomeric pinch valves isolating liquids during biochemical protocols and (ii) a silver nanoparticle microarray for fluorescent signal enhancement, using localized surface plasmon resonance. The microfluidic structures facilitated the sample and reagent to be loaded and moved between chambers with external heaters implementing thermal steps for nucleic acid amplification and detection. In a technique allowing probe DNA to be immobilised within a microfluidic system using (3D) hydrogel structures a prepolymer solution containing probe DNA was formulated and introduced into the microfluidic channel. Photo-polymerisation was undertaken forming 3D hydrogel structures attached to the microfluidic channel surface. The prepolymer material, poly-ethyleneglycol (PEG), was used to form hydrogel structures containing probe DNA. This hydrogel formulation process was fast compared to conventional biomolecule immobilization techniques and was also biocompatible with the immobilised biomolecules, as verified by on-chip hybridisation assays. This process allowed control over hydrogel height growth at the micron scale.

Investigation of p75 neurotrophin receptor and the role of its post-translational modifications in tumour cell motility and invasiveness

The p75 neurotrophin receptor (p75NTR) is a member of the tumour necrosis factor superfamily, which relies on the recruitment of cytosolic protein partners - including the TNF receptor associated factor 6 (TRAF6) E3 ubiquitin ligase - to produce cellular responses such as apoptosis, survival, and inhibition of neurite outgrowth. Recently,p75NTR was also shown to undergo γ-secretase-mediated regulated intramembrane proteolysis, and the receptor ICD was found to migrate to the nucleus where it regulates gene transcription. Moreover, γ-secretase-mediated proteolysis was shown to be involved in glioblastoma cell migration and invasion. In this study we report that TRAF6-mediated K63-linked polyubiquitination at multiple or alternative lysine residues influences p75NTR-ICD stability in vitro. In addition, we found that TRAF6-mediated ubiquitination of p75NTR is not influenced by inhibition of dynamin. Moreover, we report beta-transducin repeats-containing protein (β-TrCP) as a novel E3- ligase that ubiquitinates p75NTR, which is independent of serine phosphorylation of the p75NTR destruction motif. In contrast to its influence on other substrates, co-expression of β-TrCP did not reduce p75NTR stability. We created U87-MG glioblastoma cell lines stably expressing wild type, γ-secretaseresistant and constitutively cleaved receptor, as well as the ICD-stabilized mutant K301R. Interestingly, only wild-type p75NTR induces increased glioblastoma cell migration, which could be reversed by application of γ-secretase inhibitor. Microarray and qRT-PCR analysis of mRNA transcripts in these cell lines yielded several promising genes that might be involved in glioblastoma cell migration and invasion, such as cadherin 11 and matrix metalloproteinase 12. Analysis of potential transcription factor binding sites revealed that transcription of these genes might be regulated by well known p75NTR signalling cascades such as NF-κB or JNK signalling, which are independent of γ-secretase-mediated cleavage of the receptor. In contrast, while p75NTR overexpression was confirmed in melanoma cell lines and a patient sample of melanoma metastasis to the brain, inhibition of γ-secretase did not influence melanoma cell migration. Collectively, this study provides several avenues to better understand the physiological importance of posttranslational modifications of p75NTR and the significance of the receptor in glioblastoma cell migration and invasion.

Early biomarkers to predict grade of encephalopathy following hypoxic ischaemic injury

The standard early markers for identifying and grading HIE severity, are not sufficient to ensure all children who would benefit from treatment are identified in a timely fashion. The aim of this thesis was to explore potential early biomarkers of HIE. Methods: To achieve this a cohort of infants with perinatal depression was prospectively recruited. All infants had cord blood samples drawn and biobanked, and were assessed with standardised neurological examination, and early continuous multi-channel EEG. Cord samples from a control cohort of healthy infants were used for comparison. Biomarkers studied included; multiple inflammatory proteins using multiplex assay; the metabolomics profile using LC/MS; and the miRNA profile using microarray. Results: Eighty five infants with perinatal depression were recruited. Analysis of inflammatory proteins consisted of exploratory analysis of 37 analytes conducted in a sub-population, followed by validation of all significantly altered analytes in the remaining population. IL-6 and IL-6 differed significantly in infants with a moderate/severely abnormal vs. a normal-mildly abnormal EEG in both cohorts (Exploratory: p=0.016, p=0.005: Validation: p=0.024, p=0.039; respectively). Metabolomic analysis demonstrated a perturbation in 29 metabolites. A Cross- validated Partial Least Square Discriminant Analysis model was developed, which accurately predicted HIE with an AUC of 0.92 (95% CI: 0.84-0.97). Analysis of the miRNA profile found 70 miRNA significantly altered between moderate/severely encephalopathic infants and controls. miRNA target prediction databases identified potential targets for the altered miRNA in pathways involved in cellular metabolism, cell cycle and apoptosis, cell signaling, and the inflammatory cascade. Conclusion: This thesis has demonstrated that the recruitment of a large cohortof asphyxiated infants, with cord blood carefully biobanked, and detailed early neurophysiological and clinical assessment recorded, is feasible. Additionally the results described, provide potential alternate and novel blood based biomarkers for the identification and assessment of HIE.

Assessment of the immunomodulatory effects of raw/farm milk and probiotic bacteria

The overall aims of this study were to investigate the differences between raw/farm milk and pasteurised milk with respect to potential immune modifying effects following consumption and investigate the bacterial composition of raw milk compared to pasteurised milk. Furthermore, in this thesis, panels of potential probiotic bacteria from the Bifidobacterium and Lactobacillus genera were investigated. The overall bacterial composition of raw milk was compared with pasteurised milk using samples obtained from commercial milk producers around Ireland using next generation sequencing technology (454 pyrosequencing). Here the presence of previously unrecognised and diverse bacterial populations in unpasteurised cow’s milk was identified. Futhermore the bacterial content of pasteurised milk was found to be more diverse than previously thought. The global response of the adenocarcinoma cell line HT-29 to raw milk and pasteurised milk exposures were also characterised using whole genome microarray technology. Over one thousand differentially expressed genes were identified which were found to be involved in a plethora of cellular functions. Interestingly a reduction in immune related activity (e.g. Major histocompatability complex class II signalling and T and B cell proliferation) was identified in cells exposed to pasteurised milk compared with raw milk exposures. Further studies comparing human cell response to raw versus pasteurised milk was performed using peripheral blood mononuclear cells (PBMC) from healthy donors. A reduction in CD14 was identified following raw milk exposures compared with pasteurised milk and the pattern of cytokine production may indicate that gram positive bacteria in the raw milk were contributing to the differences in the cellular response to raw versus pasteurised milk. Panels of potentially probiotic bacteria (comprising of lactobacilli and bifidobacteria) were further assessed for immunomodulatory capabilities using cell culture based models. Gene expression and cytokine production were used to evaluate stimulated and unstimulated (LPS) cellular responses as well as interaction mechanisms

Unfolded protein response genes regulated by CED-1 are required for Caenorhabditis elegans innate immunity.

The endoplasmic reticulum stress response, also known as the unfolded protein response (UPR), has been implicated in the normal physiology of immune defense and in several disorders, including diabetes, cancer, and neurodegenerative disease. Here, we show that the apoptotic receptor CED-1 and a network of PQN/ABU proteins involved in a noncanonical UPR response are required for proper defense to pathogen infection in Caenorhabditis elegans. A full-genome microarray analysis indicates that CED-1 functions to activate the expression of pqn/abu genes. We also show that ced-1 and pqn/abu genes are required for the survival of C. elegans exposed to live Salmonella enterica, and that overexpression of pqn/abu genes confers protection against pathogen-mediated killing. The results indicate that unfolded protein response genes, regulated in a CED-1-dependent manner, are involved in the C. elegans immune response to live bacteria.

Ferrochelatase is a conserved downstream target of the blue light-sensing White collar complex in fungi.

Light is a universal signal perceived by organisms, including fungi, in which light regulates common and unique biological processes depending on the species. Previous research has established that conserved proteins, originally called White collar 1 and 2 from the ascomycete Neurospora crassa, regulate UV/blue light sensing. Homologous proteins function in distant relatives of N. crassa, including the basidiomycetes and zygomycetes, which diverged as long as a billion years ago. Here we conducted microarray experiments on the basidiomycete fungus Cryptococcus neoformans to identify light-regulated genes. Surprisingly, only a single gene was induced by light above the commonly used twofold threshold. This gene, HEM15, is predicted to encode a ferrochelatase that catalyses the final step in haem biosynthesis from highly photoreactive porphyrins. The C. neoformans gene complements a Saccharomyces cerevisiae hem15Delta strain and is essential for viability, and the Hem15 protein localizes to mitochondria, three lines of evidence that the gene encodes ferrochelatase. Regulation of HEM15 by light suggests a mechanism by which bwc1/bwc2 mutants are photosensitive and exhibit reduced virulence. We show that ferrochelatase is also light-regulated in a white collar-dependent fashion in N. crassa and the zygomycete Phycomyces blakesleeanus, indicating that ferrochelatase is an ancient target of photoregulation in the fungal kingdom.

Co-orientation of replication and transcription preserves genome integrity.

In many bacteria, there is a genome-wide bias towards co-orientation of replication and transcription, with essential and/or highly-expressed genes further enriched co-directionally. We previously found that reversing this bias in the bacterium Bacillus subtilis slows replication elongation, and we proposed that this effect contributes to the evolutionary pressure selecting the transcription-replication co-orientation bias. This selection might have been based purely on selection for speedy replication; alternatively, the slowed replication might actually represent an average of individual replication-disruption events, each of which is counter-selected independently because genome integrity is selected. To differentiate these possibilities and define the precise forces driving this aspect of genome organization, we generated new strains with inversions either over approximately 1/4 of the chromosome or at ribosomal RNA (rRNA) operons. Applying mathematical analysis to genomic microarray snapshots, we found that replication rates vary dramatically within the inverted genome. Replication is moderately impeded throughout the inverted region, which results in a small but significant competitive disadvantage in minimal medium. Importantly, replication is strongly obstructed at inverted rRNA loci in rich medium. This obstruction results in disruption of DNA replication, activation of DNA damage responses, loss of genome integrity, and cell death. Our results strongly suggest that preservation of genome integrity drives the evolution of co-orientation of replication and transcription, a conserved feature of genome organization.

Age-specific differences in oncogenic pathway deregulation seen in human breast tumors.

PURPOSE: To define the biology driving the aggressive nature of breast cancer arising in young women. EXPERIMENTAL DESIGN: Among 784 patients with early stage breast cancer, using prospectively-defined, age-specific cohorts (young or=65 years), 411 eligible patients (n = 200or=65 years) with clinically-annotated Affymetrix microarray data were identified. GSEA, signatures of oncogenic pathway deregulation and predictors of chemotherapy sensitivity were evaluated within the two age-defined cohorts. RESULTS: In comparing deregulation of oncogenic pathways between age groups, a higher probability of PI3K (p = 0.006) and Myc (p = 0.03) pathway deregulation was observed in breast tumors arising in younger women. When evaluating unique patterns of pathway deregulation, a low probability of Src and E2F deregulation in tumors of younger women, concurrent with a higher probability of PI3K, Myc, and beta-catenin, conferred a worse prognosis (HR = 4.15). In contrast, a higher probability of Src and E2F pathway activation in tumors of older women, with concurrent low probability of PI3K, Myc and beta-catenin deregulation, was associated with poorer outcome (HR = 2.7). In multivariate analyses, genomic clusters of pathway deregulation illustrate prognostic value. CONCLUSION: Results demonstrate that breast cancer arising in young women represents a distinct biologic entity characterized by unique patterns of deregulated signaling pathways that are prognostic, independent of currently available clinico-pathologic variables. These results should enable refinement of targeted treatment strategies in this clinically challenging situation.

The genomic analysis of erythrocyte microRNA expression in sickle cell diseases.

BACKGROUND: Since mature erythrocytes are terminally differentiated cells without nuclei and organelles, it is commonly thought that they do not contain nucleic acids. In this study, we have re-examined this issue by analyzing the transcriptome of a purified population of human mature erythrocytes from individuals with normal hemoglobin (HbAA) and homozygous sickle cell disease (HbSS). METHODS AND FINDINGS: Using a combination of microarray analysis, real-time RT-PCR and Northern blots, we found that mature erythrocytes, while lacking ribosomal and large-sized RNAs, contain abundant and diverse microRNAs. MicroRNA expression of erythrocytes was different from that of reticulocytes and leukocytes, and contributed the majority of the microRNA expression in whole blood. When we used microRNA microarrays to analyze erythrocytes from HbAA and HbSS individuals, we noted a dramatic difference in their microRNA expression pattern. We found that miR-320 played an important role for the down-regulation of its target gene, CD71 during reticulocyte terminal differentiation. Further investigation revealed that poor expression of miR-320 in HbSS cells was associated with their defective downregulation CD71 during terminal differentiation. CONCLUSIONS: In summary, we have discovered significant microRNA expression in human mature erythrocytes, which is dramatically altered in HbSS erythrocytes and their defect in terminal differentiation. Thus, the global analysis of microRNA expression in circulating erythrocytes can provide mechanistic insights into the disease phenotypes of erythrocyte diseases.