Ryan S. McClure

Computational analysis of bacterial RNA-Seq data

Recent advances in high-throughput RNA sequencing (RNA-seq) have enabled tremendous leaps forward in our understanding of bacterial transcriptomes. However, computational methods for analysis of bacterial transcriptome data have not kept pace with the large and growing data sets generated by RNA-seq technology. Here, we present new algorithms, specific to bacterial gene structures and transcriptomes, for analysis of RNA-seq data. The algorithms are implemented in an open source software system called Rockhopper that supports various stages of bacterial RNA-seq data analysis, including aligning sequencing reads to a genome, constructing transcriptome maps, quantifying transcript abundance, testing for differential gene expression, determining operon structures and visualizing results. We demonstrate the performance of Rockhopper using 2.1 billion sequenced reads from 75 RNA-seq experiments conducted with Escherichia coli, Neisseria gonorrhoeae, Salmonella enterica, Streptococcus pyogenes and Xenorhabdus nematophila. We find that the transcriptome maps generated by our algorithms are highly accurate when compared with focused experimental data from E. coli and N. gonorrhoeae, and we validate our system’s ability to identify novel small RNAs, operons and transcription start sites. Our results suggest that Rockhopper can be used for efficient and accurate analysis of bacterial RNA-seq data, and that it can aid with elucidation of bacterial transcriptomes.

TLR-Dependent Human Mucosal Epithelial Cell Responses to Microbial Pathogens

Toll-like receptor (TLR) signaling represents one of the best studied pathways to implement defense mechanisms against invading microbes in human being as well as in animals. TLRs respond to specific microbial ligands and to danger signals produced by the host during infection, and initiate downstream cascades that activate both innate and adaptive immunity. TLRs are expressed by professional immune cells and by the large majority of non-hematopoietic cells, including epithelial cells. In epithelial tissues, TLR functions are particularly important because these sites are constantly exposed to microorganisms, due to their location at the host interface with the environment. While at these sites specific defense mechanisms and inflammatory responses are initiated via TLR signaling against pathogens, suppression or lack of TLR activation is also observed in response to the commensal microbiota. The mechanisms by which TLR signaling is regulated in mucosal epithelial cells include differential expression and levels of TLRs (and their signaling partners), their cellular localization and positioning within the tissue in a fashion that favors responses to pathogens while dampening responses to commensals and maintaining tissue homeostasis in physiologic conditions. In this review, the expression and activation of TLRs in mucosal epithelial cells of several sites of the human body are examined. Specifically, the oral cavity, the ear canal and eye, the airways, the gut, and the reproductive tract are discussed, along with how site-specific host defense mechanisms are implemented via TLR signaling.

Role of Hfq in iron-dependent and -independent gene regulation in Neisseria meningitidis

In Neisseria meningitidis, iron-responsive gene regulation is mediated primarily by the ferric uptake regulator (Fur) protein. When complexed with iron, Fur represses gene expression by preventing transcription initiation. Fur can also indirectly activate gene expression via the repression of regulatory small RNAs (sRNA). One such Fur- and iron-regulated sRNA, NrrF, was previously identified in N. meningitidis and shown to repress expression of the sdhA and sdhC genes encoding subunits of the succinate dehydrogenase complex. In the majority of Gram-negative bacteria, sRNA-mediated regulation requires a cofactor RNA-binding protein (Hfq) for proper gene regulation and stabilization. In this study, we examined the role of Hfq in NrrF-mediated regulation of the succinate dehydrogenase genes in N. meningitidis and the effect of an hfq mutation on iron-responsive gene regulation more broadly. We first demonstrated that the stability of NrrF, as well as the regulation of sdhC and sdhA in vivo, was unaltered in the hfq mutant. Secondly, we established that iron-responsive gene regulation of the Fur-regulated sodB gene was dependent on Hfq. Finally, we demonstrated that in N. meningitidis, Hfq functions in a global manner to control expression of many ORFs and intergenic regions via iron-independent mechanisms. Collectively these studies demonstrate that in N. meningitidis, iron- and NrrF-mediated regulation of sdhC and sdhA can occur independently of Hfq, although Hfq functions more globally to control regulation of other N. meningitidis genes primarily by iron-independent mechanisms.

