Integrated functional genomics for improved manufacture of recombinant enveloped virus

Dissertation presented to obtain the Ph.D degree in Engineering and Technology Sciences-Biotechnology

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.

Comparative and functional genomics provide insights into the pathogenicity of dermatophytic fungi.

ABSTRACT: BACKGROUND: Millions of humans and animals suffer from superficial infections caused by a group of highly specialized filamentous fungi, the dermatophytes, which exclusively infect keratinized host structures. To provide broad insights into the molecular basis of the pathogenicity-associated traits, we report the first genome sequences of two closely phylogenetically related dermatophytes, Arthroderma benhamiae and Trichophyton verrucosum, both of which induce highly inflammatory infections in humans. RESULTS: 97% of the 22.5 megabase genome sequences of A. benhamiae and T. verrucosum are unambiguously alignable and collinear. To unravel dermatophyte-specific virulence-associated traits, we compared sets of potentially pathogenicity-associated proteins, such as secreted proteases and enzymes involved in secondary metabolite production, with those of closely related onygenales (Coccidioides species) and the mould Aspergillus fumigatus. The comparisons revealed expansion of several gene families in dermatophytes and disclosed the peculiarities of the dermatophyte secondary metabolite gene sets. Secretion of proteases and other hydrolytic enzymes by A. benhamiae was proven experimentally by a global secretome analysis during keratin degradation. Molecular insights into the interaction of A. benhamiae with human keratinocytes were obtained for the first time by global transcriptome profiling. Given that A. benhamiae is able to undergo mating, a detailed comparison of the genomes further unraveled the genetic basis of sexual reproduction in this species. CONCLUSIONS: Our results enlighten the genetic basis of fundamental and putatively virulence-related traits of dermatophytes, advancing future research on these medically important pathogens.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Functional genomics of intracellular bacteria.

During the genomic era, a large amount of whole-genome sequences accumulated, which identified many hypothetical proteins of unknown function. Rapidly, functional genomics, which is the research domain that assign a function to a given gene product, has thus been developed. Functional genomics of intracellular pathogenic bacteria exhibit specific peculiarities due to the fastidious growth of most of these intracellular micro-organisms, due to the close interaction with the host cell, due to the risk of contamination of experiments with host cell proteins and, for some strict intracellular bacteria such as Chlamydia, due to the absence of simple genetic system to manipulate the bacterial genome. To identify virulence factors of intracellular pathogenic bacteria, functional genomics often rely on bioinformatic analyses compared with model organisms such as Escherichia coli and Bacillus subtilis. The use of heterologous expression is another common approach. Given the intracellular lifestyle and the many effectors that are used by the intracellular bacteria to corrupt host cell functions, functional genomics is also often targeting the identification of new effectors such as those of the T4SS of Brucella and Legionella.

Revisiting sodium and water reabsorption with functional genomics tools.

PURPOSE OF REVIEW: The kidney plays an essential role in maintaining sodium and water balance, thereby controlling the volume and osmolarity of the extracellular body fluids, the blood volume and the blood pressure. The final adjustment of sodium and water reabsorption in the kidney takes place in cells of the distal part of the nephron in which a set of apical and basolateral transporters participate in vectorial sodium and water transport from the tubular lumen to the interstitium and, finally, to the general circulation. According to a current model, the activity and/or cell-surface expression of these transporters is/are under the control of a gene network composed of the hormonally regulated, as well as constitutively expressed, genes. It is proposed that this gene network may include new candidate genes for salt- and water-losing syndromes and for salt-sensitive hypertension. A new generation of functional genomics techniques have recently been applied to the characterization of this gene network. The purpose of this review is to summarize these studies and to discuss the potential of the different techniques for characterization of the renal transcriptome. RECENT FINDINGS: Recently, DNA microarrays and serial analysis of gene expression have been applied to characterize the kidney transcriptome in different in-vivo and in-vitro models. In these studies, a set of new interesting genes potentially involved in the regulation of sodium and water reabsorption by the kidney have been identified and are currently under detailed investigation. SUMMARY: Characterization of the kidney transcriptome is greatly expanding our knowledge of the gene networks involved in multiple kidney functions, including the maintenance of sodium and water homeostasis.

