Environmental DNA sequencing dataset from Lake Erie algal blooms using Oxford Nanopore MinION.

Here we describe a publicly available environmental DNA (eDNA) sequence dataset, consisting of samples collected from a National Oceanic and Atmospheric Administration (NOAA) Great Lakes Environmental Research Laboratory (GLERL) on Lake Erie. We sequenced samples drawn from before, during, and after a 2019 Microcystis harmful algal bloom (HAB) using 3rd generation sequencing with the Oxford Nanopore MinION device. We classified the eDNA reads taxonomically, and estimated the abundances of all taxa in each sample. While the taxonomic data showed evidence of significant human and E. coli contamination, we found abundant Mycrocystis, especially in the samples drawn from bloom environments. The raw sequence data are available in the Sequence Read Archive (SRA) under accession number PRJNA812770. HABs pose a significant and increasing risk, both to human health and to the Blue Economy, and genomic approaches to early detection promise to help mitigate these risks. As such, this dataset could be of interest to freshwater ecology research teams, or any stakeholders interested in the detection and mitigation of HABs.

Characterization and comparison of innate and adaptive immune responses at vaccine sites in melanoma vaccine clinical trials.

The strength and durability of systemic anti-tumor immune responses induced by cancer vaccines depends on adjuvants to support an immunogenic vaccine site microenvironment (VSME). Adjuvants include water-in-oil emulsions with incomplete Freund's adjuvant (IFA) and combinations of toll-like receptor (TLR) agonists, including a preparation containing TLR4 and TLR9 agonists with QS-21 (AS15). IFA-containing vaccines can promote immune cell accumulation at the VSME, whereas effects of AS15 are largely unexplored. Therefore, we assessed innate and adaptive immune cell accumulation and gene expression at the VSME after vaccination with AS15 and compared to effects with IFA. We hypothesized that AS15 would promote less accumulation of innate and adaptive immune cells at the VSME than IFA vaccines. In two clinical trials, patients with resected high-risk melanoma received either a multipeptide vaccine with IFA or a recombinant MAGE-A3 protein vaccine with AS15. Vaccine site biopsies were obtained after one or multiple vaccines. T cells accumulated early after vaccines with AS15, but this was not durable or of the same magnitude as vaccination in IFA. Vaccines with AS15 increased durable expression of DC- and T cell-related genes, as well as PD-L1 and IDO1, suggesting complex activation and regulation of innate and adaptive immune function with AS15. These changes were generally greater with vaccines containing IFA, but IFA induced reduction in myeloid suppressor cells markers. Evidence of tertiary lymphoid structure (TLS) formation was observed with both adjuvants. Our findings highlight adjuvant-dependent changes in immune features at the VSME that may impact systemic immune responses.

Sexually Dimorphic Crosstalk at the Maternal-Fetal Interface.

CONTEXT: Crosstalk through receptor ligand interactions at the maternal-fetal interface is impacted by fetal sex. This affects placentation in the first trimester and differences in outcomes. Sexually dimorphic signaling at early stages of placentation are not defined. OBJECTIVE: Investigate the impact of fetal sex on maternal-fetal crosstalk. DESIGN: Receptors/ligands at the maternal-fetal surface were identified from sexually dimorphic genes between fetal sexes in the first trimester placenta and defined in each cell type using single-cell RNA-Sequencing (scRNA-Seq). SETTING: Academic institution. SAMPLES: Late first trimester (~10-13 weeks) placenta (fetal) and decidua (maternal) from uncomplicated ongoing pregnancies. MAIN OUTCOME MEASURES: Transcriptomic profiling at tissue and single-cell level; immunohistochemistry of select proteins. RESULTS: We identified 91 sexually dimorphic receptor-ligand pairs across the maternal-fetal interface. We examined fetal sex differences in 5 major cell types (trophoblasts, stromal cells, Hofbauer cells, antigen-presenting cells, and endothelial cells). Ligands from the CC family chemokine ligand (CCL) family were most highly representative in females, with their receptors present on the maternal surface. Sexually dimorphic trophoblast transcripts, Mucin-15 (MUC15) and notum, palmitoleoyl-protein carboxylesterase (NOTUM) were also most highly expressed in syncytiotrophoblasts and extra-villous trophoblasts respectively. Gene Ontology (GO) analysis using sexually dimorphic genes in individual cell types identified cytokine mediated signaling pathways to be most representative in female trophoblasts. Upstream analysis demonstrated TGFB1 and estradiol to affect all cell types, but dihydrotestosterone, produced by the male fetus, was an upstream regulator most significant for the trophoblast population. CONCLUSIONS: Maternal-fetal crosstalk exhibits sexual dimorphism during placentation early in gestation.

