Cpx-signalling facilitates Hms-dependent biofilm formation by Yersinia pseudotuberculosis

Bacteria often reside in sessile communities called biofilms, where they adhere to a variety of surfaces and exist as aggregates in a viscous polymeric matrix. Biofilms are resistant to antimicrobial treatments, and are a major contributor to the persistence and chronicity of many bacterial infections. Herein, we determined that the CpxA-CpxR two-component system influenced the ability of enteropathogenic Yersinia pseudotuberculosis to develop biofilms. Mutant bacteria that accumulated the active CpxR~P isoform failed to form biofilms on plastic or on the surface of the Caenorhabditis elegans nematode. A failure to form biofilms on the worm surface prompted their survival when grown on the lawns of Y. pseudotuberculosis. Exopolysaccharide production by the hms loci is the major driver of biofilms formed by Yersinia. We used a number of molecular genetic approaches to demonstrate that active CpxR~P binds directly to the promoter regulatory elements of the hms loci to activate the repressors of hms expression and to repress the activators of hms expression. Consequently, active Cpx-signalling culminated in a loss of exopolysaccharide production. Hence, the development of Y. pseudotuberculosis biofilms on multiple surfaces is controlled by the Cpx-signalling, and at least in part this occurs through repressive effects on the Hms-dependent exopolysaccharide production.

Whole-genome sequencing analysis of uncommon Shiga toxin-producing Escherichia coli from cattle: Virulence gene profiles, antimicrobial resistance predictions, and identification of novel O-serogroups

Shiga toxin-producing E. coli (STEC) are major foodborne pathogens. While many studies have focused on the "top-7 STEC", little is known for minor serogroups. A total of 284 non-top-7 STEC strains isolated from cattle feces were subjected to whole-genome sequencing (WGS) to determine the serotypes, the presence of virulence genes and antimicrobial resistance (AMR) determinants. Nineteen typeable and three non-typeable serotypes with novel O-antigen loci were identified. Twenty-one AMR genes and point mutations in another six genes that conferred resistance to 10 antimicrobial classes were detected, as well as 46 virulence genes. The distribution of 33 virulence genes and 15 AMR determinants exhibited significant differences among serotypes (p < 0.05). Among all strains, 81.7% (n = 232) and 14.1% (n = 40) carried stx2 and stx1 only, respectively; only 4.2% (n = 12) carried both. Subtypes stx1a, stx1c, stx2a, stx2c, stx2d, and stx2g were identified. Forty-six strains carried eae and stx2a and therefore had the potential cause severe diseases; 47 strains were genetically related to human clinical strains inferred from a pan-genome phylogenetic tree. We were able to demonstrate the utility of WGS as a surveillance tool to characterize the novel serotypes, as well as AMR and virulence profiles of uncommon STEC that could potentially cause human illness.

HPLC-UV, Metabarcoding and Genome Skims of Botanical Dietary Supplements: A Case Study in Echinacea

The use of DNA-based methods to authenticate botanical dietary supplements has been vigorously debated for a variety of reasons. More comparisons of DNA-based and chemical methods are needed, and concordant evaluation of orthogonal approaches on the same products will provide data to better understand the strengths and weaknesses of both approaches. The overall application of DNA-based methods is already firmly integrated into a wide array of continually modernizing stand alone and complementary authentication protocols. Recently, the use of full-length chloroplast genome sequences provided enhanced discriminatory capacity for closely related species of Echinacea compared to traditional DNA barcoding approaches (matK and rbcL). Here, two next-generation sequencing approaches were used: (1) genome skimming and (2) PCR amplicon (metabarcoding). The two genetic approaches were then combined with HPLC-UV to evaluate 20 commercially available dietary supplements of Echinacea representing "finished" products. The trade-offs involved in different DNA approaches were discussed, with a focus on how DNA methods support existing, accepted chemical methods. In most of the products (19/20), HPLC-UV suggested the presence of Echinacea spp. While metabarcoding was not useful with this genus and instead only resolved 7 products to the family level, genome skimming was able to resolve to species (9) or genus (1) with the 10/20 products where it was successful. Additional ingredients that HPLC-UV was unable to identify were also found in four products along with the relative sequence proportion of the constituents. Additionally, genome skimming was able to identify one product that was a different Echinacea species entirely.

Polishing the Oxford Nanopore long-read assemblies of bacterial pathogens with Illumina short reads to improve genomic analyses

Oxford Nanopore sequencing has been widely used to achieve complete genomes of bacterial pathogens. However, the error rates of Oxford Nanopore long reads are high. Various polishing algorithms using Illumina short reads to correct the errors in Oxford Nanopore long-read assemblies have been developed. The impact of polishing the Oxford Nanopore long-read assemblies of bacterial pathogens with Illumina short reads on improving genomic analyses was evaluated using both simulated and real reads. Ten species (10 strains) were selected for simulated reads, while real reads were tested on 11 species (11 strains). Oxford Nanopore long reads were assembled with Unicycler to produce a draft assembly, followed by three rounds of polishing with Illumina short reads using two polishing tools, Pilon and NextPolish. One round of NextPolish polishing generated genome completeness and accuracy parameters similar to the reference genomes, whereas two or three rounds of Pilon polishing were needed, though contiguity remained unchanged after polishing. The polished assemblies of Escherichia coli O157:H7, Salmonella Typhimurium, and Cronobacter sakazakii with simulated reads did not provide accurate plasmid identifications. One round of NextPolish polishing was needed for accurately identifying plasmids in Staphylococcus aureus and E. coli O26:H11 with real reads, whereas one and two rounds of Pilon polishing were necessary for these two strains, respectively. Polishing failed to provide an accurate antimicrobial resistance (AMR) genotype for S. aureus with real reads. One round of polishing recovered an accurate AMR genotype for Klebsiella pneumoniae with real reads. The reference genome and draft assembly of Citrobacter braakii with real reads differed, which carried blaCMY-83 and fosA6, respectively, while both genes were present after one round of polishing. However, polishing did not improve the assembly of E. coli O26:H11 with real reads to achieve numbers of virulence genes similar to the reference genome. The draft and polished assemblies showed a phylogenetic tree topology comparable with the reference genomes. For multilocus sequence typing and pan-genome analyses, one round of NextPolish polishing was sufficient to obtain accurate results, while two or three rounds of Pilon polishing were needed. Overall, NextPolish outperformed Pilon for polishing the Oxford Nanopore long-read assemblies of bacterial pathogens, though both polishing strategies improved genomic analyses compared to the draft assemblies.

Suitability of DNA Sequencing Tools for Identifying Edible Seaweeds Sold in the United States

