Tightly Regulated Expression of Autographa californica Multicapsid Nucleopolyhedrovirus Immediate Early Genes Emerges from Their Interactions and Possible Collective Behaviors

To infect their hosts, DNA viruses must successfully initiate the expression of viral genes that control subsequent viral gene expression and manipulate the host environment. Viral genes that are immediately expressed upon infection play critical roles in the early infection process. In this study, we investigated the expression and regulation of five canonical regulatory immediate-early (IE) genes of Autographa californica multicapsid nucleopolyhedrovirus: ie0, ie1, ie2, me53, and pe38. A systematic transient gene-expression analysis revealed that these IE genes are generally transactivators, suggesting the existence of a highly interactive regulatory network. A genetic analysis using gene knockout viruses demonstrated that the expression of these IE genes was tolerant to the single deletions of activator IE genes in the early stage of infection. A network graph analysis on the regulatory relationships observed in the transient expression analysis suggested that the robustness of IE gene expression is due to the organization of the IE gene regulatory network and how each IE gene is activated. However, some regulatory relationships detected by the genetic analysis were contradictory to those observed in the transient expression analysis, especially for IE0-mediated regulation. Statistical modeling, combined with genetic analysis using knockout alleles for ie0 and ie1, showed that the repressor function of ie0 was due to the interaction between ie0 and ie1, not ie0 itself. Taken together, these systematic approaches provided insight into the topology and nature of the IE gene regulatory network.     

Controlling gene expression in mycobacteria with anhydrotetracycline and Tet repressor

Gene expression systems that allow the regulation of bacterial genes during an infection are valuable molecular tools but are lacking for mycobacterial pathogens. We report the development of mycobacterial gene regulation systems that allow controlling gene expression in fast and slow-growing mycobacteria, including Mycobacterium tuberculosis, using anhydrotetracycline (ATc) as inducer. The systems are based on the Escherichia coli Tn10-derived tet regulatory system and consist of a strong tet operator (tetO)-containing mycobacterial promoter, expression cassettes for the repressor TetR and the chemical inducer ATc. These systems allow gene regulation over two orders of magnitude in Mycobacterium smegmatis and M.tuberculosis. TetR-controlled gene expression was inducer concentration-dependent and maximal with ATc concentrations at least 10- and 20-fold below the minimal inhibitory concentration for M.smegmatis and M.tuberculosis, respectively. Using the essential mycobacterial gene ftsZ, we showed that these expression systems can be used to construct conditional knockouts and to analyze the function of essential mycobacterial genes. Finally, we demonstrated that these systems allow gene regulation in M.tuberculosis within the macrophage phagosome.     

L5 luciferase reporter mycobacteriophages: a sensitive tool for the detection and assay of live mycobacteria

Recombinant bacteriophages provide efficient delivery systems for introducing reporter genes into specific bacterial hosts. We have constructed mycobacteriophage L5 recombinants carrying the firefly luciferase gene inserted into the tRNA region of the phage genome. Infection of Mycobacterium smegmatis by these phages results in expression of the luciferase gene and light emission. Fortuitously, the luciferase gene is expressed continuously in lysogens surviving infection. Synthesis of luciferase from a mycobacterial promoter created by cloning enables the detection of extremely small numbers of M. smegmatis cells. These reporter phages can be used to discriminate between drug-sensitive and drug-resistant strains of M. smegmatis, and may provide tools for the rapid identification and classification of antimycobacterial agents.     

The Impact of Gene Expression Variation on the Robustness and Evolvability of a Developmental Gene Regulatory Network

