Modeling Pulsed Evolution and Time-Independent Variation Improves the Confidence Level of Ancestral and Hidden State Predictions.

Ancestral state reconstruction is not only a fundamental tool for studying trait evolution, but also very useful for predicting the unknown trait values (hidden states) of extant species. A well-known problem in ancestral and hidden state predictions is that the uncertainty associated with predictions can be so large that predictions themselves are of little use. Therefore, for meaningful interpretation of predicted traits and hypothesis testing, it is prudent to accurately assess the uncertainty of the predictions. Commonly used constant-rate Brownian motion (BM) model fails to capture the complexity of tempo and mode of trait evolution in nature, making predictions under the BM model vulnerable to lack-of-fit errors from model misspecification. Using empirical data (mammalian body size and bacterial genome size), we show that the distribution of residual Z-scores under the BM model is neither homoscedastic nor normal as expected. Consequently, the 95% confidence intervals of predicted traits are so unreliable that the actual coverage probability ranges from 33% (strongly permissive) to 100% (strongly conservative). Alternative methods such as BayesTraits and StableTraits that allow variable rates in evolution improve the predictions but are computationally expensive. Here, we develop Reconstructing Ancestral State under Pulsed Evolution in R by Gaussian Decomposition (RasperGade), a method of ancestral and hidden state prediction that uses the Levy process to explicitly model gradual evolution, pulsed evolution, and time-independent variation. Using the same empirical data, we show that RasperGade outperforms both BayesTraits and StableTraits in providing reliable confidence estimates and is orders-of-magnitude faster. Our results suggest that, when predicting the ancestral and hidden states of continuous traits, the rate variation should always be assessed and the quality of confidence estimates should always be examined. [Bacterial genomic traits; model misspecification; trait evolution.].

Influence of Binary Toxin Gene Detection and Decreased Susceptibility to Antibiotics among Clostridioides difficile Strains on Disease Severity: a Single-Center Study.

Clostridioides difficile infection (CDI) is the fifth leading cause of death from nonmalignant gastrointestinal disease in the United States. The contribution of resistance to C. difficile-active antibiotics to the outcomes of CDI is unclear. We evaluated the antimicrobial susceptibility of C. difficile isolates in a U.S. hospital and determined associations of clinical variables and binary toxin positivity with antibiotic resistance. C. difficile spores were cultured from fecal specimens of adult patients with CDI for genotyping and antimicrobial susceptibility assay (for clindamycin [CLI], fidaxomicin [FDX], metronidazole [MTZ], moxifloxacin [MXF], tigecycline [TGC], and vancomycin [VAN]). Electronic medical records were reviewed for clinical data extraction. Ninety-seven of 130 (75%) fecal samples grew toxigenic C. difficile in culture. Most of the isolates were tcdA+ tcdB+ cdtB- (80.4%), and 18.6% and 1% were tcdA+ tcdB+ cdtB+ and tcdA-tcdB+ cdtB+, respectively. Susceptibility to VAN, MTZ, FDX, TGC, MXF, and CLI was 96%, 94%, 100%, 100%, 8%, and 79%, respectively. Six isolates, all cdtB positive and belonging to the 027 ribotype, were resistant to VAN and/or MTZ. Higher MICs were found in isolates with a mutation in the VAN-related resistance gene vanR, but not vanS. In addition, cdtB+ isolates exhibited higher MICs of VAN, MTZ, TGC, CLI, and MXF compared to cdtB- strains. Patients with greater intestinal inflammation or severe disease were more likely to be infected with cdtB+ strains. Decreased susceptibility to antibiotics is not directly associated with either severe or recurrent CDI. However, antimicrobial susceptibility of C. difficile is decreased in strains positive for the binary toxin gene.

Rationale, design and baseline characteristics of the Microbiome and Insulin Longitudinal Evaluation Study (MILES).