Seaweeds have been consumed by billions of people around the world and are increasingly popular in United States (US) diets. Some seaweed species have been associated with adverse health effects-such as heavy metal toxicity-and higher priced seaweeds may be more prone to adulteration. Knowing which species of seaweeds are being marketed in the US is important for protecting human health and preventing economic adulteration. Therefore, the United States Food and Drug Administration is developing new DNA-based species identification tools to complement established chemical methods for verifying the accurate labeling of products. Here, seaweed products available in the United States were surveyed using a tiered approach to evaluate a variety of DNA extraction techniques followed by traditional DNA barcoding via Sanger sequencing; if needed, genome skimming of total extracted nuclear DNA via next-generation sequencing was performed. This two-tiered approach of DNA barcoding and genome skimming could identify most seaweed samples (41/46), even those in blends (2/2, 1 out of 3 labeled species in each). Only two commercial samples appeared to be mislabeled or to contain unintended algal species. Five samples, labeled as "hijiki" or "arame", could not be confirmed by these DNA-based identification methods.

Benchmarking Long-Read Assemblers for Genomic Analyses of Bacterial Pathogens Using Oxford Nanopore Sequencing

Oxford Nanopore sequencing can be used to achieve complete bacterial genomes. However, the error rates of Oxford Nanopore long reads are greater compared to Illumina short reads. Long-read assemblers using a variety of assembly algorithms have been developed to overcome this deficiency, which have not been benchmarked for genomic analyses of bacterial pathogens using Oxford Nanopore long reads. In this study, long-read assemblers, namely Canu, Flye, Miniasm/Racon, Raven, Redbean, and Shasta, were thus benchmarked using Oxford Nanopore long reads of bacterial pathogens. Ten species were tested for mediocre- and low-quality simulated reads, and 10 species were tested for real reads. Raven was the most robust assembler, obtaining complete and accurate genomes. All Miniasm/Racon and Raven assemblies of mediocre-quality reads provided accurate antimicrobial resistance (AMR) profiles, while the Raven assembly of Klebsiella variicola with low-quality reads was the only assembly with an accurate AMR profile among all assemblers and species. All assemblers functioned well for predicting virulence genes using mediocre-quality and real reads, whereas only the Raven assemblies of low-quality reads had accurate numbers of virulence genes. Regarding multilocus sequence typing (MLST), Miniasm/Racon was the most effective assembler for mediocre-quality reads, while only the Raven assemblies of Escherichia coli O157:H7 and K. variicola with low-quality reads showed positive MLST results. Miniasm/Racon and Raven were the best performers for MLST using real reads. The Miniasm/Racon and Raven assemblies showed accurate phylogenetic inference. For the pan-genome analyses, Raven was the strongest assembler for simulated reads, whereas Miniasm/Racon and Raven performed the best for real reads. Overall, the most robust and accurate assembler was Raven, closely followed by Miniasm/Racon.

Benchmarking hybrid assembly approaches for genomic analyses of bacterial pathogens using Illumina and Oxford Nanopore sequencing

Background

We benchmarked the hybrid assembly approaches of MaSuRCA, SPAdes, and Unicycler for bacterial pathogens using Illumina and Oxford Nanopore sequencing by determining genome completeness and accuracy, antimicrobial resistance (AMR), virulence potential, multilocus sequence typing (MLST), phylogeny, and pan genome. Ten bacterial species (10 strains) were tested for simulated reads of both mediocre- and low-quality, whereas 11 bacterial species (12 strains) were tested for real reads.

Results

Unicycler performed the best for achieving contiguous genomes, closely followed by MaSuRCA, while all SPAdes assemblies were incomplete. MaSuRCA was less tolerant of low-quality long reads than SPAdes and Unicycler. The hybrid assemblies of five antimicrobial-resistant strains with simulated reads provided consistent AMR genotypes with the reference genomes. The MaSuRCA assembly of Staphylococcus aureus with real reads contained msr(A) and tet(K), while the reference genome and SPAdes and Unicycler assemblies harbored blaZ. The AMR genotypes of the reference genomes and hybrid assemblies were consistent for the other five antimicrobial-resistant strains with real reads. The numbers of virulence genes in all hybrid assemblies were similar to those of the reference genomes, irrespective of simulated or real reads. Only one exception existed that the reference genome and hybrid assemblies of Pseudomonas aeruginosa with mediocre-quality long reads carried 241 virulence genes, whereas 184 virulence genes were identified in the hybrid assemblies of low-quality long reads. The MaSuRCA assemblies of Escherichia coli O157:H7 and Salmonella Typhimurium with mediocre-quality long reads contained 126 and 118 virulence genes, respectively, while 110 and 107 virulence genes were detected in their MaSuRCA assemblies of low-quality long reads, respectively. All approaches performed well in our MLST and phylogenetic analyses. The pan genomes of the hybrid assemblies of S. Typhimurium with mediocre-quality long reads were similar to that of the reference genome, while SPAdes and Unicycler were more tolerant of low-quality long reads than MaSuRCA for the pan-genome analysis. All approaches functioned well in the pan-genome analysis of Campylobacter jejuni with real reads.

Conclusions

Our research demonstrates the hybrid assembly pipeline of Unicycler as a superior approach for genomic analyses of bacterial pathogens using Illumina and Oxford Nanopore sequencing.

Genomic analyses of multidrug-resistant Salmonella Indiana, Typhimurium, and Enteritidis isolates using MinION and MiSeq sequencing technologies

We sequenced 25 isolates of phenotypically multidrug-resistant Salmonella Indiana (n = 11), Typhimurium (n = 8), and Enteritidis (n = 6) using both MinION long-read [SQK-LSK109 and flow cell (R9.4.1)] and MiSeq short-read (Nextera XT and MiSeq Reagent Kit v2) sequencing technologies to determine the advantages of each approach in terms of the characteristics of genome structure, antimicrobial resistance (AMR), virulence potential, whole-genome phylogeny, and pan-genome. The MinION reads were base-called in real-time using MinKnow 3.4.8 integrated with Guppy 3.0.7. The long-read-only assembly, Illumina-only assembly, and hybrid assembly pipelines of Unicycler 0.4.8 were used to generate the MinION, MiSeq, and hybrid assemblies, respectively. The MinION assemblies were highly contiguous compared to the MiSeq assemblies but lacked accuracy, a deficiency that was mitigated by adding the MiSeq short reads through the Unicycler hybrid assembly which corrected erroneous single nucleotide polymorphisms (SNPs). The MinION assemblies provided similar predictions of AMR and virulence potential compared to the MiSeq and hybrid assemblies, although they produced more total false negatives of AMR genotypes, primarily due to failure in identifying tetracycline resistance genes in 11 of the 19 MinION assemblies of tetracycline-resistant isolates. The MinION assemblies displayed a large genetic distance from their corresponding MiSeq and hybrid assemblies on the whole-genome phylogenetic tree, indicating that the lower read accuracy of MinION sequencing caused incorrect clustering. The pan-genome of the MinION assemblies contained significantly more accessory genes and less core genes compared to the MiSeq and hybrid assemblies, suggesting that although these assemblies were more contiguous, their sequencing errors reduced accurate genome annotations. Our research demonstrates that MinION sequencing by itself provides an efficient assessment of the genome structure, antimicrobial resistance, and virulence potential of Salmonella; however, it is not sufficient for whole-genome phylogenetic and pan-genome analyses. MinION in combination with MiSeq facilitated the most accurate genomic analyses.

