Identification of multivariable Boolean patterns in microbiome and microbial gene composition data

Virtually every biological system is governed by complex relations among its components. Identifying such relations requires a rigorous or heuristics-based search for patterns among variables/features of a system. Various algorithms have been developed to identify two-dimensional (involving two variables) patterns employing correlation, covariation, mutual information, etc. It seems obvious, however, that comprehensive descriptions of complex biological systems need also to include more complicated multivariable relations, which can only be described using patterns that simultaneously embrace 3, 4, and more variables. The goal of this manuscript is to (a) introduce a novel type of associations (multivariable Boolean patterns) that can be manifested between features of complex systems but cannot be identified (described) by traditional pair-vise metrics; (b) propose patterns classification method, and (c) provide a novel definition of the pattern's strength (pattern's score) able to accommodate heterogeneous multi-omics data. To demonstrate the presence of such patterns, we performed a search for all possible 2-, 3-, and 4-dimensional patterns in historical data from the Human Microbiome Project (15 body sites) and collection of H. pylori genomes associated with gastric ulcers, gastritis, and duodenal ulcers. In all datasets under consideration, we were able to identify hundreds of statistically significant multivariable patterns. These results suggest that such patterns can be common in microbial genomics/microbiomics systems.

Impairment of Tissue-Resident Mesenchymal Stem Cells in Chronic Ulcerative Colitis and Crohn's Disease

BACKGROUND AND AIMS

Little is known about the presence and function of tissue-resident mesenchymal stem cells [MtSCs] within the gastrointestinal mucosa in health and inflammatory bowel disease [IBD]. The contribution of MtSCs to the generation of inflammatory fibroblasts during IBD is also poorly understood. We hypothesized that IBD-MtSCs are impaired and contribute to the generation of the pathological myofibroblasts in IBD.

METHODS

In a cohort of clinically and endoscopically active IBD patients and normal controls, we used quantitative RT-PCR and stem cell differentiation assays, as well as confocal microscopy, to characterize MtSCs.

RESULTS

Expression of two stem cell markers, Oct4 and ALDH1A, was increased in the inflamed IBD colonic mucosa and correlated with an increase of the mesenchymal lineage marker Grem1 in ulcerative colitis [UC], but not Crohn's disease [CD]. Increased proliferation and aberrant differentiation of Oct4+Grem1+ MtSC-like cells was observed in UC, but not in CD colonic mucosa. In contrast to normal and UC-derived MtSCs, CD-MtSCs lose their clonogenic and most of their differentiation capacities. Our data also suggest that severe damage to these cells in CD may account for the pathological PD-L1low phenotype of CD myofibroblasts. In contrast, aberrant differentiation of MtSCs appears to be involved in the appearance of pathological partially differentiated PD-L1high myofibroblasts within the inflammed colonic mucosa in UC.

CONCLUSION

Our data show, for the first time, that the progenitor functions of MtSCs are differentially impaired in CD vs UC, providing a scientific rationale for the use of allogeneic MSC therapy in IBD, and particularly in CD.

Comparative transcriptomic analysis of Rickettsia conorii during in vitro infection of human and tick host cells

BACKGROUND

Pathogenic Rickettsia species belonging to the spotted fever group are arthropod-borne, obligate intracellular bacteria which exhibit preferential tropism for host microvascular endothelium in the mammalian hosts, resulting in disease manifestations attributed primarily to endothelial damage or dysfunction. Although rickettsiae are known to undergo evolution through genomic reduction, the mechanisms by which these pathogens regulate their transcriptome to ensure survival in tick vectors and maintenance by transovarial/transstadial transmission, in contrast to their ability to cause debilitating infections in human hosts remain unknown. In this study, we compare the expression profiles of rickettsial sRNAome/transcriptome and determine the transcriptional start sites (TSSs) of R. conorii transcripts during in vitro infection of human and tick host cells.

RESULTS