Incomplete Freund's adjuvant reduces arginase and enhances Th1 dominance, TLR signaling and CD40 ligand expression in the vaccine site microenvironment.

BACKGROUND: Immunogenicity of cancer vaccines is impacted by adjuvants and schedule, but systematic assessments of their effects have not been performed. Montanide ISA-51, an incomplete Freund's adjuvant (IFA), is used in many vaccine trials, but concerns have been raised about negative effects in murine studies. We found in humans that IFA enhances systemic immune responses and that repeat vaccination at one site (same site vaccination (SSV)) creates tertiary lymphoid structures (TLS) in the vaccine site microenvironment (VSME). We hypothesized that vaccination with peptides+IFA+pICLC or SSV×3 with peptides in IFA would create an immunogenic milieu locally at the VSME, with activated dendritic cells (DC), TLS-associated chemokines and a Th1-dominant VSME. METHODS: Biopsies of the VSME were obtained from participants on two clinical trials who were immunized with multiple melanoma peptides (MELITAC 12.1) in adjuvants comprising IFA and/or the TLR3-agonist pICLC. Biopsies were obtained either a week after one vaccine or a week after SSV×3. Controls included normal skin and skin injected with IFA without peptides. Gene expression analysis was performed by RNAseq. RESULTS: VSME samples were evaluated from 27 patients. One vaccine with peptides in pICLC+IFA enhanced expression of CD80, CD83, CD86 (p<0.01), CD40 and CD40L (p<0.0001) over normal skin; these effects were significantly enhanced for SSV with peptides+IFA. Vaccines containing pICLC increased expression of TBX21 (T-bet) but did not decrease GATA3 over normal skin, whereas SSV with peptides in IFA dramatically enhanced TBX21 and decreased GATA3, with high expression of IFNγ and STAT1. SSV with peptides in IFA also reduced arginase-1 (ARG1) expression and enhanced expression of TLR adapter molecules TICAM-1 (TRIF) and MYD88. Furthermore, SSV with IFA and peptides also enhanced expression of chemokines associated with TLS formation. CONCLUSIONS: These findings suggest that SSV with peptides in IFA enhances CD40L expression by CD4 T cells, supports a Th1 microenvironment, with accumulation of activated and mature DC. Increased expression of TLR adaptor proteins after SSV with peptides in IFA might implicate effects of the skin microbiome. Reduced ARG1 may reflect diminished suppressive myeloid activity in the VSME. TRIAL REGISTRATION NUMBER: (NCT00705640, NCT01585350).

A deletion in Eml1 leads to bilateral subcortical heterotopia in the tish rat.

Children with malformations of cortical development (MCD) are at risk for epilepsy, developmental delays, behavioral disorders, and intellectual disabilities. For a subset of these children, antiseizure medications or epilepsy surgery may result in seizure freedom. However, there are limited options for treating or curing the other conditions, and epilepsy surgery is not an option in all cases of pharmacoresistant epilepsy. Understanding the genetic and neurobiological mechanisms underlying MCD is a necessary step in elucidating novel therapeutic targets. The tish (telencephalic internal structural heterotopia) rat is a unique model of MCD with spontaneous seizures, but the underlying genetic mutation(s) have remained unknown. DNA and RNA-sequencing revealed that a deletion encompassing a previously unannotated first exon markedly diminished Eml1 transcript and protein abundance in the tish brain. Developmental electrographic characterization of the tish rat revealed early-onset of spontaneous spike-wave discharge (SWD) bursts beginning at postnatal day (P) 17. A dihybrid cross demonstrated that the mutant Eml1 allele segregates with the observed dysplastic cortex and the early-onset SWD bursts in monogenic autosomal recessive frequencies. Our data link the development of the bilateral, heterotopic dysplastic cortex of the tish rat to a deletion in Eml1.

