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

Association of candidate genes with flowering time and water-soluble carbohydrate content in Lolium perenne (L.)

Sk?t, L., Humphreys, J., Humphreys, M. O., Thorogood, D., Gallagher, J. A., Sanderson, R., Armstead, I. P., Thomas, I. D. (2007). Association of candidate genes with flowering time and water-soluble carbohydrate content in Lolium perenne (L.). Genetics, 177 (1), 535-547. Sponsorship: BBSRC 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.

Quantitative Analysis of Single Nucleotide Polymorphisms within Copy Number Variation

Background Single nucleotide polymorphisms (SNPs) have been used extensively in genetics and epidemiology studies. Traditionally, SNPs that did not pass the Hardy-Weinberg equilibrium (HWE) test were excluded from these analyses. Many investigators have addressed possible causes for departure from HWE, including genotyping errors, population admixture and segmental duplication. Recent large-scale surveys have revealed abundant structural variations in the human genome, including copy number variations (CNVs). This suggests that a significant number of SNPs must be within these regions, which may cause deviation from HWE. Results We performed a Bayesian analysis on the potential effect of copy number variation, segmental duplication and genotyping errors on the behavior of SNPs. Our results suggest that copy number variation is a major factor of HWE violation for SNPs with a small minor allele frequency, when the sample size is large and the genotyping error rate is 0~1%. Conclusions Our study provides the posterior probability that a SNP falls in a CNV or a segmental duplication, given the observed allele frequency of the SNP, sample size and the significance level of HWE testing.

Framingham Heart Study 100K project: genome-wide associations for cardiovascular disease outcomes

BACKGROUND:Cardiovascular disease (CVD) and its most common manifestations - including coronary heart disease (CHD), stroke, heart failure (HF), and atrial fibrillation (AF) - are major causes of morbidity and mortality. In many industrialized countries, cardiovascular disease (CVD) claims more lives each year than any other disease. Heart disease and stroke are the first and third leading causes of death in the United States. Prior investigations have reported several single gene variants associated with CHD, stroke, HF, and AF. We report a community-based genome-wide association study of major CVD outcomes.METHODS:In 1345 Framingham Heart Study participants from the largest 310 pedigrees (54% women, mean age 33 years at entry), we analyzed associations of 70,987 qualifying SNPs (Affymetrix 100K GeneChip) to four major CVD outcomes: major atherosclerotic CVD (n = 142; myocardial infarction, stroke, CHD death), major CHD (n = 118; myocardial infarction, CHD death), AF (n = 151), and HF (n = 73). Participants free of the condition at entry were included in proportional hazards models. We analyzed model-based deviance residuals using generalized estimating equations to test associations between SNP genotypes and traits in additive genetic models restricted to autosomal SNPs with minor allele frequency [greater than or equal to]0.10, genotype call rate [greater than or equal to]0.80, and Hardy-Weinberg equilibrium p-value [greater than or equal to] 0.001.RESULTS:Six associations yielded p <10-5. The lowest p-values for each CVD trait were as follows: major CVD, rs499818, p = 6.6 x 10-6; major CHD, rs2549513, p = 9.7 x 10-6; AF, rs958546, p = 4.8 x 10-6; HF: rs740363, p = 8.8 x 10-6. Of note, we found associations of a 13 Kb region on chromosome 9p21 with major CVD (p 1.7 - 1.9 x 10-5) and major CHD (p 2.5 - 3.5 x 10-4) that confirm associations with CHD in two recently reported genome-wide association studies. Also, rs10501920 in CNTN5 was associated with AF (p = 9.4 x 10-6) and HF (p = 1.2 x 10-4). Complete results for these phenotypes can be found at the dbgap website http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007.CONCLUSION:No association attained genome-wide significance, but several intriguing findings emerged. Notably, we replicated associations of chromosome 9p21 with major CVD. Additional studies are needed to validate these results. Finding genetic variants associated with CVD may point to novel disease pathways and identify potential targeted preventive therapies.

Genome-wide association with bone mass and geometry in the Framingham Heart Study

