Dynamics of an experimental microbial invasion.

The ecological dynamics underlying species invasions have been a major focus of research in macroorganisms for the last five decades. However, we still know little about the processes behind invasion by unicellular organisms. To expand our knowledge of microbial invasions, we studied the roles of propagule pressure, nutrient supply, and biotic resistance in the invasion success of a freshwater invasive alga, Prymnesium parvum, using microcosms containing natural freshwater microbial assemblages. Microcosms were subjected to a factorial design with two levels of nutrient-induced diversity and three levels of propagule pressure, and incubated for 7 d, during which P. parvum densities and microbial community composition were tracked. Successful invasion occurred in microcosms receiving high propagule pressure whereas nutrients or community diversity played no role in invasion success. Invaded communities experienced distinctive changes in composition compared with communities where the invasion was unsuccessful. Successfully invaded microbial communities had an increased abundance of fungi and ciliates, and decreased abundances of diatoms and cercozoans. Many of these changes mirrored the microbial community changes detected during a natural P. parvum bloom in the source system. This role of propagule pressure is particularly relevant for P. parvum in the reservoir-dominated southern United States because this species can form large, sustained blooms that can generate intense propagule pressures for downstream sites. Human impact and global climate change are currently causing widespread environmental changes in most southern US freshwater systems that may facilitate P. parvum establishment and, when coupled with strong propagule pressure, could put many more systems at risk for invasion.

Arhodomonas sp. strain Seminole and its genetic potential to degrade aromatic compounds under high-salinity conditions.

Arhodomonas sp. strain Seminole was isolated from a crude oil-impacted brine soil and shown to degrade benzene, toluene, phenol, 4-hydroxybenzoic acid (4-HBA), protocatechuic acid (PCA), and phenylacetic acid (PAA) as the sole sources of carbon at high salinity. Seminole is a member of the genus Arhodomonas in the class Gammaproteobacteria, sharing 96% 16S rRNA gene sequence similarity with Arhodomonas aquaeolei HA-1. Analysis of the genome predicted a number of catabolic genes for the metabolism of benzene, toluene, 4-HBA, and PAA. The predicted pathways were corroborated by identification of enzymes present in the cytosolic proteomes of cells grown on aromatic compounds using liquid chromatography-mass spectrometry. Genome analysis predicted a cluster of 19 genes necessary for the breakdown of benzene or toluene to acetyl coenzyme A (acetyl-CoA) and pyruvate. Of these, 12 enzymes were identified in the proteome of toluene-grown cells compared to lactate-grown cells. Genomic analysis predicted 11 genes required for 4-HBA degradation to form the tricarboxylic acid (TCA) cycle intermediates. Of these, proteomic analysis of 4-HBA-grown cells identified 6 key enzymes involved in the 4-HBA degradation pathway. Similarly, 15 genes needed for the degradation of PAA to the TCA cycle intermediates were predicted. Of these, 9 enzymes of the PAA degradation pathway were identified only in PAA-grown cells and not in lactate-grown cells. Overall, we were able to reconstruct catabolic steps for the breakdown of a variety of aromatic compounds in an extreme halophile, strain Seminole. Such knowledge is important for understanding the role of Arhodomonas spp. in the natural attenuation of hydrocarbon-impacted hypersaline environments.

Coral-zooxanthellae meta-transcriptomics reveals integrated response to pollutant stress.

