Shaping the Dicot Fruit: Molecular and Genomic Approaches to Fruit Development

The fruit is one of the most complex and important structures produced by flowering plants, and understanding the development and maturation process of fruits in different angiosperm species with diverse fruit structures is of immense interest. In the work presented here, molecular genetics and genomic analysis are used to explore the processes that form the fruit in two species: The model organism Arabidopsis and the diploid strawberry Fragaria vesca. One important basic question concerns the molecular genetic basis of fruit patterning. A long-standing model of Arabidopsis fruit (the gynoecium) patterning holds that auxin produced at the apex diffuses downward, forming a gradient that provides apical-basal positional information to specify different tissue types along the gynoecium’s length. The proposed gradient, however, has never been observed and the model appears inconsistent with a number of observations. I present a new, alternative model, wherein auxin acts to establish the adaxial-abaxial domains of the carpel primordia, which then ensures proper development of the final gynoecium. A second project utilizes genomics to identify genes that regulate fruit color by analyzing the genome sequences of Fragaria vesca, a species of wild strawberry. Shared and distinct SNPs among three F. vesca accessions were identified, providing a foundation for locating candidate mutations underlying phenotypic variations among different F. vesca accessions. Through systematic analysis of relevant SNP variants, a candidate SNP in FveMYB10 was identified that may underlie the fruit color in the yellow-fruited accessions, which was subsequently confirmed by functional assays. Our lab has previously generated extensive RNA-sequencing data that depict genome-scale gene expression profiles in F. vesca fruit and flower tissues at different developmental stages. To enhance the accessibility of this dataset, the web-based eFP software was adapted for this dataset, allowing visualization of gene expression in any tissues by user-initiated queries. Together, this thesis work proposes a well-supported new model of fruit patterning in Arabidopsis and provides further resources for F. vesca, including genome-wide variant lists and the ability to visualize gene expression. This work will facilitate future work linking traits of economic importance to specific genes and gaining novel insights into fruit patterning and development.

CHARACTERIZATION OF TOBACCO MOSAIC VIRUS DIRECTED REPROGRAMMING OF PHLOEM GENE EXPRESSION TO PROMOTE VIRAL PHLOEM LOADING AND SYSTMEIC MOVEMENT

Vascular phloem loading has long been recognized as an essential step in the establishment of a systemic virus infection. Yet little is known about this process and the mechanisms that control it. In this study, an interaction between the replication protein of Tobacco mosaic virus (TMV) and phloem specific auxin/indole acetic acid (Aux/IAA) transcriptional regulators was found to modulate virus phloem loading. Promoter expression studies show TMV 126/183 kDa interacting Aux/IAAs predominantly express and accumulate within the nuclei of phloem companion cells (CC). Furthermore, CC Aux/IAA nuclear localization is disrupted upon infection with an interacting virus but not during infection with a non-interacting virus. In situ analysis of virus spread shows the inability of TMV variants to disrupt Aux/IAA CC nuclear localization correlates with a reduced ability to load into the vascular tissue. Subsequent systemic movement assays also demonstrate that a virus capable of disrupting Aux/IAA localization is significantly more competitive at systemic movement than a non-interacting virus. Similarly, CC expression and over-accumulation of a degradation-resistant-interacting Aux/IAA protein was found to selectively inhibit TMV accumulation and phloem loading. Transcriptional expression studies demonstrate a role for interacting Aux/IAA proteins in the regulation of salicylic acid and jasmonic acid dependent host defense responses as well as virus specific movement factors including pectin methylesterase that are involved in regulating plasmodesmata size exclusion limits and promoting virus cell-to-cell movement. Further characterization of the phloem environment was done using two phloem specific promoters (pSUC2 and pSULTR2;2) to generate epitope-tagged polysomal-RNA complexes. Immuno-purification using the epitope tag allowed us to obtain mRNAs bound to polysomes (the translatome) specifically in phloem tissue. We found the phloem translatome is uniquely altered during TMV infection with 90% and 88% of genes down regulated in the pSUC2 and pSULTR2;2 phloem translatomes, compared to 31% of genes down regulated in the whole plant p35S translatome. Transcripts down regulated in phloem include genes involved in callose deposition at plasmodesmata, host defense responses, and RNA silencing. Combined, these findings indicate TMV reprograms gene expression within the vascular phloem as a means to enhance phloem loading and systemic spread.

