Hierarchical Multi-classification with Predictive Clustering Trees in Functional Genomics

Struyf, J., Dzeroski, S. Blockeel, H. and Clare, A. (2005) Hierarchical Multi-classification with Predictive Clustering Trees in Functional Genomics. In proceedings of the EPIA 2005 CMB Workshop

Brachypodium distachyon. A new model system for functional genomics in grasses

John Draper, Luis A.J. Mur, Glyn Jenkins, Gadab C. Ghosh-Biswas, Pauline Bablak, Robert Hasterok,and Andrew P.M. Routledge (2001). Brachypodium distachyon. A new model system for functional genomics in grasses. Plant Physiology, 127 (4), 1539-1555. Sponsorship: BBSRC / Gatsby Foundation RAE2008

Functional genomics of motile commensal intestinal bacteria

Flagella confer upon bacteria the ability to move and are therefore organelles of significant bacteriological importance. The innate immune system has evolved to recognise flagellin, (the major protein component of the bacterial flagellar filament). Flagellate microbes can potentially stimulate the immune systems of mammals, and thus have significant immunomodulatory potential. The flagellum-biogenesis genotype and phenotype of Lactobacillus ruminis, an autochthonous intestinal commensal, was studied. The flagellum-biogenesis genotypes of motile enteric Eubacterium and Roseburia species were also investigated. Flagellin proteins were recovered from these commensal species, their amino-termini were sequenced and the proteins were found to be pro-inflammatory, as assessed by measurement of interleukin-8 (IL-8) secretion from human intestinal epithelial cell lines. For L. ruminis, this IL-8 secretion required signalling through Toll Like Receptor 5. A model for the regulation of flagellum-biogenesis in L. ruminis was inferred from transcriptomics data and bioinformatics analyses. Motility gene expression in this species may be under the control of a novel regulator, LRC_15730. Potential promoters for genes encoding flagellin proteins in the Eubacterium and Roseburia genomes analysed were inferred in silico. Relative abundances of the target Eubacterium and Roseburia species in the intestinal microbiota of 25 elderly individuals were determined. These species were found to be variably abundant in these individuals. Motility genes from these species were variably detected in the shotgun metagenome databases generated by the ELDERMET project. This suggested that a greater depth of sequencing, or improved evenness of sequencing, would be required to capture the full diversity of microbial functions for specific target or low abundance species in microbial communities by metagenomics. In summary, this thesis used a functional genomics approach to describe flagellum-mediated motility in selected Gram-positive commensal bacteria. The regulation of flagellum biosynthesis in these species, and the consequences of flagella expression from a host-interaction perspective were also considered.

Comparative genomics of Bifidobacteria

Chapter 2 of this thesis describes the sequence analysis of 14 bifidobacterial genomes from various species of the genus Bifidobacterium, and the determination of their open pan-genome trend. This analysis first determined the total number of genes to be considered as the reservoir of functions available to representatives of this genus. Many identified genes are still uncharacterized, but may be involved in the adaptation to the gut environment. This comparative genomic analysis also determined a pool of ortholog functions used to infer their phylogenetic relationship, thereby providing a more robust approach compared to that based solely on the16S rRNA-encoding gene. The genome sequence of an isolate from the insect hindgut Bifidobacterium asteroides PRL2011 was fully characterized in Chapter 3, surprisingly revealing a putative respiratory metabolism, which was also found to be present in other insect isolates, suggesting that respiration was an ancient feature of this genus, but also an adaptative trait to different atmosferic oxygen levels. Chapter 4 of this thesis outlines a comparative study which focused on the analysis of representatives of the Bifidobacterium breve species, revealing that the genetic variability among members of this species principally consists of genes with a role in adaptation to host environment and gut colonization. Finally, Chapter 5 describes the analysis of the genome sequence of Bifidobacterium animalis subsp. lactis BLC1, a probiotic bacterium widely used in food industries as an ingredient of functional foods, providing information that will allow future investigations of this species.