We performed deep sequencing on total RNA from Amblyomma americanum AAE2 cells and human microvascular endothelial cells (HMECs) infected with R. conorii. Strand-specific RNA sequencing of R. conorii transcripts revealed the expression 32 small RNAs (Rc_sR's), which were preferentially expressed above the limit of detection during tick cell infection, and confirmed the expression of Rc_sR61, sR71, and sR74 by quantitative RT-PCR. Intriguingly, a total of 305 and 132 R. conorii coding genes were differentially upregulated (> 2-fold) in AAE2 cells and HMECs, respectively. Further, enrichment for primary transcripts by treatment with Terminator 5'-Phosphate-dependent Exonuclease resulted in the identification of 3903 and 2555 transcription start sites (TSSs), including 214 and 181 primary TSSs in R. conorii during the infection to tick and human host cells, respectively. Seventy-five coding genes exhibited different TSSs depending on the host environment. Finally, we also observed differential expression of 6S RNA during host-pathogen and vector-pathogen interactions in vitro, implicating an important role for this noncoding RNA in the regulation of rickettsial transcriptome depending on the supportive host niche.

CONCLUSIONS

In sum, the findings of this study authenticate the presence of novel Rc_sR's in R. conorii, reveal the first evidence for differential expression of coding transcripts and utilization of alternate transcriptional start sites depending on the host niche, and implicate a role for 6S RNA in the regulation of coding transcriptome during tripartite host-pathogen-vector interactions.

Identification of multidimensional Boolean patterns in microbial communities

BACKGROUND

Identification of complex multidimensional interaction patterns within microbial communities is the key to understand, modulate, and design beneficial microbiomes. Every community has members that fulfill an essential function affecting multiple other community members through secondary metabolism. Since microbial community members are often simultaneously involved in multiple relations, not all interaction patterns for such microorganisms are expected to exhibit a visually uninterrupted pattern. As a result, such relations cannot be detected using traditional correlation, mutual information, principal coordinate analysis, or covariation-based network inference approaches.

RESULTS

We present a novel pattern-specific method to quantify the strength and estimate the statistical significance of two-dimensional co-presence, co-exclusion, and one-way relation patterns between abundance profiles of two organisms as well as extend this approach to allow search and visualize three-, four-, and higher dimensional patterns. The proposed approach has been tested using 2380 microbiome samples from the Human Microbiome Project resulting in body site-specific networks of statistically significant 2D patterns as well as revealed the presence of 3D patterns in the Human Microbiome Project data.

CONCLUSIONS

The presented study suggested that search for Boolean patterns in the microbial abundance data needs to be pattern specific. The reported presence of multidimensional patterns (which cannot be reduced to a combination of two-dimensional patterns) suggests that multidimensional (multi-organism) relations may play important roles in the organization of microbial communities, and their detection (and appropriate visualization) may lead to a deeper understanding of the organization and dynamics of microbial communities. Video Abstract.

Matrix metalloproteinases cleave membrane-bound PD-L1 on CD90+ (myo-)fibroblasts in Crohn's disease and regulate Th1/Th17 cell responses

Increased T helper (Th)1/Th17 immune responses are a hallmark of Crohn's disease (CD) immunopathogenesis. CD90+ (myo-)fibroblasts (MFs) are abundant cells in the normal (N) intestinal mucosa contributing to mucosal tolerance via suppression of Th1 cell activity through cell surface membrane-bound PD-L1 (mPD-L1). CD-MFs have a decreased level of mPD-L1. Consequently, mPD-L1-mediated suppression of Th1 cells by CD-MFs is decreased, yet the mechanism responsible for the reduction in mPDL-1 is unknown. Increased expression of matrix metalloproteinases (MMPs) has been reported in CD. Herein we observed that when compared to N- and ulcerative colitis (UC)-MFs, CD-MFs increase in LPS-inducible levels of MMP-7 and -9 with a significant increase in both basal and inducible MMP-10. A similar pattern of MMP expression was observed in the CD-inflamed mucosa. Treatment of N-MFs with a combination of recombinant human MMP-7, -9 and -10 significantly decreased mPD-L1. In contrast, inhibition of MMP activity with MMP inhibitors or anti-MMP-10 neutralizing antibodies restores mPD-L1 on CD-MFs. CD-MFs demonstrated reduced capacity to suppress Th1 and Th17 responses from activated CD4+ T cells. By contrast, supplementation of the CD-MF:T-cell co-cultures with MMP inhibitors or anti-MMP neutralizing antibodies restored the CD-MF-mediated suppression. Our data suggest that (i) increased MMP-10 expression by CD-MFs and concomitant cleavage of PD-L1 from the surface of CD-MFs are likely to be one of the factors contributing to the decrease of mPD-L1-mediated suppression of Th1/Th17 cells in CD; and (ii) MMPs are likely to have a significant role in the intestinal mucosal immune responses.

Measuring Alphavirus Fidelity Using Non-Infectious Virus Particles

Mutations are incorporated into the genomes of RNA viruses at an optimal frequency and altering this precise frequency has been proposed as a strategy to create live-attenuated vaccines. However, determining the effect of specific mutations that alter fidelity has been difficult because of the rapid selection of the virus population during replication. By deleting residues of the structural polyprotein PE2 cleavage site, E3Δ56-59, in Venezuelan equine encephalitis virus (VEEV) TC-83 vaccine strain, non-infectious virus particles were used to assess the effect of single mutations on mutation frequency without the interference of selection that results from multiple replication cycles. Next-generation sequencing analysis revealed a significantly lower frequency of transversion mutations and overall mutation frequency for the fidelity mutants compared to VEEV TC-83 E3Δ56-59. We demonstrate that deletion of the PE2 cleavage site halts virus infection while making the virus particles available for downstream sequencing. The conservation of the site will allow the evaluation of suspected fidelity mutants across alphaviruses of medical importance.

Disturbances in Chemokine Pathways as a Mechanism for Poor Immune Containment of MTB in the setting of HIV Co-infection

Tuberculosis (TB) is one of the ten top causes of death worldwide and the leading cause in people living with Human Immunodeficiency Virus (HIV). The loss of function by diverse leukocyte populations and poor organization of the protective granuloma structure due to HIV infection is associated with poor immunity and increased proliferation of Mtb bacilli in co-infection. We employed our humanized mouse model of Mtb and HIV co-infection to investigate the HIV-mediated disturbances in the chemokines and chemokine receptors that regulate the recruitment, activation, and localization of inflammatory leukocytes at sites of pulmonary infection. Lung tissue of mice that were non-infected or infected with Mtb, HIV, or Mtb/HIV was used for multiplex ELISA of chemokine profiles in supernatants and for high throughput RNA sequencing to identify differential transcription of chemokine and chemokine receptor genes. Our results show that Mtb infection of humanized mouse lung activates human chemokine pathways important for recruitment and effector function of leukocytes at sites of infection such as CXCL1, CXCL5, and CCL2. Co-infection with HIV, however, activated differential transcription of several chemokines and chemokine receptor families, compared to Mtb mono-infection. These outcomes were associated with greater bacterial and viral burden in co-infected animals and may identify candidate mechanisms for host directed therapies to reduce aggressive TB in those living with HIV.

Broom: application for non-redundant storage of high throughput sequencing data

Motivation

The data generation capabilities of high throughput sequencing (HTS) instruments have exponentially increased over the last few years, while the cost of sequencing has dramatically decreased allowing this technology to become widely used in biomedical studies. For small labs and individual researchers, however, storage and transfer of large amounts of HTS data present a significant challenge. The recent trends in increased sequencing quality and genome coverage can be used to reconsider HTS data storage strategies.

Results

We present Broom, a stand-alone application designed to select and store only high-quality sequencing reads at extremely high compression rates. Written in C++, the application accepts single and paired-end reads in FASTQ and FASTA formats and decompresses data in FASTA format.

Availability and implementation

C++ code available at https://scsb.utmb.edu/labgroups/fofanov/broom.asp.

Supplementary information

Supplementary data are available at Bioinformatics online.

Detection of multi-dimensional co-exclusion patterns in microbial communities

Motivation

Identification of complex relationships among members of microbial communities is key to understand and control the microbiota. Co-exclusion is arguably one of the most important patterns reflecting micro-organisms' intolerance to each other's presence. Knowing these relations opens an opportunity to manipulate microbiotas, personalize anti-microbial and probiotic treatments as well as guide microbiota transplantation. The co-exclusion pattern however, cannot be appropriately described by a linear function nor its strength be estimated using covariance or (negative) Pearson and Spearman correlation coefficients. This manuscript proposes a way to quantify the strength and evaluate the statistical significance of co-exclusion patterns between two, three or more variables describing a microbiota and allows one to extend analysis beyond micro-organism abundance by including other microbiome associated measurements such as, pH, temperature etc., as well as estimate the expected numbers of false positive co-exclusion patterns in a co-exclusion network.

Results

The implemented computational pipeline (CoEx) tested against 2380 microbial profiles (samples) from The Human Microbiome Project resulted in body-site specific pairwise co-exclusion patterns.

Availability and implementation

C++ source code for calculation of the score and P-value for two, three and four dimensional co-exclusion patterns as well as source code and executable files for the CoEx pipeline are available at https://scsb.utmb.edu/labgroups/fofanov/co-exclusion_in_microbial_communities.asp.

Supplementary information

Supplementary data are available at Bioinformatics online.

Janibacter species with evidence of genomic polymorphism isolated from resected heart valve in a patient with aortic stenosis

The authors report isolation and identification of two strains of bacteria belonging to the genus Janibacter from a human patient with aortic stenosis from a rural area of the country of Georgia. The microorganisms were isolated from aortic heart valve. Two isolates with slightly distinct colony morphologies were harvested after sub-culturing from an original agar plate. Preliminary identification of the isolates is based on amplification and sequencing of a fragment of 16SrRNA. Whole genome sequencing was performed using the Illumina MiSeq instrument. Both isolates were identified as undistinguished strains of the genus Janibacter. Characterization of whole genome sequences of each culture has revealed a 15% difference in gene profile between the cultures and confirmed that both strains belong to the genus Janibacter with the closest match to J. terrae. Genomic comparison of cultures of Janibacter obtained from human cases and from environmental sources presents a promising direction for evaluating a role of these bacteria as human pathogens.

Viral RNA-dependent RNA polymerase mutants display an altered mutation spectrum resulting in attenuation in both mosquito and vertebrate hosts

The presence of bottlenecks in the transmission cycle of many RNA viruses leads to a severe reduction of number of virus particles and this occurs multiple times throughout the viral transmission cycle. Viral replication is then necessary for regeneration of a diverse mutant swarm. It is now understood that any perturbation of the mutation frequency either by increasing or decreasing the accumulation of mutations in an RNA virus results in attenuation of the virus. To determine if altering the rate at which a virus accumulates mutations decreases the probability of a successful virus infection due to issues traversing host bottlenecks, a series of mutations in the RNA-dependent RNA polymerase of Venezuelan equine encephalitis virus (VEEV), strain 68U201, were tested for mutation rate changes. All RdRp mutants were attenuated in both the mosquito and vertebrate hosts, while showing no attenuation during in vitro infections. The rescued viruses containing these mutations showed some evidence of change in fidelity, but the phenotype was not sustained following passaging. However, these mutants did exhibit changes in the frequency of specific types of mutations. Using a model of mutation production, these changes were shown to decrease the number of stop codons generated during virus replication. This suggests that the observed mutant attenuation in vivo may be due to an increase in the number of unfit genomes, which may be normally selected against by the accumulation of stop codons. Lastly, the ability of these attenuated viruses to transition through a bottleneck in vivo was measured using marked virus clones. The attenuated viruses showed an overall reduction in the number of marked clones for both the mosquito and vertebrate hosts, as well as a reduced ability to overcome the known bottlenecks in the mosquito. This study demonstrates that any perturbation of the optimal mutation frequency whether through changes in fidelity or by alterations in the mutation frequency of specific nucleotides, has significant deleterious effects on the virus, especially in the presence of host bottlenecks.

16S rRNA gene profiling and genome reconstruction reveal community metabolic interactions and prebiotic potential of medicinal herbs used in neurodegenerative disease and as nootropics

The prebiotic potential of nervine herbal medicines has been scarcely studied. We therefore used anaerobic human fecal cultivation to investigate whether medicinal herbs commonly used as treatment in neurological health and disease in Ayurveda and other traditional systems of medicine modulate gut microbiota. Profiling of fecal cultures supplemented with either Kapikacchu, Gotu Kola, Bacopa/Brahmi, Shankhapushpi, Boswellia/Frankincense, Jatamansi, Bhringaraj, Guduchi, Ashwagandha or Shatavari by 16S rRNA sequencing revealed profound changes in diverse taxa. Principal coordinate analysis highlights that each herb drives the formation of unique microbial communities predicted to display unique metabolic potential. The relative abundance of approximately one-third of the 243 enumerated species was altered by all herbs. Additional species were impacted in an herb-specific manner. In this study, we combine genome reconstruction of sugar utilization and short chain fatty acid (SCFA) pathways encoded in the genomes of 216 profiled taxa with monosaccharide composition analysis of each medicinal herb by quantitative mass spectrometry to enhance the interpretation of resulting microbial communities and discern potential drivers of microbiota restructuring. Collectively, our results indicate that gut microbiota engage in both protein and glycan catabolism, providing amino acid and sugar substrates that are consumed by fermentative species. We identified taxa that are efficient amino acid fermenters and those capable of both amino acid and sugar fermentation. Herb-induced microbial communities are predicted to alter the relative abundance of taxa encoding SCFA (butyrate and propionate) pathways. Co-occurrence network analyses identified a large number of taxa pairs in medicinal herb cultures. Some of these pairs displayed related culture growth relationships in replicate cultures highlighting potential functional interactions among medicinal herb-induced taxa.

Global Transcriptomic Profiling of Pulmonary Gene Expression in an Experimental Murine Model of Rickettsia conorii Infection

Mediterranean spotted fever develops from an infection with Rickettsia conorii, an obligate intracellular, Gram-negative, endotheliotropic, and tick-transmitted bacterial pathogen, and is an acute, febrile illness that can progress to life-threatening complications if not diagnosed and treated early with effective antibiotics. Despite significant morbidity and mortality, little is known about changes in gene expression that determine the host responses during in vivo infection. We have investigated the transcriptional landscape of host lungs as a prominently affected organ system in an established murine model of infection by RNA-sequencing. Ingenuity pathway analysis resulted in the identification of 1332 differentially expressed genes and 292 upstream regulators. Notably, genes encoding for ubiquitin D, aconitate decarboxylase, antimicrobial peptides, calgranulins, cytokines and chemokines, and guanylate binding proteins were highly up-regulated, whereas those involved in hemoglobin biosynthesis and heme homeostasis were significantly down-regulated. Amongst response regulators, nucleotide-binding oligomerization domain-containing protein 2 and killer cell lectin-like receptors were differentially expressed, and gene clustering revealed eukaryotic initiation factor-2, oxidative phosphorylation, and ubiquitination as the predominantly activated biological pathways. Collectively, this first global transcriptomic profiling has identified R. conorii-induced regulation of novel genes and pathways in the host lungs, further in-depth investigation of which will strengthen our understanding of the pathogenesis of human rickettsioses.

Enhancer Associated Long Non-coding RNA Transcription and Gene Regulation in Experimental Models of Rickettsial Infection

Recent discovery that much of the mammalian genome does not encode protein-coding genes (PCGs) has brought widespread attention to long noncoding RNAs (lncRNAs) as a novel layer of biological regulation. Enhancer lnc (elnc) RNAs from the enhancer regions of the genome carry the capacity to regulate PCGs in cis or in trans. Spotted fever rickettsioses represent the consequence of host infection with Gram-negative, obligate intracellular bacteria in the Genus Rickettsia. Despite being implicated in the pathways of infection and inflammation, the roles of lncRNAs in host response to Rickettsia species have remained a mystery. We have profiled the expression of host lncRNAs during infection of susceptible mice with R. conorii as a model closely mimicking the pathogenesis of human spotted fever rickettsioses. RNA sequencing on the lungs of infected hosts yielded reads mapping to 74,964 non-coding RNAs, 206 and 277 of which were determined to be significantly up- and down-regulated, respectively, in comparison to uninfected controls. Following removal of short non-coding RNAs and ambiguous transcripts, remaining transcripts underwent in-depth analysis of mouse lung epigenetic signatures H3K4Me1 and H3K4Me3, active transcript markers (POLR2A, p300, CTCF), and DNaseI hypersensitivity sites to identify two potentially active and highly up-regulated elncRNAs NONMMUT013718 and NONMMUT024103. Using Hi-3C sequencing resource, we further determined that genomic loci of NONMMUT013718 and NONMMUT024103 might interact with and regulate the expression of nearby PCGs, namely Id2 (inhibitor of DNA binding 2) and Apol10b (apolipoprotein 10b), respectively. Heterologous reporter assays confirmed the activity of elncRNAs as the inducers of their predicted PCGs. In the lungs of infected mice, expression of both elncRNAs and their targets was significantly higher than mock-infected controls. Induced expression of NONMMUT013718/Id2 in murine macrophages and NONMMUT024103/Apol10b in endothelial cells was also clearly evident during R. conorii infection in vitro. Finally, shRNA mediated knock-down of NONMMUT013718 and NONMMUT024103 elncRNAs resulted in reduced expression of endogenous Id2 and Apl10b, demonstrating the regulatory roles of these elncRNAs on their target PCGs. Our results provide very first experimental evidence suggesting altered expression of pulmonary lncRNAs and elncRNA-mediated regulation of PCGs involved in immunity and during host interactions with pathogenic rickettsiae.

Novel Wolbachia strains in Anopheles malaria vectors from Sub-Saharan Africa

Background:  Wolbachia, a common insect endosymbiotic bacterium that can influence pathogen transmission and manipulate host reproduction, has historically been considered absent from the  Anopheles (An.) genera, but has recently been found in  An. gambiae s.l. populations in West Africa.  As there are numerous  Anopheles species that have the capacity to transmit malaria, we analysed a range of species across five malaria endemic countries to determine  Wolbachia prevalence rates, characterise novel  Wolbachia strains and determine any correlation between the presence of  Plasmodium,  Wolbachia and the competing bacterium  Asaia. Methods:  Anopheles adult mosquitoes were collected from five malaria-endemic countries: Guinea, Democratic Republic of the Congo (DRC), Ghana, Uganda and Madagascar, between 2013 and 2017.  Molecular analysis was undertaken using quantitative PCR, Sanger sequencing,  Wolbachia multilocus sequence typing (MLST) and high-throughput amplicon sequencing of the bacterial  16S rRNA gene.  Results: Novel  Wolbachia strains were discovered in five species:  An. coluzzii,  An. gambiae s.s.,  An. arabiensis,  An. moucheti and  An. species A, increasing the number of  Anopheles species known to be naturally infected. Variable prevalence rates in different locations were observed and novel strains were phylogenetically diverse, clustering with  Wolbachia supergroup B strains.  We also provide evidence for resident strain variants within  An. species A.  Wolbachia is the dominant member of the microbiome in  An. moucheti and  An. species A but present at lower densities in  An. coluzzii.  Interestingly, no evidence of  Wolbachia/Asaia co-infections was seen and  Asaia infection densities were shown to be variable and location dependent.  Conclusions: The important discovery of novel  Wolbachia strains in  Anopheles provides greater insight into the prevalence of resident  Wolbachia strains in diverse malaria vectors.  Novel  Wolbachia strains (particularly high-density strains) are ideal candidate strains for transinfection to create stable infections in other  Anopheles mosquito species, which could be used for population replacement or suppression control strategies.

S-plot2: Rapid Visual and Statistical Analysis of Genomic Sequences

With the daily release of data from whole genome sequencing projects, tools to facilitate comparative studies are hard-pressed to keep pace. Graphical software solutions can readily recognize synteny by measuring similarities between sequences. Nevertheless, regions of dissimilarity can prove to be equally informative; these regions may harbor genes acquired via lateral gene transfer (LGT), signify gene loss or gain, or include coding regions under strong selection. Previously, we developed the software S-plot. This tool employed an alignment-free approach for comparing bacterial genomes and generated a heatmap representing the genomes’ similarities and dissimilarities in nucleotide usage. In prior studies, this tool proved valuable in identifying genome rearrangements as well as exogenous sequences acquired via LGT in several bacterial species. Herein, we present the next generation of this tool, S-plot2. Similar to its predecessor, S-plot2 creates an interactive, 2-dimensional heatmap capturing the similarities and dissimilarities in nucleotide usage between genomic sequences (partial or complete). This new version, however, includes additional metrics for analysis, new reporting options, and integrated BLAST query functionality for the user to interrogate regions of interest. Furthermore, S-plot2 can evaluate larger sequences, including whole eukaryotic chromosomes. To illustrate some of the applications of the tool, 2 case studies are presented. The first examines strain-specific variation across the Pseudomonas aeruginosa genome and strain-specific LGT events. In the second case study, corresponding human, chimpanzee, and rhesus macaque autosomes were studied and lineage specific contributions to divergence were estimated. S-plot2 provides a means to both visually and quantitatively compare nucleotide sequences, from microbial genomes to eukaryotic chromosomes. The case studies presented illustrate just 2 potential applications of the tool, highlighting its capability to identify and investigate the variation in molecular divergence rates across sequences. S-plot2 is freely available through https://bitbucket.org/lkalesinskas/splot and is supported on the Linux and MS Windows operating systems.

Novel Wolbachia strains in Anopheles malaria vectors from Sub-Saharan Africa

Background: Wolbachia, a common insect endosymbiotic bacterium that can influence pathogen transmission and manipulate host reproduction, has historically been considered absent from the Anopheles (An.) genera, but has recently been found in An. gambiae s.l. populations.  As there are numerous Anopheles species that have the capacity to transmit malaria, we analysed a range of species to determine Wolbachia prevalence rates, characterise novel Wolbachia strains and determine any correlation between the presence of Plasmodium, Wolbachia and the competing endosymbiotic bacterium Asaia. Methods: Anopheles adult mosquitoes were collected from five malaria-endemic countries: Guinea, Democratic Republic of the Congo (DRC), Ghana, Uganda and Madagascar, between 2013 and 2017.  Molecular analysis of samples was undertaken using quantitative PCR, Sanger sequencing, Wolbachia multilocus sequence typing (MLST) and high-throughput amplicon sequencing of the bacterial 16S rRNA gene.  Results: Novel Wolbachia strains were discovered in five species: An. coluzzii, An. gambiae s.s., An. arabiensis, An. moucheti and An. species ‘A’, increasing the number of Anopheles species known to be naturally infected. Variable prevalence rates in different locations were observed and novel strains were phylogenetically diverse, clustering with Wolbachia supergroup B strains.  We also provide evidence for resident strain variants within An. species ‘A’.  Wolbachia is the dominant member of the microbiome in An. moucheti and An. species ‘A’, but present at lower densities in An. coluzzii.  Interestingly, no evidence of Wolbachia/Asaia co-infections was seen and Asaia infection densities were also shown to be variable and location dependent.  Conclusions: The important discovery of novel Wolbachia strains in Anopheles provides greater insight into the prevalence of resident Wolbachia strains in diverse malaria vectors.  Novel Wolbachia strains (particularly high-density strains) are ideal candidate strains for transinfection to create stable infections in other Anopheles mosquito species, which could be used for population replacement or suppression control strategies.

Microbiome Interaction Networks and Community Structure From Laboratory-Reared and Field-Collected Aedes aegypti, Aedes albopictus, and Culex quinquefasciatus Mosquito Vectors

Microbial interactions are an underappreciated force in shaping insect microbiome communities. Although pairwise patterns of symbiont interactions have been identified, we have a poor understanding regarding the scale and the nature of co-occurrence and co-exclusion interactions within the microbiome. To characterize these patterns in mosquitoes, we sequenced the bacterial microbiome of Aedes aegypti, Ae. albopictus, and Culex quinquefasciatus caught in the field or reared in the laboratory and used these data to generate interaction networks. For collections, we used traps that attracted host-seeking or ovipositing female mosquitoes to determine how physiological state affects the microbiome under field conditions. Interestingly, we saw few differences in species richness or microbiome community structure in mosquitoes caught in either trap. Co-occurrence and co-exclusion analysis identified 116 pairwise interactions substantially increasing the list of bacterial interactions observed in mosquitoes. Networks generated from the microbiome of Ae. aegypti often included highly interconnected hub bacteria. There were several instances where co-occurring bacteria co-excluded a third taxa, suggesting the existence of tripartite relationships. Several associations were observed in multiple species or in field and laboratory-reared mosquitoes indicating these associations are robust and not influenced by environmental or host factors. To demonstrate that microbial interactions can influence colonization of the host, we administered symbionts to Ae. aegypti larvae that either possessed or lacked their resident microbiota. We found that the presence of resident microbiota can inhibit colonization of particular bacterial taxa. Our results highlight that microbial interactions in mosquitoes are complex and influence microbiome composition.

Sequential gene regulatory events leading to glucocorticoid-evoked apoptosis of CEM human leukemic cells:interactions of MAPK, MYC and glucocorticoid pathways

Gene expression responses to glucocorticoid (GC) in the hours preceding onset of apoptosis were compared in three clones of human acute lymphoblastic leukemia CEM cells. Between 2 and 20h, all three clones showed increasing numbers of responding genes. Each clone had many unique responses, but the two responsive clones showed a group of responding genes in common, different from the resistant clone. MYC levels and the balance of activities between the three major groups of MAPKs are known important regulators of glucocorticoid-driven apoptosis in several lymphoid cell systems. Common to the two sensitive clones were changed transcript levels from genes that decrease amounts or activity of anti-apoptotic ERK/MAPK1 and JNK2/MAPK9, or of genes that increase activity of pro-apoptotic p38/MAPK14. Down-regulation of MYC and several MYC-regulated genes relevant to MAPKs also occurred in both sensitive clones. Transcriptomine comparisons revealed probable NOTCH-GC crosstalk in these cells.

Abstract 5065: Cancer derived Escherichia coli induces tumor-promoting inflammatory cytokine IL-6 in cancer associated fibroblasts (CAFs) in a NF-κB/BRD4 dependent manner

Introduction: Chronic inflammatory responses within the colonic environment are critical to colorectal cancer (CRC) tumor development. While the exact causes of CRC development are unknown, recent data suggest that dysbiosis in the colonic microbiome results in the overgrowth of bacteria from Enterebacteriaceae, Fusobacteriaceae and Bacteroidaceae taxonomic families, contributing to CRC development. Among bacteria of the Enterebacteriaceae family, adherent invasive Esherichia coli has been shown to promote tumor cell growth via production of colibactin. However, the contribution of E.coli to the tumor promoting inflammatory microenvironment (TME) during sporadic CRC development remains unknown. Cancer associated fibroblasts (CAFs) are a major component of TME and contribute to tumor-promoting inflammatory responses via NF-κB dependent production of IL-6. Herein we investigate the effect of E. coli on CAFs. We hypothesize that stimulation of CAFs by cancer derived E.coli is a key processes in supporting tumor-promoting inflammation in CRC. Methods: Bacterial DNA was extracted from tumor and adjacent normal tissue of 28 CRC patients. High throughput sequencing was done using metabarcoding of 16S rDNA for bacteria and analyzed using CLC Genomics Workbench 8.0.5 Microbial Genomics Module, SILVA v119 database for 16S. E.coli PS092717 strain, isolated from CRC tumor site, was used in co-culture experiments. Real time RT-PCR and multiplex cytokine array was used to analyze fibroblasts gene expression and secretion, respectively. Results: We observed a reduction in microbial diversity in tumor vs adjacent normal tissue, with an increased prevalence of Fusobacteria, and Enterobactereacea, but not Bacteroidacea. We then tested the effect of E.coli PS092717 on colonic primary normal fibroblasts and CAFs. Exposure of N-CMFs to E.coli for 24 h only moderately increased IL-6, while CAF responded with a greater increase in IL-6 production. Analysis of the signaling mechanism demonstrated that this process was NF-κB dependent, and was abrogated in the presence of NF-kB specific inhibitor, triptolide (20 ng/mL). Bromodomain containing protein 4 (Brd4) is an epigenetic regulator that can recruit canonical NF-κB transcription factor RelA. Therefore, we analyzed its involvement in the above processes. NF-κB mediated increase in IL-6 expression and secretion induced by E.coli PS092717 was strongly reduced in CAFs than N-CMFs in the presence of BRD4 specific inhibitor (1 µM ZL0590). A similar observation was made when a major components of the E.coli cell wall, lipopolysaccharide (LPS), was used. Conclusion Taken together our data suggests that, E.coli derived from cancer tissue and its cell wall component LPS can contribute to the tumor promoting inflammation via stimulation of IL-6 production by CAFs and this process is BRD4/NF-κB-dependent.

Citation Format: Gabriela Uribe, Russel Rourke, Romain Villeger, George Golovko, Kamil Khanipov, Zhiqing Liu, Maria Pimenova, Yuriy Fofanov, Jia Zhou, Allen R. Brasier, Irina V. Pinchuk. Cancer derived Escherichia coli induces tumor-promoting inflammatory cytokine IL-6 in cancer associated fibroblasts (CAFs) in a NF-κB/BRD4 dependent manner [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5065.

GPR43 mediates microbiota metabolite SCFA regulation of antimicrobial peptide expression in intestinal epithelial cells via activation of mTOR and STAT3

The antimicrobial peptides (AMP) produced by intestinal epithelial cells (IEC) play crucial roles in the regulation of intestinal homeostasis by controlling microbiota. Gut microbiota has been shown to promote IEC expression of RegIIIγ and certain defensins. However, the mechanisms involved are still not completely understood. In this report, we found that IEC expression levels of RegIIIγ and β-defensins 1, 3, and 4 were lower in G protein-coupled receptor (GPR)43-/- mice compared to that of wild-type (WT) mice. Oral feeding with short-chain fatty acids (SCFA) promoted IEC production of RegIIIγ and defensins in mice. Furthermore, SCFA induced RegIIIγ and β-defensins in intestinal epithelial enteroids generated from WT but not GPR43-/- mice. Mechanistically, SCFA activated mTOR and STAT3 in IEC, and knockdown of mTOR and STAT3 impaired SCFA induction of AMP production. Our studies thus demonstrated that microbiota metabolites SCFA promoted IEC RegIIIγ and β-defensins in a GPR43-dependent manner. The data thereby provide a novel pathway by which microbiota regulates IEC expression of AMP and intestinal homeostasis.

Low-fidelity Venezuelan equine encephalitis virus polymerase mutants to improve live-attenuated vaccine safety and efficacy

During RNA virus replication, there is the potential to incorporate mutations that affect virulence or pathogenesis. For live-attenuated vaccines, this has implications for stability, as replication may result in mutations that either restore the wild-type phenotype via reversion or compensate for the attenuating mutations by increasing virulence (pseudoreversion). Recent studies have demonstrated that altering the mutation rate of an RNA virus is an effective attenuation tool. To validate the safety of low-fidelity mutations to increase vaccine attenuation, several mutations in the RNA-dependent RNA-polymerase (RdRp) were tested in the live-attenuated Venezuelan equine encephalitis virus vaccine strain, TC-83. Next generation sequencing after passage in the presence of mutagens revealed a mutant containing three mutations in the RdRp, TC-83 3x, to have decreased replication fidelity, while a second mutant, TC-83 4x displayed no change in fidelity, but shared many phenotypic characteristics with TC-83 3x. Both mutants exhibited increased, albeit inconsistent attenuation in an infant mouse model, as well as increased immunogenicity and complete protection against lethal challenge of an adult murine model compared with the parent TC-83. During serial passaging in a highly permissive model, the mutants increased in virulence but remained less virulent than the parent TC-83. These results suggest that the incorporation of low-fidelity mutations into the RdRp of live-attenuated vaccines for RNA viruses can confer increased immunogenicity whilst showing some evidence of increased attenuation. However, while in theory such constructs may result in more effective vaccines, the instability of the vaccine phenotype decreases the likelihood of this being an effective vaccine strategy.

Mosquito bottlenecks alter viral mutant swarm in a tissue and time-dependent manner with contraction and expansion of variant positions and diversity

Viral diversity is theorized to play a significant role during virus infections, particularly for arthropod-borne viruses (arboviruses) that must infect both vertebrate and invertebrate hosts. To determine how viral diversity influences mosquito infection and dissemination Culex taeniopus mosquitoes were infected with the Venezuelan equine encephalitis virus endemic strain 68U201. Bodies and legs/wings of the mosquitoes were collected individually and subjected to multi-parallel sequencing. Virus sequence diversity was calculated for each tissue. Greater diversity was seen in mosquitoes with successful dissemination versus those with no dissemination. Diversity across time revealed that bottlenecks influence diversity following dissemination to the legs/wings, but levels of diversity are restored by Day 12 post-dissemination. Specific minority variants were repeatedly identified across the mosquito cohort, some in nearly every tissue and time point, suggesting that certain variants are important in mosquito infection and dissemination. This study demonstrates that the interaction between the mosquito and the virus results in changes in diversity and the mutational spectrum and may be essential for successful transition of the bottlenecks associated with arbovirus infection.

Assessing the utility of metabarcoding for diet analyses of the omnivorous wild pig (Sus scrofa)

Wild pigs (Sus scrofa) are an invasive species descended from both domestic swine and Eurasian wild boar that was introduced to North America during the early 1500s. Wild pigs have since become the most abundant free‐ranging exotic ungulate in the United States. Large and ever‐increasing populations of wild pigs negatively impact agriculture, sport hunting, and native ecosystems with costs estimated to exceed $1.5 billion/year within the United States. Wild pigs are recognized as generalist feeders, able to exploit a broad array of locally available food resources, yet their feeding behaviors remain poorly understood as partially digested material is often unidentifiable through traditional stomach content analyses. To overcome the limitation of stomach content analyses, we developed a DNA sequencing‐based protocol to describe the plant and animal diet composition of wild pigs. Additionally, we developed and evaluated blocking primers to reduce the amplification and sequencing of host DNA, thus providing greater returns of sequences from diet items. We demonstrate that the use of blocking primers produces significantly more sequencing reads per sample from diet items, which increases the robustness of ascertaining animal diet composition with molecular tools. Further, we show that the overall plant and animal diet composition is significantly different between the three areas sampled, demonstrating this approach is suitable for describing differences in diet composition among the locations.

Nasopharyngeal microbiota in infants and changes during viral upper respiratory tract infection and acute otitis media

Background

Interferences between pathogenic bacteria and specific commensals are known. We determined the interactions between nasopharyngeal microbial pathogens and commensals during viral upper respiratory tract infection (URI) and acute otitis media (AOM) in infants.

Methods

We analyzed 971 specimens collected monthly and during URI and AOM episodes from 139 infants. The 16S rRNA V4 gene regions were sequenced on the Illumina MiSeq platform.

Results

Among the high abundant genus-level nasopharyngeal microbiota were Moraxella, Haemophilus, and Streptococcus (3 otopathogen genera), Corynebacterium, Dolosigranulum, Staphylococcus, Acinetobacter, Pseudomonas, and Bifidobacterium. Bacterial diversity was lower in culture-positive samples for Streptococcus pneumoniae, and Haemophilus influenzae, compared to cultured-negative samples. URI frequencies were positively associated with increasing trend in otopathogen colonization. AOM frequencies were associated with decreasing trend in Micrococcus colonization. During URI and AOM, there were increases in abundance of otopathogen genera and decreases in Pseudomonas, Myroides, Yersinia, and Sphingomonas. Otopathogen abundance was increased during symptomatic viral infection, but not during asymptomatic infection. The risk for AOM complicating URI was reduced by increased abundance of Staphylococcus and Sphingobium.

Conclusion

Otopathogen genera played the key roles in URI and AOM occurrences. Staphylococcus counteracts otopathogens thus Staphylococcal colonization may be beneficial, rather than harmful. While Sphingobium may play a role in preventing AOM complicating URI, the commonly used probiotic Bifidobacterium did not play a significant role during URI or AOM. The role of less common commensals in counteracting the deleterious effects of otopathogens requires further studies.

Transcriptional profiling of Rickettsia prowazekii coding and non-coding transcripts during in vitro host-pathogen and vector-pathogen interactions

Natural pathogen transmission of Rickettsia prowazekii, the etiologic agent of epidemic typhus, to humans is associated with arthropods, including human body lice, ticks, and ectoparasites of eastern flying squirrel. Recently, we have documented the presence of small RNAs in Rickettsia species and expression of R. prowazekii sRNAs during infection of cultured human microvascular endothelial cells (HMECs), which represent the primary target cells during human infections. Bacterial noncoding transcripts are now well established as critical post-transcriptional regulators of virulence and adaptation mechanisms in varying host environments. Despite their importance, little is known about the expression profile and regulatory activities of R. prowazekii sRNAs (Rp_sRs) in different host cells encountered as part of the natural life-cycle. To investigate the sRNA expression profile of R. prowazekii during infection of arthropod host cells, we employed an approach combining in vitro infection, bioinformatics, RNA sequencing, and PCR-based quantitation. Global analysis of R. prowazekii transcriptome by strand-specific RNA sequencing enabled us to identify 67 cis-acting (antisense) and 26 trans-acting (intergenic) Rp_sRs expressed during the infection of Amblyomma americanum (AAE2) cells. Comparative evaluation of expression during R. prowazekii infection of HMECs and AAE2 cells by quantitative RT-PCR demonstrated significantly higher expression of four selected Rp_sRs in tick AAE2 cells. Examination of the coding transcriptome revealed differential up-regulation of >150 rickettsial genes in either HMECs or AAE2 cells and yielded evidence for host cell-dependent utilization of alternative transcription start sites by 18 rickettsial genes. Our results thus suggest noticeable differences in the expression of both Rp_sRs as well as the coding transcriptome and the exploitation of multiple transcription initiation sites for select genes during the infection of human endothelium and tick vector cells as the host and yield new insights into rickettsial virulence and transmission mechanisms.

The ability of human nuclear DNA to cause false positive low-abundance heteroplasmy calls varies across the mitochondrial genome

Background

Low-abundance mutations in mitochondrial populations (mutations with minor allele frequency ≤ 1%), are associated with cancer, aging, and neurodegenerative disorders. While recent progress in high-throughput sequencing technology has significantly improved the heteroplasmy identification process, the ability of this technology to detect low-abundance mutations can be affected by the presence of similar sequences originating from nuclear DNA (nDNA). To determine to what extent nDNA can cause false positive low-abundance heteroplasmy calls, we have identified mitochondrial locations of all subsequences that are common or similar (one mismatch allowed) between nDNA and mitochondrial DNA (mtDNA).

Results

Performed analysis revealed up to a 25-fold variation in the lengths of longest common and longest similar (one mismatch allowed) subsequences across the mitochondrial genome. The size of the longest subsequences shared between nDNA and mtDNA in several regions of the mitochondrial genome were found to be as low as 11 bases, which not only allows using these regions to design new, very specific PCR primers, but also supports the hypothesis of the non-random introduction of mtDNA into the human nuclear DNA.

Conclusion

Analysis of the mitochondrial locations of the subsequences shared between nDNA and mtDNA suggested that even very short (36 bases) single-end sequencing reads can be used to identify low-abundance variation in 20.4% of the mitochondrial genome. For longer (76 and 150 bases) reads, the proportion of the mitochondrial genome where nDNA presence will not interfere found to be 44.5 and 67.9%, when low-abundance mutations at 100% of locations can be identified using 417 bases long single reads. This observation suggests that the analysis of low-abundance variations in mitochondria population can be extended to a variety of large data collections such as NCBI Sequence Read Archive, European Nucleotide Archive, The Cancer Genome Atlas, and International Cancer Genome Consortium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3375-x) contains supplementary material, which is available to authorized users.

Identification and Characterization of Novel Small RNAs in Rickettsia prowazekii

Emerging evidence implicates a critically important role for bacterial small RNAs (sRNAs) as post-transcriptional regulators of physiology, metabolism, stress/adaptive responses, and virulence, but the roles of sRNAs in pathogenic Rickettsia species remain poorly understood. Here, we report on the identification of both novel and well-known bacterial sRNAs in Rickettsia prowazekii, known to cause epidemic typhus in humans. RNA sequencing of human microvascular endothelial cells (HMECs), the preferred targets during human rickettsioses, infected with R. prowazekii revealed the presence of 35 trans-acting and 23 cis-acting sRNAs, respectively. Of these, expression of two trans-acting (Rp_sR17 and Rp_sR60) and one cis-acting (Rp_sR47) novel sRNAs and four well-characterized bacterial sRNAs (RNaseP_bact_a, α-tmRNA, 4.5S RNA, 6S RNA) was further confirmed by Northern blot or RT-PCR analyses. The transcriptional start sites of five novel rickettsial sRNAs and 6S RNA were next determined using 5′ RLM-RACE yielding evidence for their independent biogenesis in R. prowazekii. Finally, computational approaches were employed to determine the secondary structures and potential mRNA targets of novel sRNAs. Together, these results establish the presence and expression of sRNAs in R. prowazekii during host cell infection and suggest potential functional roles for these important post-transcriptional regulators in rickettsial biology and pathogenesis.

Bartonella melophagi in blood of domestic sheep (Ovis aries) and sheep keds (Melophagus ovinus) from the southwestern US: Cultures, genetic characterization, and ecological connections

Bartonella melophagi sp. nov. was isolated from domestic sheep blood and from sheep keds (Melophagus ovinus) from the southwestern United States. The sequence analyses of the reference strain performed by six molecular markers consistently demonstrated that B. melophagi relates to but differ from other Bartonella species isolated from domestic and wild ruminants. Presence of 183 genes specific for B. melophagi, being absent in genomes of other Bartonella species associated with ruminants also supports the separation of this bacterial species from species of other ruminants. Bartonella DNA was detected in all investigated sheep keds; however, culturing of these bacteria from sheep blood rejects a speculation that B. melophagi is an obligatory endosymbiont. Instead, the results support the hypothesis that the domestic sheep is a natural host reservoir for B. melophagi and the sheep ked its main vector. This bacterium was not isolated from the blood of bighorn sheep and domestic goats belonging to the same subfamily Caprinae. B. melophagi has also been shown to be zoonotic and needs to be investigated further.

Bacterial small RNAs in the Genus Rickettsia

Background

Rickettsia species are obligate intracellular Gram-negative pathogenic bacteria and the etiologic agents of diseases such as Rocky Mountain spotted fever (RMSF), Mediterranean spotted fever, epidemic typhus, and murine typhus. Genome sequencing revealed that R. prowazekii has ~25 % non-coding DNA, the majority of which is thought to be either “junk DNA” or pseudogenes resulting from genomic reduction. These characteristics also define other Rickettsia genomes. Bacterial small RNAs, whose biogenesis is predominantly attributed to either the intergenic regions (trans-acting) or to the antisense strand of an open reading frame (cis-acting), are now appreciated to be among the most important post-transcriptional regulators of bacterial virulence and growth. We hypothesize that intergenic regions in rickettsial species encode for small, non-coding RNAs (sRNAs) involved in the regulation of its transcriptome, leading to altered virulence and adaptation depending on the host niche.

Results

We employed a combination of bioinformatics and in vitro approaches to explore the presence of sRNAs in a number of species within Genus Rickettsia. Using the sRNA Identification Protocol using High-throughput Technology (SIPHT) web interface, we predicted over 1,700 small RNAs present in the intergenic regions of 16 different strains representing 13 rickettsial species. We further characterized novel sRNAs from typhus (R. prowazekii and R. typhi) and spotted fever (R. rickettsii and R. conorii) groups for their promoters and Rho-independent terminators using Bacterial Promoter Prediction Program (BPROM) and TransTermHP prediction algorithms, respectively. Strong σ70 promoters were predicted upstream of all novel small RNAs, indicating the potential for transcriptional activity. Next, we infected human microvascular endothelial cells (HMECs) with R. prowazekii for 3 h and 24 h and performed Next Generation Sequencing to experimentally validate the expression of 26 sRNA candidates predicted in R. prowazekii. Reverse transcriptase PCR was also used to further verify the expression of six putative novel sRNA candidates in R. prowazekii.

Conclusions

Our results yield clear evidence for the expression of novel R. prowazekii sRNA candidates during infection of HMECs. This is the first description of novel small RNAs for a highly pathogenic species of Rickettsia, which should lead to new insights into rickettsial virulence and adaptation mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2293-7) contains supplementary material, which is available to authorized users.

