Functional genomic studies of the structure and regulation of eukaryotic transcriptomes

The main focus of this thesis is the use of high-throughput sequencing technologies in functional genomics (in particular in the form of ChIP-seq, chromatin immunoprecipitation coupled with sequencing, and RNA-seq) and the study of the structure and regulation of transcriptomes. Some parts of it are of a more methodological nature while others describe the application of these functional genomic tools to address various biological problems. A significant part of the research presented here was conducted as part of the ENCODE (ENCyclopedia Of DNA Elements) Project.

The first part of the thesis focuses on the structure and diversity of the human transcriptome. Chapter 1 contains an analysis of the diversity of the human polyadenylated transcriptome based on RNA-seq data generated for the ENCODE Project. Chapter 2 presents a simulation-based examination of the performance of some of the most popular computational tools used to assemble and quantify transcriptomes. Chapter 3 includes a study of variation in gene expression, alternative splicing and allelic expression bias on the single-cell level and on a genome-wide scale in human lymphoblastoid cells; it also brings forward a number of critical to the practice of single-cell RNA-seq measurements methodological considerations.

The second part presents several studies applying functional genomic tools to the study of the regulatory biology of organellar genomes, primarily in mammals but also in plants. Chapter 5 contains an analysis of the occupancy of the human mitochondrial genome by TFAM, an important structural and regulatory protein in mitochondria, using ChIP-seq. In Chapter 6, the mitochondrial DNA occupancy of the TFB2M transcriptional regulator, the MTERF termination factor, and the mitochondrial RNA and DNA polymerases is characterized. Chapter 7 consists of an investigation into the curious phenomenon of the physical association of nuclear transcription factors with mitochondrial DNA, based on the diverse collections of transcription factor ChIP-seq datasets generated by the ENCODE, mouseENCODE and modENCODE consortia. In Chapter 8 this line of research is further extended to existing publicly available ChIP-seq datasets in plants and their mitochondrial and plastid genomes.

The third part is dedicated to the analytical and experimental practice of ChIP-seq. As part of the ENCODE Project, a set of metrics for assessing the quality of ChIP-seq experiments was developed, and the results of this activity are presented in Chapter 9. These metrics were later used to carry out a global analysis of ChIP-seq quality in the published literature (Chapter 10). In Chapter 11, the development and initial application of an automated robotic ChIP-seq (in which these metrics also played a major role) is presented.

The fourth part presents the results of some additional projects the author has been involved in, including the study of the role of the Piwi protein in the transcriptional regulation of transposon expression in Drosophila (Chapter 12), and the use of single-cell RNA-seq to characterize the heterogeneity of gene expression during cellular reprogramming (Chapter 13).

The last part of the thesis provides a review of the results of the ENCODE Project and the interpretation of the complexity of the biochemical activity exhibited by mammalian genomes that they have revealed (Chapters 15 and 16), an overview of the expected in the near future technical developments and their impact on the field of functional genomics (Chapter 14), and a discussion of some so far insufficiently explored research areas, the future study of which will, in the opinion of the author, provide deep insights into many fundamental but not yet completely answered questions about the transcriptional biology of eukaryotes and its regulation.

Functional genomic hypothesis generation and experimentation by a robot scientist

King R. D., Whelan, K. E., Jones, F. M., Reiser, P. G. K., Bryant, C. H., Muggleton, S., Kell, D. B. and Oliver, S. G. (2004) Functional genomic hypothesis generation and experimentation by a robot scientist. Nature 427 (6971) p247-252

Functional Genomic Screen Identifies Klebsiella pneumoniae Factors Implicated in Blocking Nuclear Factor κB (NF-κB) Signaling