Implementation and dissemination of an African American popular opinion model to improve lupus awareness: an academic-community partnership

OBJECTIVE

Lupus is a chronic, autoimmune disease that disproportionately affects African Americans. We adapted the Centers for Disease Control and Prevention's Popular Opinion Leader model to implement an intervention tailored for African American individuals that leverages an academic-community partnership and community-based social networks to disseminate culturally appropriate lupus education.

METHODS

Academic rheumatologists, social scientists, and researchers in Boston, MA and Chicago, IL partnered with local lupus support groups, community organizations, and churches in neighborhoods with higher proportions of African Americans to develop curriculum and recruit community leaders with and without lupus (Popular Opinion Leaders; POLs). POLs attended four training sessions focused on lupus education, strategies to educate others, and a review of research methods. POLs disseminated information through their social networks and recorded their impact, which was mapped using a geographic information system framework.

RESULTS

We trained 18 POLs in greater Boston and 19 in greater Chicago: 97% were African American, 97% were female; and the mean age was 57 years. Fifty-nine percent of Boston POLs and 68% of Chicago POLs had lupus. POLs at both sites engaged members of their social networks and communities in conversations about lupus, health disparities, and the importance of care. Boston POLs documented 97 encounters with 547 community members reached. Chicago POLs documented 124 encounters with 4083 community members reached.

CONCLUSIONS

An adapted, community-based POL model can be used to disseminate lupus education and increase awareness in African American communities. Further research is needed to determine the degree to which this may begin to reduce disparities in access to care and outcomes.

Alignment-free genome comparison enables accurate geographic sourcing of white oak DNA

Background

The application of genomic data and bioinformatics for the identification of restricted or illegally-sourced natural products is urgently needed. The taxonomic identity and geographic provenance of raw and processed materials have implications in sustainable-use commercial practices, and relevance to the enforcement of laws that regulate or restrict illegally harvested materials, such as timber. Improvements in genomics make it possible to capture and sequence partial-to-complete genomes from challenging tissues, such as wood and wood products.

Results

In this paper, we report the success of an alignment-free genome comparison method, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {d}_2^{\ast }, $$\end{document}d2∗, that differentiates different geographic sources of white oak (Quercus) species with a high level of accuracy with very small amount of genomic data. The method is robust to sequencing errors, different sequencing laboratories and sequencing platforms.

Conclusions

This method offers an approach based on genome-scale data, rather than panels of pre-selected markers for specific taxa. The method provides a generalizable platform for the identification and sourcing of materials using a unified next generation sequencing and analysis framework.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5253-1) contains supplementary material, which is available to authorized users.

Yersinia pseudotuberculosis BarA-UvrY Two-Component Regulatory System Represses Biofilms via CsrB

The formation of biofilms by Yersinia pseudotuberculosis (Yptb) and Y. pestis requires the hmsHFRS genes, which direct production of a polysaccharide extracellular matrix (Hms-ECM). Despite possessing identical hmsHFRS sequences, Yptb produces much less Hms-ECM than Y. pestis. The regulatory influences that control Yptb Hms-ECM production and biofilm formation are not fully understood. In this study, negative regulators of biofilm production in Yptb were identified. Inactivation of the BarA/UvrY two-component system or the CsrB regulatory RNA increased binding of Congo Red dye, which correlates with extracellular polysaccharide production. These mutants also produced biofilms that were substantially more cohesive than the wild type strain. Disruption of uvrY was not sufficient for Yptb to cause proventricular blockage during infection of Xenopsylla cheopis fleas. However, this strain was less acutely toxic toward fleas than wild type Yptb. Flow cytometry measurements of lectin binding indicated that Yptb BarA/UvrY/CsrB mutants may produce higher levels of other carbohydrates in addition to poly-GlcNAc Hms-ECM. In an effort to characterize the relevant downstream targets of the BarA/UvrY system, we conducted a proteomic analysis to identify proteins with lower abundance in the csrB::Tn5 mutant strain. Urease subunit proteins were less abundant and urease enzymatic activity was lower, which likely reduced toxicity toward fleas. Loss of CsrB impacted expression of several potential regulatory proteins that may influence biofilms, including the RcsB regulator. Overexpression of CsrB did not alter the Congo-red binding phenotype of an rcsB::Tn5 mutant, suggesting that the effect of CsrB on biofilms may require RcsB. These results underscore the regulatory and compositional differences between Yptb and Y. pestis biofilms. By activating CsrB expression, the Yptb BarA/UvrY two-component system has pleiotropic effects that impact biofilm production and stability.

Reconstructing a herbivore's diet using a novel rbcL DNA mini-barcode for plants

Next Generation Sequencing and the application of metagenomic analyses can be used to answer questions about animal diet choice and study the consequences of selective foraging by herbivores. The quantification of herbivore diet choice with respect to native versus exotic plant species is particularly relevant given concerns of invasive species establishment and their effects on ecosystems. While increased abundance of white-tailed deer (Odocoileus virginianus) appears to correlate with increased incidence of invasive plant species, data supporting a causal link is scarce. We used a metabarcoding approach (PCR amplicons of the plant rbcL gene) to survey the diet of white-tailed deer (fecal samples), from a forested site in Warren County, Virginia with a comprehensive plant species inventory and corresponding reference collection of plant barcode and chloroplast sequences. We sampled fecal pellet piles and extracted DNA from 12 individual deer in October 2014. These samples were compared to a reference DNA library of plant species collected within the study area. For 72 % of the amplicons, we were able to assign taxonomy at the species level, which provides for the first time-sufficient taxonomic resolution to quantify the relative frequency at which native and exotic plant species are being consumed by white-tailed deer. For each of the 12 individual deer we collected three subsamples from the same fecal sample, resulting in sequencing 36 total samples. Using Qiime, we quantified the plant DNA found in all 36 samples, and found that variance within samples was less than variance between samples (F = 1.73, P = 0.004), indicating additional subsamples may not be necessary. Species level diversity ranged from 60 to 93 OTUs per individual and nearly 70 % of all plant sequences recovered were from native plant species. The number of species detected did reduce significantly (range 4-12) when we excluded species whose OTU composed <1 % of each sample's total. When compared to the abundance of native and non-natives plants inventoried in the local community, our results support the observation that white-tailed deer have strong foraging preferences, but these preferences were not consistent for species in either class. Deer forage behaviour may favour some exotic species, but not all.

