BiocPkgTools: Toolkit for mining the Bioconductor package ecosystem.

Motivation: The Bioconductor project, a large collection of open source software for the comprehension of large-scale biological data, continues to grow with new packages added each week, motivating the development of software tools focused on exposing package metadata to developers and users. The resulting BiocPkgTools package facilitates access to extensive metadata in computable form covering the Bioconductor package ecosystem, facilitating downstream applications such as custom reporting, data and text mining of Bioconductor package text descriptions, graph analytics over package dependencies, and custom search approaches. Results: The BiocPkgTools package has been incorporated into the Bioconductor project, installs using standard procedures, and runs on any system supporting R. It provides functions to load detailed package metadata, longitudinal package download statistics, package dependencies, and Bioconductor build reports, all in "tidy data" form. BiocPkgTools can convert from tidy data structures to graph structures, enabling graph-based analytics and visualization. An end-user-friendly graphical package explorer aids in task-centric package discovery. Full documentation and example use cases are included. Availability: The BiocPkgTools software and complete documentation are available from Bioconductor ( https://bioconductor.org/packages/BiocPkgTools).

Orchestrating a community-developed computational workshop and accompanying training materials.

The importance of bioinformatics, computational biology, and data science in biomedical research continues to grow, driving a need for effective instruction and education. A workshop setting, with lectures and guided hands-on tutorials, is a common approach to teaching practical computational and analytical methods. Here, we detail the process we used to produce high-quality, community-authored educational materials that are available for public consumption and reuse. The coordinated efforts of 17 authors over 10 weeks resulted in 15 workshops available as a website and as a 388-page electronic book. We describe how we utilized cloud infrastructure, GitHub, and a literate programming approach to robustly deliver hands-on tutorials to participants of the annual Bioconductor conference. The scripts, raw and published workshop materials, and cloud machine image are all openly available. Our approach uses free services and software and can be adapted by workshop organizers and authors in other contests with appropriate technical backgrounds.

Relationship of Specific Bacteria in the Cervical and Vaginal Microbiotas With Cervicitis.

BACKGROUND: Cervicitis is an inflammatory condition of the cervix associated with upper genital tract infection and reproductive complications. Although cervicitis can be caused by several known pathogens, the etiology frequently remains obscure. Here we investigate vaginal bacteria associated with bacterial vaginosis as potential causes of cervicitis. METHODS: Associations between vaginal bacteria and cervicitis were assessed in a retrospective case-control study of women attending a Seattle sexually transmitted disease clinic. Individual bacterial species were detected using 2 molecular methods: quantitative polymerase chain reaction (qPCR) and broad-range 16S rRNA gene PCR with pyrosequencing. The primary finding from this initial study was evaluated using qPCR in a second cohort of Kenyan women. RESULTS: The presence of Mageeibacillus indolicus, formerly BVAB3, in the cervix was associated with cervicitis, whereas the presence of Lactobacillus jensenii was inversely associated. Quantities of these bacteria did not differ between cervicitis cases and controls, although in a model inclusive of presence and abundance, M. indolicus remained significantly associated with cervicitis after adjustment for other cervicitis-causing pathogens. M. indolicus was not associated with cervicitis in our study of Kenyan women, possibly due to differences in the clinical definition of cervicitis. CONCLUSIONS: Colonization of the endocervix with M. indolicus may contribute to the clinical manifestations of cervicitis, but further study is needed to determine whether this finding is repeatable and applicable to diverse groups of women. Colonization of the cervix with L. jensenii could be a marker of health, perhaps reducing inflammation or inhibiting pathogenic infection.

Metabolic signatures of bacterial vaginosis.

UNLABELLED: Bacterial vaginosis (BV) is characterized by shifts in the vaginal microbiota from Lactobacillus dominant to a microbiota with diverse anaerobic bacteria. Few studies have linked specific metabolites with bacteria found in the human vagina. Here, we report dramatic differences in metabolite compositions and concentrations associated with BV using a global metabolomics approach. We further validated important metabolites using samples from a second cohort of women and a different platform to measure metabolites. In the primary study, we compared metabolite profiles in cervicovaginal lavage fluid from 40 women with BV and 20 women without BV. Vaginal bacterial representation was determined using broad-range PCR with pyrosequencing and concentrations of bacteria by quantitative PCR. We detected 279 named biochemicals; levels of 62% of metabolites were significantly different in women with BV. Unsupervised clustering of metabolites separated women with and without BV. Women with BV have metabolite profiles marked by lower concentrations of amino acids and dipeptides, concomitant with higher levels of amino acid catabolites and polyamines. Higher levels of the signaling eicosanoid 12-hydroxyeicosatetraenoic acid (12-HETE), a biomarker for inflammation, were noted in BV. Lactobacillus crispatus and Lactobacillus jensenii exhibited similar metabolite correlation patterns, which were distinct from correlation patterns exhibited by BV-associated bacteria. Several metabolites were significantly associated with clinical signs and symptoms (Amsel criteria) used to diagnose BV, and no metabolite was associated with all four clinical criteria. BV has strong metabolic signatures across multiple metabolic pathways, and these signatures are associated with the presence and concentrations of particular bacteria. IMPORTANCE: Bacterial vaginosis (BV) is a common but highly enigmatic condition that is associated with adverse outcomes for women and their neonates. Small molecule metabolites in the vagina may influence host physiology, affect microbial community composition, and impact risk of adverse health outcomes, but few studies have comprehensively studied the metabolomics profile of BV. Here, we used mass spectrometry to link specific metabolites with particular bacteria detected in the human vagina by PCR. BV was associated with strong metabolic signatures across multiple pathways affecting amino acid, carbohydrate, and lipid metabolism, highlighting the profound metabolic changes in BV. These signatures were associated with the presence and concentrations of particular vaginal bacteria, including some bacteria yet to be cultivated, thereby providing clues as to the microbial origin of many metabolites. Insights from this study provide opportunities for developing new diagnostic markers of BV and novel approaches for treatment or prevention of BV.

More than meets the eye: associations of vaginal bacteria with gram stain morphotypes using molecular phylogenetic analysis.

Bacterial vaginosis (BV) is a highly prevalent condition associated with adverse health outcomes. Gram stain analysis of vaginal fluid is the standard for confirming the diagnosis of BV, wherein abundances of key bacterial morphotypes are assessed. These Lactobacillus, Gardnerella, Bacteroides, and Mobiluncus morphotypes were originally linked to particular bacterial species through cultivation studies, but no studies have systematically investigated associations between uncultivated bacteria detected by molecular methods and Gram stain findings. In this study, 16S-rRNA PCR/pyrosequencing was used to examine associations between vaginal bacteria and bacterial morphotypes in 220 women with and without BV. Species-specific quantitative PCR (qPCR) and fluorescence in Situ hybridization (FISH) methods were used to document concentrations of two bacteria with curved rod morphologies: Mobiluncus and the fastidious BV-associated bacterium-1 (BVAB1). Rank abundance of vaginal bacteria in samples with evidence of curved gram-negative rods showed that BVAB1 was dominant (26.1%), while Mobiluncus was rare (0.2% of sequence reads). BVAB1 sequence reads were associated with Mobiluncus morphotypes (p<0.001). Among women with curved rods, mean concentration of BVAB1 DNA was 2 log units greater than Mobiluncus (p<0.001) using species-specific quantitative PCR. FISH analyses revealed that mean number of BVAB1 cells was 2 log units greater than Mobiluncus cells in women with highest Nugent score (p<0.001). Prevotella and Porphyromonas spp. were significantly associated with the "Bacteroides morphotype," whereas Bacteroides species were rare. Gram-negative rods designated Mobiluncus morphotypes on Gram stain are more likely BVAB1. These findings provide a clearer picture of the bacteria associated with morphotypes on vaginal Gram stain.

Software for computing and annotating genomic ranges.

We describe Bioconductor infrastructure for representing and computing on annotated genomic ranges and integrating genomic data with the statistical computing features of R and its extensions. At the core of the infrastructure are three packages: IRanges, GenomicRanges, and GenomicFeatures. These packages provide scalable data structures for representing annotated ranges on the genome, with special support for transcript structures, read alignments and coverage vectors. Computational facilities include efficient algorithms for overlap and nearest neighbor detection, coverage calculation and other range operations. This infrastructure directly supports more than 80 other Bioconductor packages, including those for sequence analysis, differential expression analysis and visualization.

Bacterial communities in women with bacterial vaginosis: high resolution phylogenetic analyses reveal relationships of microbiota to clinical criteria.

BACKGROUND: Bacterial vaginosis (BV) is a common condition that is associated with numerous adverse health outcomes and is characterized by poorly understood changes in the vaginal microbiota. We sought to describe the composition and diversity of the vaginal bacterial biota in women with BV using deep sequencing of the 16S rRNA gene coupled with species-level taxonomic identification. We investigated the associations between the presence of individual bacterial species and clinical diagnostic characteristics of BV. METHODOLOGY/PRINCIPAL FINDINGS: Broad-range 16S rRNA gene PCR and pyrosequencing were performed on vaginal swabs from 220 women with and without BV. BV was assessed by Amsel's clinical criteria and confirmed by Gram stain. Taxonomic classification was performed using phylogenetic placement tools that assigned 99% of query sequence reads to the species level. Women with BV had heterogeneous vaginal bacterial communities that were usually not dominated by a single taxon. In the absence of BV, vaginal bacterial communities were dominated by either Lactobacillus crispatus or Lactobacillus iners. Leptotrichia amnionii and Eggerthella sp. were the only two BV-associated bacteria (BVABs) significantly associated with each of the four Amsel's criteria. Co-occurrence analysis revealed the presence of several sub-groups of BVABs suggesting metabolic co-dependencies. Greater abundance of several BVABs was observed in Black women without BV. CONCLUSIONS/SIGNIFICANCE: The human vaginal bacterial biota is heterogeneous and marked by greater species richness and diversity in women with BV; no species is universally present. Different bacterial species have different associations with the four clinical criteria, which may account for discrepancies often observed between Amsel and Nugent (Gram stain) diagnostic criteria. Several BVABs exhibited race-dependent prevalence when analyzed in separate groups by BV status which may contribute to increased incidence of BV in Black women. Tools developed in this project can be used to study microbial ecology in diverse settings at high resolution.

Genome-wide MyoD binding in skeletal muscle cells: a potential for broad cellular reprogramming.

Recent studies have demonstrated that MyoD initiates a feed-forward regulation of skeletal muscle gene expression, predicting that MyoD binds directly to many genes expressed during differentiation. We have used chromatin immunoprecipitation and high-throughput sequencing to identify genome-wide binding of MyoD in several skeletal muscle cell types. As anticipated, MyoD preferentially binds to a VCASCTG sequence that resembles the in vitro-selected site for a MyoD:E-protein heterodimer, and MyoD binding increases during differentiation at many of the regulatory regions of genes expressed in skeletal muscle. Unanticipated findings were that MyoD was constitutively bound to thousands of additional sites in both myoblasts and myotubes, and that the genome-wide binding of MyoD was associated with regional histone acetylation. Therefore, in addition to regulating muscle gene expression, MyoD binds genome wide and has the ability to broadly alter the epigenome in myoblasts and myotubes.

The effect of self-fertilization, inbreeding depression, and population size on autopolyploid establishment.

The minority cytotype exclusion principle describes how random mating between diploid and autotetraploid cytotypes hinders establishment of the rare cytotype. We present deterministic and stochastic models to ascertain how selfing, inbreeding depression, unreduced gamete production, and finite population size affect minority cytotype exclusion and the establishment of autotetraploids. Results demonstrate that higher selfing rates and lower inbreeding depression in autotetraploids facilitate establishment of autotetraploid populations. Stochastic effects due to finite population size increase the probability of polyploid establishment and decrease the mean time to tetraploid fixation. Our results extend the minority cytotype exclusion principle to include important features of plant reproduction and demonstrate that variation in mating system parameters significantly influences the conditions necessary for polyploid establishment.

Self-fertilization and the escape from pollen limitation in variable pollination environments.

Seed production in many plants is pollen limited, likely because of unpredictable variation in the pollinator environment. One way for plants to escape the consequences of pollinator variability is to evolve mating systems, such as autonomous selfing, that assure reproduction without relying on pollinators. We explore this hypothesis through the construction and analysis of heuristic models of plant population dynamics in seed- or site-limited populations. Our analysis suggests several important points: the familiar rule that inbreeding depression greater than 0.5 maintains outcrossing significantly underestimates the threshold required under pollen limited conditions with prior selfing; variability in the pollination environment erodes the ability of inbreeding depression to maintain outcrossing; and variable pollination environments can result in stable intermediate rates of prior selfing. The results reflect the importance of geometric mean fitness (which in a variable environment is less than the arithmetic mean) in the face of temporal variation.

Plant population dynamics, pollinator foraging, and the selection of self-fertilization.

Many flowering plants rely on pollinators, self-fertilization, or both for reproduction. We model the consequences of these features for plant population dynamics and mating system evolution. Our mating systems-based population dynamics model includes an Allee effect. This often leads to an extinction threshold, defined as a density below which population densities decrease. Reliance on generalist pollinators who primarily visit higher density plant species increases the extinction threshold, whereas autonomous modes of selfing decrease and can eliminate the threshold. Generalist pollinators visiting higher density plant species coupled with autonomous selfing may introduce an effect where populations decreasing in density below the extinction threshold may nonetheless persist through selfing. The extinction threshold and selfing at low density result in populations where individuals adopting a single reproductive strategy exhibit mating systems that depend on population density. The ecological and evolutionary analyses provide a mechanism where prior selfing evolves even though inbreeding depression is greater than one-half. Simultaneous consideration of ecological and evolutionary dynamics confirms unusual features (e.g., evolution into extinction or abrupt increases in population density) implicit in our separate consideration of ecological and evolutionary scenarios. Our analysis has consequences for understanding pollen limitation, reproductive assurance, and the evolution of mating systems.

Explaining phenotypic selection on plant attractive characters: male function, gender balance or ecological context?

It is widely agreed that the flowers of hermaphrodite plants evolve in response to selection acting simultaneously through male and female sexual functions, but we know very little about the pattern of gender-specific selection. We review three current hypotheses for gender-specific selection by viewing them within a single phenotypic selection framework. We compile data from phenotypic selection and manipulative studies and evaluate the fit between empirical data and the hypotheses. In this preliminary analysis, we find that neither the male-function hypothesis nor the gender-balance hypothesis is well supported. However, the context-dependence hypothesis is supported by the documented diversity of gender-specific selection and by evidence that selection through female fertility is significantly correlated with pollen limitation of seed production. Future studies contributing to our understanding of selection through male and female function in plants need to quantify and manipulate the ecological context for reproduction, as well as describe male and female fitness responses to fine-scale trait manipulation.

Coevolution in temporally variable environments.

Many potentially mutualistic interactions are conditional, with selection that varies between mutualism and antagonism over space and time. We develop a genetic model of temporally variable coevolution that incorporates stochastic fluctuations between mutualism and antagonism. We use this model to determine conditions necessary for the coevolution of matching traits between a host and a conditional mutualist. Using an analytical approximation, we show that matching traits will coevolve when the geometric mean interaction is mutualistic. When this condition does not hold, polymorphism and trait mismatching are maintained, and coevolutionary cycles may result. Numerical simulations verify this prediction and suggest that it remains robust in the presence of temporal autocorrelation. These results are compared with those from spatial models with unrestricted movement. The comparisons demonstrate that gene flow is unnecessary for generating empirical patterns predicted by the geographic mosaic theory of coevolution.

Patterns and mechanisms of selection on a family-diagnostic trait: evidence from experimental manipulation and lifetime fitness selection gradients.

Plant traits that show little variation across higher taxa are often used as diagnostic traits, but the reason for the stasis of such traits remains unclear. Wild radish, Raphanus raphanistrum, exhibits tetradynamous stamens (four long and two short, producing a dimorphism in anther height within each flower), as do the vast majority of the more than 3,000 species in the Brassicaceae. Here we examine the hypothesis that selection maintains the stasis of dimorphic anther height by investigating the effects of this trait on pollen removal, seed siring success, and seed set in R. raphanistrum using both experimental and observational methods. Observational selection gradient analysis based on lifetime seed siring success provided evidence for an optimum dimorphism that was greater than zero in one of three years. In both experimentally manipulated and unmanipulated flowers, more pollen was removed in single visits from flowers with less dimorphism. There was no significant effect of anther dimorphism on female fitness (seed set). Therefore, there is some evidence to suggest that selection is maintaining anther dimorphism in wild radish, and that higher male fitness might result from restriction of single-visit pollen removal. We discuss these results in light of pollen presentation theory.

Effects of natural rates of geitonogamy on fruit set in Asclepias speciosa (Apocynaceae): evidence favoring the plant's dilemma.

The role of geitonogamy in the evolution of inflorescence design is not well understood. The plant's dilemma hypothesis proposes that evolution of larger inflorescences is driven by selection for greater pollinator attraction, but constrained by higher rates of geitonogamy experienced by larger inflorescences. Here we investigate the role of geitonogamy on fruit set in natural populations of Asclepias speciosa. We compared fruit set from three pollination treatments: (1) inflorescences bagged before and after receiving 6 hand outcross pollinia (Bag), (2) inflorescences unbagged and receiving 6 hand outcross pollinia (Open), and (3) naturally pollinated inflorescences (Control). The Bag and Open treatments initiated significantly more fruits than the Control. Bag aborted significantly fewer fruits than Open or Control. Fruit set was significantly higher in Bag than Open, and Open had significantly higher fruit set than Control. From these results, we conclude that (1) high rates of geitonogamy significantly increase fruit abortion and reduce fruit set in natural populations of A. speciosa and (2) natural populations are compatible pollen limited. Both findings are consistent with the plant's dilemma hypothesis.

Evolution of a small-muscle polymorphism in lines of house mice selected for high activity levels.

To study the correlated evolution of locomotor behavior and exercise physiology, we conducted an artificial selection experiment. From the outbred Hsd:ICR strain of Mus domesticus, we began eight separate lines, each consisting of 10 breeding pairs. In four of the lines, we used within-family selection to increase voluntary wheel running. The remaining four lines were random-bred (within lines) to serve as controls. Various traits have been monitored to test for correlated responses. Here, we report on organ masses, with emphasis on the triceps surae muscle complex, an important extensor of the ankle. Mice from the selected lines exhibit reduced total body mass, increased relative (mass-corrected) kidney mass, and reduced relative triceps surae mass. In addition, a discrete muscle-mass polymorphism was observed: some individuals had triceps surae that were almost 50% lighter than normal for their body mass. This small-muscle phenotype was observed in only three of the eight lines: in one control line, it has fluctuated in frequency between zero and 10%, whereas in two of the selected lines it has increased in frequency to approximately 50% by generation 22. Data from a set of parents and offspring (generations 23 and 24) are consistent with inheritance as a single autosomal recessive allele. Evidence for the adaptive significance of the small-muscle allele was obtained by fitting multiple-generation data to hierarchical models that include effects of genetic drift and/or selection. The small-muscle allele is estimated to have been present at low frequency (approximately 7%) in the base population, and analysis indicates that strong selection favors the allele in the selected but not control lines. We hypothesize that the small muscles possess functional characteristics and/or that the underlying allele causes pleiotropic effects (e.g., reduced total body mass; increased relative heart, liver, and kidney mass) that facilitate high levels of wheel running. Nevertheless, at generation 22, wheel running of affected individuals did not differ significantly from those with normal-sized muscles, and the magnitude of response to selection has been similar in all four selected lines, indicating that multiple genetic "solutions" are possible in response to selection for high activity levels.

ATTRACTIVE STRUCTURES AND THE STABILITY OF HERMAPHRODITIC SEX EXPRESSION IN FLOWERING PLANTS.

Plant attributes serving both male and female sex functions are thought to favor hermaphroditic breeding systems over systems where sex expression is separated between individuals. Morphological features used by plants to attract biotic pollinators can be considered in this light, because such "attractive" structures are prerequisite for both pollen donation and receipt. A model involving allocation of a limiting resource between attractive structures benefiting both sex functions, structures benefiting male function only, and structures benefiting female function only, is constructed and analyzed using the evolutionary stable strategy (ESS) approach. Contrary to expectation, conditions resulting in large allocation to attractive structures tend to destabilize hermaphroditism, although the effect is slight. Substantial asymmetry in the benefits accruing to male and female fertility for a given investment in pollinator attraction produces more pronounced effects, reducing the parameter space in which hermaphroditism is expected. Results of this model are used to assess the importance of "accelerating gain" curves previously suggested for the evolution of dioecy in flowering plants. Accelerating gains associated with fruit production and dispersal represent female function only, and can be sufficient to initiate evolution toward dioecy. Accelerating gains associated with pollinator attraction may contribute to the evolution of separate sexes only if male function benefits from increased investment in attractive structures at a disproportionate rate compared to female function. Even in this instance, though, hermaphroditism may be evolutionarily stable.

THE DISSOLUTION OF A COMPLEX GENETIC POLYMORPHISM: THE EVOLUTION OF SELF-FERTILIZATION IN TRISTYLOUS EICHHORNIA PANICULATA (PONTEDERIACEAE).

Eichhornia paniculata (Pontederiaceae) displays a wide range of outcrossing levels as a result of the dissolution of the tristylous genetic polymorphism and the evolution of semihomostyly. Population surveys, comparison of fitness components of the style morphs, and computer simulations were used to investigate the breakdown of tristyly and the selective mechanisms responsible for the evolution of self-fertilization. Of 110 populations surveyed in northeast Brazil and Jamaica, 53% were trimorphic, 25% were dimorphic, and 22% were monomorphic for style morph. The short (S) morph was underrepresented in trimorphic populations and absent from nontrimorphic populations. The mid (M) morph predominated in dimorphic populations and was the only morph in monomorphic populations. Stamen modifications promoting selfing, associated with semihomostyle evolution, were largely confined to the M morph. They were rare in trimorphic populations, common in dimorphic populations, and often fixed in monomorphic populations. Stochastic simulations and comparisons of fruit set in natural populations indicate that founder events, population bottlenecks, and lowered fertility of the S morph due to an absence of long-tongued pollinators can each account for loss of the S morph from trimorphic populations. A reduced level of disassortative mating can accentuate the rate at which the S morph is lost by both random and deterministic processes. Nontrimorphic populations occur at the geographical margins of the region surveyed and tend to be smaller and less dense than trimorphic populations. These observations and the higher fruit set of the M morph relative to the L morph in dimorphic populations suggest that reproductive assurance, favoring selfing variants of the M morph under conditions of low pollinator service, has been of primary importance in the origin of most monomorphic populations. Where pollinator service is reliable, however, automatic selection of selfing genes, aided by mating asymmetries between the morphs, can cause the M morph to spread to fixation in dimorphic populations.