Functional genomics of commensal Lactobacilli

Catabolic flexibility affords a bacterium the ability to utilise different sugar sources as carbon for energy. This is important for commensal lactobacilli like Lactobacillus ruminis which can be exposed to a variety of carbohydrates in vivo. However, little is known about the fermentation capabilities, metabolic pathways, genetic diversity or potential survival mechanisms used by L. ruminis in vivo. A combination of in vitro and in silico techniques was used to identify the catabolic pathways of L. ruminis. I also compared 16 L. ruminis strains using a panel of biochemical and survival assays, genetically, whole genome sequencing and RNA sequencing. Multi locus sequence typing revealed that strains clustered according to their host sources. Transcriptome analysis by RNAseq of two motile strains under three growth conditions, including swarming, identified the up-regulation of carbohydrate-related genes under swarming conditions. This suggests that carbohydrate flexibility may have an uncharacterised role in L. ruminis swarming. Following on from the assessment of L. ruminis catabolic flexibility, the porcine diet was supplemented with galactooligosaccharides or L. ruminis ATCC 25644 plus galactooligosaccharides. Supplementation of the porcine diet with galactooligosaccharide had no effect on microbiota diversity. In contrast, the L. ruminis plus galactooligosaccharide treatment significantly reduced the microbiota diversity. Diet is a major factor that affects the diversity of the gut microbiota. In order to get a more thorough understanding of diet and gut health in animals such as racehorses and domesticated herbivores, I determined the core microbiota of animals consuming different feeds. Interestingly, the gut microbiota diversity correlated with the host phylogeny of the animal. The genome of Lactobacillus equi (2.19 Mb), isolated from a healthy Irish thoroughbred was also sequenced and annotated, and comprised 2,263 predicted genes. The large repertoire of predicted carbohydrate-related genes may offer L. equi an advantage in the complex and harsh hindgut environment. In summary, this thesis uses functional genomics to assess the effect that carbohydrates have on commensal lactobacilli and the microbiota as a whole.

Helicobacter pylori: comparative genomics and structure-function analysis of the flagellum biogenesis protein HP0958

Helicobacter pylori is a gastric pathogen which infects ~50% of the global population and can lead to the development of gastritis, gastric and duodenal ulcers and carcinoma. Genome sequencing of H. pylori revealed high levels of genetic variability; this pathogen is known for its adaptability due to mechanisms including phase variation, recombination and horizontal gene transfer. Motility is essential for efficient colonisation by H. pylori. The flagellum is a complex nanomachine which has been studied in detail in E. coli and Salmonella. In H. pylori, key differences have been identified in the regulation of flagellum biogenesis, warranting further investigation. In this study, the genomes of two H. pylori strains (CCUG 17874 and P79) were sequenced and published as draft genome sequences. Comparative studies identified the potential role of restriction modification systems and the comB locus in transformation efficiency differences between these strains. Core genome analysis of 43 H. pylori strains including 17874 and P79 defined a more refined core genome for the species than previously published. Comparative analysis of the genome sequences of strains isolated from individuals suffering from H. pylori related diseases resulted in the identification of “disease-specific” genes. Structure-function analysis of the essential motility protein HP0958 was performed to elucidate its role during flagellum assembly in H. pylori. The previously reported HP0958-FliH interaction could not be substantiated in this study and appears to be a false positive. Site-directed mutagenesis confirmed that the coiled-coil domain of HP0958 is involved in the interaction with RpoN (74-284), while the Zn-finger domain is required for direct interaction with the full length flaA mRNA transcript. Complementation of a non-motile hp0958-null derivative strain of P79 with site-directed mutant alleles of hp0958 resulted in cells producing flagellar-type extrusions from non-polar positions. Thus, HP0958 may have a novel function in spatial localisation of flagella in H. pylori

Bayesian variable selection in structured high-dimensional covariate spaces with applications in genomics

