Functional genomic analysis of heat stress in Vitis vinifera

Grapevine (Vitis vinifera) is one of most agro-economically important fruit crops worldwide, with a special relevance in Portugal where over 300 varieties are used for wine production. Due to global warming, temperature stress is currently a serious issue affecting crop production especially in temperate climates. Mobile genetic elements such as retrotransposons have been shown to be involved in environmental stress induced genetic and epigenetic modifications. In this study, sequences related to Grapevine Retrotransposon 1 (Gret1) were utilized to determine heat induced genomic and transcriptomic modifications in Touriga Nacional, a traditional Portuguese grapevine variety. For this purpose, growing canes were treated to 42 oC for four hours and leaf genomic DNA and RNA was utilized for various techniques to observe possible genomic alterations and variation in transcription levels of coding and non-coding sequences between non-treated plants and treated plants immediately after heat stress (HS-0 h) or after a 24 hour recovery period (HS-24 h). Heat stress was found to induce a significant decrease in Gret1 related sequences in HS-24 h leaves, indicating an effect of heat stress on genomic structure. In order to identify putative heat induced DNA modifications, genome wide approaches such as Amplified Fragment Length Polymorphism were utilized. This resulted in the identification of a polymorphic DNA fragment in HS-0 h and HS-24 h leaves whose sequence mapped to a genomic region flanking a house keeping gene (NADH) that is represented in multiple copies in the Vitis vinifera genome. Heat stress was also found to affect the transcript levels of various non-coding and gene coding sequences. Accordingly, quantitative real time PCR results established that Gret1 related sequences are up regulated immediately after heat stress whereas the level of transcript of genes involved in identification and repair of double strand breaks are significantly down regulated in HS-0 h plants. Taken together, the results of this work demonstrated heat stress affects both genomic integrity and transcription levels.

RNA interference: a new and powerful tool for functional genomic analysis.

In many species, the introduction of double-stranded RNA induces potent and specific gene silencing, referred to as RNA interference. This phenomenon, which is based on targeted degradation of mRNAs and occurs in almost any eukaryote, from trypanosomes to mice including plants and fungi, has sparked general interest from both applied and fundamental standpoints. RNA interference, which is currently used to investigate gene function in a variety of systems, is linked to natural resistance to viruses and transposon silencing, as if it were a primitive immune system involved in genome surveillance. Here, we review the mechanism of RNA interference in post-transcriptional gene silencing, its function in nature, its value for functional genomic analysis, and the modifications and improvements that may make it more efficient and inheritable. We also discuss the future directions of this versatile technique in both fundamental and applied science.

Functional Genomic Approaches for Deciphering Plant-Virus Interactions

Plant-virus interactions are very complex in nature and lead to disease and symptom formation by causing various physiological, metabolic and developmental changes in the host plants. These interactions are mainly the outcomes of viral hijacking of host components to complete their infection cycles and of host defensive responses to restrict the viral infections. Viral genomes contain only a small number of genes often encoding for multifunctional proteins, and all are essential in establishing a viral infection. Thus, it is important to understand the specific roles of individual viral genes and their contribution to the viral life cycles. Among the most important viral proteins are the suppressors of RNA silencing (VSRs). These proteins function to suppress host defenses mediated by RNA silencing and can also serve in other functions, e.g. in viral movement, transactivation of host genes, virus replication and protein processing. Thus these proteins are likely to have a significant impact on host physiology and metabolism. In the present study, I have examined the plant-virus interactions and the effects of three different VSRs on host physiology and gene expression levels by microarray analysis of transgenic plants that express these VSR genes. I also studied the gene expression changes related to the expression of the whole genome of Tobacco mosaic virus (TMV) in transgenic tobacco plants. Expression of the VSR genes in the transgenic tobacco plants causes significant changes in the gene expression profiles. HC-Pro gene derived from the Potyvirus Y (PVY) causes alteration of 748 and 332 transcripts, AC2 gene derived from the African cassava mosaic virus (ACMV) causes alteration of 1118 and 251transcripts, and P25 gene derived from the Potyvirus X (PVX) causes alterations of 1355 and 64 transcripts in leaves and flowers, respectively. All three VSRs cause similar up-regulation in defense, hormonally regulated and different stress-related genes and down-regulation in the photosynthesis and starch metabolism related genes. They also induce alterations that are specific to each viral VSR. The phenotype and transcriptome alterations of the HC-Pro expressing transgenic plants are similar to those observed in some Potyvirus-infected plants. The plants show increased protein degradation, which may be due to the HC-Pro cysteine endopeptidase and thioredoxin activities. The AC2-expressing transgenic plants show a similar phenotype and gene expression pattern as HC-Pro-expressing plants, but also alter pathways related to jasmonic acid, ethylene and retrograde signaling. In the P25 expressing transgenic plants, high numbers of genes (total of 1355) were up-regulated in the leaves, compared to a very low number of down-regulated genes (total of 5). Despite of strong induction of the transcripts, only mild growth reduction and no other distinct phenotype was observed in these plants. As an example of whole virus interactions with its host, I also studied gene expression changes caused by Tobacco mosaic virus (TMV) in tobacco host in three different conditions, i.e. in transgenic plants that are first resistant to the virus, and then become susceptible to it and in wild type plants naturally infected with this virus. The microarray analysis revealed up and down-regulation of 1362 and 1422 transcripts in the TMV resistant young transgenic plants, and up and down-regulation of a total of 1150 and 1200 transcripts, respectively, in the older plants, after the resistance break. Natural TMV infections in wild type plants caused up-regulation of 550 transcripts and down-regulation of 480 transcripts. 124 up-regulated and 29 down-regulated transcripts were commonly altered between young and old TMV transgenic plants, and only 6 up-regulated and none of the down-regulated transcripts were commonly altered in all three plants. During the resistant stage, the strong down-regulation in translation-related transcripts (total of 750 genes) was observed. Additionally, transcripts related to the hormones, protein degradation and defense pathways, cell division and stress were distinctly altered. All these alterations may contribute to the TMV resistance in the young transgenic plants, and the resistance may also be related to RNA silencing, despite of the low viral abundance and lack of viral siRNAs or TMV methylation activity in the plants.

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.

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