FUNCTIONAL CHARACTERIZATION OF THE LINK BETWEEN CARBOHYDRATE METABOLISM AND THE PATHOGENESIS OF THE INVASIVE M1T1 GROUP A STREPTOCOCCUS

The Group A Streptococcus (GAS), or Streptococcus pyogenes, is a strict human pathogen that colonizes a variety of sites within the host. Infections can vary from minor and easily treatable, to life-threatening, invasive forms of disease. In order to adapt to niches, GAS utilizes environmental cues, such as carbohydrates, to coordinate the expression of virulence factors. Research efforts to date have focused on identifying how either components of the phosphoenolpyruvate-phosphotransferase system (PTS) or global transcriptional networks affect the regulation of virulence factors, but not the synergistic relationship between the two. The present study investigates the role of a putative PTS-fructose operon encoded by fruRBA and its role in virulence in the M1T1 strain 5448. Growth in fructose resulted in induction of fruRBA. RT-PCR showed that fruRBA formed an operon, which was repressed by FruR in the absence of fructose. Growth and carbon utilization profiles revealed that although the entire fruRBA operon was required for growth in fructose, FruA was the main fructose transporter. The ability of both ΔfruR and ΔfruB mutants to survive in whole human blood or neutrophils was impaired. However, the phenotypes were not reproduced in murine whole blood or in a mouse intraperitoneal infection, indicating a human-specific mechanism. While it is known that the PTS can affect activity of the Mga virulence regulator, further characterization of the mechanism by which sugars and its protein domains affect activity have not been studied. Transcriptional studies revealed that the core Mga regulon is activated more in a glucose-rich than a glucose-poor environment. This activation correlates with the differential phosphorylation of Mga at its PTS regulatory domains (PRDs). Using a 5448 mga mutant, transcriptome studies in THY or C media established that the Mga regulon reflects the media used. Interestingly, Mga regulates phage-encoded DNases in a low glucose environment. We also show that Mga activity is dependent on C-terminal amino acid interactions that aid in the formation of homodimers. Overall, the studies presented sought to define how external environmental cues, specifically carbohydrates, control complex regulatory networks used by GAS, contribute to pathogenesis, and aid in adaptation to various nutrient conditions encountered.

Genome Wide Association Studies of Phagocytosis and the Cellular Immune Response in Drosophila melanogaster.

Phagocytosis of bacteria by specialized blood cells, known as hemocytes, is a vital component of Drosophila cellular immunity. To identify novel genes that mediate the cellular response to bacteria, we conducted three separate genetic screens using the Drosophila Genetic Reference Panel (DGRP). Adult DGRP lines were tested for the ability of their hemocytes to phagocytose the Gram-positive bacteria Staphylococcus aureus or the Gram-negative bacteria Escherichia coli. The DGRP lines were also screened for the ability of their hemocytes to clear S. aureus infection through the process of phagosome maturation. Genome-wide association analyses were performed to identify potentially relevant single nucleotide polymorphisms (SNPs) associated with the cellular immune phenotypes. The S. aureus phagosome maturation screen identified SNPs near or in 528 candidate genes, many of which have no known role in immunity. Three genes, dpr10, fred, and CG42673, were identified whose loss-of-function in blood cells significantly impaired the innate immune response to S. aureus. The DGRP S. aureus screens identified variants in the gene, Ataxin 2 Binding Protein-1 (A2bp1) as important for the cellular immune response to S. aureus. A2bp1 belongs to the highly conserved Fox-1 family of RNA-binding proteins. Genetic studies revealed that A2bp1 transcript levels must be tightly controlled for hemocytes to successfully phagocytose S. aureus. The transcriptome of infected and uninfected hemocytes from wild type and A2bp1 mutant flies was analyzed and it was found that A2bp1 negatively regulates the expression of the Immunoglobulin-superfamily member Down syndrome adhesion molecule 4 (Dscam4). Silencing of A2bp1 and Dscam4 in hemocytes rescues the fly’s immune response to S. aureus indicating that Dscam4 negatively regulates S. aureus phagocytosis. Overall, we present an examination of the cellular immune response to bacteria with the aim of identifying and characterizing roles for novel mediators of innate immunity in Drosophila. By screening panel of lines in which all genetic variants are known, we successfully identified a large set of candidate genes that could provide a basis for future studies of Drosophila cellular immunity. Finally, we describe a novel, immune-specific role for the highly conserved Fox-1 family member, A2bp1.