Functional genomics of O-glucosyltransferases from Concord grape (Vitis labrusca)

Grape (Vitis spp.) is a culturally and economically important crop plant that has been cultivated for thousands of years, primarily for the production of wine. Grape berries accumulate a myriad of phenylpropanoid secondary metabolites, many of which are glucosylated in plantae More than 90 O-glucosyltransferases have been cloned and biochemically characterized from plants, only two of which have been isolated from Vitis spp. The world-wide economic importance of grapes as a crop plant, the human health benefits associated with increased consumption of grape-derived metabolites, the biological relevance of glucosylation, and the lack of information about Vitis glucosyltransferases has inspired the identification, cloning and biochemical characterization of five novel "family 1" O-glucosyltransferases from Concord grape (Vitis labrusca cv. Concord). Protein purification and associated protein sequencIng led to the molecular cloning of UDP-glucose: resveratrollhydroxycinnamic acid O-glucosyltransferase (VLRSGT) from Vitis labrusca berry mesocarp tissue. In addition to being the first glucosyltransferase which accepts trans-resveratrol as a substrate to be characterized in vitro, the recombinant VLRSGT preferentially produces the glucose esters of hydroxycinnamic acids at pH 6.0, and the glucosides of trans-resveratrol and flavonols at 'pH 9.0; the first demonstration of pH-dependent bifunctional glucosylation for this class of enzymes. Gene expression and metabolite profiling support a role for this enzyme in the bifuncitonal glucosylation ofstilbenes and hydroxycinnamic acids in plantae A homology-based approach to cloning was used to identify three enzymes from the Vitis vinifera TIGR grape gene index which had high levels of protein sequence iii identity to previously characterized UDP-glucose: anthocyanin 5-0-glucosyltransferases. Molecular cloning and biochemical characterization demonstrated that these enzymes (rVLOGTl, rVLOGT2, rVLOGT3) glucosylate the 7-0-position of flavonols and the xenobiotic 2,4,5-trichlorophenol (TCP), but not anthocyanins. Variable gene expression throughout grape berry development and enzyme assays with native grape berry protein are consistent with a role for these enzymes in the glucosylation of flavonols; while the broad substrate specificity, the ability of these enzymes to glucosylate TCP and expression of these genes in tissues which are subject to pathogen attack (berry, flower, bud) is consistent with a role for these genes in the plant defense response. Additionally, the Vitis labrusca UDP-glucose: flavonoid 3-0-glucosyltransferase (VL3GT) was identified, cloned and characterized. VL3GT has 96 % protein sequence identity to the previously characterized Vitis vinifera flavonoid 3-0-glucosyltransferase (VV3GT); and glucosylates the 3-0-position of anthocyanidins and flavonols in vitro. Despite high levels of protein sequence identity, VL3GT has distinct biochemical characteristics (as compared to VV3GT), including a preference for B-ring methylated flavonoids and the inability to use UDP-galactose as a donor substrate. RT-PCR analysis of VL3GT gene expression and enzyme assays with native grape protein is consistent with an in planta role for this enzyme in the glucosylation of anthocyanidins,but not flavonols. These studies reveal the power of combining several biochemistry- and molecular biology-based tools to identify, clone, biochemically characterize and elucidate the in planta function of several biologically relevant O-glucosyltransferases from Vitis spp.

Functional genomics of BAHD acyltransferases in different varieties of the wine grape, Vilis vinifera

