MeFiT: merging and filtering tool for illumina paired-end reads for 16S rRNA amplicon sequencing.

BACKGROUND: Recent advances in next-generation sequencing have revolutionized genomic research. 16S rRNA amplicon sequencing using paired-end sequencing on the MiSeq platform from Illumina, Inc., is being used to characterize the composition and dynamics of extremely complex/diverse microbial communities. For this analysis on the Illumina platform, merging and quality filtering of paired-end reads are essential first steps in data analysis to ensure the accuracy and reliability of downstream analysis. RESULTS: We have developed the Merging and Filtering Tool (MeFiT) to combine these pre-processing steps into one simple, intuitive pipeline. MeFiT invokes CASPER (context-aware scheme for paired-end reads) for merging paired-end reads and provides users the option to quality filter the reads using the traditional average Q-score metric or using a maximum expected error cut-off threshold. CONCLUSIONS: MeFiT provides an open-source solution that permits users to merge and filter paired end illumina reads. The tool has been implemented in python and the source-code is freely available at https://github.com/nisheth/MeFiT .

Dynamical System Modeling to Simulate Donor T Cell Response to Whole Exome Sequencing-Derived Recipient Peptides Demonstrates Different Alloreactivity Potential in HLA-Matched and -Mismatched Donor-Recipient Pairs.

Immune reconstitution kinetics and subsequent clinical outcomes in HLA-matched recipients of allogeneic stem cell transplantation (SCT) are variable and difficult to predict. Considering SCT as a dynamical system may allow sequence differences across the exomes of the transplant donors and recipients to be used to simulate an alloreactive T cell response, which may allow better clinical outcome prediction. To accomplish this, whole exome sequencing was performed on 34 HLA-matched SCT donor-recipient pairs (DRPs) and the nucleotide sequence differences translated to peptides. The binding affinity of the peptides to the relevant HLA in each DRP was determined. The resulting array of peptide-HLA binding affinity values in each patient was considered as an operator modifying a hypothetical T cell repertoire vector, in which each T cell clone proliferates in accordance with the logistic equation of growth. Using an iterating system of matrices, each simulated T cell clone's growth was calculated with the steady-state population being proportional to the magnitude of the binding affinity of the driving HLA-peptide complex. Incorporating competition between T cell clones responding to different HLA-peptide complexes reproduces a number of features of clinically observed T cell clonal repertoire in the simulated repertoire, including sigmoidal growth kinetics of individual T cell clones and overall repertoire, Power Law clonal frequency distribution, increase in repertoire complexity over time with increasing clonal diversity, and alteration of clonal dominance when a different antigen array is encountered, such as in SCT. The simulated, alloreactive T cell repertoire was markedly different in HLA-matched DRPs. The patterns were differentiated by rate of growth and steady-state magnitude of the simulated T cell repertoire and demonstrate a possible correlation with survival. In conclusion, exome wide sequence differences in DRPs may allow simulation of donor alloreactive T cell response to recipient antigens and may provide a quantitative basis for refining donor selection and titration of immunosuppression after SCT.

Comparison of Lactobacillus crispatus isolates from Lactobacillus-dominated vaginal microbiomes with isolates from microbiomes containing bacterial vaginosis-associated bacteria.

Vaginal lactobacilli can inhibit colonization by and growth of other bacteria, thereby preventing development of bacterial vaginosis (BV). Amongst the lactobacilli, Lactobacillus crispatus appears to be particularly effective at inhibiting growth of BV-associated bacteria. Nonetheless, some women who are colonized with this species can still develop clinical BV. Therefore, we sought to determine whether strains of L. crispatus that colonize women with lactobacilli-dominated vaginal microbiomes are distinct from strains that colonize women who develop BV. The genomes of L. crispatus isolates from four women with lactobacilli-dominated vaginal microbiomes ( <1% 16S rRNA reads above threshold from genera other than Lactobacillus) and four women with microbiomes containing BV-associated bacteria (>12% 16S rRNA reads from bacterial taxa associated with BV) were sequenced and compared. Lactic acid production by the different strains was quantified. Phage induction in the strains was also analysed. There was considerable genetic diversity between strains, and several genes were exclusive to either the strains from Lactobacillus-dominated microbiomes or those containing BV-associated bacteria. Overall, strains from microbiomes dominated by lactobacilli did not differ from strains from microbiomes containing BV-associated bacteria with respect to lactic acid production. All of the strains contained multiple phage, but there was no clear distinction between the presence or absence of BV-associated bacteria with respect to phage-induced lysis. Genes found to be exclusive to the Lactobacillus-dominated versus BV-associated bacteria-containing microbiomes could play a role in the maintenance of vaginal health and the development of BV, respectively.