Survey of viral populations within Lake Michigan nearshore waters at four Chicago area beaches

In comparison to the oceans, freshwater environments represent a more diverse community of microorganisms, exhibiting comparatively high levels of variability both temporally and spatially Maranger and Bird, Microb. Ecol. 31 (1996) 141–151. This level of variability is likely to extend to the world of viruses as well, in particular bacteria-infecting viruses (bacteriophages). Phages are known to influence bacterial diversity, and therefore key processes, in environmental niches across the globe Clokie et al., Bacteriophage 1 (2011) 31–45; Jacquet et al., Adv. Ocean Limn. 1 (2010) 97–141; Wilhelm and Suttle, Bioscience 49 (1999) 781–788; Bratback et al., Microb. Ecol. 28 (1994) 209–221. Despite their prevalence and likely critical role in freshwater environments, very few viral species have been characterized. Metagenomic approaches, however, have allowed for a glimpse into phage diversity. We collected surface water samples from four Chicago area beaches – Gillson Park, Montrose Beach, 57th Street Beach, and Calumet Beach – every two weeks from May 13 through August 5, 2014. Sampling was conducted with four biological replicates for each sampling date and location, resulting in 112 samples. DNA isolated from each of the individual samples for a given collection date/location was pooled together, with one exception – Calumet Beach on August 5, 2014 – in which each biological replicate was sequenced individually. Raw sequence data is available via NCBI’s SRA database (part of BioProject PRJNA248239).