We consider the problem of variable selection in regression modeling in high-dimensional spaces where there is known structure among the covariates. This is an unconventional variable selection problem for two reasons: (1) The dimension of the covariate space is comparable, and often much larger, than the number of subjects in the study, and (2) the covariate space is highly structured, and in some cases it is desirable to incorporate this structural information in to the model building process. We approach this problem through the Bayesian variable selection framework, where we assume that the covariates lie on an undirected graph and formulate an Ising prior on the model space for incorporating structural information. Certain computational and statistical problems arise that are unique to such high-dimensional, structured settings, the most interesting being the phenomenon of phase transitions. We propose theoretical and computational schemes to mitigate these problems. We illustrate our methods on two different graph structures: the linear chain and the regular graph of degree k. Finally, we use our methods to study a specific application in genomics: the modeling of transcription factor binding sites in DNA sequences. © 2010 American Statistical Association.

Genomics research: world survey of public funding.

BACKGROUND: Over the past two decades, genomics has evolved as a scientific research discipline. Genomics research was fueled initially by government and nonprofit funding sources, later augmented by private research and development (R&D) funding. Citizens and taxpayers of many countries have funded much of the research, and have expectations about access to the resulting information and knowledge. While access to knowledge gained from all publicly funded research is desired, access is especially important for fields that have broad social impact and stimulate public dialogue. Genomics is one such field, where public concerns are raised for reasons such as health care and insurance implications, as well as personal and ancestral identification. Thus, genomics has grown rapidly as a field, and attracts considerable interest. RESULTS: One way to study the growth of a field of research is to examine its funding. This study focuses on public funding of genomics research, identifying and collecting data from major government and nonprofit organizations around the world, and updating previous estimates of world genomics research funding, including information about geographical origins. We initially identified 89 publicly funded organizations; we requested information about each organization's funding of genomics research. Of these organizations, 48 responded and 34 reported genomics research expenditures (of those that responded but did not supply information, some did not fund such research, others could not quantify it). The figures reported here include all the largest funders and we estimate that we have accounted for most of the genomics research funding from government and nonprofit sources. CONCLUSION: Aggregate spending on genomics research from 34 funding sources averaged around $2.9 billion in 2003-2006. The United States spent more than any other country on genomics research, corresponding to 35% of the overall worldwide public funding (compared to 49% US share of public health research funding for all purposes). When adjusted to genomics funding intensity, however, the United States dropped below Ireland, the United Kingdom, and Canada, as measured both by genomics research expenditure per capita and per Gross Domestic Product.

Integrative Genomics Reveals a Role for GNA13 in Lymphomagenesis

Lymphomas comprise a diverse group of malignancies derived from immune cells. High throughput sequencing has recently emerged as a powerful and versatile method for analysis of the cancer genome and transcriptome. As these data continue to emerge, the crucial work lies in sorting through the wealth of information to hone in on the critical aspects that will give us a better understanding of biology and new insight for how to treat disease. Finding the important signals within these large data sets is one of the major challenges of next generation sequencing.

In this dissertation, I have developed several complementary strategies to describe the genetic underpinnings of lymphomas. I begin with developing a better method for RNA sequencing that enables strand-specific total RNA sequencing and alternative splicing profiling in the same analysis. I then combine this RNA sequencing technique with whole exome sequencing to better understand the global landscape of aberrations in these diseases. Finally, I use traditional cell and molecular biology techniques to define the consequences of major genetic alterations in lymphoma.

Through this analysis, I find recurrent silencing mutations in the G alpha binding protein GNA13 and associated focal adhesion proteins. I aim to describe how loss-of-function mutations in GNA13 can be oncogenic in the context of germinal center B cell biology. Using in vitro techniques including liquid chromatography-mass spectrometry and knockdown and overexpression of genes in B cell lymphoma cell lines, I determine protein binding partners and downstream effectors of GNA13. I also develop a transgenic mouse model to study the role of GNA13 in the germinal center in vivo to determine effects of GNA13 deletion on germinal center structure and cell migration.