BACKGROUND: Corals represent symbiotic meta-organisms that require harmonization among the coral animal, photosynthetic zooxanthellae and associated microbes to survive environmental stresses. We investigated integrated-responses among coral and zooxanthellae in the scleractinian coral Acropora formosa in response to an emerging marine pollutant, the munitions constituent, 1,3,5-trinitro-1,3,5 triazine (RDX; 5 day exposures to 0 (control), 0.5, 0.9, 1.8, 3.7, and 7.2 mg/L, measured in seawater). RESULTS: RDX accumulated readily in coral soft tissues with bioconcentration factors ranging from 1.1 to 1.5. Next-generation sequencing of a normalized meta-transcriptomic library developed for the eukaryotic components of the A. formosa coral holobiont was leveraged to conduct microarray-based global transcript expression analysis of integrated coral/zooxanthellae responses to the RDX exposure. Total differentially expressed transcripts (DET) increased with increasing RDX exposure concentrations as did the proportion of zooxanthellae DET relative to the coral animal. Transcriptional responses in the coral demonstrated higher sensitivity to RDX compared to zooxanthellae where increased expression of gene transcripts coding xenobiotic detoxification mechanisms (i.e. cytochrome P450 and UDP glucuronosyltransferase 2 family) were initiated at the lowest exposure concentration. Increased expression of these detoxification mechanisms was sustained at higher RDX concentrations as well as production of a physical barrier to exposure through a 40% increase in mucocyte density at the maximum RDX exposure. At and above the 1.8 mg/L exposure concentration, DET coding for genes involved in central energy metabolism, including photosynthesis, glycolysis and electron-transport functions, were decreased in zooxanthellae although preliminary data indicated that zooxanthellae densities were not affected. In contrast, significantly increased transcript expression for genes involved in cellular energy production including glycolysis and electron-transport pathways was observed in the coral animal. CONCLUSIONS: Transcriptional network analysis for central energy metabolism demonstrated highly correlated responses to RDX among the coral animal and zooxanthellae indicative of potential compensatory responses to lost photosynthetic potential within the holobiont. These observations underscore the potential for complex integrated responses to RDX exposure among species comprising the coral holobiont and highlight the need to understand holobiont-species interactions to accurately assess pollutant impacts.

Fluorescence in situ hybridization and optical mapping to correct scaffold arrangement in the tomato genome.

The order and orientation (arrangement) of all 91 sequenced scaffolds in the 12 pseudomolecules of the recently published tomato (Solanum lycopersicum, 2n = 2x = 24) genome sequence were positioned based on marker order in a high-density linkage map. Here, we report the arrangement of these scaffolds determined by two independent physical methods, bacterial artificial chromosome-fluorescence in situ hybridization (BAC-FISH) and optical mapping. By localizing BACs at the ends of scaffolds to spreads of tomato synaptonemal complexes (pachytene chromosomes), we showed that 45 scaffolds, representing one-third of the tomato genome, were arranged differently than predicted by the linkage map. These scaffolds occur mostly in pericentric heterochromatin where 77% of the tomato genome is located and where linkage mapping is less accurate due to reduced crossing over. Although useful for only part of the genome, optical mapping results were in complete agreement with scaffold arrangement by FISH but often disagreed with scaffold arrangement based on the linkage map. The scaffold arrangement based on FISH and optical mapping changes the positions of hundreds of markers in the linkage map, especially in heterochromatin. These results suggest that similar errors exist in pseudomolecules from other large genomes that have been assembled using only linkage maps to predict scaffold arrangement, and these errors can be corrected using FISH and/or optical mapping. Of note, BAC-FISH also permits estimates of the sizes of gaps between scaffolds, and unanchored BACs are often visualized by FISH in gaps between scaffolds and thus represent starting points for filling these gaps.

BAC-pool sequencing and analysis of large segments of A12 and D12 homoeologous chromosomes in upland cotton.

Although new and emerging next-generation sequencing (NGS) technologies have reduced sequencing costs significantly, much work remains to implement them for de novo sequencing of complex and highly repetitive genomes such as the tetraploid genome of Upland cotton (Gossypium hirsutum L.). Herein we report the results from implementing a novel, hybrid Sanger/454-based BAC-pool sequencing strategy using minimum tiling path (MTP) BACs from Ctg-3301 and Ctg-465, two large genomic segments in A12 and D12 homoeologous chromosomes (Ctg). To enable generation of longer contig sequences in assembly, we implemented a hybrid assembly method to process ~35x data from 454 technology and 2.8-3x data from Sanger method. Hybrid assemblies offered higher sequence coverage and better sequence assemblies. Homology studies revealed the presence of retrotransposon regions like Copia and Gypsy elements in these contigs and also helped in identifying new genomic SSRs. Unigenes were anchored to the sequences in Ctg-3301 and Ctg-465 to support the physical map. Gene density, gene structure and protein sequence information derived from protein prediction programs were used to obtain the functional annotation of these genes. Comparative analysis of both contigs with Arabidopsis genome exhibited synteny and microcollinearity with a conserved gene order in both genomes. This study provides insight about use of MTP-based BAC-pool sequencing approach for sequencing complex polyploid genomes with limited constraints in generating better sequence assemblies to build reference scaffold sequences. Combining the utilities of MTP-based BAC-pool sequencing with current longer and short read NGS technologies in multiplexed format would provide a new direction to cost-effectively and precisely sequence complex plant genomes.