Secondary Analysis of the NCI-60 Whole Exome Sequencing Data Indicates Significant Presence of Propionibacterium acnes Genomic Material in Leukemia (RPMI-8226) and Central Nervous System (SF-295, SF-539, and SNB-19) Cell Lines

The NCI-60 human tumor cell line panel has been used in a broad range of cancer research over the last two decades. A landmark 2013 whole exome sequencing study of this panel added an exceptional new resource for cancer biologists. The complementary analysis of the sequencing data produced by this study suggests the presence of Propionibacterium acnes genomic sequences in almost half of the datasets, with the highest abundance in the leukemia (RPMI-8226) and central nervous system (SF-295, SF-539, and SNB-19) cell lines. While the origin of these contaminating bacterial sequences remains to be determined, observed results suggest that computational control for the presence of microbial genomic material is a necessary step in the analysis of the high throughput sequencing (HTS) data.

Mycoplasma CG- and GATC-specific DNA methyltransferases selectively and efficiently methylate the host genome and alter the epigenetic landscape in human cells

Aberrant DNA methylation is frequently observed in disease, including many cancer types, yet the underlying mechanisms remain unclear. Because germline and somatic mutations in the genes that are responsible for DNA methylation are infrequent in malignancies, additional mechanisms must be considered. Mycoplasmas spp., including Mycoplasma hyorhinis, efficiently colonize human cells and may serve as a vehicle for delivery of enzymatically active microbial proteins into the intracellular milieu. Here, we performed, for the first time, genome-wide and individual gene mapping of methylation marks generated by the M. hyorhinis CG- and GATC-specific DNA cytosine methyltransferases (MTases) in human cells. Our results demonstrated that, upon expression in human cells, MTases readily translocated to the cell nucleus. In the nucleus, MTases selectively and efficiently methylated the host genome at the DNA sequence sites free from pre-existing endogenous methylation, including those in a variety of cancer-associated genes. We also established that mycoplasma is widespread in colorectal cancers, suggesting that either the infection contributed to malignancy onset or, alternatively, that tumors provide a favorable environment for mycoplasma growth. In the human genome, ∼ 11% of GATC sites overlap with CGs (e.g., CGAT(m)CG); therefore, the methylated status of these sites can be perpetuated by human DNMT1. Based on these results, we now suggest that the GATC-specific methylation represents a novel type of infection-specific epigenetic mark that originates in human cells with a previous exposure to infection. Overall, our findings unveil an entirely new panorama of interactions between the human microbiome and epigenome with a potential impact in disease etiology.