The truth about metagenomics: quantifying and counteracting bias in 16S rRNA studies.

BACKGROUND: Characterizing microbial communities via next-generation sequencing is subject to a number of pitfalls involving sample processing. The observed community composition can be a severe distortion of the quantities of bacteria actually present in the microbiome, hampering analysis and threatening the validity of conclusions from metagenomic studies. We introduce an experimental protocol using mock communities for quantifying and characterizing bias introduced in the sample processing pipeline. We used 80 bacterial mock communities comprised of prescribed proportions of cells from seven vaginally-relevant bacterial strains to assess the bias introduced in the sample processing pipeline. We created two additional sets of 80 mock communities by mixing prescribed quantities of DNA and PCR product to quantify the relative contribution to bias of (1) DNA extraction, (2) PCR amplification, and (3) sequencing and taxonomic classification for particular choices of protocols for each step. We developed models to predict the "true" composition of environmental samples based on the observed proportions, and applied them to a set of clinical vaginal samples from a single subject during four visits. RESULTS: We observed that using different DNA extraction kits can produce dramatically different results but bias is introduced regardless of the choice of kit. We observed error rates from bias of over 85% in some samples, while technical variation was very low at less than 5% for most bacteria. The effects of DNA extraction and PCR amplification for our protocols were much larger than those due to sequencing and classification. The processing steps affected different bacteria in different ways, resulting in amplified and suppressed observed proportions of a community. When predictive models were applied to clinical samples from a subject, the predicted microbiome profiles were better reflections of the physiology and diagnosis of the subject at the visits than the observed community compositions. CONCLUSIONS: Bias in 16S studies due to DNA extraction and PCR amplification will continue to require attention despite further advances in sequencing technology. Analysis of mock communities can help assess bias and facilitate the interpretation of results from environmental samples.

Skin-to-Skin Care and the Development of the Preterm Infant Oral Microbiome.

OBJECTIVE: The oral cavity represents an initial entry way for oral and gut indigenous colonization. Skin-to-skin (STS) care, in which the mother holds the diaper clad naked preterm (PT) infant between her breasts, is associated with improved digestive function, decreased stress, and improved survival. This study evaluated the development of oral microbial colonization repertoires and health characteristics in PT infants with or without STS exposure. METHODS: Saliva from 42 PT infants (<32 weeks of gestation at birth) was collected prospectively at 1 month and/or at discharge. High-throughput 16S rRNA sequencing identified microbial diversity and prevalence of bacterial signatures correlated with clinical STS or non-STS care. RESULTS: Corrected for gestational age (CGA) at sampling, bacterial taxa demonstrated increased Streptococcus as a signature of oral repertoire maturation. STS was associated with increased Streptococcus (p < 0.024), while non-STS was associated with greater Corynebacterium (p < 0.023) and Pseudomonas (p < 0.019) in infants ≤ 32 weeks CGA. In infants > 32 weeks CGA, Neisseria and Acinetobacter were more prevalent, 50 vs. 16.7% and 40 vs. 0%, respectively. STS care was associated with shorter hospitalization (p < 0.039). CONCLUSION: STS care during earlier gestation was associated with a distinct microbial pattern and an accelerated pace of oral microbial repertoire maturity.

Whole exome sequencing to estimate alloreactivity potential between donors and recipients in stem cell transplantation.