Seasonality and disturbance: annual pattern and response of the bacterial and microbial eukaryotic assemblages in a freshwater ecosystem.

High-throughput pyrosequencing of SSU rDNA genes was used to obtain monthly snapshots of eukaryotic and bacterial diversity and community structure at two locations in Lake Texoma, a low salinity lake in the south central United States, over 1 year. The lake experienced two disturbance events (i) a localized bloom of Prymnesium parvum restricted to one of the locations that lasted from January to April, and (ii) a large (17 cm), global rain event in the beginning of May, overlaid onto seasonal environmental change. Eukaryotic species richness as well as both eukaryotic and bacterial community similarity exhibited seasonal patterns, including distinct responses to the rain event. The P. parvum bloom created a natural experiment in which to directly explore the effects of an Ecosystem Disruptive Algal Bloom (EDAB) on the microbial community separated from seasonal changes. Microbial species richness was unaffected by the bloom, however, the eukaryotic community structure (evenness) and the patterns of both eukaryotic and bacterial community similarity at bloom and non-bloom sites were statistically distinct during the 4 months of the bloom. These results indicate that physical and biological disturbances as well as seasonal environmental forces contribute to the structure of both the eukaryotic and bacterial communities.

Plant-symbiotic fungi as chemical engineers: multi-genome analysis of the clavicipitaceae reveals dynamics of alkaloid loci.

The fungal family Clavicipitaceae includes plant symbionts and parasites that produce several psychoactive and bioprotective alkaloids. The family includes grass symbionts in the epichloae clade (Epichloë and Neotyphodium species), which are extraordinarily diverse both in their host interactions and in their alkaloid profiles. Epichloae produce alkaloids of four distinct classes, all of which deter insects, and some-including the infamous ergot alkaloids-have potent effects on mammals. The exceptional chemotypic diversity of the epichloae may relate to their broad range of host interactions, whereby some are pathogenic and contagious, others are mutualistic and vertically transmitted (seed-borne), and still others vary in pathogenic or mutualistic behavior. We profiled the alkaloids and sequenced the genomes of 10 epichloae, three ergot fungi (Claviceps species), a morning-glory symbiont (Periglandula ipomoeae), and a bamboo pathogen (Aciculosporium take), and compared the gene clusters for four classes of alkaloids. Results indicated a strong tendency for alkaloid loci to have conserved cores that specify the skeleton structures and peripheral genes that determine chemical variations that are known to affect their pharmacological specificities. Generally, gene locations in cluster peripheries positioned them near to transposon-derived, AT-rich repeat blocks, which were probably involved in gene losses, duplications, and neofunctionalizations. The alkaloid loci in the epichloae had unusual structures riddled with large, complex, and dynamic repeat blocks. This feature was not reflective of overall differences in repeat contents in the genomes, nor was it characteristic of most other specialized metabolism loci. The organization and dynamics of alkaloid loci and abundant repeat blocks in the epichloae suggested that these fungi are under selection for alkaloid diversification. We suggest that such selection is related to the variable life histories of the epichloae, their protective roles as symbionts, and their associations with the highly speciose and ecologically diverse cool-season grasses.

Detection of members of the Secoviridae in the Tallgrass Prairie Preserve, Osage County, Oklahoma, USA.

Viruses are most frequently discovered because they cause disease. To expand knowledge of plant-associated viruses beyond these narrow constraints, non-cultivated plants of the Tallgrass Prairie of the United States were systematically surveyed for evidence of viruses. This report discusses putative viruses of the family Secoviridae identified by the survey. Sequence analysis suggests the presence of at least six viruses in the study site, including Bean pod mottle virus, Maize chlorotic dwarf virus, three previously undescribed viruses within the subfamily Comovirinae and one unclassifiable virus.

Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine.

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

Elevated carbon dioxide alters the structure of soil microbial communities.

Pyrosequencing analysis of 16S rRNA genes was used to examine impacts of elevated CO(2) (eCO(2)) on soil microbial communities from 12 replicates each from ambient CO(2) (aCO(2)) and eCO(2) settings. The results suggest that the soil microbial community composition and structure significantly altered under conditions of eCO(2), which was closely associated with soil and plant properties.

Molecular characterization, ecology, and epidemiology of a novel Tymovirus in Asclepias viridis from Oklahoma.