A SUBSET OF 3´ TO 5´ EXORIBONUCLEASES ARE THE PRIMARY ENZYMES RESPONSIBLE FOR THE DEGRADATION OF PGPG IN CYCLIC DI-GMP SIGNALING

Bis-(3´-5´)-cyclic dimeric guanosine monophosphate, or cyclic di-GMP (c-di-GMP) is a ubiquitous bacterial second messenger that regulates processes such biofilm formation, motility, and virulence. C-di-GMP is synthesized by diguanylate cyclases (DGCs), while phosphodiesterases (PDE-As) end signaling by linearizing c-di-GMP to 5ʹ-phosphoguanylyl-(3ʹ,5ʹ)-guanosine (pGpG), which is then hydrolyzed to two GMPs by previously unidentified enzymes termed PDE-Bs. To identify the PDE-B responsible for pGpG turnover, a screen for pGpG binding proteins in a Vibrio cholerae open reading frame library was conducted to identify potential pGpG binding proteins. This screen led to identification of oligoribonuclease (Orn). Purified Orn binds to pGpG and can cleave pGpG to GMP in vitro. A deletion mutant of orn in Pseudomonas aeruginosa was highly defective in pGpG turnover and accumulated pGpG. Deletion of orn also resulted in accumulation c-di-GMP, likely through pGpG-mediated inhibition of the PDE-As, causing an increase in c-di-GMP-governed auto-aggregation and biofilm. Thus, we found that Orn serves as the primary PDE-B enzyme in P. aeruginosa that removes pGpG, which is necessary to complete the final step in the c-di-GMP degradation pathway. However, not all bacteria that utilize c-di-GMP signaling also have an ortholog of orn, suggesting that other PDE-Bs must be present. Therefore, we asked whether RNases that cleave small oligoribonucleotides in other species could also act as PDE-Bs. NrnA, NrnB, and NrnC can rapidly degrade pGpG to GMP. Furthermore, they can reduce the elevated aggregation and biofilm formation in P. aeruginosa ∆orn. Together, these results indicate that rather than having a single dedicated PDE-B, different bacteria utilize distinct RNases to cleave pGpG and complete c-di-GMP signaling. The ∆orn strain also has a growth defect, indicating changes in other regulatory processes that could be due to pGpG accumulation, c-di-GMP accumulation, or another effect due to loss of Orn. We sought to investigate the genetic pathways responsible for these growth defect phenotypes by use of a transposon suppressor screen, and also investigated transcriptional changes using RNA-Seq. This work identifies that c-di-GMP degradation intersects with RNA degradation at the point of the Orn and the functionally related RNases.

Comparative transcriptome profiling of human foreskin fibroblasts infected with the Sylvio and Y strains of Trypanosoma cruzi

Trypanosoma cruzi, the causative agent of Chagas Disease, is phylogenetically distributed into nearly identical genetic strains which show divergent clinical presentations including differences in rates of cardiomyopathy in humans, different vector species and transmission cycles, and differential congenital transmission in a mouse model. The population structure of these strains divides into two groups, which are geographically and clinically distinct. The aim of this study was to compare the transcriptome of two strains of T. cruzi, Sylvio vs. Y to identify differences in expression that could account for clinical and biochemical differences. We collected and sequenced RNA from T. cruzi-infected and control Human Foreskin Fibroblasts at three timepoints. Differential expression analysis identified gene expression profiles at different timepoints in Sylvio infections, and between Sylvio and Y infections in both parasite and host. The Sylvio strain parasite and the host response to Sylvio infection largely mirrored the host-pathogen interaction seen in our previous Y strain work. IL-8 was more highly expressed in Sylvio-infected HFFs than in Y-infected HFFs.

Characterization of Borrelia burgdorferi Gene Product BBD18 for Its Role(s) in Pathogen Biology and Infectivity