BACKGROUND:Osteoporosis is characterized by low bone mass and compromised bone structure, heritable traits that contribute to fracture risk. There have been no genome-wide association and linkage studies for these traits using high-density genotyping platforms.METHODS:We used the Affymetrix 100K SNP GeneChip marker set in the Framingham Heart Study (FHS) to examine genetic associations with ten primary quantitative traits: bone mineral density (BMD), calcaneal ultrasound, and geometric indices of the hip. To test associations with multivariable-adjusted residual trait values, we used additive generalized estimating equation (GEE) and family-based association tests (FBAT) models within each sex as well as sexes combined. We evaluated 70,987 autosomal SNPs with genotypic call rates [greater than or equal to]80%, HWE p [greater than or equal to] 0.001, and MAF [greater than or equal to]10% in up to 1141 phenotyped individuals (495 men and 646 women, mean age 62.5 yrs). Variance component linkage analysis was performed using 11,200 markers.RESULTS:Heritability estimates for all bone phenotypes were 30-66%. LOD scores [greater than or equal to]3.0 were found on chromosomes 15 (1.5 LOD confidence interval: 51,336,679-58,934,236 bp) and 22 (35,890,398-48,603,847 bp) for femoral shaft section modulus. The ten primary phenotypes had 12 associations with 100K SNPs in GEE models at p < 0.000001 and 2 associations in FBAT models at p < 0.000001. The 25 most significant p-values for GEE and FBAT were all less than 3.5 x 10-6 and 2.5 x 10-5, respectively. Of the 40 top SNPs with the greatest numbers of significantly associated BMD traits (including femoral neck, trochanter, and lumbar spine), one half to two-thirds were in or near genes that have not previously been studied for osteoporosis. Notably, pleiotropic associations between BMD and bone geometric traits were uncommon. Evidence for association (FBAT or GEE p < 0.05) was observed for several SNPs in candidate genes for osteoporosis, such as rs1801133 in MTHFR; rs1884052 and rs3778099 in ESR1; rs4988300 in LRP5; rs2189480 in VDR; rs2075555 in COLIA1; rs10519297 and rs2008691 in CYP19, as well as SNPs in PPARG (rs10510418 and rs2938392) and ANKH (rs2454873 and rs379016). All GEE, FBAT and linkage results are provided as an open-access results resource at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007.CONCLUSION:The FHS 100K SNP project offers an unbiased genome-wide strategy to identify new candidate loci and to replicate previously suggested candidate genes for osteoporosis.

Genetic correlates of longevity and selected age-related phenotypes: a genome-wide association study in the Framingham Study

BACKGROUND: Family studies and heritability estimates provide evidence for a genetic contribution to variation in the human life span. METHODS:We conducted a genome wide association study (Affymetrix 100K SNP GeneChip) for longevity-related traits in a community-based sample. We report on 5 longevity and aging traits in up to 1345 Framingham Study participants from 330 families. Multivariable-adjusted residuals were computed using appropriate models (Cox proportional hazards, logistic, or linear regression) and the residuals from these models were used to test for association with qualifying SNPs (70, 987 autosomal SNPs with genotypic call rate [greater than or equal to]80%, minor allele frequency [greater than or equal to]10%, Hardy-Weinberg test p [greater than or equal to] 0.001).RESULTS:In family-based association test (FBAT) models, 8 SNPs in two regions approximately 500 kb apart on chromosome 1 (physical positions 73,091,610 and 73, 527,652) were associated with age at death (p-value < 10-5). The two sets of SNPs were in high linkage disequilibrium (minimum r2 = 0.58). The top 30 SNPs for generalized estimating equation (GEE) tests of association with age at death included rs10507486 (p = 0.0001) and rs4943794 (p = 0.0002), SNPs intronic to FOXO1A, a gene implicated in lifespan extension in animal models. FBAT models identified 7 SNPs and GEE models identified 9 SNPs associated with both age at death and morbidity-free survival at age 65 including rs2374983 near PON1. In the analysis of selected candidate genes, SNP associations (FBAT or GEE p-value < 0.01) were identified for age at death in or near the following genes: FOXO1A, GAPDH, KL, LEPR, PON1, PSEN1, SOD2, and WRN. Top ranked SNP associations in the GEE model for age at natural menopause included rs6910534 (p = 0.00003) near FOXO3a and rs3751591 (p = 0.00006) in CYP19A1. Results of all longevity phenotype-genotype associations for all autosomal SNPs are web posted at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007. CONCLUSION: Longevity and aging traits are associated with SNPs on the Affymetrix 100K GeneChip. None of the associations achieved genome-wide significance. These data generate hypotheses and serve as a resource for replication as more genes and biologic pathways are proposed as contributing to longevity and healthy aging.

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