Metagenomic analysis of the airborne environment in urban spaces

The organisms in aerosol microenvironments, especially densely populated urban areas, are relevant to maintenance of public health and detection of potential epidemic or biothreat agents. To examine aerosolized microorganisms in this environment, we performed sequencing on the material from an urban aerosol surveillance program. Whole metagenome sequencing was applied to DNA extracted from air filters obtained during periods from each of the four seasons. The composition of bacteria, plants, fungi, invertebrates, and viruses demonstrated distinct temporal shifts. Bacillus thuringiensis serovar kurstaki was detected in samples known to be exposed to aerosolized spores, illustrating the potential utility of this approach for identification of intentionally introduced microbial agents. Together, these data demonstrate the temporally dependent metagenomic complexity of urban aerosols and the potential of genomic analytical techniques for biosurveillance and monitoring of threats to public health.

The presence of nitrate dramatically changed the predominant microbial community in perchlorate degrading cultures under saline conditions

Background

Perchlorate contamination has been detected in both ground water and drinking water. An attractive treatment option is the use of ion-exchange to remove and concentrate perchlorate in brine. Biological treatment can subsequently remove the perchlorate from the brine. When nitrate is present, it will also be concentrated in the brine and must also be removed by biological treatment. The primary objective was to obtain an in-depth characterization of the microbial populations of two salt-tolerant cultures each of which is capable of metabolizing perchlorate. The cultures were derived from a single ancestral culture and have been maintained in the laboratory for more than 10 years. One culture was fed perchlorate only, while the other was fed both perchlorate and nitrate.