AIM: To investigate the role of the gut microbiome in regulating key insulin homeostasis traits (insulin sensitivity, insulin secretion and insulin clearance) whose dysfunction leads to type 2 diabetes (T2D). MATERIALS AND METHODS: The Microbiome and Insulin Longitudinal Evaluation Study (MILES) focuses on African American and non-Hispanic white participants aged 40-80 years without diabetes. Three study visits are planned (at baseline, 15 and 30 months). Baseline measurements include assessment of the stool microbiome and administration of an oral glucose tolerance test, which will yield indexes of insulin sensitivity, insulin secretion and insulin clearance. The gut microbiome profile (composition and function) will be determined using whole metagenome shotgun sequencing along with analyses of plasma short chain fatty acids. Additional data collected include dietary history, sociodemographic factors, health habits, anthropometry, medical history, medications and family history. Most assessments are repeated 15 and 30 months following baseline. RESULTS: After screening 875 individuals, 129 African American and 224 non-Hispanic white participants were enrolled. At baseline, African American participants have higher blood pressure, weight, body mass index, waist and hip circumferences but similar waist-hip ratio compared with the non-Hispanic white participants. On average, African American participants are less insulin-sensitive and have higher acute insulin secretion and lower insulin clearance. CONCLUSIONS: The longitudinal design and robust characterization of potential mediators will allow for the assessment of glucose and insulin homeostasis and gut microbiota as they change over time, improving our ability to discern causal relationships between the microbiome and the insulin homeostasis traits whose deterioration determines T2D, setting the stage for future microbiome-directed therapies to prevent and treat T2D.

Innate Immune Response and Outcome of Clostridium difficile Infection Are Dependent on Fecal Bacterial Composition in the Aged Host.

Background: Clostridium difficile infection (CDI) is a serious threat for an aging population. Using an aged mouse model, we evaluated the effect of age and the roles of innate immunity and intestinal microbiota. Methods: Aged (18 months) and young (8 weeks) mice were infected with C difficile, and disease severity, immune response, and intestinal microbiome were compared. The same experiment was repeated with intestinal microbiota exchange between aged and young mice before infection. Results: Higher mortality was observed in aged mice with weaker neutrophilic mobilization in blood and intestinal tissue and depressed proinflammatory cytokines in early infection. Microbiota exchange improved survival and early immune response in aged mice. Microbiome analysis revealed that aged mice have significant deficiencies in Bacteroidetes phylum and, specifically, Bacteroides, Alistipes, and rc4-4 genera, which were replenished by cage switching. Conclusions: Microbiota-dependent alteration in innate immune response early on during infection may explain poor outcome in aged host with CDI.

Evidence for the Role of the Cecal Microbiome in Maintenance of Immune Regulation and Homeostasis.

OBJECTIVE (S): Our objective was to investigate alterations in the cecal microbial composition during the development of type 1 diabetes (T1D) with or without IgM therapy, and correlate these alterations with the corresponding immune profile. METHODS: (1) Female nonobese diabetic (NOD) mice treated with IgM or saline (n = 20/group) were divided into 5-week-old nondiabetic; 9 to 12-week-old prehyperglycemic stage-1; ≥13-week-old prehyperglycemic stage-2; and diabetic groups. 16S rRNA libraries were prepared from bacterial DNA and deep-sequenced. (2) New-onset diabetic mice were treated with IgM (200 µg on Days 1, 3, and 5) and their blood glucose monitored for 2 months. RESULTS: Significant dysbiosis was observed in the cecal microbiome with the progression of T1D development. The alteration in microbiome composition was characterized by an increase in the bacteroidetes:firmicutes ratio. In contrast, IgM conserved normal bacteroidetes:firmicutes ratio and this effect was long-lasting. Furthermore, oral gavage using cecal content from IgM-treated mice significantly diminished the incidence of diabetes compared with controls, indicating that IgM specifically affected mucosa-associated microbes, and that the affect was causal and not an epiphenomenon. Also, regulatory immune cell populations (myeloid-derived suppressor cells and regulatory T cells) were expanded and insulin autoantibody production diminished in the IgM-treated mice. In addition, IgM therapy reversed hyperglycemia in 70% of new-onset diabetic mice (n = 10) and the mice remained normoglycemic for the entire post-treatment observation period. CONCLUSIONS: The cecal microbiome appears to be important in maintaining immune homeostasis and normal immune responses.

