Cis-acting elements and developmental regulators govern toxin gene expression in Bacillus anthracis

Expression of the structural genes for the anthrax toxin proteins is coordinately controlled by host-related signals such as elevated CO2 , and the trans-acting positive regulator, AtxA. No specific binding of AtxA to the toxin gene promoters has been demonstrated and no sequence-based similarities are apparent in the promoter regions of toxin genes. We hypothesized that the toxin genes possess common structural features that are required for positive regulation. To test this hypothesis, I performed an extensive characterization of the toxin gene promoters. I determined the minimal sequences required for atxA-mediated toxin gene expression and compared these sequences for structural similarities. In silico modeling and in vitro experiments indicated significant curvature within these regions. Random mutagenesis revealed that point mutations associated with reduced transcriptional activity, mostly mapped to areas of high curvature. This work enabled the identification of two potential cis-acting elements implicated in AtxA-mediated regulation of the toxin genes. In addition to the growth condition requirements and AtxA, toxin gene expression is under growth phase regulation. The transition state regulator AbrB represses atxA expression to influence toxin synthesis. Here I report that toxin gene expression also requires sigH, a gene encoding the RNA polymerase sigma factor associated with development in B. subtilis. In the well-studied B. subtilis system, σH is part of a feedback control pathway that involves AbrB and the major response regulator of sporulation initiation, Spo0A. My data indicate that in B. anthracis, regulatory relationships exist between these developmental regulators and atxA . Interestingly, during growth in toxin-inducing conditions, sigH and abrB expression deviates from that described for B. subtilis, affecting expression of the atxA gene. These findings, combined with previous observations, suggest that the steady state level of atxA expression is critical for optimal toxin gene transcription. I propose a model whereby, under toxin-inducing conditions, control of toxin gene expression is fine-tuned by the independent effects of the developmental regulators on the expression of atxA . The growth condition-dependent changes in expression of these regulators may be crucial for the correct timing and uninterrupted expression of the toxin genes during infection. ^

Systems Biology Approaches to Probe Gene Regulation in Bacteria

Mechanisms that allow pathogens to colonize the host are not the product of isolated genes, but instead emerge from the concerted operation of regulatory networks. Therefore, identifying components and the systemic behavior of networks is necessary to a better understanding of gene regulation and pathogenesis. To this end, I have developed systems biology approaches to study transcriptional and post-transcriptional gene regulation in bacteria, with an emphasis in the human pathogen Mycobacterium tuberculosis (Mtb). First, I developed a network response method to identify parts of the Mtb global transcriptional regulatory network utilized by the pathogen to counteract phagosomal stresses and survive within resting macrophages. As a result, the method unveiled transcriptional regulators and associated regulons utilized by Mtb to establish a successful infection of macrophages throughout the first 14 days of infection. Additionally, this network-based analysis identified the production of Fe-S proteins coupled to lipid metabolism through the alkane hydroxylase complex as a possible strategy employed by Mtb to survive in the host. Second, I developed a network inference method to infer the small non-coding RNA (sRNA) regulatory network in Mtb. The method identifies sRNA-mRNA interactions by integrating a priori knowledge of possible binding sites with structure-driven identification of binding sites. The reconstructed network was useful to predict functional roles for the multitude of sRNAs recently discovered in the pathogen, being that several sRNAs were postulated to be involved in virulence-related processes. Finally, I applied a combined experimental and computational approach to study post-transcriptional repression mediated by small non-coding RNAs in bacteria. Specifically, a probabilistic ranking methodology termed rank-conciliation was developed to infer sRNA-mRNA interactions based on multiple types of data. The method was shown to improve target prediction in Escherichia coli, and therefore is useful to prioritize candidate targets for experimental validation.

Cross -packing among retroviruses and characterization of the feline immunodeficiency virus (FIV) packaging signal: Implications for the development of FIV-based gene transfer systems