Regulatory interactions buffer development against genetic and environmental perturbations, but adaptation requires phenotypes to change. We investigated the relationship between robustness and evolvability within the gene regulatory network underlying development of the larval skeleton in the sea urchin Strongylocentrotus purpuratus. We find extensive variation in gene expression in this network throughout development in a natural population, some of which has a heritable genetic basis. Switch-like regulatory interactions predominate during early development, buffer expression variation, and may promote the accumulation of cryptic genetic variation affecting early stages. Regulatory interactions during later development are typically more sensitive (linear), allowing variation in expression to affect downstream target genes. Variation in skeletal morphology is associated primarily with expression variation of a few, primarily structural, genes at terminal positions within the network. These results indicate that the position and properties of gene interactions within a network can have important evolutionary consequences independent of their immediate regulatory role.     

Inferring Regulatory Networks from Expression Data Using Tree-Based Methods

One of the pressing open problems of computational systems biology is the elucidation of the topology of genetic regulatory networks (GRNs) using high throughput genomic data, in particular microarray gene expression data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) challenge aims to evaluate the success of GRN inference algorithms on benchmarks of simulated data. In this article, we present GENIE3, a new algorithm for the inference of GRNs that was best performer in the DREAM4 In Silico Multifactorial challenge. GENIE3 decomposes the prediction of a regulatory network between p genes into p different regression problems. In each of the regression problems, the expression pattern of one of the genes (target gene) is predicted from the expression patterns of all the other genes (input genes), using tree-based ensemble methods Random Forests or Extra-Trees. The importance of an input gene in the prediction of the target gene expression pattern is taken as an indication of a putative regulatory link. Putative regulatory links are then aggregated over all genes to provide a ranking of interactions from which the whole network is reconstructed. In addition to performing well on the DREAM4 In Silico Multifactorial challenge simulated data, we show that GENIE3 compares favorably with existing algorithms to decipher the genetic regulatory network of Escherichia coli. It doesn''t make any assumption about the nature of gene regulation, can deal with combinatorial and non-linear interactions, produces directed GRNs, and is fast and scalable. In conclusion, we propose a new algorithm for GRN inference that performs well on both synthetic and real gene expression data. The algorithm, based on feature selection with tree-based ensemble methods, is simple and generic, making it adaptable to other types of genomic data and interactions.     

Host-Parasite Interactions and Purifying Selection in a Microsporidian Parasite of Honey Bees

To clarify the mechanisms of Nosema ceranae parasitism, we deep-sequenced both honey bee host and parasite mRNAs throughout a complete 6-day infection cycle. By time-series analysis, 1122 parasite genes were significantly differently expressed during the reproduction cycle, clustering into 4 expression patterns. We found reactive mitochondrial oxygen species modulator 1 of the host to be significantly down regulated during the entire infection period. Our data support the hypothesis that apoptosis of honey bee cells was suppressed during infection. We further analyzed genome-wide genetic diversity of this parasite by comparing samples collected from the same site in 2007 and 2013. The number of SNP positions per gene and the proportion of non-synonymous substitutions per gene were significantly reduced over this time period, suggesting purifying selection on the parasite genome and supporting the hypothesis that a subset of N. ceranae strains might be dominating infection.     

Characterization of the mIHF Gene of Mycobacterium smegmatis

Integration of mycobacteriophage L5 requires the mycobacterial integration host factor (mIHF) in vitro. mIHF is a 105-residue heat-stable polypeptide that is not obviously related to HU or any other small DNA-binding proteins. mIHF is most abundant just prior to entry into stationary phase and is essential for the viability of Mycobacterium smegmatis.     

Characterization and classification of Bo4 as a cluster G mycobacteriophage that can infect and lyse M. tuberculosis