An analysis of Echinacea chloroplast genomes: Implications for future botanical identification

Echinacea is a common botanical used in dietary supplements, primarily to treat upper respiratory tract infections and to support immune function. There are currently thought to be nine species in the genus Echinacea. Due to very low molecular divergence among sister species, traditional DNA barcoding has not been successful for differentiation of Echinacea species. Here, we present the use of full chloroplast genomes to distinguish between all 9 reported species. Total DNA was extracted from specimens stored at the National Museum of Natural History, Smithsonian Institution, which had been collected from the wild with species identification documented by experts in the field. We used Next Generation Sequencing (NGS) and CLC Genomics Workbench to assemble complete chloroplast genomes for all nine species. Full chloroplasts unambiguously differentiated all nine species, compared with the very few single nucleotide polymorphisms (SNPs) available with core DNA barcoding markers. SNPs for any two Echinacea chloroplast genomes ranged from 181 to 910, and provided robust data for unambiguous species delimitation. Implications for DNA-based species identification assays derived from chloroplast genome sequences are discussed in light of product safety, adulteration and quality issues.

Landscape Change in an Agricultural Watershed: The Effect of Parcelization on Riparian Forest Cover

In this research we address factors contributing to landscape change in a midwestern (USA) watershed. Specifically, the relationship between the parcelization of agricultural land (ownership subdivisions) and changes in amounts of riparian forest cover are explored. The study area is the River Raisin Watershed in southeastern Michigan, which is typical of rural watersheds in the lower Great Lakes region. Two townships within the watershed were sampled from data spanning a 20-year time period. Using land-use and land-cover data for 1968, 1978, and 1988, combined with ownership records for the same years, we determined a relationship between rates of parcelization and changes in forest cover. These findings indicate that increasing trends in riparian forest area follow increases in land parcelization at the township scale. For these two townships, agricultural decline and land-use diversification seem positively related to larger areas and widths of riparian forests. As the number of acres being farmed and the number of farms have decreased, the land has been divided into more and smaller lots and the forests along rivers like the Raisin and its tributaries have increased in area. There are several possible factors at the township and farm scale which may influence this landscape phenomenon, including land-use policies, governmental programs, agricultural mechanization, and agricultural economics.

Lipopolysaccharide Biosynthesis Genes of Yersinia pseudotuberculosis Promote Resistance to Antimicrobial Chemokines

Antimicrobial chemokines (AMCs) are a recently described family of host defense peptides that play an important role in protecting a wide variety of organisms from bacterial infection. Very little is known about the bacterial targets of AMCs or factors that influence bacterial susceptibility to AMCs. In an effort to understand how bacterial pathogens resist killing by AMCs, we screened Yersinia pseudotuberculosis transposon mutants for those with increased binding to the AMCs CCL28 and CCL25. Mutants exhibiting increased binding to AMCs were subjected to AMC killing assays, which revealed their increased sensitivity to chemokine-mediated cell death. The majority of the mutants exhibiting increased binding to AMCs contained transposon insertions in genes related to lipopolysaccharide biosynthesis. A particularly strong effect on susceptibility to AMC mediated killing was observed by disruption of the hldD/waaF/waaC operon, necessary for ADP-L-glycero-D-manno-heptose synthesis and a complete lipopolysaccharide core oligosaccharide. Periodate oxidation of surface carbohydrates also enhanced AMC binding, whereas enzymatic removal of surface proteins significantly reduced binding. These results suggest that the structure of Y. pseudotuberculosis LPS greatly affects the antimicrobial activity of AMCs by shielding a protein ligand on the bacterial cell surface.

Closely-related taxa influence woody species discrimination via DNA barcoding: evidence from global forest dynamics plots

To determine how well DNA barcodes from the chloroplast region perform in forest dynamics plots (FDPs) from global CTFS-ForestGEO network, we analyzed DNA barcoding sequences of 1277 plant species from a wide phylogenetic range (3 FDPs in tropics, 5 in subtropics and 5 in temperate zone) and compared the rates of species discrimination (RSD). We quantified RSD by two DNA barcode combinations (rbcL + matK and rbcL + matK + trnH-psbA) using a monophyly-based method (GARLI). We defined two indexes of closely-related taxa (Gm/Gt and S/G ratios) and correlated these ratios with RSD. The combination of rbcL + matK averagely discriminated 88.65%, 83.84% and 72.51% at the local, regional and global scales, respectively. An additional locus trnH-psbA increased RSD by 2.87%, 1.49% and 3.58% correspondingly. RSD varied along a latitudinal gradient and were negatively correlated with ratios of closely-related taxa. Successes of species discrimination generally depend on scales in global FDPs. We suggested that the combination of rbcL + matK + trnH-psbA is currently applicable for DNA barcoding-based phylogenetic studies on forest communities.

A well-resolved phylogeny of the trees of Puerto Rico based on DNA barcode sequence data

Background

The use of phylogenetic information in community ecology and conservation has grown in recent years. Two key issues for community phylogenetics studies, however, are (i) low terminal phylogenetic resolution and (ii) arbitrarily defined species pools.

Methodology/principal findings

We used three DNA barcodes (plastid DNA regions rbcL, matK, and trnH-psbA) to infer a phylogeny for 527 native and naturalized trees of Puerto Rico, representing the vast majority of the entire tree flora of the island (89%). We used a maximum likelihood (ML) approach with and without a constraint tree that enforced monophyly of recognized plant orders. Based on 50% consensus trees, the ML analyses improved phylogenetic resolution relative to a comparable phylogeny generated with Phylomatic (proportion of internal nodes resolved: constrained ML = 74%, unconstrained ML = 68%, Phylomatic = 52%). We quantified the phylogenetic composition of 15 protected forests in Puerto Rico using the constrained ML and Phylomatic phylogenies. We found some evidence that tree communities in areas of high water stress were relatively phylogenetically clustered. Reducing the scale at which the species pool was defined (from island to soil types) changed some of our results depending on which phylogeny (ML vs. Phylomatic) was used. Overall, the increased terminal resolution provided by the ML phylogeny revealed additional patterns that were not observed with a less-resolved phylogeny.

Conclusions/significance

With the DNA barcode phylogeny presented here (based on an island-wide species pool), we show that a more fully resolved phylogeny increases power to detect nonrandom patterns of community composition in several Puerto Rican tree communities. Especially if combined with additional information on species functional traits and geographic distributions, this phylogeny will (i) facilitate stronger inferences about the role of historical processes in governing the assembly and composition of Puerto Rican forests, (ii) provide insight into Caribbean biogeography, and (iii) aid in incorporating evolutionary history into conservation planning.