Feline immunodeficiency virus (FIV)-based gene transfer systems are being seriously considered for human gene therapy as an alternative to vectors based on primate lentiviruses, a genetically complex group of retroviruses capable of infecting non-dividing cells. The greater phylogenetic distance between the feline and primate lentiviruses is thought to reduce chances of the generation of recombinant viruses. However, safety of FIV-based vector systems has not been tested experimentally. Since primate lentiviruses such as human and simian immunodeficiency viruses (HIV/SIV) can cross-package each other's genomes, we tested this trait with respect to FIV. Unexpectedly, both feline and primate lentiviruses were reciprocally able to both cross-package and propagate each other's RNA genomes. This was largely due to the recognition of viral packaging signals by the heterologous proteins. However, a simple retrovirus such as Mason-Pfizer monkey virus (MPMV) was unable to package FIV RNA. Interestingly, FIV could package MPMV RNA, but not propagate it for further steps of replication. These findings suggest that upon co-infection of the same host, cross-packaging may allow distinct retroviruses to generate chimeric variants with unknown pathogenic potential. ^ In order to understand the packaging determinants in FIV, we conducted a detailed mutational analysis of the region thought to contain FIV packaging signal. We show that the first 90–120 nt of the 5′ untranslated region (UTR) and the first 90 nt of gag were simultaneously required for efficient FIV RNA packaging. These results suggest that the primary FIV packaging signal is multipartite and discontinuous, composed of two core elements separated by 150 nt of the 5 ′UTR. ^ The above studies are being used towards the development of safer FIV-based self-inactivating (SIN) vectors. These vectors are being designed to eliminate the ability of FIV transfer vector RNAs to be mobilized by primate lentiviral proteins that may be present in the target cells. Preliminary test of the first generation of these vectors has revealed that they are incapable of being propagated by feline proteins. The inability of FIV transfer vectors to express packageable vector RNA after integration should greatly increase the safety of FIV vectors for human gene therapy. ^

Physiological response to elevated temperature and pCO2 varies across four Pacific coral species: Understanding the unique host + symbiont response

The physiological response to individual and combined stressors of elevated temperature and pCO2 were measured over a 24-day period in four Pacific corals and their respective symbionts (Acropora millepora/Symbiodinium C21a, Pocillopora damicornis/Symbiodinium C1c-d-t, Montipora monasteriata/Symbiodinium C15, and Turbinaria reniformis/Symbiodinium trenchii). Multivariate analyses indicated that elevated temperature played a greater role in altering physiological response, with the greatest degree of change occurring within M. monasteriata and T. reniformis. Algal cellular volume, protein, and lipid content all increased for M. monasteriata. Likewise, S. trenchii volume and protein content in T. reniformis also increased with temperature. Despite decreases in maximal photochemical efficiency, few changes in biochemical composition (i.e. lipids, proteins, and carbohydrates) or cellular volume occurred at high temperature in the two thermally sensitive symbionts C21a and C1c-d-t. Intracellular carbonic anhydrase transcript abundance increased with temperature in A. millepora but not in P. damicornis, possibly reflecting differences in host mitigated carbon supply during thermal stress. Importantly, our results show that the host and symbiont response to climate change differs considerably across species and that greater physiological plasticity in response to elevated temperature may be an important strategy distinguishing thermally tolerant vs. thermally sensitive species.

Populations genomics analysis of legume host preference for specific rhizobial genotypes in the Rhizobium leguminosarum bv. viciae symbioses

Rhizobium leguminosarum bv. viciae establishes root nodule symbioses with several legume genera. Although most isolates are equally effective in establishing symbioses with all host genera, previous evidence suggests that hosts select specific rhizobial genotypes among those present in the soil. We have used population genomics to further investigate this observation. P. sativum, L. culinaris, V. sativa, and V. faba plants were used to trap rhizobia from a well-characterized soil, and pooled genomic DNAs from one-hundred isolates from each plant were sequenced. Sequence reads were aligned to the R. leguminosarum bv. viciae 3841 reference genome. High overall conservation of sequences was observed in all subpopulations, although several multigenic regions were absent from the soil population. A large fraction (16-22%) of sequence reads could not be recruited to the reference genome, suggesting that they represent sequences specific to that particular soil population. Although highly conserved, the 16S-23S rRNA gene region presented single nucleotide polymorphisms (SNPs) regarding the reference genome, but no striking differences could be found among plant-selected subpopulations. Plant-specific SNP patterns were, however, clearly observed within the nod gene cluster, supporting the existence of a plant preference for specific rhizobial genotypes. This was also shown after genome-wide analysis of SNP patterns.

Metagenomic Anlaysis of microsymbiont selection by the legume plant host