Thus, I have developed complementary approaches that span the spectrum from discovery to context-dependent gene models that afford a better understanding of the biological function of aberrant events and ultimately result in a better understanding of disease.

The past, present, and future of human centromere genomics.

The centromere is the chromosomal locus essential for chromosome inheritance and genome stability. Human centromeres are located at repetitive alpha satellite DNA arrays that compose approximately 5% of the genome. Contiguous alpha satellite DNA sequence is absent from the assembled reference genome, limiting current understanding of centromere organization and function. Here, we review the progress in centromere genomics spanning the discovery of the sequence to its molecular characterization and the work done during the Human Genome Project era to elucidate alpha satellite structure and sequence variation. We discuss exciting recent advances in alpha satellite sequence assembly that have provided important insight into the abundance and complex organization of this sequence on human chromosomes. In light of these new findings, we offer perspectives for future studies of human centromere assembly and function.

RCN4GSC Workshop Report: Managing Data at the Interface of Biodiversity and (Meta)Genomics, March 2011.

Building on the planning efforts of the RCN4GSC project, a workshop was convened in San Diego to bring together experts from genomics and metagenomics, biodiversity, ecology, and bioinformatics with the charge to identify potential for positive interactions and progress, especially building on successes at establishing data standards by the GSC and by the biodiversity and ecological communities. Until recently, the contribution of microbial life to the biomass and biodiversity of the biosphere was largely overlooked (because it was resistant to systematic study). Now, emerging genomic and metagenomic tools are making investigation possible. Initial research findings suggest that major advances are in the offing. Although different research communities share some overlapping concepts and traditions, they differ significantly in sampling approaches, vocabularies and workflows. Likewise, their definitions of 'fitness for use' for data differ significantly, as this concept stems from the specific research questions of most importance in the different fields. Nevertheless, there is little doubt that there is much to be gained from greater coordination and integration. As a first step toward interoperability of the information systems used by the different communities, participants agreed to conduct a case study on two of the leading data standards from the two formerly disparate fields: (a) GSC's standard checklists for genomics and metagenomics and (b) TDWG's Darwin Core standard, used primarily in taxonomy and systematic biology.

Bordetella pertussis from functional genomics to intranasal vaccination.

Whooping cough still represents a major health problem, despite the use of effective vaccines for several decades. Being classically a typical childhood disease, whooping cough in young adults is now more common than it used to be, suggesting that protection after vaccination wanes during adolescence. As an alternative to the current vaccines, we wish to develop live attenuated vaccines to be delivered by the nasal route, such as to mimic the natural route of infection and to induce long lasting immunity. Bordetella pertussis, the etiological agent of whooping cough, produces a number of virulence factors, including toxins. Its recently determined genome sequence makes it now possible to apply functional genomics, such as transcriptomics and systematic knock-out mutagenesis. The expression of most known B. pertussis virulence genes is controlled by the two-component system BvgA/S. DNA microarray analyses have led to the identification of novel genes in the BvgA/S regulon, some of which are activated by BvgA/S and others are repressed by BvgA/S. In addition, some genes appear to be differentially modulated by nicotinic acid and MgSO4, both known to modulate the expression of BvgA/S-regulated genes. Among others, the functional genomics approach has uncovered two strongly BvgA/S-activated genes, named hotA and hotB (for 'homolog of toxin'), the products of which show high sequence similarities to pertussis toxin subunits. The identification of the full array of virulence factors, as well as an integrated understanding of the bacterial physiology should allow us to design attenuated B. pertussis strains useful for intranasal vaccination. A first generation of attenuated strains has already shown full protection in mice after a single intranasal administration. Such strains may also serve as vaccine carriers for heterologous antigens, in order to vaccinate against several different pathogens simultaneously.