Comparative evolutionary diversity and phylogenetic structure across multiple forest dynamics plots: a mega-phylogeny approach

Forest dynamics plots, which now span longitudes, latitudes, and habitat types across the globe, offer unparalleled insights into the ecological and evolutionary processes that determine how species are assembled into communities. Understanding phylogenetic relationships among species in a community has become an important component of assessing assembly processes. However, the application of evolutionary information to questions in community ecology has been limited in large part by the lack of accurate estimates of phylogenetic relationships among individual species found within communities, and is particularly limiting in comparisons between communities. Therefore, streamlining and maximizing the information content of these community phylogenies is a priority. To test the viability and advantage of a multi-community phylogeny, we constructed a multi-plot mega-phylogeny of 1347 species of trees across 15 forest dynamics plots in the ForestGEO network using DNA barcode sequence data (rbcL, matK, and psbA-trnH) and compared community phylogenies for each individual plot with respect to support for topology and branch lengths, which affect evolutionary inference of community processes. The levels of taxonomic differentiation across the phylogeny were examined by quantifying the frequency of resolved nodes throughout. In addition, three phylogenetic distance (PD) metrics that are commonly used to infer assembly processes were estimated for each plot [PD, Mean Phylogenetic Distance (MPD), and Mean Nearest Taxon Distance (MNTD)]. Lastly, we examine the partitioning of phylogenetic diversity among community plots through quantification of inter-community MPD and MNTD. Overall, evolutionary relationships were highly resolved across the DNA barcode-based mega-phylogeny, and phylogenetic resolution for each community plot was improved when estimated within the context of the mega-phylogeny. Likewise, when compared with phylogenies for individual plots, estimates of phylogenetic diversity in the mega-phylogeny were more consistent, thereby removing a potential source of bias at the plot-level, and demonstrating the value of assessing phylogenetic relationships simultaneously within a mega-phylogeny. An unexpected result of the comparisons among plots based on the mega-phylogeny was that the communities in the ForestGEO plots in general appear to be assemblages of more closely related species than expected by chance, and that differentiation among communities is very low, suggesting deep floristic connections among communities and new avenues for future analyses in community ecology.

Wood nitrogen concentrations in tropical trees: phylogenetic patterns and ecological correlates

In tropical and temperate trees, wood chemical traits are hypothesized to covary with species' life-history strategy along a 'wood economics spectrum' (WES), but evidence supporting these expected patterns remains scarce. Due to its role in nutrient storage, we hypothesize that wood nitrogen (N) concentration will covary along the WES, being higher in slow-growing species with high wood density (WD), and lower in fast-growing species with low WD. In order to test this hypothesis we quantified wood N concentrations in 59 Panamanian hardwood species, and used this dataset to examine ecological correlates and phylogenetic patterns of wood N. Wood N varied > 14-fold among species between 0.04 and 0.59%; closely related species were more similar in wood N than expected by chance. Wood N was positively correlated with WD, and negatively correlated with log-transformed relative growth rates, although these relationships were relatively weak. We found evidence for co-evolution between wood N and both WD and log-transformed mortality rates. Our study provides evidence that wood N covaries with tree life-history parameters, and that these patterns consistently co-evolve in tropical hardwoods. These results provide some support for the hypothesized WES, and suggest that wood is an increasingly important N pool through tropical forest succession.

The coding region of the UFGT gene is a source of diagnostic SNP markers that allow single-locus DNA genotyping for the assessment of cultivar identity and ancestry in grapevine (Vitis vinifera L.)

BACKGROUND

Vitis vinifera L. is one of society's most important agricultural crops with a broad genetic variability. The difficulty in recognizing grapevine genotypes based on ampelographic traits and secondary metabolites prompted the development of molecular markers suitable for achieving variety genetic identification.

FINDINGS

Here, we propose a comparison between a multi-locus barcoding approach based on six chloroplast markers and a single-copy nuclear gene sequencing method using five coding regions combined with a character-based system with the aim of reconstructing cultivar-specific haplotypes and genotypes to be exploited for the molecular characterization of 157 V. vinifera accessions. The analysis of the chloroplast target regions proved the inadequacy of the DNA barcoding approach at the subspecies level, and hence further DNA genotyping analyses were targeted on the sequences of five nuclear single-copy genes amplified across all of the accessions. The sequencing of the coding region of the UFGT nuclear gene (UDP-glucose: flavonoid 3-0-glucosyltransferase, the key enzyme for the accumulation of anthocyanins in berry skins) enabled the discovery of discriminant SNPs (1/34 bp) and the reconstruction of 130 V. vinifera distinct genotypes. Most of the genotypes proved to be cultivar-specific, and only few genotypes were shared by more, although strictly related, cultivars.

CONCLUSION

On the whole, this technique was successful for inferring SNP-based genotypes of grapevine accessions suitable for assessing the genetic identity and ancestry of international cultivars and also useful for corroborating some hypotheses regarding the origin of local varieties, suggesting several issues of misidentification (synonymy/homonymy).

Tropical plant-herbivore networks: reconstructing species interactions using DNA barcodes

Plants and their associated insect herbivores, represent more than 50% of all known species on earth. The first step in understanding the mechanisms generating and maintaining this important component of biodiversity is to identify plant-herbivore associations. In this study we determined insect-host plant associations for an entire guild of insect herbivores using plant DNA extracted from insect gut contents. Over two years, in a tropical rain forest in Costa Rica (La Selva Biological Station), we recorded the full diet breadth of rolled-leaf beetles, a group of herbivores that feed on plants in the order Zingiberales. Field observations were used to determine the accuracy of diet identifications using a three-locus DNA barcode (rbcL, trnH-psbA and ITS2). Using extraction techniques for ancient DNA, we obtained high-quality sequences for two of these loci from gut contents (rbcL and ITS2). Sequences were then compared to a comprehensive DNA barcode library of the Zingiberales. The rbcL locus identified host plants to family (success/sequence = 58.8%) and genus (success/sequence = 47%). For all Zingiberales except Heliconiaceae, ITS2 successfully identified host plants to genus (success/sequence = 67.1%) and species (success/sequence = 61.6%). Kindt’s sampling estimates suggest that by collecting ca. four individuals representing each plant-herbivore interaction, 99% of all host associations included in this study can be identified to genus. For plants that amplified ITS2, 99% of the hosts can be identified to species after collecting at least four individuals representing each interaction. Our study demonstrates that host plant identifications at the species-level using DNA barcodes are feasible, cost-effective, and reliable, and that reconstructing plant-herbivore networks with these methods will become the standard for a detailed understanding of these interactions.