Native virus-plant interactions require more understanding and their study will provide a basis from which to identify potential sources of emerging destructive viruses in crops. A novel tymovirus sequence was detected in Asclepias viridis (green milkweed), a perennial growing in a natural setting in the Tallgrass Prairie Preserve (TGPP) of Oklahoma. It was abundant within and frequent among A. viridis plants and, to varying extents, within other dicotyledonous and one grass (Panicum virgatum) species obtained from the TGPP. Extracts from A. viridis containing the sequence were infectious to a limited number of species. The virus genome was cloned and determined to be closely related to Kennedya yellow mosaic virus. The persistence of the virus within the Oklahoma A. viridis population was monitored for five successive years. Virus was present in a high percentage of plants within representative areas of the TGPP in all years and was spreading to additional plants. Virus was present in regions adjacent to the TGPP but not in plants sampled from central and south-central Oklahoma. Virus was present in the underground caudex of the plant during the winter, suggesting overwintering in this tissue. The RNA sequence encoding the virus coat protein varied considerably between individual plants (≈3%), likely due to drift rather than selection. An infectious clone was constructed and the virus was named Asclepias asymptomatic virus (AsAV) due to the absence of obvious symptoms on A. viridis.

Detection of members of the Tombusviridae in the Tallgrass Prairie Preserve, Osage County, Oklahoma, USA.

Viruses are most frequently discovered because they cause disease in organisms of importance to humans. To expand knowledge of plant-associated viruses beyond these narrow constraints, non-cultivated plants of the Tallgrass Prairie Preserve, Osage County, Oklahoma, USA were systematically surveyed for evidence of the presence of viruses. This report discusses viruses of the family Tombusviridae putatively identified by the survey. Evidence of two carmoviruses, a tombusvirus, a panicovirus and an unclassifiable tombusvirid was found. The complete genome sequence was obtained for putative TGP carmovirus 1 from the legume Lespedeza procumbens, and the virus was detected in several other plant species including the fern Pellaea atropurpurea. Phylogenetic analysis of the sequence and partial sequence of a related virus supported strongly the placement of these viruses in the genus Carmovirus. Polymorphisms in the sequences suggested existence of two populations of TGP carmovirus 1 in the study area and year-to-year variations in infection by TGP carmovirus 3.

Genotype and phenotypes of an intestine-adapted Escherichia coli K-12 mutant selected by animal passage for superior colonization.

We previously isolated a spontaneous mutant of Escherichia coli K-12, strain MG1655, following passage through the streptomycin-treated mouse intestine, that has colonization traits superior to the wild-type parent strain (M. P. Leatham et al., Infect. Immun. 73:8039-8049, 2005). This intestine-adapted strain (E. coli MG1655*) grew faster on several different carbon sources than the wild type and was nonmotile due to deletion of the flhD gene. We now report the results of several high-throughput genomic analysis approaches to further characterize E. coli MG1655*. Whole-genome pyrosequencing did not reveal any changes on its genome, aside from the deletion at the flhDC locus, that could explain the colonization advantage of E. coli MG1655*. Microarray analysis revealed modest yet significant induction of catabolic gene systems across the genome in both E. coli MG1655* and an isogenic flhD mutant constructed in the laboratory. Catabolome analysis with Biolog GN2 microplates revealed an enhanced ability of both E. coli MG1655* and the isogenic flhD mutant to oxidize a variety of carbon sources. The results show that intestine-adapted E. coli MG1655* is more fit than the wild type for intestinal colonization, because loss of FlhD results in elevated expression of genes involved in carbon and energy metabolism, resulting in more efficient carbon source utilization and a higher intestinal population. Hence, mutations that enhance metabolic efficiency confer a colonization advantage.

Genome-wide DNA methylation maps in follicular lymphoma cells determined by methylation-enriched bisulfite sequencing.