The characteristic "foxy" aroma of Vilis labrusca Concord grapes is due in large part to methyl anthranilate, a volatile ester formed by the enzyme anthraniloyl- CoA:methanol anthraniloyltransferase (VIAMAT) of the superfamily of BARD acyltransferases. The publication of the genome ofthe closely related wine grape Vilis vinifera, which does not accumulate methyl anthranilate, permitted the searching for any putative VlAU4T-like genes, with the result of 5 highly homologous candidates being found, with one candidate sharing 95% identity to VlAU4T. Probing the gene expression of 18 different cultivars of V. vinifora ripe berries by RT -PCR showed that many varieties do indeed express VlAU4T-like genes. Subsequent cloning of the full-length open reading frame of one of these genes from eDNA prepared from the cultivar Sauvignon Blanc permitted preliminary biochemical characterization of the enzyme after heterologous expression in E. coli. It was determined that this alcohol acyltransferase (named VvsbAATl) catalyzes the formation of cis-3-hexenyl acetate, a "green-leaf' volatile. Although the cloned gene from Sauvignon Blanc had 95% identity at the amino acid level to VIAMAT, it displayed an altered substrate specificity and expression pattern. These results highlight the difficulty in predicting substrate specificity and function of enzymes through the basis of sequence homology, which is a common finding in the study of BARD acyltransferases. Also, the determination of function of VvsbAATl and other BARD acyltransferases in V. vinifera could be used as a genetic marker for certain aroma characteristics in grape breeding programs.

Functional genomics of the Drosophila melanogaster X chromosome and the role of DWnt5 during development

The collection of X chromosome insertions (PX) lethal lines, which was isolated from a screen for essential genes on the X chromosome, was characterized by means of cloning the insertion sites, mapping the sites within genomic DNA and determination of the associated reporter gene expresssion patterns. The established STS flanking the P element insertion sites were submitted to EMBL nucleotide databases and their in situ data together with the enhancer trap expression patterns have been deposited in the FlyView database. The characterized lines are now available to be used by the scientific community for a detailed analysis of the newly established lethal gene functions. One of the isolated genes on the X chromosome was the Drosophila gene Wnt5 (DWnt5). From two independent screens, one lethal and three homozygous viable alleles were recovered, allowing the identification of two distinct functions for DWnt5 in the fly. Observations on the developing nervous system of mutant embryos suggest that DWnt5 activity affects axon projection pattern. Elevated levels of DWNT5 activity in the midline cells of the central nervous system causes improper establishment and maintenance of the axonal pathways. Our analysis of the expression and mutant phenotype indicates that DWnt5 function in a process needed for proper organization of the nervous system. A second and novel function of DWnt5 is the control of the body size by regulation of the cell number rather than affecting the size of cells. Moreover, experimentally increased DWnt5 levels in a post-mitotic region of the eye imaginal disc causes abnormal cell cycle progression, resulting in additional ommatidia in the adult eye when compared to wild type. The increased cell number and the effects on the cell cycle after exposure to high DWNT5 levels is the result of a failure to downregulate cyclin B and therefore the unsuccessful establishment of a G1 arrest.

A functional genomics approach reveals novel quantitative trait loci associated with platelet signalling pathways

Platelet response to activation varies widely between individuals but shows interindividual consistency and strong heritability. The genetic basis of this variation has not been properly explored. We therefore systematically measured the effect on function of sequence variation in 97 candidate genes in the collagen and adenosine-diphosphate (ADP) signaling pathways. Resequencing of the genes in 48 European DNA samples nearly doubled the number of known single nucleotide polymorphisms (SNPs) and informed the selection of 1327 SNPs for genotyping in 500 healthy Northern European subjects with known platelet responses to collagen-related peptide (CRP-XL) and ADP. This identified 17 novel associations with platelet function (P < .005) accounting for approximately 46% of the variation in response. Further investigations with platelets of known genotype explored the mechanisms behind some of the associations. SNPs in PEAR1 associated with increased platelet response to CRP-XL and increased PEAR1 protein expression after platelet degranulation. The minor allele of a 3' untranslated region (UTR) SNP (rs2769668) in VAV3 was associated with higher protein expression (P = .03) and increased P-selectin exposure after ADP activation (P = .004). Furthermore the minor allele of the intronic SNP rs17786144 in ITPR1 modified Ca2+ levels after activation with ADP (P < .004). These data provide novel insights into key hubs within platelet signaling networks.

Functional genomics in zebrafish permits rapid characterization of novel platelet membrane proteins