Temporal turnover in the composition of tropical tree communities: functional determinism and phylogenetic stochasticity

The degree to which turnover in biological communities is structured by deterministic or stochastic factors and the identities of influential deterministic factors are fundamental, yet unresolved, questions in ecology. Answers to these questions are particularly important for projecting the fate of forests with diverse disturbance histories worldwide. To uncover the processes governing turnover we use species-level molecular phylogenies and functional trait data sets for two long-term tropical forest plots with contrasting disturbance histories: one forest is older-growth, and one was recently disturbed. Having both phylogenetic and functional information further allows us to parse out the deterministic influences of different ecological filters. With the use of null models we find that compositional turnover was random with respect to phylogeny on average, but highly nonrandom with respect to measured functional traits. Furthermore, as predicted by a deterministic assembly process, the older-growth and disturbed forests were characterized by less than and greater than expected functional turnover, respectively. These results suggest that the abiotic environment, which changes due to succession in the disturbed forest, strongly governs the temporal dynamics of disturbed and undisturbed tropical forests. Predicting future changes in the composition of disturbed and undisturbed forests may therefore be tractable when using a functional-trait-based approach.

Gene expression analysis of Xenopsylla cheopis (Siphonaptera: Pulicidae) suggests a role for reactive oxygen species in response to Yersinia pestis infection

Fleas are vectors for a number of pathogens including Yersinia pestis, yet factors that govern interactions between fleas and Y. pestis are not well understood. Examining gene expression changes in infected fleas could reveal pathways that affect Y. pestis survival in fleas and subsequent transmission. We used suppression subtractive hybridization to identify genes that are induced in Xenopsylla cheopis (Rothschild) (Siphonaptera: Pulicidae) in response to oral or hemocoel infection with Y. pestis. Overall, the transcriptional changes we detected were very limited. We identified several genes that are likely involved in the production or removal of reactive oxygen species (ROS). Midgut ROS levels were higher in infected fleas and antioxidant treatment before infection reduced ROS levels and resulted in higher bacterial loads. An ROS-sensitive mutant strain of Y. pestis lacking the OxyR transcriptional regulator showed reduced growth early after infection. Our results indicate that ROS may limit Y. pestis early colonization of fleas and that bacterial strategies to overcome ROS may enhance transmission.

Development and characterization of microsatellite loci for two Caribbean Heliconia (Heliconiaceae: H. bihai and H. caribaea)

PREMISE OF THE STUDY

Microsatellite loci were developed to characterize genetic variation and population subdivision in Heliconia bihai and H. caribaea from the Caribbean Islands.

METHODS AND RESULTS

A total of 13 new microsatellite markers were developed and characterized in the two Caribbean heliconias. Di-, tri-, and tetranucleotide repeats were identified with one to 17 alleles per locus, and the observed heterozygosity ranged from 0.13 to 0.87. Additionally, cross-species amplification was successful in eight out of 13 loci.

CONCLUSIONS

The microsatellite loci developed have discriminatory potential to be used in genetic characterizations of Caribbean Heliconia. Both H. bihai and H. caribaea are known to have adaptive interactions with their hummingbird pollinators, and the characterized microsatellite markers will be used to study mating system, genetic structure, and phylogeographic patterns in Caribbean Heliconia.

Generating plant DNA barcodes for trees in long-term forest dynamics plots

Long-term forest dynamics plots, such as those maintained and coordinated by the Center for Tropical Forest Science and Smithsonian Institution Global Earth Observatories (CTFS/SIGEO), are a rich source of biological data that describe the demographics, ecology, and evolution of pristine and disturbed forest habitats across ecosystems. As molecular techniques for plant systematic and ecological studies, including DNA barcodes, have improved so have the methods for collecting tissue samples, generating DNA sequences, and managing genetic data. Tissue samples can be processed at the point of collection and stored in silica gel for extended periods of time or samples can be taken from historical museum collections with sufficient DNA yields for study. In this chapter, we provide a workflow that includes the tracking of data from field collection of tissue samples to the DNA barcode sequence laboratory to final analyses for forensic and phylogenetic investigations.

DNA barcodes: methods and protocols

DNA barcoding, a new method for the quick identification of any species based on extracting a DNA sequence from a tiny tissue sample of any organism, is now being applied to taxa across the tree of life. As a research tool for taxonomists, DNA barcoding assists in identification by expanding the ability to diagnose species by including all life history stages of an organism. As a biodiversity discovery tool, DNA barcoding helps to flag species that are potentially new to science. As a biological tool, DNA barcoding is being used to address fundamental ecological and evolutionary questions, such as how species in plant communities are assembled. The process of DNA barcoding entails two basic steps: (1) building the DNA barcode library of known species and (2) matching the barcode sequence of the unknown sample against the barcode library for identification. Although DNA barcoding as a methodology has been in use for less than a decade, it has grown exponentially in terms of the number of sequences generated as barcodes as well as its applications. This volume provides the latest information on generating, applying, and analyzing DNA barcodes across the Tree of Life from animals and fungi to protists, algae, and plants.

Future directions

It is a risky task to attempt to predict the direction that DNA barcoding and its applications may take in the future. In a very short time, the endeavor of DNA barcoding has gone from being a tool to facilitate taxonomy in difficult to identify species, to an ambitious, global initiative that seeks to tackle such pertinent and challenging issues as quantifying global biodiversity, revolutionizing the forensic identifications of species, advancing the study of interactions among species, and promoting the reconstruction of evolutionary relationships within communities. The core element of DNA barcoding will always remain the same: the generation of a set of well-identified samples collected and genotyped at one or more genetic barcode markers and assembled into a properly curated database. But the application of this body of data will depend on the creativity and need of the research community in using a "gold standard" of annotated DNA sequence data at the species level. We foresee several areas where the application of DNA barcode data is likely to yield important evolutionary, ecological, and societal insights, and while far from exclusive, provide examples of how DNA barcode data will continue to empower scientists to address hypothesis-driven research. Three areas of immediate and obvious concern are (1) biodiversity inventories, (2) phylogenetic applications, and (3) species interactions.

Construction and analysis of phylogenetic trees using DNA barcode data

The assembly of sequence data obtained from DNA barcodes into phylogenies or NJ trees has proven highly useful in estimating relatedness among species as well as providing a framework in which hypotheses regarding the evolution of traits or species distributions may be investigated. In this chapter, we outline the process by which DNA sequence data is assembled into a phylogenetically informative matrix, and then provide details on the methods to reconstruct NJ or phylogenetic trees that employ DNA barcode data, using only barcode data or combining barcodes with other data.

The roots of diversity: below ground species richness and rooting distributions in a tropical forest revealed by DNA barcodes and inverse modeling

BACKGROUND