The Gonococcal Transcriptome during Infection of the Lower Genital Tract in Women.

Gonorrhea is a highly prevalent disease resulting in significant morbidity worldwide, with an estimated 106 cases reported annually. Neisseria gonorrhoeae, the causative agent of gonorrhea, colonizes and infects the human genital tract and often evades host immune mechanisms until successful antibiotic treatment is used. The alarming increase in antibiotic-resistant strains of N. gonorrhoeae, the often asymptomatic nature of this disease in women and the lack of a vaccine directed at crucial virulence determinants have prompted us to perform transcriptome analysis to understand gonococcal gene expression patterns during natural infection. We sequenced RNA extracted from cervico-vaginal lavage samples collected from women recently exposed to infected male partners and determined the complete N. gonorrhoeae transcriptome during infection of the lower genital tract in women. On average, 3.19% of total RNA isolated from female samples aligned to the N. gonorrhoeae NCCP11945 genome and 1750 gonococcal ORFs (65% of all protein-coding genes) were transcribed. High expression in vivo was observed in genes encoding antimicrobial efflux pumps, iron response, phage production, pilin structure, outer membrane structures and hypothetical proteins. A parallel analysis was performed using the same strains grown in vitro in a chemically defined media (CDM). A total of 140 genes were increased in expression during natural infection compared to growth in CDM, and 165 genes were decreased in expression. Large differences were found in gene expression profiles under each condition, particularly with genes involved in DNA and RNA processing, iron, transposase, pilin and lipoproteins. We specifically interrogated genes encoding DNA binding regulators and iron-scavenging proteins, and identified increased expression of several iron-regulated genes, including tbpAB and fbpAB, during infection in women as compared to growth in vitro, suggesting that during infection of the genital tract in women, the gonococcus is exposed to an iron deplete environment. Collectively, we demonstrate that a large portion of the gonococcal genome is expressed and regulated during mucosal infection including genes involved in regulatory functions and iron scavenging.

Identification of sRNAs expressed by the human pathogen Neisseria gonorrhoeae under disparate growth conditions

In the last several years, bacterial gene regulation via small RNAs (sRNAs) has been recognized as an important mechanism controlling expression of essential proteins that are critical to bacterial growth and metabolism. Technologies such as RNA-seq are rapidly expanding the field of sRNAs and are enabling a global view of the “sRNAome” of several bacterial species. While numerous sRNAs have been identified in a variety of both Gram-negative and Gram-positive bacteria, only a very small number have been fully characterized in the human pathogen Neisseria gonorrhoeae, the etiological agent of the STD gonorrhea. Here we present the first analysis of N. gonorrhoeae specifically focused on the identification of sRNAs through RNA-seq analysis of the organism cultured under different in vitro growth conditions. Using a new computational program, Rockhopper, to analyze prokaryotic RNA-seq data obtained from N. gonorrhoeae we identified several putative sRNAs and confirmed their expression and size through Northern blot analysis. In addition, RNA was collected from four different growth conditions (iron replete and deplete, as well as with and without co-culture with human endocervical cells). Many of the putative sRNAs identified shoed varying expression levels relative to the different growth conditions examine or were detected only under certain conditions but not others. Comparisons of identified sRNAs with the regulatory pattern of putative mRNA targets revealed possible functional roles for these sRNAs. These studies are the first to carry out a global analysis of N. gonorrhoeae specifically focused on sRNAs and show that RNA-mediated regulation may be an important mechanism of gene control in this human pathogen.

Unlocking the Constraints of Cyanobacterial Productivity: Acclimations Enabling Ultrafast Growth