Klebsiella pneumoniae is etiologic agent of community-acquired and nosocomial pneumonia. It has been shown that K. pneumoniae infections are characterized by reduced early inflammatory response. Recently our group have shown that K. pneumoniae dampens the activation of inflammatory responses by antagonizing the activation of the NF-κB canonical pathway. Our results revealed that K. pneumoniae capsule (CPS) was necessary but not sufficient to attenuate inflammation. To identify additional Klebsiella factors required to dampen inflammation, we standardized and applied a high-throughput gain-on-function screen to examine a Klebsiella transposon mutant library. We identified 114 mutants that triggered the activation of NF-κB. Two gene ontology categories accounted for half of the loci identified in the screening, that of metabolism and transport (32% of the mutants), and of enveloperelated genes (17%). Characterization of the mutants revealed that the lack of the enterobactin siderophore was linked to a reduced CPS expression which in turn underlined the NF- κB activation induced by the mutant. The lipopolysaccharide (LPS) O-polysaccharide and the pullulanase (PulA) type 2 secretion system (T2SS) are required for full effectiveness of immune evasion. Importantly, these factors do not play a redundant role. The fact that LPS Opolysaccharide and T2SS mutants-induced responses were dependent on TLR2-TLR4- MyD88 activation suggested that LPS Opolysaccharide and PulA perturbed TLRdependent recognition of K. pneumoniae. Finally, we demonstrate that LPS O-polysaccharide and pulA mutants are attenuated in the pneumonia mouse model. We propose that LPS Opolysaccharide and PulA T2SS could be new targets for designing new antimicrobials. Increasing TLR-governed defence responses might provide also selective alternatives for the management of K. pneumoniae pneumonia.

Mapping the platelet profile for functional genomic studies and demonstration of the effect size of the GP6 locus

BACKGROUND: Evidence suggests the wide variation in platelet response within the population is genetically controlled. Unraveling the complex relationship between sequence variation and platelet phenotype requires accurate and reproducible measurement of platelet response. OBJECTIVE: To develop a methodology suitable for measuring signaling pathway-specific platelet phenotype, to use this to measure platelet response in a large cohort, and to demonstrate the effect size of sequence variation in a relevant model gene. METHODS: Three established platelet assays were evaluated: mobilization of [Ca(2+)](i), aggregometry and flow cytometry, each in response to adenosine 5'-diphosphate (ADP) or the glycoprotein (GP) VI-specific crosslinked collagen-related peptide (CRP). Flow cytometric measurement of fibrinogen binding and P-selectin expression in response to a single, intermediate dose of each agonist gave the best combination of reproducibility and inter-individual variability and was used to measure the platelet response in 506 healthy volunteers. Pathway specificity was ensured by blocking the main subsidiary signaling pathways. RESULTS: Individuals were identified who were hypo- or hyper-responders for both pathways, or who had differential responses to the two agonists, or between outcomes. 89 individuals, retested three months later using the same methodology, showed high concordance between the two visits in all four assays (r(2) = 0.872, 0.868, 0.766 and 0.549); all subjects retaining their phenotype at recall. The effect of sequence variation at the GP6 locus accounted for approximately 35% of the variation in the CRP-XL response. CONCLUSION: Genotyping-phenotype association studies in a well-characterized, large cohort provides a powerful strategy to measure the effect of sequence variation in genes regulating the platelet response.

A graph theoretic approach to testing associations between disparate sources of functional genomic data

The last few years have seen the advent of high-throughput technologies to analyze various properties of the transcriptome and proteome of several organisms. The congruency of these different data sources, or lack thereof, can shed light on the mechanisms that govern cellular function. A central challenge for bioinformatics research is to develop a unified framework for combining the multiple sources of functional genomics information and testing associations between them, thus obtaining a robust and integrated view of the underlying biology. We present a graph theoretic approach to test the significance of the association between multiple disparate sources of functional genomics data by proposing two statistical tests, namely edge permutation and node label permutation tests. We demonstrate the use of the proposed tests by finding significant association between a Gene Ontology-derived "predictome" and data obtained from mRNA expression and phenotypic experiments for Saccharomyces cerevisiae. Moreover, we employ the graph theoretic framework to recast a surprising discrepancy presented in Giaever et al. (2002) between gene expression and knockout phenotype, using expression data from a different set of experiments.

Functional genomic and molecular analyses of transcripts related to glycerol synthesis and trafficking in rainbow smelt (Osmerus mordax) using cold temperature-induced models of glycerol production