Plants interact with each other, nutrients, and microbial communities in soils through extensive root networks. Understanding these below ground interactions has been difficult in natural systems, particularly those with high plant species diversity where morphological identification of fine roots is difficult. We combine DNA-based root identification with a DNA barcode database and above ground stem locations in a floristically diverse lowland tropical wet forest on Barro Colorado Island, Panama, where all trees and lianas >1 cm diameter have been mapped to investigate richness patterns below ground and model rooting distributions.

METHODOLOGY/PRINCIPAL FINDINGS

DNA barcode loci, particularly the cpDNA locus trnH-psba, can be used to identify fine and small coarse roots to species. We recovered 33 species of roots from 117 fragments sequenced from 12 soil cores. Despite limited sampling, we recovered a high proportion of the known species in the focal hectare, representing approximately 14% of the measured woody plant richness. This high value is emphasized by the fact that we would need to sample on average 13 m(2) at the seedling layer and 45 m(2) for woody plants >1 cm diameter to obtain the same number of species above ground. Results from inverse models parameterized with the locations and sizes of adults and the species identifications of roots and sampling locations indicates a high potential for distal underground interactions among plants.

CONCLUSIONS

DNA barcoding techniques coupled with modeling approaches should be broadly applicable to studying root distributions in any mapped vegetation plot. We discuss the implications of our results and outline how second-generation sequencing technology and environmental sampling can be combined to increase our understanding of how root distributions influence the potential for plant interactions in natural ecosystems.

Biodiversity studies in Phaseolus species by DNA barcoding

The potential of DNA barcoding was tested as a system for studying genetic diversity and genetic traceability in bean germplasm. This technique was applied to several pure lines of Phaseolus vulgaris L. belonging to wild, domesticated, and cultivated common beans, along with some accessions of Phaseolus coccineus L., Phaseolus lunatus L., and Vigna unguiculata (L.) Walp. A multilocus approach was exploited using three chloroplast genic regions (rbcL, trnL, and matK), four intergenic spacers (rpoB-trnC, atpBrbcL, trnT-trnL, and psbA-trnH), and nuclear ITS1 and ITS2 rDNA sequences. Our main goals were to identify the markers and SNPs that show the best discriminant power at the variety level in common bean germplasm, to examine two methods (tree based versus character based) for biodiversity analysis and traceability assays, and to evaluate the overall utility of chloroplast DNA barcodes for reconstructing the origins of modern Italian varieties. Our results indicate that the neighbor-joining method is a powerful approach for comparing genetic diversity within plant species, but it is relatively uninformative for the genetic traceability of plant varieties. In contrast, the character-based method was able to identify several distinct haplotypes over all target regions corresponding to Mesoamerican or Andean accessions; Italian accessions originated from both gene pools. On the whole, our findings raise some concerns about the use of DNA barcoding for intraspecific genetic diversity studies in common beans and highlights its limitations for resolving genetic relationships between landraces and varieties.

Exploring tree-habitat associations in a Chinese subtropical forest plot using a molecular phylogeny generated from DNA barcode loci

Elucidating the ecological mechanisms underlying community assembly in subtropical forests remains a central challenge for ecologists. The assembly of species into communities can be due to interspecific differences in habitat associations, and there is increasing evidence that these associations may have an underlying phylogenetic structure in contemporary terrestrial communities. In other words, by examining the degree to which closely related species prefer similar habitats and the degree to which they co-occur, ecologists are able to infer the mechanisms underlying community assembly. Here we implement this approach in a diverse subtropical tree community in China using a long-term forest dynamics plot and a molecular phylogeny generated from three DNA barcode loci. We find that there is phylogenetic signal in plant-habitat associations (i.e. closely related species tend to prefer similar habitats) and that patterns of co-occurrence within habitats are typically non-random with respect to phylogeny. In particular, we found phylogenetic clustering in valley and low-slope habitats in this forest, indicating a filtering of lineages plays a dominant role in structuring communities in these habitats and we found evidence of phylogenetic overdispersion in high-slope, ridge-top and high-gully habitats, indicating that distantly related species tended to co-occur in these high elevation habitats and that lineage filtering is less important in structuring these communities. Thus we infer that non-neutral niche-based processes acting upon evolutionarily conserved habitat preferences explain the assembly of local scale communities in the forest studied.

Advances in the use of DNA barcodes to build a community phylogeny for tropical trees in a Puerto Rican forest dynamics plot

BACKGROUND

Species number, functional traits, and phylogenetic history all contribute to characterizing the biological diversity in plant communities. The phylogenetic component of diversity has been particularly difficult to quantify in species-rich tropical tree assemblages. The compilation of previously published (and often incomplete) data on evolutionary relationships of species into a composite phylogeny of the taxa in a forest, through such programs as Phylomatic, has proven useful in building community phylogenies although often of limited resolution. Recently, DNA barcodes have been used to construct a robust community phylogeny for nearly 300 tree species in a forest dynamics plot in Panama using a supermatrix method. In that study sequence data from three barcode loci were used to generate a well-resolved species-level phylogeny.

METHODOLOGY/PRINCIPAL FINDINGS

Here we expand upon this earlier investigation and present results on the use of a phylogenetic constraint tree to generate a community phylogeny for a diverse, tropical forest dynamics plot in Puerto Rico. This enhanced method of phylogenetic reconstruction insures the congruence of the barcode phylogeny with broadly accepted hypotheses on the phylogeny of flowering plants (i.e., APG III) regardless of the number and taxonomic breadth of the taxa sampled. We also compare maximum parsimony versus maximum likelihood estimates of community phylogenetic relationships as well as evaluate the effectiveness of one- versus two- versus three-gene barcodes in resolving community evolutionary history.

CONCLUSIONS/SIGNIFICANCE

As first demonstrated in the Panamanian forest dynamics plot, the results for the Puerto Rican plot illustrate that highly resolved phylogenies derived from DNA barcode sequence data combined with a constraint tree based on APG III are particularly useful in comparative analysis of phylogenetic diversity and will enhance research on the interface between community ecology and evolution.

Plant DNA barcodes and a community phylogeny of a tropical forest dynamics plot in Panama

The assembly of DNA barcode libraries is particularly relevant within species-rich natural communities for which accurate species identifications will enable detailed ecological forensic studies. In addition, well-resolved molecular phylogenies derived from these DNA barcode sequences have the potential to improve investigations of the mechanisms underlying community assembly and functional trait evolution. To date, no studies have effectively applied DNA barcodes sensu strictu in this manner. In this report, we demonstrate that a three-locus DNA barcode when applied to 296 species of woody trees, shrubs, and palms found within the 50-ha Forest Dynamics Plot on Barro Colorado Island (BCI), Panama, resulted in >98% correct identifications. These DNA barcode sequences are also used to reconstruct a robust community phylogeny employing a supermatrix method for 281 of the 296 plant species in the plot. The three-locus barcode data were sufficient to reliably reconstruct evolutionary relationships among the plant taxa in the plot that are congruent with the broadly accepted phylogeny of flowering plants (APG II). Earlier work on the phylogenetic structure of the BCI forest dynamics plot employing less resolved phylogenies reveals significant differences in evolutionary and ecological inferences compared with our data and suggests that unresolved community phylogenies may have increased type I and type II errors. These results illustrate how highly resolved phylogenies based on DNA barcode sequence data will enhance research focused on the interface between community ecology and evolution.