bbd18 is a differentially expressed Borrelia burgdorferi gene that is transcribed at almost undetectable levels in spirochetes grown in vitro but dramatically upregulated during tick infection. The gene also displays low yet detectable expression at various times in tissues of murine hosts. As the gene product bears no homology to known proteins, its biological significance remains enigmatic. To understand the gene function, we created isogenic bbd18-deletion mutants as well as genetically-complemented isolates from an infectious wild-type B. burgdorferi strain. Compared to parental isolates, bbd18 mutants - but not complemented spirochetes - displayed slower in vitro growth. The bbd18 mutants also reflect significantly reduced ability to persist or remain undetectable both in immunocompetent and SCID mice, yet were able to survive in ticks. This suggests BBD18 function is essential in mammalian hosts but redundant in the arthropod vector. Notably, although bbd18 expression and in vitro growth defects are restored in the complemented isolates, their phenotype is similar to the mutants - being unable to persist in mice but able to survive in ticks. Despite low expression in cultured wild-type B. burgdorferi, bbd18 deletion downregulated several genes. Interestingly, expression of some, including ospD and bbi39, could be complemented, while that of others could not be restored via bbd18 re-expression. Correspondingly, bbd18 mutants displayed altered production of several proteins, and similar to RNA levels, some were restored in the bbd18 complement and others not. To understand how bbd18 deletion results in apparently permanent and noncomplementable phenotypic defects, we sought to genetically disturb the DNA topology surrounding the bbd18 locus without deleting the gene. Spirochetes with an antibiotic cassette inserted downstream of the gene, between bbd17 and bbd18, were significantly attenuated in mice, while a similar upstream insertion, between bbd18 and bbd19, did not affect infectivity, suggesting that an unidentified cis element downstream of bbd18 may encode a virulence-associated factor critical for infection.

HIV Reverse Transcriptase Fidelity And Inhibition Are Modulated By Divalent Cations In A Concentration-Dependent Manner In Vitro

Human immunodeficiency virus (HIV) rapidly evolves through generation and selection of mutants that can escape drug therapy. This process is fueled, in part, by the presumably highly error prone polymerase reverse transcriptase (RT). Fidelity of polymerases can be influenced by cation co-factors. Physiologically, magnesium (Mg2+) is used as a co-factor by RT to perform catalysis, however, alternative cations including manganese (Mn2+), cobalt (Co2+), and zinc (Zn2+) can also be used. I demonstrate here that fidelity and inhibition of HIV RT can be influenced differently, in vitro, by divalent cations depending on their concentration. The reported mutation frequency for purified HIV RT in vitro is typically in the 10-4 range (per nucleotide addition), making the enzyme several-fold less accurate than most polymerases. Paradoxically, results examining HIV replication in cells indicate an error frequency that is ~10 times lower than the error rate obtained in the test tube. Here, I reconcile, at least in part, these discrepancies by showing that HIV RT fidelity in vitro is in the same range as cellular results, in physiological concentrations of free Mg2+ (~0.25 mM). At low Mg2+, mutation rates were 5-10 times lower compared to high Mg2+ conditions (5-10 mM). Alternative divalent cations also have a concentration-dependent effect on RT fidelity. Presumed promutagenic cations Mn2+ and Co2+ decreases the fidelity of RT only at elevated concentrations, and Zn2+, when present in low concentration, increases the fidelity of HIV-1 RT by ~2.5 fold compared to Mg2+. HIV-1 and HIV-2 RT inhibition by nucleoside (NRTIs) and non-nucleoside RT inhibitors (NNRTIs) in vitro is also affected by the Mg2+ concentration. NRTIs lacking 3'-OH group inhibited both enzymes less efficiently in low Mg2+ than in high Mg2+; whereas inhibition by the “translocation defective RT inhibitor”, which retains the 3ʹ-OH, was unaffected by Mg2+ concentration, suggesting that NRTIs with a 3ʹ-OH group may be more potent than other NRTIs. In contrast, NNRTIs were more effective in low vs. high Mg2+ conditions. Overall, the studies presented reveal strategies for designing novel RT inhibitors and strongly emphasize the need for studying HIV RT and RT inhibitors in physiologically relevant low Mg2+ conditions.

KINESIN MOTOR PROTEINS ARE ESSENTIAL FOR MALE GAMETOPHYTE DEVELOPMENT IN MARSILEA VESTITA