Development, diet and dynamism: longitudinal and cross-sectional predictors of gut microbial communities in wild baboons.

Gut bacterial communities play essential roles in host biology, but to date we lack information on the forces that shape gut microbiota between hosts and over time in natural populations. Understanding these forces in wild primates provides a valuable comparative context that enriches scientific perspectives on human gut microbiota. To this end, we tested predictors of gut microbial composition in a well-studied population of wild baboons. Using cross-sectional and longitudinal samples collected over 13 years, we found that baboons harbour gut microbiota typical of other omnivorous primates, albeit with an especially high abundance of Bifidobacterium. Similar to previous work in humans and other primates, we found strong effects of both developmental transitions and diet on gut microbial composition. Strikingly, baboon gut microbiota appeared to be highly dynamic such that samples collected from the same individual only a few days apart were as different from each other as samples collected over 10 years apart. Despite the dynamic nature of baboon gut microbiota, we identified a set of core taxa that is common among primates, supporting the hypothesis that microbiota codiversify with their host species. Our analysis identified two tentative enterotypes in adult baboons that differ from those of humans and chimpanzees.

Does adaptive radiation of a host lineage promote ecological diversity of its bacterial communities? A test using gut microbiota of Anolis lizards.

Adaptive radiations provide unique opportunities to test whether and how recent ecological and evolutionary diversification of host species structures the composition of entire bacterial communities. We used 16S rRNA gene sequencing of faecal samples to test for differences in the gut microbiota of six species of Puerto Rican Anolis lizards characterized by the evolution of distinct 'ecomorphs' related to differences in habitat use. We found substantial variation in the composition of the microbiota within each species and ecomorph (trunk-crown, trunk-ground, grass-bush), but no differences in bacterial alpha diversity among species or ecomorphs. Beta diversity analyses revealed subtle but significant differences in bacterial composition related to host phylogeny and species, but these differences were not consistently associated with Anolis ecomorph. Comparison of a trunk-ground species from this clade (A. cristatellus) with a distantly related member of the same ecomorph class (A. sagrei) where the two species have been introduced and are now sympatric in Florida revealed pronounced differences in the alpha diversity and beta diversity of their microbiota despite their ecological similarity. Comparisons of these populations with allopatric conspecifics also revealed geographic differences in bacterial alpha diversity and beta diversity within each species. Finally, we observed high intraindividual variation over time and strong effects of a simplified laboratory diet on the microbiota of A. sagrei. Collectively, our results indicate that bacterial communities are only weakly shaped by the diversification of their lizard hosts due to the strikingly high levels of bacterial diversity and variation observed within Anolis species.

Cytonuclear interactions and relaxed selection accelerate sequence evolution in organelle ribosomes.

Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus (Silene) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.

Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates.

Genome size and complexity vary tremendously among eukaryotic species and their organelles. Comparisons across deeply divergent eukaryotic lineages have suggested that variation in mutation rates may explain this diversity, with increased mutational burdens favoring reduced genome size and complexity. The discovery that mitochondrial mutation rates can differ by orders of magnitude among closely related angiosperm species presents a unique opportunity to test this hypothesis. We sequenced the mitochondrial genomes from two species in the angiosperm genus Silene with recent and dramatic accelerations in their mitochondrial mutation rates. Contrary to theoretical predictions, these genomes have experienced a massive proliferation of noncoding content. At 6.7 and 11.3 Mb, they are by far the largest known mitochondrial genomes, larger than most bacterial genomes and even some nuclear genomes. In contrast, two slowly evolving Silene mitochondrial genomes are smaller than average for angiosperms. Consequently, this genus captures approximately 98% of known variation in organelle genome size. The expanded genomes reveal several architectural changes, including the evolution of complex multichromosomal structures (with 59 and 128 circular-mapping chromosomes, ranging in size from 44 to 192 kb). They also exhibit a substantial reduction in recombination and gene conversion activity as measured by the relative frequency of alternative genome conformations and the level of sequence divergence between repeat copies. The evolution of mutation rate, genome size, and chromosome structure can therefore be extremely rapid and interrelated in ways not predicted by current evolutionary theories. Our results raise the hypothesis that changes in recombinational processes, including gene conversion, may be a central force driving the evolution of both mutation rate and genome structure.