Mycobacteriophage therapy is a potential alternative treatment for Mycobacterium tuberculosis infection. Here, we further characterized a mycobacteriophage, Bo4, and evaluated its ability to infect and kill M. tuberculosis. We first found that Bo4 can infect M. tuberculosis and Mycobacterium smegmatis. The observed clear plaques created by Bo4 infection indicated that Bo4 might be a lytic phage able to lyse mycobacterial strains, which was confirmed by phage antimicrobial activity. Bo4 formed clear zones in a medium with pH values of 7.4 or 5.0, suggesting the possibility that Bo4 could lyse mycobacteria, such as M. tuberculosis, in blood as well as in lysosomal macrophages. Further investigation into the Bo4 genome revealed that Bo4 had a dsDNA genome. Moreover, Bo4 contained ~39,318 bp comprised of 66.76 % G+C content. Complete genome sequencing showed high nucleotide identity with cluster G mycobacteriophages, thus classifying Bo4 as a member of the cluster G family. Additionally, annotation of the Bo4 genome indicated that it was a lytic bacteriophage and did not contain any harmful genes that increased mycobacterial virulence or decreased human immunity. Overall, the results of investigation indicate that the Bo4 possesses the potential to destroy M. tuberculosis, making it a potentially useful tool for diagnosing and treating tuberculosis.     

LANA-Mediated Recruitment of Host Polycomb Repressive Complexes onto the KSHV Genome during De Novo Infection

One of the hallmarks of the latent phase of Kaposi’s sarcoma-associated herpesvirus (KSHV) infection is the global repression of lytic viral gene expression. Following de novo KSHV infection, the establishment of latency involves the chromatinization of the incoming viral genomes and recruitment of the host Polycomb repressive complexes (PRC1 and PRC2) to the promoters of lytic genes, which is accompanied by the inhibition of lytic genes. However, the mechanism of how PRCs are recruited to the KSHV episome is still unknown. Utilizing a genetic screen of latent genes in the context of KSHV genome, we identified the latency-associated nuclear antigen (LANA) to be responsible for the genome-wide recruitment of PRCs onto the lytic promoters following infection. We found that LANA initially bound to the KSHV genome right after infection and subsequently recruited PRCs onto the viral lytic promoters, thereby repressing lytic gene expression. Furthermore, both the DNA and chromatin binding activities of LANA were required for the binding of LANA to the KSHV promoters, which was necessary for the recruitment of PRC2 to the lytic promoters during de novo KSHV infection. Consequently, the LANA-knockout KSHV could not recruit PRCs to its viral genome upon de novo infection, resulting in aberrant lytic gene expression and dysregulation of expression of host genes involved in cell cycle and proliferation pathways. In this report, we demonstrate that KSHV LANA recruits host PRCs onto the lytic promoters to suppress lytic gene expression following de novo infection.     

Identification and characterization of genes controlled by the sporulation-regulatory gene spo0H in Bacillus subtilis.

We describe a general strategy for the identification of genes that are controlled by a specific regulatory factor in vivo and the use of this strategy to identify genes in Bacillus subtilis that are controlled by spo0H, a regulatory gene required for the initiation of sporulation. The general strategy makes use of a cloned regulatory gene fused to an inducible promoter to control expression of the regulatory gene and random gene fusions to a reporter gene to monitor expression in the presence and absence of the regulatory gene product. spo0H encodes a sigma factor of RNA polymerase, sigma H, and is required for the extensive reprograming of gene expression during the transition from growth to stationary phase and during the initiation of sporulation. We identified 18 genes that are controlled by sigma H (csh genes) in vivo by monitoring expression of random gene fusions to lacZ, made by insertion mutagenesis with the transposon Tn917lac, in the presence and absence of sigma H. These genes had lower levels of expression in the absence of sigma H than in the presence of sigma H. Patterns of expression of the csh genes during growth and sporulation in wild-type and spo0H mutant cells indicated that other regulatory factors are probably involved in controlling expression of some of these genes. Three of the csh::Tn917lac insertion mutations caused noticeable phenotypes. One caused a defect in vegetative growth, but only in combination with a spo0H mutation. Two others caused a partial defect in sporulation. One of these also caused a defect in the development of genetic competence. Detailed characterization of some of the csh genes and their regulatory regions should help define the role of spo0H in the regulation of gene expression during the transition from growth to stationary phase and during the initiation of sporulation.     

DNA sequence, structure and gene expression of mycobacteriophage L5: a phage system for mycobacterial genetics