SAGE transcript profiling of the juvenile cambial region of Eucalyptus grandis

Despite the importance of Eucalyptus spp. in the pulp and paper industry, functional genomic approaches have only recently been applied to understand wood formation in this genus. We attempted to establish a global view of gene expression in the juvenile cambial region of Eucalyptus grandis Hill ex Maiden. The expression profile was obtained from serial analysis of gene expression (SAGE) library data produced from 3- and 6-year-old trees. Fourteen-base expressed sequence tags (ESTs) were searched against public Eucalyptus ESTs and annotated with GenBank. Altogether 43,304 tags were generated producing 3066 unigenes with three or more copies each, 445 with a putative identity, 215 with unknown function and 2406 without an EST match. The expression profile of the juvenile cambial region revealed the presence of highly frequent transcripts related to general metabolism and energy metabolism, cellular processes, transport, structural components and information pathways. We made a quantitative analysis of a large number of genes involved in the biosynthesis of cellulose, pectin, hemicellulose and lignin. Our findings provide insight into the expression of functionally related genes involved in juvenile wood formation in young fast-growing E. grandis trees.

Silencing of HTLV-1 gag and env genes by small interfering RNAs in HEK 293 cells

Since the discovery of RNAi technology, several functional genomic and disease therapy studies have been conducted using this technique in the field of oncology and virology. RNAi-based antiviral therapies are being studied for the treatment of retroviruses such as HIV-1. These studies include the silencing of regulatory, infectivity and structural genes. The HTLV-1 structural genes are responsible for the synthesis of proteins involved in the entry, assembly and release of particles during viral infection. To examine the possibility of silencing HTLV-1 genes gag and env by RNA interference technology, these genes were cloned into reporter plasmids. These vectors expressed the target mRNAs fused to EGFP reporter genes. Three small interference RNAs (siRNAs) corresponding to gag and three corresponding to env were designed to analyze the effect of silencing by RNAi technology. The plasmids and siRNAs were co-transfected into HEK 293 cells. The results demonstrated that the expression of the HTLV-1 gag and env genes decreased significantly in vitro. Thus, siRNAs can be used to inhibit HTLV-1 structural genes in transformed cells, which could provide a tool for clarifying the roles of HTLV-1 structural genes, as well as a therapy for this infection. (C) 2011 Elsevier B.V. All rights reserved.

Progress with schistosome transgenesis