A phylogeny-driven genomic encyclopaedia of Bacteria and Archaea.

Sequencing of bacterial and archaeal genomes has revolutionized our understanding of the many roles played by microorganisms. There are now nearly 1,000 completed bacterial and archaeal genomes available, most of which were chosen for sequencing on the basis of their physiology. As a result, the perspective provided by the currently available genomes is limited by a highly biased phylogenetic distribution. To explore the value added by choosing microbial genomes for sequencing on the basis of their evolutionary relationships, we have sequenced and analysed the genomes of 56 culturable species of Bacteria and Archaea selected to maximize phylogenetic coverage. Analysis of these genomes demonstrated pronounced benefits (compared to an equivalent set of genomes randomly selected from the existing database) in diverse areas including the reconstruction of phylogenetic history, the discovery of new protein families and biological properties, and the prediction of functions for known genes from other organisms. Our results strongly support the need for systematic 'phylogenomic' efforts to compile a phylogeny-driven 'Genomic Encyclopedia of Bacteria and Archaea' in order to derive maximum knowledge from existing microbial genome data as well as from genome sequences to come.

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.

NtERF32: a non-NIC2 locus AP2/ERF transcription factor required in jasmonate-inducible nicotine biosynthesis in tobacco.

Nicotine biosynthesis in tobacco (Nicotiana tabacum L.) is highly regulated by jasmonic acid (JA). Two nuclear loci, A and B (renamed NIC1 and NIC2) have been identified that mediate JA-inducible nicotine formation and total alkaloid accumulation. NIC2 was recently shown to be a cluster of seven genes encoding Apetala2/Ethylene-Response Factor (AP2/ERF)-domain transcription factors (TFs) in Group IX of the tobacco AP2/ERF family. Here we report the characterization of several NtERF TF genes that are not within the NIC2 locus, but required for methyl JA (MeJA)-induced nicotine biosynthesis. Expression of NtERF1, NtERF32, and NtERF121 is rapidly induced (<30 min) by MeJA treatment. All three of these TFs specifically bind the GCC box-like element of the GAG motif required for MeJA-induced transcription of NtPMT1a, a gene encoding putrescine N-methyltransferase, the first committed step in the synthesis of the nicotine pyrrolidine ring. Ectopic overexpression of NtERF32 increases expression of NtPMT1a in vivo and elevates total alkaloid contents, whereas RNAi-mediated knockdown of NtERF32 reduces the mRNA levels of multiple genes in the nicotine biosynthetic pathway including NtPMT1a and quinolinate phosphoribosyltransferase (NtQPT2), and lowers nicotine and total alkaloid levels. We conclude that NtERF32 and related ERF genes are important non-NIC2 locus associated transcriptional regulators of nicotine and total alkaloid formation.

A phylum-level bacterial phylogenetic marker database.

Large-scale, genome-level molecular phylogenetic analyses present both opportunities and challenges for bacterial evolutionary and ecological studies. We constructed a phylum-level bacterial phylogenetic marker database by surveying all complete bacterial genomes and identifying single-copy genes that were widely distributed in each of the 20 bacterial phyla. We showed that phylum trees made using these markers were highly resolved and were more robust than the bacterial genome tree based on 31 universal bacterial marker genes. In addition, using the Global Ocean Sampling data set as an example, we demonstrated that the expanded marker database greatly increased the power of metagenomic phylotyping. We incorporated the database into an automated phylogenomic inference application (Phyla-AMPHORA) and made it publicly available. We believe that this centralized resource should have broad applicability in bacterial systematics, phylogenetics, and metagenomic studies.