Inflammation Triggers Liver X Receptor-Dependent Lipogenesis.

Immune cell function can be modulated by changes in lipid metabolism. Our studies indicate that cholesterol and fatty acid synthesis increases in macrophages between 12 and 18 h after the activation of Toll-like receptors with proinflammatory stimuli and that the upregulation of lipogenesis may contribute to the resolution of inflammation. The inflammation-dependent increase in lipogenesis requires the induction of the liver X receptors, members of the nuclear receptor superfamily of transcription factors, by type I interferons in response to inflammatory signals. Instead of the well-established role for liver X receptors in stimulating cholesterol efflux, we demonstrate that liver X receptors are necessary for the proper resumption of cholesterol synthesis in response to inflammatory signals. Thus, liver X receptors function as bidirectional regulators of cholesterol homeostasis, driving efflux when cholesterol levels are high and facilitating synthesis in response to inflammatory signals. Liver X receptor activity is also required for the proper shutdown of a subset of type I interferon-stimulated genes as inflammation subsides, placing the receptors in a negative-feedback loop that may contribute to the resolution of the inflammatory response.

Comparative Transcriptomic Analysis Identifies a Range of Immunologically Related Functional Elaborations of Lymph Node Associated Lymphatic and Blood Endothelial Cells.

Lymphatic and blood vessels are formed by specialized lymphatic endothelial cells (LEC) and blood endothelial cells (BEC), respectively. These endothelial populations not only form peripheral tissue vessels, but also critical supporting structures in secondary lymphoid organs, particularly the lymph node (LN). Lymph node LEC (LN-LEC) also have been shown to have important immunological functions that are not observed in LEC from tissue lymphatics. LN-LEC can maintain peripheral tolerance through direct presentation of self-antigen via MHC-I, leading to CD8 T cell deletion; and through transfer of self-antigen to dendritic cells for presentation via MHC-II, resulting in CD4 T cell anergy. LN-LEC also can capture and archive foreign antigens, transferring them to dendritic cells for maintenance of memory CD8 T cells. The molecular basis for these functional elaborations in LN-LEC remain largely unexplored, and it is also unclear whether blood endothelial cells in LN (LN-BEC) might express similar enhanced immunologic functionality. Here, we used RNA-Seq to compare the transcriptomic profiles of freshly isolated murine LEC and BEC from LN with one another and with freshly isolated LEC from the periphery (diaphragm). We show that LN-LEC, LN-BEC, and diaphragm LEC (D-LEC) are transcriptionally distinct from one another, demonstrating both lineage and tissue-specific functional specializations. Surprisingly, tissue microenvironment differences in gene expression profiles were more numerous than those determined by endothelial cell lineage specification. In this regard, both LN-localized endothelial cell populations show a variety of functional elaborations that suggest how they may function as antigen presenting cells, and also point to as yet unexplored roles in both positive and negative regulation of innate and adaptive immune responses. The present work has defined in depth gene expression differences that point to functional specializations of endothelial cell populations in different anatomical locations, but especially the LN. Beyond the analyses provided here, these data are a resource for future work to uncover mechanisms of endothelial cell functionality.

In Vivo Gene Essentiality and Metabolism in Bordetella pertussis.

