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.

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.

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.

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, E3D56-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 E3D56-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.

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.

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.

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.

Interleukin-33 Promotes REG3γ Expression in Intestinal Epithelial Cells and Regulates Gut Microbiota.

BACKGROUND & AIMS: Regenerating islet-derived protein (REG3γ), an antimicrobial peptide, typically expressed by intestinal epithelial cells (IEC), plays crucial roles in intestinal homeostasis and controlling gut microbiota. However, the mechanisms that regulate IEC expression of REG3γ are still largely unclear. In this study, we investigated whether and how interleukin (IL) 33, an alarmin produced by IEC in response to injury, regulates REG3γ expression in IEC, thus contributing to intestinal homeostasis. METHODS: IEC were isolated from wild-type and IL33-/- mice to determine expression of REG3γ and other antimicrobial peptides by quantitative real-time polymerase chain reaction and Western blot. IEC cell lines were used for mechanistic studies. 16S rRNA pyrosequencing analysis was used for measuring gut microbiota. Citrobacter rodentium was used for enteric infections. RESULTS: The expression of REG3γ, but not ß-defensins, in IECs of IL33-/- mice was significantly lower than wild-type mice. IL33 treatment induced IEC expression of REG3γ in both mice and human cell lines. Mechanistically, IL33 activated STAT3, mTOR, and ERK1/2 in IEC. Inhibition of these pathways abrogated IL33-induction of REG3γ. IL33-/- mice demonstrated higher bacteria loads and altered microbiota composition. IL33 did not directly inhibit bacterial growth, but promoted wild-type, not REG3γKO, IECs to kill bacteria in vitro. Consistently, C rodentium infection induced IEC IL33 expression, and IL33-/- mice demonstrated an impaired bacterial clearance with C rodentium infection. CONCLUSIONS: Our study demonstrated that IL33, which is produced by IEC in response to injury and inflammatory stimulation, in turn promotes IEC expression of REG3γ, and controls the gut microbiota of the host.

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.

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.

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.

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.

Microbiota dysbiosis and its pathophysiological significance in bowel obstruction.

Bowel obstruction (OB) causes local and systemic dysfunctions. Here we investigated whether obstruction leads to alterations in microbiota community composition and total abundance, and if so whether these changes contribute to dysfunctions in OB. Partial colon obstruction was maintained in rats for 7 days. The mid colon and its intraluminal feces - proximal to the obstruction - were studied. OB did not cause bacterial overgrowth or mucosa inflammation, but induced profound changes in fecal microbiota composition and diversity. At the phylum level, the 16S rRNA sequencing showed a significant decrease in the relative abundance of Firmicutes with corresponding increases in Proteobacteria and Bacteroidetes in OB compared with sham controls. Daily treatment using broad spectrum antibiotics dramatically reduced total bacterial abundance, but increased the relative presence of Proteobacteria. Antibiotics eliminated viable bacteria in the spleen and liver, but not in the mesentery lymph node in OB. Although antibiotic treatment decreased muscle contractility in sham rats, it had little effect on OB-associated suppression of muscle contractility or inflammatory changes in the muscle layer. In conclusion, obstruction leads to marked dysbiosis in the colon. Antibiotic eradication of microbiota had limited effects on obstruction-associated changes in inflammation, motility, or bacterial translocation.

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.

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.

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.

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.