Whole exome sequencing (WES) was performed on stem cell transplant donor-recipient (D-R) pairs to determine the extent of potential antigenic variation at a molecular level. In a small cohort of D-R pairs, a high frequency of sequence variation was observed between the donor and recipient exomes independent of human leucocyte antigen (HLA) matching. Nonsynonymous, nonconservative single nucleotide polymorphisms were approximately twice as frequent in HLA-matched unrelated, compared with related D-R pairs. When mapped to individual chromosomes, these polymorphic nucleotides were uniformly distributed across the entire exome. In conclusion, WES reveals extensive nucleotide sequence variation in the exomes of HLA-matched donors and recipients.

Differences in vaginal microbiome in African American women versus women of European ancestry.

Women of European ancestry are more likely to harbour a Lactobacillus-dominated microbiome, whereas African American women are more likely to exhibit a diverse microbial profile. African American women are also twice as likely to be diagnosed with bacterial vaginosis and are twice as likely to experience preterm birth. The objective of this study was to further characterize and contrast the vaginal microbial profiles in African American versus European ancestry women. Through the Vaginal Human Microbiome Project at Virginia Commonwealth University, 16S rRNA gene sequence analysis was used to compare the microbiomes of vaginal samples from 1268 African American women and 416 women of European ancestry. The results confirmed significant differences in the vaginal microbiomes of the two groups and identified several taxa relevant to these differences. Major community types were dominated by Gardnerella vaginalis and the uncultivated bacterial vaginosis-associated bacterium-1 (BVAB1) that were common among African Americans. Moreover, the prevalence of multiple bacterial taxa that are associated with microbial invasion of the amniotic cavity and preterm birth, including Mycoplasma, Gardnerella, Prevotella and Sneathia, differed between the two ethnic groups. We investigated the contributions of intrinsic and extrinsic factors, including pregnancy, body mass index, diet, smoking and alcohol use, number of sexual partners, and household income, to vaginal community composition. Ethnicity, pregnancy and alcohol use correlated significantly with the relative abundance of bacterial vaginosis-associated species. Trends between microbial profiles and smoking and number of sexual partners were observed; however, these associations were not statistically significant. These results support and extend previous findings that there are significant differences in the vaginal microbiome related to ethnicity and demonstrate that these differences are pronounced even in healthy women.

Second-generation shRNA libraries covering the mouse and human genomes.

Loss-of-function phenotypes often hold the key to understanding the connections and biological functions of biochemical pathways. We and others previously constructed libraries of short hairpin RNAs that allow systematic analysis of RNA interference-induced phenotypes in mammalian cells. Here we report the construction and validation of second-generation short hairpin RNA expression libraries designed using an increased knowledge of RNA interference biochemistry. These constructs include silencing triggers designed to mimic a natural microRNA primary transcript, and each target sequence was selected on the basis of thermodynamic criteria for optimal small RNA performance. Biochemical and phenotypic assays indicate that the new libraries are substantially improved over first-generation reagents. We generated large-scale-arrayed, sequence-verified libraries comprising more than 140,000 second-generation short hairpin RNA expression plasmids, covering a substantial fraction of all predicted genes in the human and mouse genomes. These libraries are available to the scientific community.

Preliminary Evaluation of a Novel Neural Network-Based Hybrid Simulator for Surgical Training and Performance Assessment of Neonatal Thoracentesis.

INTRODUCTION: Tension pneumothorax is a rare and life-threatening situation in neonates requiring immediate intervention through thoracentesis. Significant complications can arise while performing thoracentesis in the case of inadequate skill level or exposure to the condition. Although simulation-based training (SBT) has proven to be effective in learning surgical skills, training sessions are long, subjective, and expensive, because of which they cannot be held regularly. This article attempts to improve traditional SBT for neonatal thoracentesis through an autonomous simulator that can provide real-time objective feedback during surgical training and assessment. METHODS: The simulator incorporates a custom manikin and virtual reality software interfaced through electromagnetic sensors that track the motion of surgical instruments. The software application reads and stores instrument motion information to replicate physical actions in the virtual environment, play back previously stored surgical performances and analyze data through a pretrained neural network. The simulator encapsulates the experience of SBT by allowing trainees to watch and replicate an ideal method of conducting the procedure, providing simplified, real-time autonomous guidance during practice and an objective taskwise assessment of the performance during testing. RESULTS: The preliminary trial held at the University of Illinois Hospital in the presence of 1 neonatologist and 4 fellows revealed that all the participants used the autonomous guidance more than once, and all found simulation experience to be accurate and overall effective in learning thoracentesis. CONCLUSION: Although the sample size is small, the simulator shows potential in being a viable alternative approach for training and assessment for thoracentesis.