Bordetella pertussis is the causative agent of whooping cough, a serious respiratory illness affecting children and adults, associated with prolonged cough and potential mortality. Whooping cough has reemerged in recent years, emphasizing a need for increased knowledge of basic mechanisms of B. pertussis growth and pathogenicity. While previous studies have provided insight into in vitro gene essentiality of this organism, very little is known about in vivo gene essentiality, a critical gap in knowledge, since B. pertussis has no previously identified environmental reservoir and is isolated from human respiratory tract samples. We hypothesize that the metabolic capabilities of B. pertussis are especially tailored to the respiratory tract and that many of the genes involved in B. pertussis metabolism would be required to establish infection in vivo In this study, we generated a diverse library of transposon mutants and then used it to probe gene essentiality in vivo in a murine model of infection. Using the CON-ARTIST pipeline, 117 genes were identified as conditionally essential at 1 day postinfection, and 169 genes were identified as conditionally essential at 3 days postinfection. Most of the identified genes were associated with metabolism, and we utilized two existing genome-scale metabolic network reconstructions to probe the effects of individual essential genes on biomass synthesis. This analysis suggested a critical role for glucose metabolism and lipooligosaccharide biosynthesis in vivo This is the first genome-wide evaluation of in vivo gene essentiality in B. pertussis and provides tools for future exploration.IMPORTANCE Our study describes the first in vivo transposon sequencing (Tn-seq) analysis of B. pertussis and identifies genes predicted to be essential for in vivo growth in a murine model of intranasal infection, generating key resources for future investigations into B. pertussis pathogenesis and vaccine design.

Sex differences in the late first trimester human placenta transcriptome.

BACKGROUND: Development of the placenta during the late first trimester is critical to ensure normal growth and development of the fetus. Developmental differences in this window such as sex-specific variation are implicated in later placental disease states, yet gene expression at this time is poorly understood. METHODS: RNA-sequencing was performed to characterize the transcriptome of 39 first trimester human placentas using chorionic villi following genetic testing (17 females, 22 males). Gene enrichment analysis was performed to find enriched canonical pathways and gene ontologies in the first trimester. DESeq2 was used to find sexually dimorphic gene expression. Patient demographics were analyzed for sex differences in fetal weight at time of chorionic villus sampling and birth. RESULTS: RNA-sequencing analyses detected 14,250 expressed genes, with chromosome 19 contributing the greatest proportion (973/2852, 34.1% of chromosome 19 genes) and Y chromosome contributing the least (16/568, 2.8%). Several placenta-enriched genes as well as histone-coding genes were identified to be unique to the first trimester and common to both sexes. Further, we identified 58 genes with significantly different expression between males and females: 25 X-linked, 15 Y-linked, and 18 autosomal genes. Genes that escape X inactivation were highly represented (59.1%) among X-linked genes upregulated in females. Many genes differentially expressed by sex consisted of X/Y gene pairs, suggesting that dosage compensation plays a role in sex differences. These X/Y pairs had roles in parallel, ancient canonical pathways important for eukaryotic cell growth and survival: chromatin modification, transcription, splicing, and translation. CONCLUSIONS: This study is the first characterization of the late first trimester placenta transcriptome, highlighting similarities and differences among the sexes in ongoing human pregnancies resulting in live births. Sexual dimorphism may contribute to pregnancy outcomes, including fetal growth and birth weight, which was seen in our cohort, with males significantly heavier than females at birth. This transcriptome provides a basis for development of early diagnostic tests of placental function that can indicate overall pregnancy heath, fetal-maternal health, and long-term adult health.

Correction: Prevalence and extent of heteroresistance by next generation sequencing of multidrug-resistant tuberculosis.

[This corrects the article DOI: 10.1371/journal.pone.0176522.].

Prevalence and extent of heteroresistance by next generation sequencing of multidrug-resistant tuberculosis.