The rainbow smelt (Osmerus mordax) is an anadromous teleost that produces type II antifreeze protein (AFP) and accumulates modest urea and high glycerol levels in plasma and tissues as adaptive cryoprotectant mechanisms in sub-zero temperatures. It is known that glyceroneogenesis occurs in liver via a branch in glycolysis and gluconeogenesis and is activated by low temperature; however, the precise mechanisms of glycerol synthesis and trafficking in smelt remain to be elucidated. The objective of this thesis was to provide further insight using functional genomic techniques [e.g. suppression subtractive hybridization (SSH) cDNA library construction, microarray analyses] and molecular analyses [e.g. cloning, quantitative reverse transcription - polymerase chain reaction (QPCR)]. Novel molecular mechanisms related to glyceroneogenesis were deciphered by comparing the transcript expression profiles of glycerol (cold temperature) and non-glycerol (warm temperature) accumulating hepatocytes (Chapter 2) and livers from intact smelt (Chapter 3). Briefly, glycerol synthesis can be initiated from both amino acids and carbohydrate; however carbohydrate appears to be the preferred source when it is readily available. In glycerol accumulating hepatocytes, levels of the hepatic glucose transporter (GLUT2) plummeted and transcript levels of a suite of genes (PEPCK, MDH2, AAT2, GDH and AQP9) associated with the mobilization of amino acids to fuel glycerol synthesis were all transiently higher. In contrast, in glycerol accumulating livers from intact smelt, glycerol synthesis was primarily fuelled by glycogen degradation with higher PGM and PFK (glycolysis) transcript levels. Whether initiated from amino acids or carbohydrate, there were common metabolic underpinnings. Increased PDK2 (an inhibitor of PDH) transcript levels would direct pyruvate derived from amino acids and / or DHAP derived from G6P to glycerol as opposed to oxidation via the citric acid cycle. Robust LIPL (triglyceride catabolism) transcript levels would provide free fatty acids that could be oxidized to fuel ATP synthesis. Increased cGPDH (glyceroneogenesis) transcript levels were not required for increased glycerol production, suggesting that regulation is more likely by post-translational modification. Finally, levels of a transcript potentially encoding glycerol-3-phosphatase, an enzyme not yet characterized in any vertebrate species, were transiently higher. These comparisons also led to the novel discoveries that increased G6Pase (glucose synthesis) and increased GS (glutamine synthesis) transcript levels were part of the low temperature response in smelt. Glucose may provide increased colligative protection against freezing; whereas glutamine could serve to store nitrogen released from amino acid catabolism in a non-toxic form and / or be used to synthesize urea via purine synthesis-uricolysis. Novel key aspects of cryoprotectant osmolyte (glycerol and urea) trafficking were elucidated by cloning and characterizing three aquaglyceroporin (GLP)-encoding genes from smelt at the gene and cDNA levels in Chapter 4. GLPs are integral membrane proteins that facilitate passive movement of water, glycerol and urea across cellular membranes. The highlight was the discovery that AQP10ba transcript levels always increase in posterior kidney only at low temperature. This AQP10b gene paralogue may have evolved to aid in the reabsorption of urea from the proximal tubule. This research has contributed significantly to a general understanding of the cold adaptation response in smelt, and more specifically to the development of a working scenario for the mechanisms involved in glycerol synthesis and trafficking in this species.

Rhipicephalus (Boophilus) microplus tick in vitro feeding methods for functional (dsRNA) and vaccine candidate (antibody) screening

Rhipicephalus (Boophilus) microplus (Acari: Ixodidae) ticks cause economic losses for cattle industries throughout tropical and subtropical regions of the world estimated at $US2.5 billion annually. Lack of access to efficacious long-lasting vaccination regimes and increases in tick acaricide resistance have led to the investigation of targets for the development of novel tick vaccines and treatments. In vitro tick feeding has been used for many tick species to study the effect of new acaricides on the transmission of tick-borne pathogens. Few studies have reported the use of in vitro feeding for functional genomic studies using RNA interference and/or the effect of specific anti-tick antibodies. In particular, in vitro feeding reports for the cattle tick are limited due to its relatively short hypostome. Previously published methods were further modified to broaden optimal tick sizes/weights, feeding sources including bovine and ovine serum, optimisation of commercially available blood anti-coagulant tubes, and IgG concentrations for effective antibody delivery. Ticks are fed overnight and monitored for ∼5–6 weeks to determine egg output and success of larval emergence using a humidified incubator. Lithium-heparin blood tubes provided the most reliable anti-coagulant for bovine blood feeding compared with commercial citrated (CPDA) and EDTA tubes. Although >30 mg semi-engorged ticks fed more reliably, ticks as small as 15 mg also fed to repletion to lay viable eggs. Ticks which gained less than ∼10 mg during in vitro feeding typically did not lay eggs. One mg/ml IgG from Bm86-vaccinated cattle produced a potent anti-tick effect in vitro (83% efficacy) similar to that observed in vivo. Alternatively, feeding of dsRNA targeting Bm86 did not demonstrate anti-tick effects (11% efficacy) compared with the potent effects of ubiquitin dsRNA. This study optimises R. microplus tick in vitro feeding methods which support the development of cattle tick vaccines and treatments.