Genetic studies of Mycobacterium tuberculosis and other mycobacterial pathogens have suffered from the lack of a sophisticated genetic system. To address this issue we have developed a viral system through a detailed characterization of mycobacteriophage L5, a temperate phage that infects both fast- and slow-growing mycobacteria. We describe here the complete DNA sequence of the L5 genome and initial characterization of L5 virion structure and gene expression. In addition to providing a genetic‘toolbox’for the mycobacteria we find that L5 offers a new paradigm for dsDNA phages, being phenotypically temperate but employing genetic strategies for phage growth usually associated with lytic bacteriophages.     

Phage Therapy Is Effective against Infection by Mycobacterium ulcerans in a Murine Footpad Model

Background

Buruli Ulcer (BU) is a neglected, necrotizing skin disease caused by Mycobacterium ulcerans. Currently, there is no vaccine against M. ulcerans infection. Although the World Health Organization recommends a combination of rifampicin and streptomycin for the treatment of BU, clinical management of advanced stages is still based on the surgical resection of infected skin. The use of bacteriophages for the control of bacterial infections has been considered as an alternative or to be used in association with antibiotherapy. Additionally, the mycobacteriophage D29 has previously been shown to display lytic activity against M. ulcerans isolates.

Methodology/Principal findings

We used the mouse footpad model of M. ulcerans infection to evaluate the therapeutic efficacy of treatment with mycobacteriophage D29. Analyses of macroscopic lesions, bacterial burdens, histology and cytokine production were performed in both M. ulcerans-infected footpads and draining lymph nodes (DLN). We have demonstrated that a single subcutaneous injection of the mycobacteriophage D29, administered 33 days after bacterial challenge, was sufficient to decrease pathology and to prevent ulceration. This protection resulted in a significant reduction of M. ulcerans numbers accompanied by an increase of cytokine levels (including IFN-γ), both in footpads and DLN. Additionally, mycobacteriophage D29 treatment induced a cellular infiltrate of a lymphocytic/macrophagic profile.

Conclusions/Significance

Our observations demonstrate the potential of phage therapy against M. ulcerans infection, paving the way for future studies aiming at the development of novel phage-related therapeutic approaches against BU.     

Genome-wide identification and expression profiling of DNA methyltransferase gene family in maize

Key message

In this study, we identified eight DNA MTase genes in maize and the diversity of expression patterns of them was presented by EST mining, microarray and semi-quantitative expression profile analyses.

Abstract

DNA methylation plays a pivotal role in promoting genomic stability through diverse biological processes including regulation of gene expression during development and chromatin organization. Although this important biological process is mainly regulated by several conserved Cytosine-5 DNA methyltransferases encoded by a smaller multigene family in plants, investigation of the plant C5-MTase-encoding gene family will serve to elucidate the epigenetic mechanism diversity in plants. Recently, genome-wide identification and evolutionary analyses of the C5-MTase-encoding gene family have been characterized in multiple plant species including Arabidopsis, rice, carrot and wheat. However, little is known regarding the C5-MTase-encoding genes in the entire maize genome. Here, genome-wide identification and expression profile analyses of maize C5-MTase-encoding genes (ZmMETs) were performed from the latest version of the maize (B73) genome. Phylogenetic analysis indicated that the orthologs from the three species (maize, Arabidopsis and rice) were categorized into four classes. Chromosomal location of these genes revealed that they are unevenly distributed on 6 of all 10 chromosomes with three chromosomal/segmental duplication events, suggesting that gene duplication played a key role in expansion of the maize C5-MTase-encoding gene family. Furthermore, EST expression data mining, microarray data and semi-quantitative expression profile analyses detected in the leaves by two different abiotic stress treatments have demonstrated that these genes had temporal and spatial expression pattern and exhibited different expression levels in stress treatments, suggesting that functional diversification of ZmMET genes family. Overall, our study will serve to present signification insights to explore the plant C5-MTase-encoding gene expression and function and also be beneficial for future experimental research to further unravel the mechanisms of epigenetic regulation in plants.