Amplicon-based Next Generation Sequencing (NGS) is an emerging method for Mycobacterium tuberculosis drug susceptibility testing (DST) but has not been well described. We examined 158 clinical multidrug-resistant M. tuberculosis isolates via NGS of 11 resistance-associated gene regions covering 3519 nucleotides. Across these gene regions, complete resistance or heteroresistance (defined as 1%-99% mutation) was present in at least one isolate in 6.3% of loci. The number of isolates with heteroresistance was highest for gyrA codon 94, rpoB codons 526 and 531, and embB codons 306, 372 and 406 (range 11-26% of isolates exhibited heteroresistance). 57% of MDR strains had heteroresistance of one or more recognized resistance-associated mutation. Heteroresistant loci generally exhibited high or low degrees of mutation (>90% or <10%). The deep sensitivity of NGS for detecting low level pncA heteroresistance appeared to improve genotypic-phenotypic PZA susceptibility correlations over that of Sanger. NGS demonstrates that heteroresistance in TB in the regions of key genes is common and will need to be bioinformatically managed. The clinical significance of such heteroresistance is unclear, and further study of pncA should be pursued.

Targetable T-type Calcium Channels Drive Glioblastoma.

Glioblastoma (GBM) stem-like cells (GSC) promote tumor initiation, progression, and therapeutic resistance. Here, we show how GSCs can be targeted by the FDA-approved drug mibefradil, which inhibits the T-type calcium channel Cav3.2. This calcium channel was highly expressed in human GBM specimens and enriched in GSCs. Analyses of the The Cancer Genome Atlas and REMBRANDT databases confirmed upregulation of Cav3.2 in a subset of tumors and showed that overexpression associated with worse prognosis. Mibefradil treatment or RNAi-mediated attenuation of Cav3.2 was sufficient to inhibit the growth, survival, and stemness of GSCs and also sensitized them to temozolomide chemotherapy. Proteomic and transcriptomic analyses revealed that Cav3.2 inhibition altered cancer signaling pathways and gene transcription. Cav3.2 inhibition suppressed GSC growth in part by inhibiting prosurvival AKT/mTOR pathways and stimulating proapoptotic survivin and BAX pathways. Furthermore, Cav3.2 inhibition decreased expression of oncogenes (PDGFA, PDGFB, and TGFB1) and increased expression of tumor suppressor genes (TNFRSF14 and HSD17B14). Oral administration of mibefradil inhibited growth of GSC-derived GBM murine xenografts, prolonged host survival, and sensitized tumors to temozolomide treatment. Our results offer a comprehensive characterization of Cav3.2 in GBM tumors and GSCs and provide a preclinical proof of concept for repurposing mibefradil as a mechanism-based treatment strategy for GBM. Cancer Res; 77(13); 3479-90. ©2017 AACR.

Investigation of Genetic Variants Associated with Alzheimer Disease in Parkinson Disease Cognition.

BACKGROUND: Meta-analysis of genome-wide association studies have implicated multiple single nucleotide polymorphisms (SNPs) and associated genes with Alzheimer disease. The role of these SNPs in cognitive impairment in Parkinson disease (PD) remains incompletely evaluated. OBJECTIVE: The objective of this study was to test alleles associated with risk of Alzheimer disease for association with cognitive impairment in Parkinson disease (PD). METHODS: Two datasets with PD subjects accessed through the NIH database of Genotypes and Phenotypes contained both single nucleotide polymorphism (SNP) arrays and mini-mental state exam (MMSE) scores. Genetic data underwent rigorous quality control and we selected SNPs for genes associated with AD other than APOE. We constructed logistic regression and ordinal regression models, adjusted for sex, age at MMSE, and duration of PD, to assess the association between selected SNPs and MMSE score. RESULTS: In one dataset, PICALM rs3851179 was associated with cognitive impairment (MMSE <  24) in PD subjects > 70 years old (OR = 2.3; adjusted p-value = 0.017; n = 250) but not in PD subjects ≤ 70 years old. CONCLUSIONS: Our finding suggests that PICALM rs3851179 could contribute to cognitive impairment in older patients with PD. It is important that future studies consider the interaction of age and genetic risk factors in the development of cognitive impairment in PD.

Species matter: the role of competition in the assembly of congeneric bacteria.