The male gametophyte of the semi-aquatic fern, Marsilea vestita, produces multiciliated spermatozoids in a rapid developmental sequence that is controlled post-transcriptionally when dry microspores are placed in water. Development can be divided into two phases, mitosis and differentiation. During the mitotic phase, a series of nine successive division cycles produce 7 sterile cells and 32 spermatids in 4.5-5 hours. During the next 5-6 hours, each spermatid differentiates into a corkscrew-shaped motile spermatozoid with ~140 cilia. This document focuses on the role of motor proteins in the regulation of male gametophyte development and during ciliogenesis. In order to study the mechanisms that regulate spermatogenesis, RNAseq was used to generate a reference transcriptome that allowed us to assess the abundance of transcripts at different stages of development. Over 120 kinesin-like sequences were identified in the transcriptome that represent 56 unique kinesin transcripts. Members of the kinesin-2, -4, -5, -7, -8, -9, -12, -13, and -14 families, in addition to several plant specific and ‘orphan’ kinesins are present. Most (91%) of these kinesin transcripts change in abundance throughout gametophyte development, with 52% of kinesin mRNAs enriched during the mitotic phase and 39% enriched during differentiation. Functional analyses show that the temporal regulation of kinesin transcripts during gametogenesis directly correlates with kinesin protein function. Specifically, Marsilea makes one kinesin-2 (MvKinesin-2) and two kinesin-9 (MvKinesin-9A and MvKinesin-9B) transcripts, which are present during spermatid differentiation and ciliogenesis. Silencing experiments showed that MvKinesin-2 and MvKinesin-9A are required for ciliogenesis and motility in the Marsilea male gametophyte; however, these kinesins display atypical roles during these processes. In contrast, spermatozoids produced after the silencing of MvKinesin-9B exhibit normal morphology. MvKinesin-2 is necessary for cytokinesis as well as for regulating ciliary length and MvKinesin-9A is needed for the correct orientation of basal bodies, events not typically associated with these proteins. In addition, Marsilea makes motile, ciliated gametophytes without the help of IFT dynein, outer arm dynein, or the BBsome. These results are the first to investigate the kinesin-linked mechanisms that regulate ciliogenesis in a land plant.

Alternate Conformations Regulate Ribosomal Recoding in a Positive-sense RNA Virus

Positive-sense RNA viruses are important animal, plant, insect and bacteria pathogens and constitute the largest group of RNA viruses. Due to the relatively small size of their genomes, these viruses have evolved a variety of non-canonical translation mechanisms to optimize coding capacity expanding their proteome diversity. One such strategy is codon redefinition or recoding. First described in viruses, recoding is a programmed translation event in which codon alterations are context dependent. Recoding takes place in a subset of messenger RNA (mRNAs) with some products reflecting new, and some reflecting standard, meanings. The ratio between the two is both critical and highly regulated. While a variety of recoding mechanisms have been documented, (ribosome shunting, stop-carry on, termination-reinitiation, and translational bypassing), the two most extensively employed by RNA viruses are Programmed Ribosomal Frameshifting (PRF) and Programmed Ribosomal Readthrough (PRT). While both PRT and PRF subvert normal decoding for expression of C-terminal extension products, the former involves an alteration of reading frame, and the latter requires decoding of a non-sense codon. Both processes occur at a low but defined frequency, and both require Recoding Stimulatory Elements (RSE) for regulation and optimum functionality. These stimulatory signals can be embedded in the RNA in the form of sequence or secondary structure, or trans-acting factors outside the mRNA such as proteins or micro RNAs (miRNA). Despite 40+ years of study, the precise mechanisms by which viral RSE mediate ribosome recoding for the synthesis of their proteins, or how the ratio of these products is maintained, is poorly defined. This study reveals that in addition to a long distance RNA:RNA interaction, three alternate conformations and a phylogenetically conserved pseudoknot regulate PRT in the carmovirus Turnip crinkle virus (TCV).

Dinoflagellate Genomic Organization and Phylogenetic Marker Discovery Utilizing Deep Sequencing Data