Functional genomics to identify therapeutic prophylactic targets

Infectious diseases are a leading cause of global human mortality. The use of antimicrobials remains the most common strategy for treatment. However, the isolation of pathogens resistant to virtually all antimicrobials makes it urgent to develop effective therapeutics based on new targets. Here we review a new drug discovery paradigm focusing on identifying and targeting host factors important for infection as well as pathogen determinants involved in disease progression. We summarize innovative strategies which by combining bioinformatics with transcriptomics and chemical genetics have already identified host factors essential for pathogen entry, survival and replication. We describe how the discovery of RNA interference which allows loss-of-function studies has facilitated functional genomic studies in human cells. It is expected that these studies will identify targets to be used as host-directed drug therapy which, together with antimicrobials targeting microbial virulence factors, will efficiently eliminate the invading pathogen.

Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens

Background: Pseudomonas fluorescens are common soil bacteria that can improve plant health through nutrient cycling, pathogen antagonism and induction of plant defenses. The genome sequences of strains SBW25 and Pf0-1 were determined and compared to each other and with P. fluorescens Pf-5. A functional genomic in vivo expression technology (IVET) screen provided insight into genes used by P. fluorescens in its natural environment and an improved understanding of the ecological significance of diversity within this species. Results: Comparisons of three P. fluorescens genomes (SBW25, Pf0-1, Pf-5) revealed considerable divergence: 61% of genes are shared, the majority located near the replication origin. Phylogenetic and average amino acid identity analyses showed a low overall relationship. A functional screen of SBW25 defined 125 plant-induced genes including a range of functions specific to the plant environment. Orthologues of 83 of these exist in Pf0-1 and Pf-5, with 73 shared by both strains. The P. fluorescens genomes carry numerous complex repetitive DNA sequences, some resembling Miniature Inverted-repeat Transposable Elements (MITEs). In SBW25, repeat density and distribution revealed 'repeat deserts' lacking repeats, covering approximately 40% of the genome. Conclusions: P. fluorescens genomes are highly diverse. Strain-specific regions around the replication terminus suggest genome compartmentalization. The genomic heterogeneity among the three strains is reminiscent of a species complex rather than a single species. That 42% of plant-inducible genes were not shared by all strains reinforces this conclusion and shows that ecological success requires specialized and core functions. The diversity also indicates the significant size of genetic information within the Pseudomonas pan genome.

Genomic characterization of the italian wolf (Canis lupus): the genes involved in black coat colour determination and application of microarray technique for snps detection.

This study provides a comprehensive genetic overview on the endangered Italian wolf population. In particular, it focuses on two research lines. On one hand, we focalised on melanism in wolf in order to isolate a mutation related with black coat colour in canids. With several reported black individuals (an exception at European level), the Italian wolf population constituted a challenging research field posing many unanswered questions. As found in North American wolf, we reported that melanism in the Italian population is caused by a different melanocortin pathway component, the K locus, in which a beta-defensin protein acts as an alternative ligand for the Mc1r. This research project was conducted in collaboration with Prof. Gregory Barsh, Department of Genetics and Paediatrics, Stanford University. On the other hand, we performed analysis on a high number of SNPs thanks to a customized Canine microarray useful to integrate or substitute the STR markers for genotyping individuals and detecting wolf-dog hybrids. Thanks to DNA microchip technology, we obtained an impressive amount of genetic data which provides a solid base for future functional genomic studies. This study was undertaken in collaboration with Prof. Robert K. Wayne, Department of Ecology and Evolutionary Biology, University of California, Los Angeles (UCLA).

Flexible strategies for association analysis with genomic phenotypes

Thesis (Ph.D.)--University of Washington, 2016-08