BACKGROUND: Follicular lymphoma (FL) is a form of non-Hodgkin's lymphoma (NHL) that arises from germinal center (GC) B-cells. Despite the significant advances in immunotherapy, FL is still not curable. Beyond transcriptional profiling and genomics datasets, there currently is no epigenome-scale dataset or integrative biology approach that can adequately model this disease and therefore identify novel mechanisms and targets for successful prevention and treatment of FL. METHODOLOGY/PRINCIPAL FINDINGS: We performed methylation-enriched genome-wide bisulfite sequencing of FL cells and normal CD19(+) B-cells using 454 sequencing technology. The methylated DNA fragments were enriched with methyl-binding proteins, treated with bisulfite, and sequenced using the Roche-454 GS FLX sequencer. The total number of bases covered in the human genome was 18.2 and 49.3 million including 726,003 and 1.3 million CpGs in FL and CD19(+) B-cells, respectively. 11,971 and 7,882 methylated regions of interest (MRIs) were identified respectively. The genome-wide distribution of these MRIs displayed significant differences between FL and normal B-cells. A reverse trend in the distribution of MRIs between the promoter and the gene body was observed in FL and CD19(+) B-cells. The MRIs identified in FL cells also correlated well with transcriptomic data and ChIP-on-Chip analyses of genome-wide histone modifications such as tri-methyl-H3K27, and tri-methyl-H3K4, indicating a concerted epigenetic alteration in FL cells. CONCLUSIONS/SIGNIFICANCE: This study is the first to provide a large scale and comprehensive analysis of the DNA methylation sequence composition and distribution in the FL epigenome. These integrated approaches have led to the discovery of novel and frequent targets of aberrant epigenetic alterations. The genome-wide bisulfite sequencing approach developed here can be a useful tool for profiling DNA methylation in clinical samples.

Insights into evolution of multicellular fungi from the assembled chromosomes of the mushroom Coprinopsis cinerea (Coprinus cinereus).

The mushroom Coprinopsis cinerea is a classic experimental model for multicellular development in fungi because it grows on defined media, completes its life cycle in 2 weeks, produces some 10(8) synchronized meiocytes, and can be manipulated at all stages in development by mutation and transformation. The 37-megabase genome of C. cinerea was sequenced and assembled into 13 chromosomes. Meiotic recombination rates vary greatly along the chromosomes, and retrotransposons are absent in large regions of the genome with low levels of meiotic recombination. Single-copy genes with identifiable orthologs in other basidiomycetes are predominant in low-recombination regions of the chromosome. In contrast, paralogous multicopy genes are found in the highly recombining regions, including a large family of protein kinases (FunK1) unique to multicellular fungi. Analyses of P450 and hydrophobin gene families confirmed that local gene duplications drive the expansions of paralogous copies and the expansions occur in independent lineages of Agaricomycotina fungi. Gene-expression patterns from microarrays were used to dissect the transcriptional program of dikaryon formation (mating). Several members of the FunK1 kinase family are differentially regulated during sexual morphogenesis, and coordinate regulation of adjacent duplications is rare. The genomes of C. cinerea and Laccaria bicolor, a symbiotic basidiomycete, share extensive regions of synteny. The largest syntenic blocks occur in regions with low meiotic recombination rates, no transposable elements, and tight gene spacing, where orthologous single-copy genes are overrepresented. The chromosome assembly of C. cinerea is an essential resource in understanding the evolution of multicellularity in the fungi.

High-throughput sequence analysis of Ciona intestinalis SL trans-spliced mRNAs: alternative expression modes and gene function correlates.

Pre-mRNA 5' spliced-leader (SL) trans-splicing occurs in some metazoan groups but not in others. Genome-wide characterization of the trans-spliced mRNA subpopulation has not yet been reported for any metazoan. We carried out a high-throughput analysis of the SL trans-spliced mRNA population of the ascidian tunicate Ciona intestinalis by 454 Life Sciences (Roche) pyrosequencing of SL-PCR-amplified random-primed reverse transcripts of tailbud embryo RNA. We obtained approximately 250,000 high-quality reads corresponding to 8790 genes, approximately 58% of the Ciona total gene number. The great depth of this data revealed new aspects of trans-splicing, including the existence of a significant class of "infrequently trans-spliced" genes, accounting for approximately 28% of represented genes, that generate largely non-trans-spliced mRNAs, but also produce trans-spliced mRNAs, in part through alternative promoter use. Thus, the conventional qualitative dichotomy of trans-spliced versus non-trans-spliced genes should be supplanted by a more accurate quantitative view recognizing frequently and infrequently trans-spliced gene categories. Our data include reads representing approximately 80% of Ciona frequently trans-spliced genes. Our analysis also revealed significant use of closely spaced alternative trans-splice acceptor sites which further underscores the mechanistic similarity of cis- and trans-splicing and indicates that the prevalence of +/-3-nt alternative splicing events at tandem acceptor sites, NAGNAG, is driven by spliceosomal mechanisms, and not nonsense-mediated decay, or selection at the protein level. The breadth of gene representation data enabled us to find new correlations between trans-splicing status and gene function, namely the overrepresentation in the frequently trans-spliced gene class of genes associated with plasma/endomembrane system, Ca(2+) homeostasis, and actin cytoskeleton.