Dinoflagellates possess large genomes in which most genes are present in many copies. This has made studies of their genomic organization and phylogenetics challenging. Recent advances in sequencing technology have made deep sequencing of dinoflagellate transcriptomes feasible. This dissertation investigates the genomic organization of dinoflagellates to better understand the challenges of assembling dinoflagellate transcriptomic and genomic data from short read sequencing methods, and develops new techniques that utilize deep sequencing data to identify orthologous genes across a diverse set of taxa. To better understand the genomic organization of dinoflagellates, a genomic cosmid clone of the tandemly repeated gene Alchohol Dehydrogenase (AHD) was sequenced and analyzed. The organization of this clone was found to be counter to prevailing hypotheses of genomic organization in dinoflagellates. Further, a new non-canonical splicing motif was described that could greatly improve the automated modeling and annotation of genomic data. A custom phylogenetic marker discovery pipeline, incorporating methods that leverage the statistical power of large data sets was written. A case study on Stramenopiles was undertaken to test the utility in resolving relationships between known groups as well as the phylogenetic affinity of seven unknown taxa. The pipeline generated a set of 373 genes useful as phylogenetic markers that successfully resolved relationships among the major groups of Stramenopiles, and placed all unknown taxa on the tree with strong bootstrap support. This pipeline was then used to discover 668 genes useful as phylogenetic markers in dinoflagellates. Phylogenetic analysis of 58 dinoflagellates, using this set of markers, produced a phylogeny with good support of all branches. The Suessiales were found to be sister to the Peridinales. The Prorocentrales formed a monophyletic group with the Dinophysiales that was sister to the Gonyaulacales. The Gymnodinales was found to be paraphyletic, forming three monophyletic groups. While this pipeline was used to find phylogenetic markers, it will likely also be useful for finding orthologs of interest for other purposes, for the discovery of horizontally transferred genes, and for the separation of sequences in metagenomic data sets.

Exploring the Function and Regulation of Arabidopsis EIN2 in Ethylene Signaling

Ethylene is an essential plant hormone involved in nearly all stages of plant growth and development. EIN2 (ETHYLENE INSENSITIVE2) is a master positive regulator in the ethylene signaling pathway, consisting of an N-terminal domain and a C-terminal domain. The EIN2 N-terminal domain localizes to the endoplasmic reticulum (ER) membrane and shows sequence similarity to Nramp metal ion transporters. The cytosolic C-terminal domain is unique to plants and signals downstream. There have been several major gaps in our knowledge of EIN2 function. It was unknown how the ethylene signal gets relayed from the known upstream component CTR1 (CONSTITUTIVE RESPONSE1) a Ser/Thr kinase at the ER, to EIN2. How the ethylene signal was transduced from EIN2 to the next downstream component transcription factor EIN3 (ETHYLENE INSENSITIVE3) in the nucleus was also unknown. The N-terminal domain of EIN2 shows homology to Nramp metal ion transporters and whether EIN2 can also function as a metal transporter has been a question plaguing the ethylene field for almost two decades. Here, EIN2 was found to interact with the CTR1 protein kinase, leading to the discovery that CTR1 phosphorylates the C-terminal domain of EIN2 in Arabidopsis thaliana. Using tags at the termini of EIN2, it was deduced that in the presence of ethylene, the EIN2 C-terminal domain is cleaved and translocates into the nucleus, where it could somehow activate downstream ethylene responses. The EIN2 C-terminal domain interacts with nuclear proteins, RTE3 and EER5, which are components of the TREX-2 mRNA export complex, although the role of these interactions remains unclear. The EIN2 N-terminal domain was found to be capable of divalent metal transport when expressed in E. coli and S. cerevisiae leading to the hypothesis that metal transport plays a role in ethylene signaling. This hypothesis was tested using a novel missense allele, ein2 G36E, substituting a highly conserved residue that is required for metal transport in Nramp proteins. This G36E substitution did not disrupt metal ion transport of EIN2, but the ethylene insensitive phenotype of this mutant indicates that the EIN2 N-terminal domain is important for positively regulating the C-terminal domain. The defect of the ein2 G36E mutant does not prevent proper expression or subcellular localization, but might affect protein modifications. The ein2 G36E allele is partially dominant, mostly likely displaying haploinsufficiency. Overexpression of the EIN2 N-terminal domain in the ein2 G36E mutant did not rescue ethylene insensitivity, suggesting the N-terminal domain functions in cis to regulate the C-terminal domain. These findings advance our knowledge of EIN2, which is critical to understanding ethylene signaling.