Bacterial communities associated with flea vectors of plague

The microbial flora associated with fleas may affect their ability to transmit specific pathogens, including Yersinia pestis, and also could be used to develop paratransgenesis-based approaches to interfere with transmission. To begin addressing this hypothesis, the microbial flora associated with the relatively efficient Y. pestis vectors Xenopsylla cheopis (Rothschild) (Siphonaptera: Pulicidae) and Oropsyllamontana (Baker) (Siphonaptera: Ceratophyllidae), and the inefficient vector Ctenocephalides felis felis (Bouché) (Siphonaptera: Pulicidae) were investigated using polymerase chain reaction amplification of 16S rDNA genes. DNA sequencing revealed that these species harbor distinct communities of microbial flora and suggest that Acinetobacter sp. might be used in developing anti-transmission strategies.

DNA BARCODING IN LAND PLANTS: DEVELOPING STANDARDS TO QUANTIFY AND MAXIMIZE SUCCESS

The selection of a DNA barcode in plants has been impeded in part due to the relatively low rates of nucleotide substitution observed at the most accessible plastid markers. However, the absence of consensus also reflects a lack of standards for comparing potential barcode markers. While many publications have suggested a host of plant DNA barcodes, the studies cannot be readily compared with each other through any quantitative or statistical parameter, partly because they put forward no single compelling rationale relevant to the adoption of a DNA barcode in plants. Here, we argue that the efficacy of any particular plant DNA barcode selection should reflect the anticipated performance of the resulting barcode database in assignment of a query sequence to species. While legitimate scientific disagreement exists over the criteria relevant to "database performance", the notion gives a unifying rationale for prioritizing selection criteria. Accordingly, we suggest a measure of barcode efficacy based on the rationale of database performance, "the probability of correct identification" (PCI). Moreover, the definition of PCI is left flexible enough to handle most of the scientific disagreement over how to best evaluate DNA barcodes. Finally, we consider how different types of barcodes might require different methods of analysis and database design and indicate how the analysis might affect the selection of the most broadly effective barcode for land plants.

Yersinia pestis biofilm in the flea vector and its role in the transmission of plague

Transmission by fleabite is a relatively recent evolutionary adaptation of Yersinia pestis, the bacterial agent of bubonic plague. To produce a transmissible infection, Y. pestis grows as an attached biofilm in the foregut of the flea vector. Biofilm formation both in the flea foregut and in vitro is dependent on an extracellular matrix (ECM) synthesized by the Yersinia hms gene products. The hms genes are similar to the pga and ica genes of Escherichia coli and Staphylococcus epidermidis, respectively, that act to synthesize a poly-beta-1,6-N-acetyl-d-glucosamine ECM required for biofilm formation. As with extracellular polysaccharide production in many other bacteria, synthesis of the Hms-dependent ECM is controlled by intracellular levels of cyclic-di-GMP. Yersinia pseudotuberculosis, the food- and water-borne enteric pathogen from which Y. pestis evolved recently, possesses identical hms genes and can form biofilm in vitro but not in the flea. The genetic changes in Y. pestis that resulted in adapting biofilm-forming capability to the flea gut environment, a critical step in the evolution of vector-borne transmission, have yet to be identified. During a flea bite, Y. pestis is regurgitated into the dermis in a unique biofilm phenotype, and this has implications for the initial interaction with the mammalian innate immune response.

Acute oral toxicity of Yersinia pseudotuberculosis to fleas: implications for the evolution of vector-borne transmission of plague

Yersinia pestis diverged from Yersinia pseudotuberculosis</= 20 000 years ago, during which time it evolved to be transmitted by fleas. In comparing the ability of these closely related species to infect the rat flea Xenopsylla cheopis, we found that Y. pseudotuberculosis, unlike Y. pestis, is orally toxic to fleas. Fleas showed signs of acute toxicity, including diarrhoea, immediately after feeding on blood containing Y. pseudotuberculosis in response to protein toxin(s) produced by the bacteria. Adherence of Y. pseudotuberculosis to the midgut and large intracellular vacuoles in midgut epithelial cells were detected during the first 24 h after infection. The insect pathogen Photorhabdus luminescens and its TcdA1 and TcdB1-TccC1 insecticidal toxin complexes were similarly toxic to fleas, implicating the toxin complex (tc) genes also present in Yersinia species. However, the Y. pestis and Y. pseudotuberculosis TcaAB and TcaC-TccC proteins were non-toxic to fleas, and Y. pseudotuberculosis mutants deleted of tc genes retained acute toxicity. Our results indicate that loss of one or more insect gut toxins was a critical step in the recent evolution of flea-borne transmission in the genus Yersinia. Changes in the tc insecticidal genes do not appear to have been responsible, but may have had other effects on Yersinia-flea interactions.

A two-locus global DNA barcode for land plants: the coding rbcL gene complements the non-coding trnH-psbA spacer region

BACKGROUND

A useful DNA barcode requires sufficient sequence variation to distinguish between species and ease of application across a broad range of taxa. Discovery of a DNA barcode for land plants has been limited by intrinsically lower rates of sequence evolution in plant genomes than that observed in animals. This low rate has complicated the trade-off in finding a locus that is universal and readily sequenced and has sufficiently high sequence divergence at the species-level.

METHODOLOGY/PRINCIPAL FINDINGS

Here, a global plant DNA barcode system is evaluated by comparing universal application and degree of sequence divergence for nine putative barcode loci, including coding and non-coding regions, singly and in pairs across a phylogenetically diverse set of 48 genera (two species per genus). No single locus could discriminate among species in a pair in more than 79% of genera, whereas discrimination increased to nearly 88% when the non-coding trnH-psbA spacer was paired with one of three coding loci, including rbcL. In silico trials were conducted in which DNA sequences from GenBank were used to further evaluate the discriminatory power of a subset of these loci. These trials supported the earlier observation that trnH-psbA coupled with rbcL can correctly identify and discriminate among related species.

CONCLUSIONS/SIGNIFICANCE

A combination of the non-coding trnH-psbA spacer region and a portion of the coding rbcL gene is recommended as a two-locus global land plant barcode that provides the necessary universality and species discrimination.