Results

A metagenomic characterization was performed using Illumina DNA sequencing technology, and the 16S rDNA of several pure strains isolated from the mixed cultures were sequenced. In the absence of nitrate, members of the Rhodobacteraceae constituted the prevailing taxonomic group. Second in abundance were the Rhodocyclaceae. In the nitrate fed culture, the Rhodobacteraceae are essentially absent. They are replaced by a major expansion of the Rhodocyclaceae and the emergence of the Alteromonadaceae as a significant community member. Gene sequences exhibiting significant homology to known perchlorate and nitrate reduction enzymes were found in both cultures.

Conclusions

The structure of the two microbial ecosystems of interest has been established and some representative strains obtained in pure culture. The results illustrate that under favorable conditions a group of organisms can readily dominate an ecosystem and yet be effectively eliminated when their advantage is lost. Almost all known perchlorate-reducing organisms can also effectively reduce nitrate. This is certainly not the case for the Rhodobacteraceae that were found to dominate in the absence of nitrate, but effectively disappeared in its presence. This study is significant in that it reveals the existence of a novel group of organisms that play a role in the reduction of perchlorate under saline conditions. These Rhodobacteraceae especially, as well as other organisms present in these communities may be a promising source of unique salt-tolerant enzymes for perchlorate reduction.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-014-0225-3) contains supplementary material, which is available to authorized users.