Genomic heterogeneity and ecological speciation within one subspecies of Bacillus subtilis.

Closely related bacterial genomes usually differ in gene content, suggesting that nearly every strain in nature may be ecologically unique. We have tested this hypothesis by sequencing the genomes of extremely close relatives within a recognized taxon and analyzing the genomes for evidence of ecological distinctness. We compared the genomes of four Death Valley isolates plus the laboratory strain W23, all previously classified as Bacillus subtilis subsp. spizizenii and hypothesized through multilocus analysis to be members of the same ecotype (an ecologically homogeneous population), named putative ecotype 15 (PE15). These strains showed a history of positive selection on amino acid sequences in 38 genes. Each of the strains was under a different regimen of positive selection, suggesting that each strain is ecologically unique and represents a distinct ecological speciation event. The rate of speciation appears to be much faster than can be resolved with multilocus sequencing. Each PE15 strain contained unique genes known to confer a function for bacteria. Remarkably, no unique gene conferred a metabolic system or subsystem function that was not already present in all the PE15 strains sampled. Thus, the origin of ecotypes within this clade shows no evidence of qualitative divergence in the set of resources utilized. Ecotype formation within this clade is consistent with the nanoniche model of bacterial speciation, in which ecotypes use the same set of resources but in different proportions, and genetic cohesion extends beyond a single ecotype to the set of ecotypes utilizing the same resources.

A recurring syndrome of accelerated plastid genome evolution in the angiosperm tribe Sileneae (Caryophyllaceae).

In flowering plants, plastid genomes are generally conserved, exhibiting slower rates of sequence evolution than the nucleus and little or no change in structural organization. However, accelerated plastid genome evolution has occurred in scattered angiosperm lineages. For example, some species within the genus Silene have experienced a suite of recent changes to their plastid genomes, including inversions, shifts in inverted repeat boundaries, large indels, intron losses, and rapid rates of amino acid sequence evolution in a subset of protein genes, with the most extreme divergence occurring in the protease gene clpP. To investigate the relationship between the rates of sequence and structural evolution, we sequenced complete plastid genomes from three species (Silene conoidea, S. paradoxa, and Lychnis chalcedonica), representing independent lineages within the tribe Sileneae that were previously shown to have accelerated rates of clpP evolution. We found a high degree of parallel evolution. Elevated rates of amino acid substitution have occurred repeatedly in the same subset of plastid genes and have been accompanied by a recurring pattern of structural change, including cases of identical inversions and intron loss. This "syndrome" of changes was not observed in the closely related outgroup Agrostemma githago or in the more slowly evolving Silene species that were sequenced previously. Although no single mechanism has yet been identified to explain the correlated suite of changes in plastid genome sequence and structure that has occurred repeatedly in angiosperm evolution, we discuss a possible mixture of adaptive and non-adaptive forces that may be responsible.

Consistent spatial patterns in microbial taxa of red squirrel gut microbiomes.

Gut microbiomes are diverse ecosystems whose drivers of variation remain largely unknown, especially in time and space. We analysed a dataset with over 900 red squirrel (Tamiasciurus hudsonicus) gut microbiome samples to identify the drivers of gut microbiome composition in this territorial rodent. The large-scale spatiotemporal replication in the data analysed was an essential component of understanding the assembly of these microbial communities. We identified that the spatial location of the sampled squirrels in their local environment is a key contributor to gut microbial community composition. The non-core gut microbiome (present in less than 75% of gut microbiome samples) had highly localised spatial patterns throughout different seasons and different study areas in the host squirrel population. The core gut microbiome, on the other hand, showed some spatial patterns, though fewer than in the non-core gut microbiome. Environmental transmission of microbiota is the likely contributor to the spatiotemporal distribution observed in the North American red squirrel gut microbiome.