ABSTRACT Harnessing the metabolic potential of photosynthetic microbes for next-generation biotechnology objectives requires detailed scientific understanding of the physiological constraints and regulatory controls affecting carbon partitioning between biomass, metabolite storage pools, and bioproduct synthesis. We dissected the cellular mechanisms underlying the remarkable physiological robustness of the euryhaline unicellular cyanobacterium Synechococcus sp. strain PCC 7002 (Synechococcus 7002) and identify key mechanisms that allow cyanobacteria to achieve unprecedented photoautotrophic productivities (~2.5-h doubling time). Ultrafast growth of Synechococcus 7002 was supported by high rates of photosynthetic electron transfer and linked to significantly elevated transcription of precursor biosynthesis and protein translation machinery. Notably, no growth or photosynthesis inhibition signatures were observed under any of the tested experimental conditions. Finally, the ultrafast growth in Synechococcus 7002 was also linked to a 300% expansion of average cell volume. We hypothesize that this cellular adaptation is required at high irradiances to support higher cell division rates and reduce deleterious effects, corresponding to high light, through increased carbon and reductant sequestration. IMPORTANCE Efficient coupling between photosynthesis and productivity is central to the development of biotechnology based on solar energy. Therefore, understanding the factors constraining maximum rates of carbon processing is necessary to identify regulatory mechanisms and devise strategies to overcome productivity constraints. Here, we interrogate the molecular mechanisms that operate at a systems level to allow cyanobacteria to achieve ultrafast growth. This was done by considering growth and photosynthetic kinetics with global transcription patterns. We have delineated putative biological principles that allow unicellular cyanobacteria to achieve ultrahigh growth rates through photophysiological acclimation and effective management of cellular resource under different growth regimes. Efficient coupling between photosynthesis and productivity is central to the development of biotechnology based on solar energy. Therefore, understanding the factors constraining maximum rates of carbon processing is necessary to identify regulatory mechanisms and devise strategies to overcome productivity constraints. Here, we interrogate the molecular mechanisms that operate at a systems level to allow cyanobacteria to achieve ultrafast growth. This was done by considering growth and photosynthetic kinetics with global transcription patterns. We have delineated putative biological principles that allow unicellular cyanobacteria to achieve ultrahigh growth rates through photophysiological acclimation and effective management of cellular resource under different growth regimes.

Characterization of the Neisseria gonorrhoeae Iron and Fur Regulatory Network

UNLABELLED The Neisseria gonorrhoeae ferric uptake regulator (Fur) protein controls expression of iron homeostasis genes in response to intracellular iron levels. In this study, using transcriptome sequencing (RNA-seq) analysis of an N. gonorrhoeae fur strain, we defined the gonococcal Fur and iron regulons and characterized Fur-controlled expression of an ArsR-like DNA binding protein. We observed that 158 genes (8% of the genome) showed differential expression in response to iron in an N. gonorrhoeae wild-type or fur strain, while 54 genes exhibited differential expression in response to Fur. The Fur regulon was extended to additional regulators, including NrrF and 13 other small RNAs (sRNAs), and two transcriptional factors. One transcriptional factor, coding for an ArsR-like regulator (ArsR), exhibited increased expression under iron-replete conditions in the wild-type strain but showed decreased expression across iron conditions in the fur strain, an effect that was reversed in a fur-complemented strain. Fur was shown to bind to the promoter region of the arsR gene downstream of a predicted σ(70) promoter region. Electrophoretic mobility shift assay (EMSA) analysis confirmed binding of the ArsR protein to the norB promoter region, and sequence analysis identified two additional putative targets, NGO1411 and NGO1646. A gonococcal arsR strain demonstrated decreased survival in human endocervical epithelial cells compared to that of the wild-type and arsR-complemented strains, suggesting that the ArsR regulon includes genes required for survival in host cells. Collectively, these results demonstrate that the N. gonorrhoeae Fur functions as a global regulatory protein to repress or activate expression of a large repertoire of genes, including additional transcriptional regulatory proteins. IMPORTANCE Gene regulation in bacteria in response to environmental stimuli, including iron, is of paramount importance to both bacterial replication and, in the case of pathogenic bacteria, successful infection. Bacterial DNA binding proteins are a common mechanism utilized by pathogens to control gene expression under various environmental conditions. Here, we show that the DNA binding protein Fur, expressed by the human pathogen Neisseria gonorrhoeae, controls the expression of a large repertoire of genes and extends this regulon by controlling expression of additional DNA binding proteins. One of these proteins, an ArsR-like regulator, was required for N. gonorrhoeae survival within host cells. These results show that the Fur regulon extends to additional regulatory proteins, which together contribute to gonococcal mechanisms of pathogenesis.