Genome sequences for Schistosoma japonicum and Schistosoma mansoni are now available. The schistosome genome encodes ~13,000 protein encoding genes for which the function of only a minority is understood. There is a valuable role for transgenesis in functional genomic investigations of these new schistosome gene sequences. In gain-of-function approaches, transgenesis can lead to integration of transgenes into the schistosome genome which can facilitate insertional mutagenesis screens. By contrast, transgene driven, vector-based RNA interference (RNAi) offers powerful loss-of-function manipulations. Our laboratory has focused on development of tools to facilitate schistosome transgenesis. We have investigated the utility of retroviruses and transposons to transduce schistosomes. Vesicular stomatitis virus glycoprotein (VSVG) pseudotyped murine leukemia virus (MLV) can transduce developmental stages of S. mansoni including eggs. We have also observed that the piggyBac transposon is transpositionally active in schistosomes. Approaches with both VSVG-MLV and piggyBac have resulted in somatic transgenesis and have lead to integration of active reporter transgenes into schistosome chromosomes. These findings provided the first reports of integration of reporter transgenes into schistosome chromosomes. Experience with these systems is reviewed herewith, along with findings with transgene mediated RNAi and germ line transgenesis, in addition to pioneering and earlier reports of gene manipulation for schistosomes.

zyg-11and cul-2 promote the metaphase to anaphase transition and M phase exit of meiosis II and prevent polarity establishment in C.elegans

Résumé Les mécanismes qui coordonnent la progression du cycle cellulaire lors de la méiose avec les événements du développement embryonnaire précoce, y compris la formation des axes de polarité embryonnaire, sont peu compris. Dans le zygote du vers Caenorhabditis elegans, les premiers signes de polarité Antéro-Postérieur (A-P) embryonnaire apparaissent après que la méiose soit terminée. La nature des protéines et des mécanismes moléculaires qui cassent la symétrie du zygote n'est pas connue. Nous démontrons que zyg-11 et cul-2 promeuvent la transition métaphase - anaphase et la sortie de la phase M lors de la seconde division méiotique. Nos résultats indiquent que ZYG-11 agit comme unité recrutant le substrat d'une ligase E3 comprennant CUL-2. Nos résultats montrent aussi que le délai de sortie de la phase M dépend de l'accumulation de la Cyclin B, CYB-3. Nous démontrons que dans des embryons zyg-11(RNAi) ou cul-2(RNAi), une polarité inversée est établie lors du délai de méiosis II. Enfin nous montrons que les défauts de cycle cellulaire et ceux de polarité peuvent être séparés. De plus, nous faisons apparaitre que l'établissement d'une polarité inversée pendant le délai de méiose II des embryons zyg-11(RNAi), comme l'établissement de la A-P polarité des embryons sauvage ne semblent pas requérir les microtubules. Nous montrons également les premiers résultats d'un crible deux hybrides ainsi qu'un crible génomique qui vise à identifier des gènes dont l'inactivation augmente ou supprime les défauts de mutants pour le gène zyg-11, afin d'identifier les gènes qui intéragissent avec ZYG-11 pour assumer ses deux fonctions séparables. Par conséquent, nos trouvailles suggèrent un modèle selon lequel ZYG-11 est une sous-unité qui recrute les substrats d'une ligase E3 basée sur CUL-2 qui promeut la progression du cycle cellulaire et empêche l'établissement de la polarité pendant la méiose II, et où le centrosome agit comme la clé qui polarise l'embryon à la fin de la méiose. Summary The mechanisms that couple meiotic cell cycle progression to subsequent developmental events, including specification of embryonic axes, are poorly understood. In the one cell stage embryos of Caenorhabditis elegans, the first signs of Antero-Posterior (A-P) polarity appear after meiosis completion. A centrosome ¬derived component breaks symmetry of the embryo, but the molecular nature of this polarity signal is not known. We established that zyg-11 and cul-2 promote the metaphase to anaphase transition and M phase exit at meiosis II. Our results indicate that ZYG-11 acts as a substrate recruitment subunit of a CUL-2-based E3 ligase. Moreover, we find that the delayed meiosis II exit of embryos lacking zyg-11 is caused by accumulation of the B-type cyclin, CYB-3. We demonstrate that inverted A-P polarity is established during the meiosis II delay in zyg-11(RNAi) and cul¬2(RNAi) embryos. Importantly, we demonstrate that the polarity defects following zyg-11 or cul-2 inactivation can be uncoupled from the cell cycle defects. Furthermore, we found that microtubules appear dispensable for inverted polarity during the meiosis II delay in zyg-11(RNAi) embryos, as well as for A-P polarity during the first mitotic cell cycle in wild-type embryos. We also show the initial results from a comprehensive yeast two hybrid, as well as an RNAi-based functional genomic enhancer and suppressor screen, that may lead to identification of proteins that interact with zyg-11 to ensure the two functions. Our findings suggest a model in which ZYG-11 is a substrate recruitment subunit of an CUL-2-based E3 ligase that promotes cell cycle progression and prevents polarity establishment during meiosis II, and in which the centrosome acts as a cue to polarize the embryo along the AP axis after exit from the meiotic cell cycle.