16S rRNA survey revealed complex bacterial communities and evidence of bacterial interference on human adenoids.

Adenoid microbiota plays an important role in the development of various infectious and non-infectious diseases of the upper airways, such as otitis media, adenotonsillitis, rhinosinusitis and adenoid hypertrophy. Studies have suggested that adenoids could act as a potential reservoir of opportunistic pathogens. However, previous bacterial surveys of adenoids were mainly culture based and therefore might only provide an incomplete and potentially biased assessment of the microbial diversity. To develop an in-depth and comprehensive understanding of the adenoid microbial communities and test the 'pathogen reservoir hypothesis', we carried out a 16S rRNA based, culture-independent survey of bacterial communities on 67 human adenoids removed by surgery. Our survey revealed highly diverse adenoid bacterial communities distinct from those of other body habitats. Despite large interpersonal variations, adenoid microbiota shared a core set of taxa and can be classified into at least five major types based on its bacterial species composition. Our results support the 'pathogen reservoir hypothesis' as we found common pathogens of otitis media to be both prevalent and abundant. Co-occurrence analyses revealed evidence consistent with the bacterial interference theory in that multiple common pathogens showed 'non-coexistence' relationships with non-pathogenic members of the commensal microflora.

Phylogenomic analysis of bacterial and archaeal sequences with AMPHORA2.

SUMMARY: With the explosive growth of bacterial and archaeal sequence data, large-scale phylogenetic analyses present both opportunities and challenges. Here we describe AMPHORA2, an automated phylogenomic inference tool that can be used for high-throughput, high-quality genome tree reconstruction and metagenomic phylotyping. Compared with its predecessor, AMPHORA2 has several major enhancements and new functions: it has a greatly expanded phylogenetic marker database and can analyze both bacterial and archaeal sequences; it incorporates probability-based sequence alignment masks that improve the phylogenetic accuracy; it can analyze DNA as well as protein sequences and is more sensitive in marker identification; finally, it is over 100× faster in metagenomic phylotyping. AVAILABILITY: http://wolbachia.biology.virginia.edu/WuLab/Software.html. CONTACT: mw4yv@virginia.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Surface-exposed proteins of pathogenic mycobacteria and the role of cu-zn superoxide dismutase in macrophages and neutrophil survival.

Pathogenic mycobacteria are important agents causing human disease. Mycobacterium avium subsp. hominissuis (M. avium) is a species of recalcitrant environmental pathogen. The bacterium forms robust biofilms that allow it to colonize and persist in austere environments, such as residential and commercial water systems. M. avium is also an opportunistic pathogen that is a significant source of mortality for immune-compromised individuals. Proteins exposed at the bacterial surface play a central role in mediating the relationship between the bacterium and its environment. The processes underlying both biofilm formation and pathogenesis are directly dependent on this essential subset of the bacterial proteome. Therefore, the characterization of the surface-exposed proteome is an important step towards an improved understanding of the mycobacterial biology and pathogenesis. Here we examined the complement of surface exposed proteins from Mycobacterium avium 104, a clinical isolate and reference strain of Mycobacterium avium subsp. hominissuis. To profile the surface-exposed proteins of viable M. avium 104, bacteria were covalently labeled with a membrane impermeable biotinylation reagent and labeled proteins were affinity purified via the biotin-streptavidin interaction. The results provide a helpful snapshot of the surface-exposed proteome of this frequently utilized reference strain of M. avium. A Cu-Zn SOD knockout mutant, MAV_2043, a surface identified protein, was evaluated regarding its role in the survival in both macrophages and neutrophils.

Incorporating 16S gene copy number information improves estimates of microbial diversity and abundance.

The abundance of different SSU rRNA ("16S") gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly - from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.