In this study, we demonstrate the suitability of the vertebrate Danio rerio (zebrafish) for functional screening of novel platelet genes in vivo by reverse genetics. Comparative transcript analysis of platelets and their precursor cell, the megakaryocyte, together with nucleated blood cell elements, endothelial cells, and erythroblasts, identified novel platelet membrane proteins with hitherto unknown roles in thrombus formation. We determined the phenotype induced by antisense morpholino oligonucleotide (MO)–based knockdown of 5 of these genes in a laser-induced arterial thrombosis model. To validate the model, the genes for platelet glycoprotein (GP) IIb and the coagulation protein factor VIII were targeted. MO-injected fish showed normal thrombus initiation but severely impaired thrombus growth, consistent with the mouse knockout phenotypes, and concomitant knockdown of both resulted in spontaneous bleeding. Knockdown of 4 of the 5 novel platelet proteins altered arterial thrombosis, as demonstrated by modified kinetics of thrombus initiation and/or development. We identified a putative role for BAMBI and LRRC32 in promotion and DCBLD2 and ESAM in inhibition of thrombus formation. We conclude that phenotypic analysis of MO-injected zebrafish is a fast and powerful method for initial screening of novel platelet proteins for function in thrombosis.

Coupling Virus-Induced Gene Silencing to Exogenous Green Fluorescence Protein Expression Provides a Highly Efficient System for Functional Genomics in Arabidopsis and across All Stages of Tomato Fruit Development

Since the advent of the postgenomic era, efforts have focused on the development of rapid strategies for annotating plant genes of unknown function. Given its simplicity and rapidity, virus-induced gene silencing (VIGS) has become one of the preeminent approaches for functional analyses. However, several problems remain intrinsic to the use of such a strategy in the study of both metabolic and developmental processes. The most prominent of these is the commonly observed phenomenon of ""sectoring"" the tissue regions that are not effectively targeted by VIGS. To better discriminate these sectors, an effective marker system displaying minimal secondary effects is a prerequisite. Utilizing a VIGS system based on the tobacco rattle virus vector, we here studied the effect of silencing the endogenous phytoene desaturase gene (pds) and the expression and subsequent silencing of the exogenous green fluorescence protein (gfp) on the metabolism of Arabidopsis (Arabidopsis thaliana) leaves and tomato (Solanum lycopersicum) fruits. In leaves, we observed dramatic effects on primary carbon and pigment metabolism associated with the photobleached phenotype following the silencing of the endogenous pds gene. However, relatively few pleiotropic effects on carbon metabolism were observed in tomato fruits when pds expression was inhibited. VIGS coupled to gfp constitutive expression revealed no significant metabolic alterations after triggering of silencing in Arabidopsis leaves and a mild effect in mature green tomato fruits. By contrast, a wider impact on metabolism was observed in ripe fruits. Silencing experiments with an endogenous target gene of interest clearly demonstrated the feasibility of cosilencing in this system; however, carefully constructed control experiments are a prerequisite to prevent erroneous interpretation.

Functional Genomics and Cell Biology of the Dolphin (Tursiops runcatus): Establishment of Novel Molecular Tools to Study Marine Mammals in Changing Environments

The dolphin (Tursiops truncatus) is a mammal that is adapted to life in a totally aquatic environment. Despite the popularity and even iconic status of the dolphin, our knowledge of its physiology, its unique adaptations and the effects on it of environmental stressors are limited. One approach to improve this limited understanding is the implementation of established cellular and molecular methods to provide sensitive and insightful information for dolphin biology. We initiated our studies with the analysis of wild dolphin peripheral blood leukocytes, which have the potential to be informative of the animal’s global immune status. Transcriptomic profiles from almost 200 individual samples were analyzed using a newly developed species-specific microarray to assess its value as a prognostic and diagnostic tool. Functional genomics analyses were informative of stress-induced gene expression profiles and also of geographical location specific transcriptomic signatures, determined by the interaction of genetic, disease and environmental factors. We have developed quantitative metrics to unambiguously characterize the phenotypic properties of dolphin cells in culture. These quantitative metrics can provide identifiable characteristics and baseline data which will enable identification of changes in the cells due to time in culture. We have also developed a novel protocol to isolate primary cultures from cryopreserved tissue of stranded marine mammals, establishing a tissue (and cell) biorepository, a new approach that can provide a solution to the limited availability of samples. The work presented represents the development and application of tools for the study of the biology, health and physiology of the dolphin, and establishes their relevance for future studies of the impact on the dolphin of environmental infection and stress.