Rhizobium leguminosarum bv.viciae is able to establish nitrogen-fixing symbioses with legumes of the genera Pisum, Lens, Lathyrus and Vicia. Classic studies using trap plants (Laguerre et al., Young et al.) provided evidence that different plant hosts are able to select different rhizobial genotypes among those available in a given soil. However, these studies were necessarily limited by the paucity of relevant biodiversity markers. We have now reappraised this problem with the help of genomic tools. A well-characterized agricultural soil (INRA Bretennieres) was used as source of rhizobia. Plants of Pisum sativum, Lens culinaris, Vicia sativa and V. faba were used as traps. Isolates from 100 nodules were pooled, and DNA from each pool was sequenced (BGI-Hong Kong; Illumina Hiseq 2000, 500 bp PE libraries, 100 bp reads, 12 Mreads). Reads were quality filtered (FastQC, Trimmomatic), mapped against reference R. leguminosarum genomes (Bowtie2, Samtools), and visualized (IGV). An important fraction of the filtered reads were not recruited by reference genomes, suggesting that plant isolates contain genes that are not present in the reference genomes. For this study, we focused on three conserved genomic regions: 16S-23S rDNA, atpD and nodDABC, and a Single Nucleotide Polymorphism (SNP) analysis was carried out with meta / multigenomes from each plant. Although the level of polymorphism varied (lowest in the rRNA region), polymorphic sites could be identified that define the specific soil population vs. reference genomes. More importantly, a plant-specific SNP distribution was observed. This could be confirmed with many other regions extracted from the reference genomes (data not shown). Our results confirm at the genomic level previous observations regarding plant selection of specific genotypes. We expect that further, ongoing comparative studies on differential meta / multigenomic sequences will identify specific gene components of the plant-selected genotypes

Effect of the Legume Plant Host on Rhizobial Soil Population Dynamics

Rhizobium leguminosarum bv viciae (Rlv) is a soil bacterium able to establish specific root-nodule symbioses with legumes of four different genera: Pisum, Vicia, Lens and Lathyrus. Rlv isolates from nodules of any of these legumes can nodulate any of them; however, it has been shown that plants select specific rhizobial genotypes from those present in the soil (1,2). We have previously shown this at the genomic level by following a population genomics approach. Pool genomic sequences from 100 isolates from each of four plant species: P. sativum, L. culinaris, V. faba and V. sativa, show different, specific profiles at the single nucleotide polymorphism (SNP) level for relevant genes. In this work, the extent of Rlv selection from a well-characterized soil population by different legume plant hosts: P. sativum, L. culinaris, V. faba and V. sativa, after a medium-term mesocosm study is described. Direct soil isolates from each of these mesocosm studies have been tested for specific rhizobial genes (glnII and fnrN) and symbiotic genes (nodC and nifH). Different populations were characterized further by Sanger sequencing of both the rpoB phylogenetic marker gene and the symbiotic genes nodC and nifH. The distribution and size of the rhizobial population for each legume host showed changes during the medium-term mesocosm study. Particularly, a non-symbiotic group of rhizobia was enriched by all four hosts, in contrast to the symbiotic rhizobia profile, which was specific for each legume plant host.

Extracellular enveloped vaccinia virus is resistant to complement because of incorporation of host complement control proteins into its envelope

Vaccinia virus (VV) produces two antigenically and structurally distinct infectious virions, intracellular mature virus (IMV) and extracellular enveloped virus (EEV). Here we have investigated the resistance of EEV and IMV to neutralization by complement in the absence of immune antibodies. When EEV is challenged with complement from the same species as the cells used to grow the virus, EEV is resistant to neutralization by complement, whereas IMV is not. EEV resistance was not a result of EEV protein B5R, despite its similarity to proteins of the regulators of complement activation (RCA) family, or to any of the other EEV proteins tested (A34R, A36R, and A56R gene products). EEV was sensitive to complement when the virus was grown in one species and challenged with complement from a different species, suggesting that complement resistance might be mediated by host RCA incorporated into the EEV outer envelope. This hypothesis was confirmed by several observations: (i) immunoblot analysis revealed that cellular membrane proteins CD46, CD55, CD59, CD71, CD81, and major histocompatibility complex class I antigen were detected in purified EEV but not IMV; (ii) immunoelectron microscopy revealed cellular RCA on the surface of EEV retained on the cell surface; and (iii) EEV derived from rat cells expressing the human RCA CD55 or CD55 and CD59 were more resistant to human complement than EEV derived from control rat cells that expressed neither CD55 nor CD59. These data justify further analysis of the roles of these (and possible other) cellular proteins in EEV biology.

Transgenic DNA integrated into the oat genome is frequently interspersed by host DNA