Species-level classification of the vaginal microbiome.

BACKGROUND: The application of next-generation sequencing to the study of the vaginal microbiome is revealing the spectrum of microbial communities that inhabit the human vagina. High-resolution identification of bacterial taxa, minimally to the species level, is necessary to fully understand the association of the vaginal microbiome with bacterial vaginosis, sexually transmitted infections, pregnancy complications, menopause, and other physiological and infectious conditions. However, most current taxonomic assignment strategies based on metagenomic 16S rDNA sequence analysis provide at best a genus-level resolution. While surveys of 16S rRNA gene sequences are common in microbiome studies, few well-curated, body-site-specific reference databases of 16S rRNA gene sequences are available, and no such resource is available for vaginal microbiome studies. RESULTS: We constructed the Vaginal 16S rDNA Reference Database, a comprehensive and non-redundant database of 16S rDNA reference sequences for bacterial taxa likely to be associated with vaginal health, and we developed STIRRUPS, a new method that employs the USEARCH algorithm with a curated reference database for rapid species-level classification of 16S rDNA partial sequences. The method was applied to two datasets of V1-V3 16S rDNA reads: one generated from a mock community containing DNA from six bacterial strains associated with vaginal health, and a second generated from over 1,000 mid-vaginal samples collected as part of the Vaginal Human Microbiome Project at Virginia Commonwealth University. In both datasets, STIRRUPS, used in conjunction with the Vaginal 16S rDNA Reference Database, classified more than 95% of processed reads to a species-level taxon using a 97% global identity threshold for assignment. CONCLUSIONS: This database and method provide accurate species-level classifications of metagenomic 16S rDNA sequence reads that will be useful for analysis and comparison of microbiome profiles from vaginal samples. STIRRUPS can be used to classify 16S rDNA sequence reads from other ecological niches if an appropriate reference database of 16S rDNA sequences is available.

Racioethnic diversity in the dynamics of the vaginal microbiome during pregnancy.

The microbiome of the female reproductive tract has implications for women's reproductive health. We examined the vaginal microbiome in two cohorts of women who experienced normal term births: a cross-sectionally sampled cohort of 613 pregnant and 1,969 non-pregnant women, focusing on 300 pregnant and 300 non-pregnant women of African, Hispanic or European ancestry case-matched for race, gestational age and household income; and a longitudinally sampled cohort of 90 pregnant women of African or non-African ancestry. In these women, the vaginal microbiome shifted during pregnancy toward Lactobacillus-dominated profiles at the expense of taxa often associated with vaginal dysbiosis. The shifts occurred early in pregnancy, followed predictable patterns, were associated with simplification of the metabolic capacity of the microbiome and were significant only in women of African or Hispanic ancestry. Both genomic and environmental factors are likely contributors to these trends, with socioeconomic status as a likely environmental influence.

The vaginal microbiome and preterm birth.

The incidence of preterm birth exceeds 10% worldwide. There are significant disparities in the frequency of preterm birth among populations within countries, and women of African ancestry disproportionately bear the burden of risk in the United States. In the present study, we report a community resource that includes 'omics' data from approximately 12,000 samples as part of the integrative Human Microbiome Project. Longitudinal analyses of 16S ribosomal RNA, metagenomic, metatranscriptomic and cytokine profiles from 45 preterm and 90 term birth controls identified harbingers of preterm birth in this cohort of women predominantly of African ancestry. Women who delivered preterm exhibited significantly lower vaginal levels of Lactobacillus crispatus and higher levels of BVAB1, Sneathia amnii, TM7-H1, a group of Prevotella species and nine additional taxa. The first representative genomes of BVAB1 and TM7-H1 are described. Preterm-birth-associated taxa were correlated with proinflammatory cytokines in vaginal fluid. These findings highlight new opportunities for assessment of the risk of preterm birth.

Comprehensive splice-site analysis using comparative genomics.