Interspecific competition is an important driver of community assembly in plants and animals, but phylogenetic evidence for interspecific competition in bacterial communities has been elusive. This could indicate that other processes such as habitat filtering or neutral processes are more important in bacterial community assembly. Alternatively, this could be a consequence of the lack of a consistent and meaningful species definition in bacteria. We hypothesize that competition in bacterial community assembly has gone undetected at least partly because overly broad measures of bacterial diversity units were used in previous studies. First, we tested our hypothesis in a simulation where we showed that how species are defined can dramatically affect whether phylogenetic overdispersion (a signal consistent with competitive exclusion) will be detected. Second, we demonstrated that using finer-scale Operational Taxonomic Units (OTUs) (with more stringent 16S rRNA sequence identity cutoffs or based on fast-evolving protein coding genes) in natural populations revealed previously undetected overdispersion. Finally, we argue that bacterial ecotypes, diversity units incorporating ecological and evolutionary theory, are superior to OTUs for the purpose of studying community assembly.

Surprisingly extensive mixed phylogenetic and ecological signals among bacterial Operational Taxonomic Units.

The lack of a consensus bacterial species concept greatly hampers our ability to understand and organize bacterial diversity. Operational taxonomic units (OTUs), which are clustered on the basis of DNA sequence identity alone, are the most commonly used microbial diversity unit. Although it is understood that OTUs can be phylogenetically incoherent, the degree and the extent of the phylogenetic inconsistency have not been explicitly studied. Here, we tested the phylogenetic signal of OTUs in a broad range of bacterial genera from various phyla. Strikingly, we found that very few OTUs were monophyletic, and many showed evidence of multiple independent origins. Using previously established bacterial habitats as benchmarks, we showed that OTUs frequently spanned multiple ecological habitats. We demonstrated that ecological heterogeneity within OTUs is caused by their phylogenetic inconsistency, and not merely due to 'lumping' of taxa resulting from using relaxed identity cut-offs. We argue that ecotypes, as described by the Stable Ecotype Model, are phylogenetically and ecologically more consistent than OTUs and therefore could serve as an alternative unit for bacterial diversity studies. In addition, we introduce QuickES, a new wrapper program for the Ecotype Simulation algorithm, which is capable of demarcating ecotypes in data sets with tens of thousands of sequences.

Global marine bacterial diversity peaks at high latitudes in winter.

Genomic approaches to characterizing bacterial communities are revealing significant differences in diversity and composition between environments. But bacterial distributions have not been mapped at a global scale. Although current community surveys are way too sparse to map global diversity patterns directly, there is now sufficient data to fit accurate models of how bacterial distributions vary across different environments and to make global scale maps from these models. We apply this approach to map the global distributions of bacteria in marine surface waters. Our spatially and temporally explicit predictions suggest that bacterial diversity peaks in temperate latitudes across the world's oceans. These global peaks are seasonal, occurring 6 months apart in the two hemispheres, in the boreal and austral winters. This pattern is quite different from the tropical, seasonally consistent diversity patterns observed for most macroorganisms. However, like other marine organisms, surface water bacteria are particularly diverse in regions of high human environmental impacts on the oceans. Our maps provide the first picture of bacterial distributions at a global scale and suggest important differences between the diversity patterns of bacteria compared with other organisms.

Speedy speciation in a bacterial microcosm: new species can arise as frequently as adaptations within a species.

Microbiologists are challenged to explain the origins of enormous numbers of bacterial species worldwide. Contributing to this extreme diversity may be a simpler process of speciation in bacteria than in animals and plants, requiring neither sexual nor geographical isolation between nascent species. Here, we propose and test a novel hypothesis for the extreme diversity of bacterial species-that splitting of one population into multiple ecologically distinct populations (cladogenesis) may be as frequent as adaptive improvements within a single population's lineage (anagenesis). We employed a set of experimental microcosms to address the relative rates of adaptive cladogenesis and anagenesis among the descendants of a Bacillus subtilis clone, in the absence of competing species. Analysis of the evolutionary trajectories of genetic markers indicated that in at least 7 of 10 replicate microcosm communities, the original population founded one or more new, ecologically distinct populations (ecotypes) before a single anagenetic event occurred within the original population. We were able to support this inference by identifying putative ecotypes formed in these communities through differences in genetic marker association, colony morphology and microhabitat association; we then confirmed the ecological distinctness of these putative ecotypes in competition experiments. Adaptive mutations leading to new ecotypes appeared to be about as common as those improving fitness within an existing ecotype. These results suggest near parity of anagenesis and cladogenesis rates in natural populations that are depauperate of bacterial diversity.