Integration of transgenic DNA into the plant genome was investigated in 13 transgenic oat (Avena sativa L.) lines produced using microprojectile bombardment with one or two cotransformed plasmids. In all transformation events, the transgenic DNA integrated into the plant genome consisted of intact transgene copies that were accompanied by multiple, rearranged, and/or truncated transgene fragments. All fragments of transgenic DNA cosegregated, indicating that they were integrated at single gene loci. Analysis of the structure of the transgenic loci indicated that the transgenic DNA was interspersed by the host genomic DNA. The number of insertions of transgenic DNA within the transgene loci varied from 2 to 12 among the 13 lines. Restriction endonucleases that do not cleave the introduced plasmids produced restriction fragments ranging from 3.6 to about 60 kb in length hybridizing to a probe comprising the introduced plasmids. Although the size of the interspersing host DNA within the transgene locus is unknown, the sizes of the transgene-hybridizing restriction fragments indicated that the entire transgene locus must be at least from 35–280 kb. The observation that all transgenic lines analyzed exhibited genomic interspersion of multiple clustered transgenes suggests a predominating integration mechanism. We propose that transgene integration at multiple clustered DNA replication forks could account for the observed interspersion of transgenic DNA with host genomic DNA within transgenic loci.

Altered 3′-terminal RNA structure in phage Qβ adapted to host factor-less Escherichia coli

The RNA phage Qβ requires for the replication of its genome an RNA binding protein called Qβ host factor or Hfq protein. Our previous results suggested that this protein mediates the access of replicase to the 3′-end of the Qβ plus strand RNA. Here we report the results of an evolutionary experiment in which phage Qβ was adapted to an Escherichia coli Q13 host strain with an inactivated host factor (hfq) gene. This strain initially produced phage at a titer ≈10,000-fold lower than the wild-type strain and with minute plaque morphology, but after 12 growth cycles, phage titer and plaque size had evolved to levels near those of the wild-type host. RNAs isolated from adapted Qβ mutants were efficient templates for replicase without host factor in vitro. Electron microscopy showed that mutant RNAs, in contrast to wild-type RNA, efficiently interacted with replicase at the 3′-end in the absence of host factor. The same set of four mutations in the 3′-terminal third of the genome was found in several independently evolved phage clones. One mutation disrupts the base pairing of the 3′-terminal CCCoh sequence, suggesting that the host factor stimulates activity of the wild-type RNA template by melting out its 3′-end.

Selective maintenance of allozyme differences among sympatric host races of the apple maggot fly

Whether phytophagous insects can speciate in sympatry when they shift and adapt to new host plants is a controversial question. One essential requirement for sympatric speciation is that disruptive selection outweighs gene flow between insect populations using different host plants. Empirical support for host-related selection (i.e., fitness trade-offs) is scant, however. Here, we test for host-dependent selection acting on apple (Malus pumila)- and hawthorn (Crataegus spp.)-infesting races of Rhagoletis pomonella (Diptera: Tephritidae). In particular, we examine whether the earlier fruiting phenology of apple trees favors pupae in deeper states of diapause (or with slower metabolisms/development rates) in the apple fly race. By experimentally lengthening the time period preceding winter, we exposed hawthorn race pupae to environmental conditions typically faced by apple flies. This exposure induced a significant genetic response at six allozyme loci in surviving hawthorn fly adults toward allele frequencies found in the apple race. The sensitivity of hawthorn fly pupae to extended periods of warm weather therefore selects against hawthorn flies that infest apples and helps to maintain the genetic integrity of the apple race by counteracting gene flow from sympatric hawthorn populations. Our findings confirm that postzygotic reproductive isolation can evolve as a pleiotropic consequence of host-associated adaptation, a central tenet of nonallopatric speciation. They also suggest that one reason for the paucity of reported fitness trade-offs is a failure to consider adequately costs associated with coordinating an insect’s life cycle with the phenology of its host plant.

Noncytopathic Sindbis virus RNA vectors for heterologous gene expression

Infection of vertebrate cells with alphaviruses normally leads to prodigious expression of virus-encoded genes and a dramatic inhibition of host protein synthesis. Recombinant Sindbis viruses and replicons have been useful as vectors for high level foreign gene expression, but the cytopathic effects of viral replication have limited their use to transient studies. We recently selected Sindbis replicons capable of persistent, noncytopathic growth in BHK cells and describe here a new generation of Sindbis vectors useful for long-term foreign gene expression based on such replicons. Foreign genes of interest as well as the dominant selectable marker puromycin N-acteyltransferase, which confers resistance to the drug puromycin, were expressed as subgenomic transcripts of noncytopathic replicons or defective-interfering genomes complemented in trans by a replicon. Based on these strategies, we developed vectors that can be initiated via either RNA or DNA transfection and analyzed them for their level and stability of foreign gene expression. Noncytopathic Sindbis vectors express reasonably high levels of protein in nearly every cell. These vectors should prove to be flexible tools for the rapid expression of heterologous genes under conditions in which cellular metabolism is not perturbed, and we illustrate their utility with a number of foreign proteins.