We have collected over half a million splice sites from five species-Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans and Arabidopsis thaliana-and classified them into four subtypes: U2-type GT-AG and GC-AG and U12-type GT-AG and AT-AC. We have also found new examples of rare splice-site categories, such as U12-type introns without canonical borders, and U2-dependent AT-AC introns. The splice-site sequences and several tools to explore them are available on a public website (SpliceRack). For the U12-type introns, we find several features conserved across species, as well as a clustering of these introns on genes. Using the information content of the splice-site motifs, and the phylogenetic distance between them, we identify: (i) a higher degree of conservation in the exonic portion of the U2-type splice sites in more complex organisms; (ii) conservation of exonic nucleotides for U12-type splice sites; (iii) divergent evolution of C.elegans 3' splice sites (3'ss) and (iv) distinct evolutionary histories of 5' and 3'ss. Our study proves that the identification of broad patterns in naturally-occurring splice sites, through the analysis of genomic datasets, provides mechanistic and evolutionary insights into pre-mRNA splicing.

Ethanol-Induced Behavioral Sensitization Alters the Synaptic Transcriptome and Exon Utilization in DBA/2J Mice.

Alcoholism is a complex behavioral disorder characterized by loss of control in limiting intake, and progressive compulsion to seek and consume ethanol. Prior studies have suggested that the characteristic behaviors associated with escalation of drug use are caused, at least in part, by ethanol-evoked changes in gene expression affecting synaptic plasticity. Implicit in this hypothesis is a dependence on new protein synthesis and remodeling at the synapse. It is well established that mRNA can be transported to distal dendritic processes, where it can undergo localized translation. It is unknown whether such modulation of the synaptic transcriptome might contribute to ethanol-induced synaptic plasticity. Using ethanol-induced behavioral sensitization as a model of neuroplasticity, we investigated whether repeated exposure to ethanol altered the synaptic transcriptome, contributing to mechanisms underlying subsequent increases in ethanol-evoked locomotor activity. RNAseq profiling of DBA/2J mice subjected to acute ethanol or ethanol-induced behavioral sensitization was performed on frontal pole synaptoneurosomes to enrich for synaptic mRNA. Genomic profiling showed distinct functional classes of mRNA enriched in the synaptic vs. cytosolic fractions, consistent with their role in synaptic function. Ethanol sensitization regulated more than twice the number of synaptic localized genes compared to acute ethanol exposure. Synaptic biological processes selectively perturbed by ethanol sensitization included protein folding and modification as well as and mitochondrial respiratory function, suggesting repeated ethanol exposure alters synaptic energy production and the processing of newly translated proteins. Additionally, marked differential exon usage followed ethanol sensitization in both synaptic and non-synaptic cellular fractions, with little to no perturbation following acute ethanol exposure. Altered synaptic exon usage following ethanol sensitization strongly affected genes related to RNA processing and stability, translational regulation, and synaptic function. These genes were also enriched for targets of the FMRP RNA-binding protein and contained consensus sequence motifs related to other known RNA binding proteins, suggesting that ethanol sensitization altered selective mRNA trafficking mechanisms. This study provides a foundation for investigating the role of ethanol in modifying the synaptic transcriptome and inducing changes in synaptic plasticity.

Effects of combined oral contraceptives, depot medroxyprogesterone acetate and the levonorgestrel-releasing intrauterine system on the vaginal microbiome.