Epigenetic alteration by DNA-demethylating treatment restores apoptotic response to glucocorticoids in dexamethasone-resistant human malignant lymphoid cells

Background

Glucocorticoids (GCs) are often included in the therapy of lymphoid malignancies because they kill several types of malignant lymphoid cells. GCs activate the glucocorticoid receptor (GR), to regulate a complex genetic network, culminating in apoptosis. Normal lymphoblasts and many lymphoid malignancies are sensitive to GC-driven apoptosis. Resistance to GCs can be a significant clinical problem, however, and correlates with resistance to several other major chemotherapeutic agents.

Methods

We analyzed the effect of treatment with the cytosine analogue 5 aza-2’ deoxycytidine (AZA) on GC resistance in two acute lymphoblastic leukemia (T or pre-T ALL) cell lines- CEM and Molt-4- and a (B-cell) myeloma cell line, RPMI 8226. Methods employed included tissue culture, flow cytometry, and assays for clonogenicity, cytosine extension, immunochemical identification of proteins, and gene transactivation. High throughput DNA sequencing was used to confirm DNA methylation status.

Conclusions

Treatment of these cells with AZA resulted in altered DNA methylation and restored GC-evoked apoptosis in all 3 cell lines. In CEM cells the altered epigenetic state resulted in site-specific phosphorylation of the GR, increased GR potency, and GC-driven induction of the GR from promoters that lie in CpG islands. In RPMI 8226 cells, expression of relevant coregulators of GR function was altered. Activation of p38 mitogen-activated protein kinase (MAPK), which is central to a feed-forward mechanism of site-specific GR phosphorylation and ultimately, apoptosis, occurred in all 3 cell lines. These data show that in certain malignant hematologic B- and T-cell types, epigenetically controlled GC resistance can be reversed by cell exposure to a compound that causes DNA demethylation. The results encourage studies of application to in vivo systems, looking towards eventual clinical applications.

Where does Neisseria acquire foreign DNA from: an examination of the source of genomic and pathogenic islands and the evolution of the Neisseria genus

Background

Pathogenicity islands (PAIs) or genomic islands (GEIs) are considered to be the result of a recent horizontal transfer. Detecting PAIs/GEIs as well as their putative source can provide insight into the organism’s pathogenicity within its host. Previously we introduced a tool called S-plot which provides a visual representation of the variation in compositional properties across and between genomic sequences. Utilizing S-plot and new functionality developed here, we examined 18 publicly available Neisseria genomes, including strains of both pathogenic and non-pathogenic species, in order to identify regions of unusual compositional properties (RUCPs) using both a sliding window as well as a gene-by-gene approach.

Results

Numerous GEIs and PAIs were identified including virulence genes previously found within the pathogenic Neisseria species. While some genes were conserved amongst all species, only pathogenic species, or an individual species, a number of genes were detected that are unique to an individual strain. While the majority of such genes have an origin unknown, a number of putative sources including pathogenic and capsule-containing bacteria were determined, indicative of gene exchange between Neisseria spp. and other bacteria within their microhabitat. Furthermore, we uncovered evidence that both N. meningitidis and N. gonorrhoeae have separately acquired DNA from their human host. Data suggests that all three Neisseria species have received horizontally transferred elements post-speciation.

Conclusions

Using this approach, we were able to not only find previously identified regions of virulence but also new regions which may be contributing to the virulence of the species. This comparative analysis provides a means for tracing the evolutionary history of the acquisition of foreign DNA within this genus. Looking specifically at the RUCPs present within the 18 genomes considered, a stronger similarity between N. meningitidis and N. lactamica is observed, suggesting that N. meningitidis arose before N. gonorrhoeae.

Whole genome sequencing of mutation accumulation lines reveals a low mutation rate in the social amoeba Dictyostelium discoideum

Spontaneous mutations play a central role in evolution. Despite their importance, mutation rates are some of the most elusive parameters to measure in evolutionary biology. The combination of mutation accumulation (MA) experiments and whole-genome sequencing now makes it possible to estimate mutation rates by directly observing new mutations at the molecular level across the whole genome. We performed an MA experiment with the social amoeba Dictyostelium discoideum and sequenced the genomes of three randomly chosen lines using high-throughput sequencing to estimate the spontaneous mutation rate in this model organism. The mitochondrial mutation rate of 6.76×10(-9), with a Poisson confidence interval of 4.1×10(-9) - 9.5×10(-9), per nucleotide per generation is slightly lower than estimates for other taxa. The mutation rate estimate for the nuclear DNA of 2.9×10(-11), with a Poisson confidence interval ranging from 7.4×10(-13) to 1.6×10(-10), is the lowest reported for any eukaryote. These results are consistent with low microsatellite mutation rates previously observed in D. discoideum and low levels of genetic variation observed in wild D. discoideum populations. In addition, D. discoideum has been shown to be quite resistant to DNA damage, which suggests an efficient DNA-repair mechanism that could be an adaptation to life in soil and frequent exposure to intracellular and extracellular mutagenic compounds. The social aspect of the life cycle of D. discoideum and a large portion of the genome under relaxed selection during vegetative growth could also select for a low mutation rate. This hypothesis is supported by a significantly lower mutation rate per cell division in multicellular eukaryotes compared with unicellular eukaryotes.

Slim-filter: an interactive Windows-based application for illumina genome analyzer data assessment and manipulation

Background

The emergence of Next Generation Sequencing technologies has made it possible for individual investigators to generate gigabases of sequencing data per week. Effective analysis and manipulation of these data is limited due to large file sizes, so even simple tasks such as data filtration and quality assessment have to be performed in several steps. This requires (potentially problematic) interaction between the investigator and a bioinformatics/computational service provider. Furthermore, such services are often performed using specialized computational facilities.

Results