Identification of microRNAs and their targets in switchgrass, a model biofuel plant species.

In recent years, several plant species such as switchgrass, Miscanthus and Brachypodium have been recognized as potential model plant species for cellulosic bioenergy production. Of these, switchgrass has attracted much attention in the United States and worldwide because it can grow well on marginal lands and tolerate frequent drought spells. However, little is known about the basic biology of the traits that control these important characteristics in switchgrass. Genome-encoded approximately 21-24nt microRNAs (miRNAs) have emerged as critical regulators of gene expression important for normal growth and development and adaptation to abiotic stress, including nutrient-deprived conditions. To understand miRNA-guided post-transcriptional gene regulatory networks in this plant species, we sought to identify miRNAs in switchgrass. Using computational and experimental approaches, we identified approximately 20 conserved miRNA families. Temporal expression analysis indicated that some miRNAs have distinct tissue-specific expression, although most are ubiquitously expressed. Unlike in Arabidopsis and other plants, miR395 and miR399 were detected in plants grown on optimal levels of sulfate or phosphate in switchgrass, and were only slightly altered when exposed to sulfate or phosphate deficit conditions. Thirty-seven genes were predicted as targets for miRNAs, and 4 target mRNAs (Squamosa promoter binding-like factor, apetala 2-like, NAC domain containing transcription factor and HD-Zip homologs) were validated by 5'-RACE assays. These findings provide a snapshot of the miRNA component and possible targets in switchgrass.

A novel conjugative plasmid from Enterococcus faecalis E99 enhances resistance to ultraviolet radiation.

Enterococcus faecalis has emerged as a prominent healthcare-associated pathogen frequently encountered in bacteremia, endocarditis, urinary tract infection, and as a leading cause of antibiotic-resistant infections. We recently demonstrated a capacity for high-level biofilm formation by a clinical E. faecalis isolate, E99. This high biofilm-forming phenotype was attributable to a novel locus, designated bee, specifying a pilus at the bacterial cell surface and localized to a large approximately 80 kb conjugative plasmid. To better understand the origin of the bee locus, as well as to potentially identify additional factors important to the biology and pathogenesis of strain E99, we sequenced the entire plasmid. The nucleotide sequence of the plasmid, designated pBEE99, revealed large regions of identity to the previously characterized conjugative plasmid pCF10. In addition to the bee locus, pBEE99 possesses an open reading frame potentially encoding aggregation substance, as well as open reading frames putatively encoding polypeptides with 60% to 99% identity at the amino acid level to proteins involved in regulation of the pheromone response and conjugal transfer of pCF10. However, strain E99 did not respond to the cCF10 pheromone in clumping assays. While pBEE99 was found to be devoid of any readily recognizable antibiotic resistance determinants, it carries two non-identical impB/mucB/samB-type genes, as well as genes potentially encoding a two-component bacteriocin similar to that encoded on pYI14. Although no bacteriocin activity was detected from an OG1RF transconjugant carrying pBEE99 against strain FA2-2, it was approximately an order of magnitude more resistant to ultraviolet radiation. Moreover, curing strain E99 of this plasmid significantly reduced its ability to survive UV exposure. Therefore, pBEE99 represents a novel conjugative plasmid that confers biofilm-forming and enhanced UV resistance traits that might potentially impact the virulence and/or fitness of E. faecalis.

Ecogenomics: using massively parallel pyrosequencing to understand virus ecology.

Environmental samples have been analysed for viruses in metagenomic studies, but these studies have not linked individual viruses to their hosts. We designed a strategy to isolate double-stranded RNA, a hallmark of RNA virus infection, from individual plants and convert this to cDNA with a unique four nucleotide Tag at each end. Using 96 different Tags allowed us to pool samples and still retain the link to the original sample. We then analysed the sequence of pooled samples using massively parallel sequencing with Roche 454 pyrosequencing such that 384 samples could be assessed per picotiter plate. Using this method we have been able to analyse thousands of plants, and we have discovered several thousand new plant viruses, all linked to their specific plant hosts. Here we describe the method in detail, including the results and analysis for eight pools of samples. This technology will be extremely useful in understanding the full scope of plant virus biodiversity.