Up-regulation of specific tyrosinase mRNAs in mouse melanomas with the c2j gene substituted for the wild-type tyrosinase allele: Utilization in design of syngeneic immunotherapy models

The expression of cell-specialization genes is likely to be changing in tumor cells as their differentiation declines. Functional changes in these genes might yield unusual peptide epitopes with anti-tumor potential and could occur without modification in the DNA sequence of the gene. Melanomas undergo a characteristic decline in melanization that may reflect altered contributions of key melanocytic genes such as tyrosinase. Quantitative reverse transcriptase–PCR of the wild-type (C) tyrosinase gene in transgenic (C57BL/6 strain) mouse melanomas has revealed a shift toward alternative splicing of the pre-mRNA that generated increased levels of the Δ1b and Δ1d mRNA splice variants. The spontaneous c2j albino mutation of tyrosinase (in the C57BL/6 strain) changes the pre-mRNA splicing pattern. In c2j/c2j melanomas, alternative splicing was again increased. However, while some mRNAs (notably Δ1b) present in C/C were obligatorily absent, others (Δ3 and Δ1d) were elevated. In c2j/c2j melanomas, the percentage of total tyrosinase transcripts attributable to Δ3 reached approximately 2-fold the incidence in c2j/c2j or C/C skin melanocytes. The percentage attributable to Δ1d rose to approximately 2-fold the incidence in c2j/c2j skin, and to 10-fold that in C/C skin. These differences provide a basis for unique mouse models in which the melanoma arises in skin grafted from a C/C or c2j/c2j transgenic donor to a transgenic host of the same or opposite tyrosinase genotype. Immunotherapy designs then could be based on augmenting those antigenic peptides that are novel or overrepresented in a tumor relative to the syngeneic host.

Recurrent invasion and extinction of a selfish gene

Homing endonuclease genes show super-Mendelian inheritance, which allows them to spread in populations even when they are of no benefit to the host organism. To test the idea that regular horizontal transmission is necessary for the long-term persistence of these genes, we surveyed 20 species of yeasts for the ω-homing endonuclease gene and associated group I intron. The status of ω could be categorized into three states (functional, nonfunctional, or absent), and status was not clustered on the host phylogeny. Moreover, the phylogeny of ω differed significantly from that of the host, strong evidence of horizontal transmission. Further analyses indicate that horizontal transmission is more common than transposition, and that it occurs preferentially between closely related species. Parsimony analysis and coalescent theory suggest that there have been 15 horizontal transmission events in the ancestry of our yeast species, through simulations indicate that this value is probably an underestimate. Overall, the data support a cyclical model of invasion, degeneration, and loss, followed by reinvasion, and each of these transitions is estimated to occur about once every 2 million years. The data are thus consistent with the idea that frequent horizontal transmission is necessary for the long-term persistence of homing endonuclease genes, and further, that this requirement limits these genes to organisms with easily accessible germ lines. The data also show that mitochondrial DNA sequences are transferred intact between yeast species; if other genes do not show such high levels of horizontal transmission, it would be due to lack of selection, rather than lack of opportunity.

Host factor I, Hfq, binds to Escherichia coli ompA mRNA in a growth rate-dependent fashion and regulates its stability

The stability of the ompA mRNA depends on the bacterial growth rate. The 5′ untranslated region is the stability determinant of this transcript and the target of the endoribonuclease, RNase E, the key player of mRNA degradation. An RNA-binding protein with affinity for the 5′ untranslated region ompA was purified and identified as Hfq, a host factor initially recognized for its function in phage Qβ replication. The ompA RNA-binding activity parallels the amount of Hfq, which is elevated in bacteria cultured at slow growth rate, a condition leading to facilitated degradation of the ompA mRNA. In hfq mutant cells with a deficient Hfq gene product, the RNA-binding activity is missing, and analysis of the ompA mRNA showed that the growth-rate dependence of degradation is lost. Furthermore, the half-life of the ompA mRNA is prolonged in the mutant cells, irrespective of growth rate. Hfq has no affinity for the lpp transcript whose degradation, like that of bulk mRNA, is not affected by bacterial growth rate. Compatible with our results, we found that the intracellular concentration of RNase E and its associated degradosome components is independent of bacterial growth rate. Thus our results suggest a regulatory role for Hfq that specifically facilitates the ompA mRNA degradation in a growth rate-dependent manner.