We present a Windows-based application, Slim-Filter designed to interactively examine the statistical properties of sequencing reads produced by Illumina Genome Analyzer and to perform a broad spectrum of data manipulation tasks including: filtration of low quality and low complexity reads; filtration of reads containing undesired subsequences (such as parts of adapters and PCR primers used during the sample and sequencing libraries preparation steps); excluding duplicated reads (while keeping each read’s copy number information in a specialized data format); and sorting reads by copy numbers allowing for easy access and manual editing of the resulting files. Slim-Filter is organized as a sequence of windows summarizing the statistical properties of the reads. Each data manipulation step has roll-back abilities, allowing for return to previous steps of the data analysis process. Slim-Filter is written in C++ and is compatible with fasta, fastq, and specialized AS file formats presented in this manuscript. Setup files and a user’s manual are available for download at the supplementary web site ( https://www.bioinfo.uh.edu/Slim_Filter/).

Conclusion

The presented Windows-based application has been developed with the goal of providing individual investigators with integrated sequencing reads analysis, curation, and manipulation capabilities.

Detection and typing of viruses using broadly sensitive cocktail-PCR and mass spectrometric cataloging: demonstration with dengue virus

Virus detection and taxonomic identification of serotypes, strains, or genotypes provide important information relevant for diagnosis, and for the epidemiological characterization and tracking of new strains in an endemic region. In the specific case of dengue virus, rapid serotype identification can also be useful in the treatment of secondary infections that may cause the more severe dengue hemorrhagic fever and dengue shock syndrome. In this work, dengue virus was used as a model to test a new approach of combining broadly sensitive RT-PCR amplification of nearly any virus strain with subsequent serotype- and finer-level identification by mass spectrometry. PCR primers were appended with promoter sequences, such that the resulting PCR products could be transcribed into RNA. RNA fragments generated by guanosine-specific RNase T(1) digestion were analyzed by matrix-assisted laser desorption/ionization-time of flight mass spectrometry. Viral serotypes were identified by comparing the pattern of observed fragment masses to a mass database. The database was created by computationally fragmenting 2517 dengue strains after each guanosine residue using the same primers. Computationally, all 2517 strains in the mass database were correctly identified at the serotype level from the predicted PCR product. The methodology was successfully demonstrated experimentally by identifying the serotypes of eight test strains using mosquito cell cultures infected with strains of all four serotypes and with full-length cDNA clones.

PCR detection of nearly any dengue virus strain using a highly sensitive primer 'cocktail'

PCR detection of viral pathogens is extremely useful, but suffers from the challenge of detecting the many variant strains of a given virus that arise over time. Here, we report the computational derivation and initial experimental testing of a combination of 10 PCR primers to be used in a single high-sensitivity mixed PCR reaction for the detection of dengue virus. Primer sequences were computed such that their probability of mispriming with human DNA is extremely low. A 'cocktail' of 10 primers was shown experimentally to be able to detect cDNA clones representing the four serotypes and dengue virus RNA spiked into total human whole blood RNA. Computationally, the primers are predicted to detect 95% of the 1688 dengue strains analyzed (with perfect primer match). Allowing up to one mismatch and one insertion per primer, the primer set detects 99% of strains. Primer sets from three previous studies have been compared with the present set of primers and their relative sensitivity for dengue virus is discussed. These results provide the formulation and demonstration of a mixed primer PCR reagent that may enable the detection of nearly any dengue strain irrespective of serotype, in a single PCR reaction, and illustrate an approach to the broad problem of detecting highly mutable RNA viruses.

PIQA: pipeline for Illumina G1 genome analyzer data quality assessment

Summary: PIQA is a quality analysis pipeline designed to examine genomic reads produced by Next Generation Sequencing technology (Illumina G1 Genome Analyzer). A short statistical summary, as well as tile-by-tile and cycle-by-cycle graphical representation of clusters density, quality scores and nucleotide frequencies allow easy identification of various technical problems including defective tiles, mistakes in sample/library preparations and abnormalities in the frequencies of appearance of sequenced genomic reads. PIQA is written in the R statistical programming language and is compatible with bustard, fastq and scarf Illumina G1 Genome Analyzer data formats.

Availability: The PIQA pipeline, installation instructions and examples are available at the supplementary web site (http://bioinfo.uh.edu/PIQA).

Contact:yfofanov@bioinfo.uh.edu

Ultraspecific probes for high throughput HLA typing

Background

The variations within an individual's HLA (Human Leukocyte Antigen) genes have been linked to many immunological events, e.g. susceptibility to disease, response to vaccines, and the success of blood, tissue, and organ transplants. Although the microarray format has the potential to achieve high-resolution typing, this has yet to be attained due to inefficiencies of current probe design strategies.

Results

We present a novel three-step approach for the design of high-throughput microarray assays for HLA typing. This approach first selects sequences containing the SNPs present in all alleles of the locus of interest and next calculates the number of base changes necessary to convert a candidate probe sequences to the closest subsequence within the set of sequences that are likely to be present in the sample including the remainder of the human genome in order to identify those candidate probes which are "ultraspecific" for the allele of interest. Due to the high specificity of these sequences, it is possible that preliminary steps such as PCR amplification are no longer necessary. Lastly, the minimum number of these ultraspecific probes is selected such that the highest resolution typing can be achieved for the minimal cost of production. As an example, an array was designed and in silico results were obtained for typing of the HLA-B locus.

Conclusion

The assay presented here provides a higher resolution than has previously been developed and includes more alleles than previously considered. Based upon the in silico and preliminary experimental results, we believe that the proposed approach can be readily applied to any highly polymorphic gene system.

Effect of the mutation rate and background size on the quality of pathogen identification

MOTIVATION

Genomic-based methods have significant potential for fast and accurate identification of organisms or even genes of interest in complex environmental samples (air, water, soil, food, etc.), especially when isolation of the target organism cannot be performed by a variety of reasons. Despite this potential, the presence of the unknown, variable and usually large quantities of background DNA can cause interference resulting in false positive outcomes.

RESULTS

In order to estimate how the genomic diversity of the background (total length of all of the different genomes present in the background), target length and target mutation rate affect the probability of misidentifications, we introduce a mathematical definition for the quality of an individual signature in the presence of a background based on its length and number of mismatches needed to transform the signature into the closest subsequence present in the background. This definition, in conjunction with a probabilistic framework, allows one to predict the minimal signature length required to identify the target in the presence of different sizes of backgrounds and the effect of the target's mutation rate on the quality of its identification. The model assumptions and predictions were validated using both Monte Carlo simulations and real genomic data examples. The proposed model can be used to determine appropriate signature lengths for various combinations of target and background genome sizes. It also predicted that any genomic signatures will be unable to identify target if its mutation rate is >5%.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exhaustive computational identification of pathogen sequences far-distant from background genomes: Identification and experimental verification of human-blind dengue PCR primers

We recently developed novel algorithms for exhaustive identification of all nucleotide subsequences present in a pathogen genome which differ by at least a chosen number of mismatches from the sequences of host/background organisms. This type of exhaustive computational analysis will be useful in reducing false positives and cross-reactivity in PCR and hybridization assays. We present the first experimental test of the method by showing that the subsequences identified when used as 18-mer PCR primers can detect the presence of dengue virus (DENV) even in the presence of a large excess of complex human genomic DNA. From our computations, 715 serotype-specific primer pairs were identified for three different DENV serotypes in which each primer sequence lies at least two mismatches from the nearest human sequence. DNA clones of representative strains of DENV-1, DENV-2, and DENV-4 viruses were subjected to real-time PCR testing using eight primer pairs each. Efficiencies were uniformly very high (mean+/-S.D.=99.6+/-3%), and amplification of human DNA was never observed within 35 cycles, even at a 5.5-fold molar excess of human DNA. Exhaustive primer/probe screening can potentially produce more selective and sensitive diagnostic assays for pathogens, especially in the presence of complex backgrounds.

PIDA:A new algorithm for pattern identification

Algorithms for motif identification in sequence space have predominately been focused on recognizing patterns of a fixed length containing regions of perfect conservation with possible regions of unconstrained sequence. Such motifs can be found in everything from proteins with distinct active sites to non-coding RNAs with specific structural elements that are necessary to maintain functionality. In the event that an insertion/deletion has occurred within an unconstrained portion of the pattern, it is possible that the pattern retains its functionality. In such a case the length of the pattern is now variable and may be overlooked when utilizing existing motif detection methods. The Pattern Island Detection Algorithm (PIDA) presented here has been developed to recognize patterns that have occurrences of varying length within sequences of any size alphabet. PIDA works by identifying all regions of perfect conservation (for lengths longer than a user-specified threshold), and then builds those conservation "islands" into fixed-length patterns. Next the algorithm modifies these fixed-length patterns by identifying additional (and different) islands that can be incorporated into each pattern through insertions/deletions within the "water" separating the islands. To provide some benchmarks for this analysis, PIDA was used to search for patterns within randomly generated sequences as well as sequences known to contain conserved patterns. For each of the patterns found, the statistical significance is calculated based upon the pattern's likelihood to appear by chance, thus providing a means to determine those patterns which are likely to have a functional role. The PIDA approach to motif finding is designed to perform best when searching for patterns of variable length although it is also able to identify patterns of a fixed length. PIDA has been created to be as generally applicable as possible since there are a variety of sequence problems of this type. The algorithm was implemented in C++ and is freely available upon request from the authors.

A computational tool for the genomic identification of regions of unusual compositional properties and its utilization in the detection of horizontally transferred sequences

Similarity Plot (S-plot) is a Windows-based application for large-scale comparisons and 2-dimensional visualization of compositional similarities between genomic sequences. This application combines 2 approaches widely used in genomics: window analysis of statistical characteristics along genomes and dot-plot visual representation. S-plot is effective in identifying highly similar regions between genomes as well as regions with unusual compositional properties (RUCPs) within a single genome, which may be indicative of horizontal gene transfer or of locus-specific selective forces. We use S-plot to identify regions that may have originated through horizontal gene transfer through a 2-step approach, by first comparing a genomic sequence to itself and, subsequently, comparing it to the genomic sequence of a closely related taxon. Moreover, by comparing these suspect sequences to one another, we can estimate a minimum number of sources for these putative xenologous sequences. We illustrate the uses of S-plot in a comparison involving Escherichia coli K12 and E. coli O157:H7. In O157:H7, we found 145 regions that have most probably originated through horizontal gene transfer. By using S-plot to compare each of these regions with 277 completely sequenced prokaryotic genomes, 1 sequence was found to have similar compositional properties to the Yersinia pseudotuberculosis genome, indicating a transfer from a Yersinia or Yersinia relative. Based upon our analysis of RUCPs in O157:H7, we infer that there were at least 53 sources of horizontally transferred sequences.

Using Mutual Information to Discover Temporal Patterns in Gene Expression Data

Finding relations among gene expressions involves the definition of the similarity between experimental data. A simplest similarity measure is the Correlation Coefficient. It is able to identify linear dependences only; moreover, is sensitive to experimental errors. An alternative measure, the Shannon Mutual Information (MI), is free from the above mentioned weaknesses. However, the calculation of MI for continuous variables from the finite number of experimental points, N, involves an ambiguity arising when one divides the range of values of the continuous variable into boxes. Then the distribution of experimental points among the boxes (and, therefore, MI) depends on the box size. An algorithm for the calculation of MI for continuous variables is proposed. We find the optimum box sizes for a given N from the condition of minimum entropy variation with respect to the change of the box sizes. We have applied this technique to the gene expression dataset from Stanford, containing microarray data at 18 time points from yeast Saccharomyces cerevisiae cultures (Spellman et al.,[3]). We calculated MI for all of the pairs of time points. The MI analysis allowed us to identify time patterns related to different biological processes in the cell.