OBJECTIVES: Prior studies suggest that the composition of the vaginal microbiome may positively or negatively affect susceptibility to sexually transmitted infections (STIs) and bacterial vaginosis (BV). Some female hormonal contraceptive methods also appear to positively or negatively influence STI transmission and BV. Therefore, changes in the vaginal microbiome that are associated with different contraceptive methods may explain, in part, effects on STI transmission and BV. STUDY DESIGN: We performed a retrospective study of 16S rRNA gene survey data of vaginal samples from a subset of participants from the Human Vaginal Microbiome Project at Virginia Commonwealth University. The subset included 682 women who reported using a single form of birth control that was condoms, combined oral contraceptives (COCs), depot medroxyprogesterone acetate (DMPA) or the levonorgestrel-releasing intrauterine system (LNG-IUS). RESULTS: Women using COCs [adjusted odds ratio (aOR) 0.29, 95% confidence interval (CI) 0.13-0.64] and DMPA (aOR 0.34, 95% CI 0.13-0.89), but not LNG-IUS (aOR 1.55, 95% CI 0.72-3.35), were less likely to be colonized by BV-associated bacteria relative to women who used condoms. Women using COCs (aOR 1.94, 95% CI 1.25-3.02) were more likely to be colonized by beneficial H2O2-producing Lactobacillus species compared with women using condoms, while women using DMPA (aOR 1.09, 95% CI 0.63-1.86) and LNG-IUS (aOR 0.74, 95% CI 0.48-1.15) were not. CONCLUSIONS: Use of COCs is significantly associated with increased vaginal colonization by healthy lactobacilli and reduced BV-associated taxa. IMPLICATIONS: COC use may positively influence gynecologic health through an increase in healthy lactobacilli and a decrease in BV-associated bacterial taxa.

Rare mutations and potentially damaging missense variants in genes encoding fibrillar collagens and proteins involved in their production are candidates for risk for preterm premature rupture of membranes.

Preterm premature rupture of membranes (PPROM) is the leading identifiable cause of preterm birth with ~ 40% of preterm births being associated with PPROM and occurs in 1% - 2% of all pregnancies. We hypothesized that multiple rare variants in fetal genes involved in extracellular matrix synthesis would associate with PPROM, based on the assumption that impaired elaboration of matrix proteins would reduce fetal membrane tensile strength, predisposing to unscheduled rupture. We performed whole exome sequencing (WES) on neonatal DNA derived from pregnancies complicated by PPROM (49 cases) and healthy term deliveries (20 controls) to identify candidate mutations/variants. Genotyping for selected variants from the WES study was carried out on an additional 188 PPROM cases and 175 controls. All mothers were self-reported African Americans, and a panel of ancestry informative markers was used to control for genetic ancestry in all genetic association tests. In support of the primary hypothesis, a statistically significant genetic burden (all samples combined, SKAT-O p-value = 0.0225) of damaging/potentially damaging rare variants was identified in the genes of interest-fibrillar collagen genes, which contribute to fetal membrane strength and integrity. These findings suggest that the fetal contribution to PPROM is polygenic, and driven by an increased burden of rare variants that may also contribute to the disparities in rates of preterm birth among African Americans.

Association between statin use, the vaginal microbiome, and Gardnerella vaginalis vaginolysin-mediated cytotoxicity.

BACKGROUND: Bacterial vaginosis (BV) is the leading dysbiosis of the vaginal microbiome. The pathways leading towards the development of BV are not well understood. Gardnerella vaginalis is frequently associated with BV. G. vaginalis produces the cholesterol-dependent cytolysin (CDC), vaginolysin, which can lyse a variety of human cells and is thought to play a role in pathogenesis. Because membrane cholesterol is required for vaginolysin to function, and because HMG-CoA reductase inhibitors (statins) affect not only serum levels of cholesterol but membrane levels as well, we hypothesized that statins might affect the vaginal microbiome. METHODS: To investigate the relationship between use of the statins and the vaginal microbiome, we analyzed 16S rRNA gene taxonomic surveys performed on vaginal samples from 133 women who participated in the Vaginal Human Microbiome Project and who were taking statins at the time of sampling, 152 women who reported high cholesterol levels but were not taking statins, and 316 women who did not report high cholesterol. To examine the effect of statins on the cytolytic effect of vaginolysin, the cholesterol-dependent cytolysin (CDC) produced by Gardnerella vaginalis, we assessed the effect of simvastatin pretreatment of VK2E6/E7 vaginal epithelial cells on vaginolysin-mediated cytotoxicity. RESULTS: The mean proportion of G. vaginalis among women taking statins was significantly lower relative to women not using statins. Women using statins had higher mean proportions of Lactobacillus crispatus relative to women with normal cholesterol levels, and higher levels of Lactobacillus jensenii relative to women with high cholesterol but not taking statins. In vitro, vaginal epithelial cells pretreated with simvastatin were relatively resistant to vaginolysin and this effect was inhibited by cholesterol. CONCLUSIONS: In this cross-sectional study, statin use was associated with reduced proportions of G. vaginalis and greater proportions of beneficial lactobacilli within the vaginal microbiome. The negative association between statin use and G. vaginalis may be related to inhibition of vaginolysin function.