Lineage-dependent ecological coherence in bacteria.

Bacteria comprise an essential element of all ecosystems, including those present on and within the human body. Understanding bacterial diversity therefore offers enormous scientific and medical benefit, but significant questions remain regarding how best to characterize that diversity and organize it into biologically meaningful units. Bacterial communities are routinely characterized based on the relative abundances of taxa at the genus or even the phylum level, but the ecological coherence of these high-level taxonomic units is uncertain. Using human microbiota from the skin and gut as our model systems, we tested the ecological coherence of bacteria by investigating the habitat associations of bacteria at all levels of the taxonomic hierarchy. We observed four distinct taxonomic patterns of habitat association, reflecting different levels of ecological coherence among taxa. Our results support the hypothesis that deep-branch bacterial clades could be ecologically coherent and suggest that the phylogenetic depth of ecological coherence varies among the bacterial lineages and is an important factor to consider in studies of human microbiome associations.

Ecology of speciation in the genus Bacillus.

Microbial ecologists and systematists are challenged to discover the early ecological changes that drive the splitting of one bacterial population into two ecologically distinct populations. We have aimed to identify newly divergent lineages ("ecotypes") bearing the dynamic properties attributed to species, with the rationale that discovering their ecological differences would reveal the ecological dimensions of speciation. To this end, we have sampled bacteria from the Bacillus subtilis-Bacillus licheniformis clade from sites differing in solar exposure and soil texture within a Death Valley canyon. Within this clade, we hypothesized ecotype demarcations based on DNA sequence diversity, through analysis of the clade's evolutionary history by Ecotype Simulation (ES) and AdaptML. Ecotypes so demarcated were found to be significantly different in their associations with solar exposure and soil texture, suggesting that these and covarying environmental parameters are among the dimensions of ecological divergence for newly divergent Bacillus ecotypes. Fatty acid composition appeared to contribute to ecotype differences in temperature adaptation, since those ecotypes with more warm-adapting fatty acids were isolated more frequently from sites with greater solar exposure. The recognized species and subspecies of the B. subtilis-B. licheniformis clade were found to be nearly identical to the ecotypes demarcated by ES, with a few exceptions where a recognized taxon is split at most into three putative ecotypes. Nevertheless, the taxa recognized do not appear to encompass the full ecological diversity of the B. subtilis-B. licheniformis clade: ES and AdaptML identified several newly discovered clades as ecotypes that are distinct from any recognized taxon.

The origins of ecological diversity in prokaryotes.

The urkingdoms and major divisions of prokaryotes are enormously diverse in their metabolic capabilities and membrane architectures. These ancient differences likely have a strong influence on the kinds of ecological adaptations that may evolve today. Some ecological transitions have been identified as having occurred primarily in the distant past, including transitions between saline and non-saline habitats. At the microevolutionary level, the likely existence of a billion prokaryotic species challenges microbiologists to determine what might promote rapid speciation in prokaryotes, and to identify the ecological dimensions upon which new species diverge and by which they may coexist. Rapid speciation in prokaryotes is fostered by several unique properties of prokaryotic genetic exchange, including their propensity to acquire novel gene loci by horizontal genetic transfer, as well as the rarity of their genetic exchange, which allows speciation by ecological divergence alone, without a requirement for sexual isolation. The ecological dimensions of prokaryotic speciation may be identified by comparing the ecology of the most newly divergent, ecologically distinct populations (ecotypes). This program is challenged by our ignorance of the physiological and ecological features most likely responsible for adaptive divergence between closely related ecotypes in any given clade. This effort will require development of universal approaches to hypothesize demarcations of ecotypes, and to confirm and characterize their ecological distinctness, without prior knowledge of a given clade's ecology.