Neisseria Prophage Repressor Implicated in Gonococcal Pathogenesis

ABSTRACT Neisseria gonorrhoeae, the causative agent of the sexually transmitted disease gonorrhea, can infect and colonize multiple mucosal sites in both men and women. The ability to cope with different environmental conditions requires tight regulation of gene expression. In this study, we identified and characterized a gonococcal transcriptional regulatory protein (Neisseria phage repressor [Npr]) that was previously annotated as a putative gonococcal phage repressor protein. Npr was found to repress transcription of NGNG_00460 to NGNG_00463 (NGNG_00460-00463), an operon present within the phage locus NgoΦ4. Npr binding sites within the NGNG_00460-00463 promoter region were found to overlap the −10 and −35 promoter motifs. A gonococcal npr mutant demonstrated increased adherence to and invasion of human endocervical epithelial cells compared to a wild-type gonococcal strain. Likewise, the gonococcal npr mutant exhibited enhanced colonization in a gonococcal mouse model of mucosal infection. Analysis of the gonococcal npr mutant using RNA sequence (RNA-seq) analysis demonstrated that the Npr regulon is limited to the operon present within the phage locus. Collectively, our studies have defined a new gonococcal phage repressor protein that controls the transcription of genes implicated in gonococcal pathogenesis.

Network analysis of transcriptomics expands regulatory landscapes in Synechococcus sp. PCC 7002

Cyanobacterial regulation of gene expression must contend with a genome organization that lacks apparent functional context, as the majority of cellular processes and metabolic pathways are encoded by genes found at disparate locations across the genome and relatively few transcription factors exist. In this study, global transcript abundance data from the model cyanobacterium Synechococcus sp. PCC 7002 grown under 42 different conditions was analyzed using Context-Likelihood of Relatedness (CLR). The resulting network, organized into 11 modules, provided insight into transcriptional network topology as well as grouping genes by function and linking their response to specific environmental variables. When used in conjunction with genome sequences, the network allowed identification and expansion of novel potential targets of both DNA binding proteins and sRNA regulators. These results offer a new perspective into the multi-level regulation that governs cellular adaptations of the fast-growing physiologically robust cyanobacterium Synechococcus sp. PCC 7002 to changing environmental variables. It also provides a methodological high-throughput approach to studying multi-scale regulatory mechanisms that operate in cyanobacteria. Finally, it provides valuable context for integrating systems-level data to enhance gene grouping based on annotated function, especially in organisms where traditional context analyses cannot be implemented due to lack of operon-based functional organization.

Multi-Omic Dynamics Associate Oxygenic Photosynthesis with Nitrogenase-Mediated H2 Production in Cyanothece sp. ATCC 51142

To date, the proposed mechanisms of nitrogenase-driven photosynthetic H2 production by the diazotrophic unicellular cyanobacterium Cyanothece sp. ATCC 51142 have assumed that reductant and ATP requirements are derived solely from glycogen oxidation and cyclic-electron flow around photosystem I. Through genome-scale transcript and protein profiling, this study presents and tests a new hypothesis on the metabolic relationship between oxygenic photosynthesis and nitrogenase-mediated H2 production in Cyanothece 51142. Our results show that net-positive rates of oxygenic photosynthesis and increased expression of photosystem II reaction centers correspond and are synchronized with nitrogenase expression and H2 production. These findings provide a new and more complete view on the metabolic processes contributing to the energy budget of photosynthetic H2 production and highlight the role of concurrent photocatalytic H2O oxidation as a participating process.