Dynamical system modeling to simulate donor T cell response to whole exome sequencing-derived recipient peptides: Understanding randomness in alloreactivity incidence following stem cell transplantation.

Quantitative relationship between the magnitude of variation in minor histocompatibility antigens (mHA) and graft versus host disease (GVHD) pathophysiology in stem cell transplant (SCT) donor-recipient pairs (DRP) is not established. In order to elucidate this relationship, whole exome sequencing (WES) was performed on 27 HLA matched related (MRD), & 50 unrelated donors (URD), to identify nonsynonymous single nucleotide polymorphisms (SNPs). An average 2,463 SNPs were identified in MRD, and 4,287 in URD DRP (p<0.01); resulting peptide antigens that may be presented on HLA class I molecules in each DRP were derived in silico (NetMHCpan ver2.0) and the tissue expression of proteins these were derived from determined (GTex). MRD DRP had an average 3,670 HLA-binding-alloreactive peptides, putative mHA (pmHA) with an IC50 of <500 nM, and URD, had 5,386 (p<0.01). To simulate an alloreactive donor cytotoxic T cell response, the array of pmHA in each patient was considered as an operator matrix modifying a hypothetical cytotoxic T cell clonal vector matrix; each responding T cell clone's proliferation was determined by the logistic equation of growth, accounting for HLA binding affinity and tissue expression of each alloreactive peptide. The resulting simulated organ-specific alloreactive T cell clonal growth revealed marked variability, with the T cell count differences spanning orders of magnitude between different DRP. Despite an estimated, uniform set of constants used in the model for all DRP, and a heterogeneously treated group of patients, higher total and organ-specific T cell counts were associated with cumulative incidence of moderate to severe GVHD in recipients. In conclusion, exome wide sequence differences and the variable alloreactive peptide binding to HLA in each DRP yields a large range of possible alloreactive donor T cell responses. Our findings also help understand the apparent randomness observed in the development of alloimmune responses.

BOTUX: bayesian-like operational taxonomic unit examiner.

Bayesian-like operational taxonomic unit examiner (BOTUX) is a new tool for the classification of 16S rRNA gene sequences into operational taxonomic units (OTUs) that addresses the problem of overestimation caused by errors introduced during PCR amplification and DNA sequencing steps. BOTUX utilises a grammar-based assignment strategy, where Bayesian models are built from each word of a given length (e.g., 8-mers). de novo analysis is possible with BOTUX as it does not require a training set, and updates probabilistic models as new sequences are recruited to an OTU. In benchmarking tests performed with real and simulated datasets of 16S rDNA sequences, BOTUX accurately identifies OTUs with comparable or better clustering efficiency and lower execution times than other OTU algorithms tested. BOTUX is the only OTU classifier, which allows incremental analysis of large datasets, and is also adept in clustering both 454 and Illumina datasets in a reasonable timeframe.

An emerging mycoplasma associated with trichomoniasis, vaginal infection and disease.

Humans are colonized by thousands of bacterial species, but it is difficult to assess the metabolic and pathogenic potential of the majority of these because they have yet to be cultured. Here, we characterize an uncultivated vaginal mycoplasma tightly associated with trichomoniasis that was previously known by its 16S rRNA sequence as "Mnola." In this study, the mycoplasma was found almost exclusively in women infected with the sexually transmitted pathogen Trichomonas vaginalis, but rarely observed in women with no diagnosed disease. The genomes of four strains of this species were reconstructed using metagenome sequencing and assembly of DNA from four discrete mid-vaginal samples, one of which was obtained from a pregnant woman with trichomoniasis who delivered prematurely. These bacteria harbor several putative virulence factors and display unique metabolic strategies. Genes encoding proteins with high similarity to potential virulence factors include two collagenases, a hemolysin, an O-sialoglycoprotein endopeptidase and a feoB-type ferrous iron transport system. We propose the name "Candidatus Mycoplasma girerdii" for this potential new pathogen.