Sheila Podell

De novo metagenomic assembly reveals abundant novel major lineage of Archaea in hypersaline microbial communities

This study describes reconstruction of two highly unusual archaeal genomes by de novo metagenomic assembly of multiple, deeply sequenced libraries from surface waters of Lake Tyrrell (LT), a hypersaline lake in NW Victoria, Australia. Lineage-specific probes were designed using the assembled genomes to visualize these novel archaea, which were highly abundant in the 0.1–0.8 μm size fraction of lake water samples. Gene content and inferred metabolic capabilities were highly dissimilar to all previously identified hypersaline microbial species. Distinctive characteristics included unique amino acid composition, absence of Gvp gas vesicle proteins, atypical archaeal metabolic pathways and unusually small cell size (approximately 0.6 μm diameter). Multi-locus phylogenetic analyses demonstrated that these organisms belong to a new major euryarchaeal lineage, distantly related to halophilic archaea of class Halobacteria. Consistent with these findings, we propose creation of a new archaeal class, provisionally named ‘Nanohaloarchaea’. In addition to their high abundance in LT surface waters, we report the prevalence of Nanohaloarchaea in other hypersaline environments worldwide. The simultaneous discovery and genome sequencing of a novel yet ubiquitous lineage of uncultivated microorganisms demonstrates that even historically well-characterized environments can reveal unexpected diversity when analyzed by metagenomics, and advances our understanding of the ecology of hypersaline environments and the evolutionary history of the archaea.

The Natural Product Domain Seeker NaPDoS: A Phylogeny Based Bioinformatic Tool to Classify Secondary Metabolite Gene Diversity

New bioinformatic tools are needed to analyze the growing volume of DNA sequence data. This is especially true in the case of secondary metabolite biosynthesis, where the highly repetitive nature of the associated genes creates major challenges for accurate sequence assembly and analysis. Here we introduce the web tool Natural Product Domain Seeker (NaPDoS), which provides an automated method to assess the secondary metabolite biosynthetic gene diversity and novelty of strains or environments. NaPDoS analyses are based on the phylogenetic relationships of sequence tags derived from polyketide synthase (PKS) and non-ribosomal peptide synthetase (NRPS) genes, respectively. The sequence tags correspond to PKS-derived ketosynthase domains and NRPS-derived condensation domains and are compared to an internal database of experimentally characterized biosynthetic genes. NaPDoS provides a rapid mechanism to extract and classify ketosynthase and condensation domains from PCR products, genomes, and metagenomic datasets. Close database matches provide a mechanism to infer the generalized structures of secondary metabolites while new phylogenetic lineages provide targets for the discovery of new enzyme architectures or mechanisms of secondary metabolite assembly. Here we outline the main features of NaPDoS and test it on four draft genome sequences and two metagenomic datasets. The results provide a rapid method to assess secondary metabolite biosynthetic gene diversity and richness in organisms or environments and a mechanism to identify genes that may be associated with uncharacterized biochemistry.

Genomic islands link secondary metabolism to functional adaptation in marine Actinobacteria

Genomic islands have been shown to harbor functional traits that differentiate ecologically distinct populations of environmental bacteria. A comparative analysis of the complete genome sequences of the marine Actinobacteria Salinispora tropica and Salinispora arenicola reveals that 75% of the species-specific genes are located in 21 genomic islands. These islands are enriched in genes associated with secondary metabolite biosynthesis providing evidence that secondary metabolism is linked to functional adaptation. Secondary metabolism accounts for 8.8% and 10.9% of the genes in the S. tropica and S. arenicola genomes, respectively, and represents the major functional category of annotated genes that differentiates the two species. Genomic islands harbor all 25 of the species-specific biosynthetic pathways, the majority of which occur in S. arenicola and may contribute to the cosmopolitan distribution of this species. Genome evolution is dominated by gene duplication and acquisition, which in the case of secondary metabolism provide immediate opportunities for the production of new bioactive products. Evidence that secondary metabolic pathways are exchanged horizontally, coupled with earlier evidence for fixation among globally distributed populations, supports a functional role and suggests that the acquisition of natural product biosynthetic gene clusters represents a previously unrecognized force driving bacterial diversification. Species-specific differences observed in clustered regularly interspaced short palindromic repeat sequences suggest that S. arenicola may possess a higher level of phage immunity, whereas a highly duplicated family of polymorphic membrane proteins provides evidence for a new mechanism of marine adaptation in Gram-positive bacteria.

Predicting N-terminal myristoylation sites in plant proteins

BackgroundN-terminal myristoylation plays a vital role in membrane targeting and signal transduction in plant responses to environmental stress. Although N-myristoyltransferase enzymatic function is conserved across plant, animal, and fungal kingdoms, exact substrate specificities vary, making it difficult to predict protein myristoylation accurately within specific taxonomic groups.ResultsA new method for predicting N-terminal myristoylation sites specifically in plants has been developed and statistically tested for sensitivity, specificity, and robustness. Compared to previously available methods, the new model is both more sensitive in detecting known positives, and more selective in avoiding false positives. Scores of myristoylated and non-myristoylated proteins are more widely separated than with other methods, greatly reducing ambiguity and the number of sequences giving intermediate, uninformative results. The prediction model is available at http://plantsp.sdsc.edu/myrist.html.ConclusionSuperior performance of the new model is due to the selection of a plant-specific training set, covering 266 unique sequence examples from 40 different species, the use of a probability-based hidden Markov model to obtain predictive scores, and a threshold cutoff value chosen to provide maximum positive-negative discrimination. The new model has been used to predict 589 plant proteins likely to contain N-terminal myristoylation signals, and to analyze the functional families in which these proteins occur.

DarkHorse: a method for genome-wide prediction of horizontal gene transfer

A new approach to rapid, genome-wide identification and ranking of horizontal transfer candidate proteins is presented. The method is quantitative, reproducible, and computationally undemanding. It can be combined with genomic signature and/or phylogenetic tree-building procedures to improve accuracy and efficiency. The method is also useful for retrospective assessments of horizontal transfer prediction reliability, recognizing orthologous sequences that may have been previously overlooked or unavailable. These features are demonstrated in bacterial, archaeal, and eukaryotic examples.

Metabolic versatility of the Riftia pachyptila endosymbiont revealed through metagenomics

The facultative symbiont of Riftia pachyptila, named here Candidatus Endoriftia persephone, has evaded culture to date, but much has been learned regarding this symbiosis over the past three decades since its discovery. The symbiont population metagenome was sequenced in order to gain insight into its physiology. The population genome indicates that the symbionts use a partial Calvin-Benson Cycle for carbon fixation and the reverse TCA cycle (an alternative pathway for carbon fixation) that contains an unusual ATP citrate lyase. The presence of all genes necessary for heterotrophic metabolism, a phosphotransferase system, and dicarboxylate and ABC transporters indicate that the symbiont can live mixotrophically. The metagenome has a large suite of signal transduction, defence (both biological and environmental) and chemotaxis mechanisms. The physiology of Candidatus Endoriftia persephone is explored with respect to functionality while associated with a eukaryotic host, versus free-living in the hydrothermal environment.

Genomic insights into the physiology and ecology of the marine filamentous cyanobacterium Lyngbya majuscula.

Filamentous cyanobacteria of the genus Lyngbya are important contributors to coral reef ecosystems, occasionally forming dominant cover and impacting the health of many other co-occurring organisms. Moreover, they are extraordinarily rich sources of bioactive secondary metabolites, with 35% of all reported cyanobacterial natural products deriving from this single pantropical genus. However, the true natural product potential and life strategies of Lyngbya strains are poorly understood because of phylogenetic ambiguity, lack of genomic information, and their close associations with heterotrophic bacteria and other cyanobacteria. To gauge the natural product potential of Lyngbya and gain insights into potential microbial interactions, we sequenced the genome of Lyngbya majuscula 3L, a Caribbean strain that produces the tubulin polymerization inhibitor curacin A and the molluscicide barbamide, using a combination of Sanger and 454 sequencing approaches. Whereas ∼293,000 nucleotides of the draft genome are putatively dedicated to secondary metabolism, this is far too few to encode a large suite of Lyngbya metabolites, suggesting Lyngbya metabolites are strain specific and may be useful in species delineation. Our analysis revealed a complex gene regulatory network, including a large number of sigma factors and other regulatory proteins, indicating an enhanced ability for environmental adaptation or microbial associations. Although Lyngbya species are reported to fix nitrogen, nitrogenase genes were not found in the genome or by PCR of genomic DNA. Subsequent growth experiments confirmed that L. majuscula 3L is unable to fix atmospheric nitrogen. These unanticipated life history characteristics challenge current views of the genus Lyngbya.

Characterization of cyanobacterial hydrocarbon composition and distribution of biosynthetic pathways

Cyanobacteria possess the unique capacity to naturally produce hydrocarbons from fatty acids. Hydrocarbon compositions of thirty-two strains of cyanobacteria were characterized to reveal novel structural features and insights into hydrocarbon biosynthesis in cyanobacteria. This investigation revealed new double bond (2- and 3-heptadecene) and methyl group positions (3-, 4- and 5-methylheptadecane) for a variety of strains. Additionally, results from this study and literature reports indicate that hydrocarbon production is a universal phenomenon in cyanobacteria. All cyanobacteria possess the capacity to produce hydrocarbons from fatty acids yet not all accomplish this through the same metabolic pathway. One pathway comprises a two-step conversion of fatty acids first to fatty aldehydes and then alkanes that involves a fatty acyl ACP reductase (FAAR) and aldehyde deformylating oxygenase (ADO). The second involves a polyketide synthase (PKS) pathway that first elongates the acyl chain followed by decarboxylation to produce a terminal alkene (olefin synthase, OLS). Sixty-one strains possessing the FAAR/ADO pathway and twelve strains possessing the OLS pathway were newly identified through bioinformatic analyses. Strains possessing the OLS pathway formed a cohesive phylogenetic clade with the exception of three Moorea strains and Leptolyngbya sp. PCC 6406 which may have acquired the OLS pathway via horizontal gene transfer. Hydrocarbon pathways were identified in one-hundred-forty-two strains of cyanobacteria over a broad phylogenetic range and there were no instances where both the FAAR/ADO and the OLS pathways were found together in the same genome, suggesting an unknown selective pressure maintains one or the other pathway, but not both.

Genomic Insights into Mn(II) Oxidation by the Marine Alphaproteobacterium Aurantimonas sp. Strain SI85-9A1

ABSTRACT Microbial Mn(II) oxidation has important biogeochemical consequences in marine, freshwater, and terrestrial environments, but many aspects of the physiology and biochemistry of this process remain obscure. Here, we report genomic insights into Mn(II) oxidation by the marine alphaproteobacterium Aurantimonas sp. strain SI85-9A1, isolated from the oxic/anoxic interface of a stratified fjord. The SI85-9A1 genome harbors the genetic potential for metabolic versatility, with genes for organoheterotrophy, methylotrophy, oxidation of sulfur and carbon monoxide, the ability to grow over a wide range of O2 concentrations (including microaerobic conditions), and the complete Calvin cycle for carbon fixation. Although no growth could be detected under autotrophic conditions with Mn(II) as the sole electron donor, cultures of SI85-9A1 grown on glycerol are dramatically stimulated by addition of Mn(II), suggesting an energetic benefit from Mn(II) oxidation. A putative Mn(II) oxidase is encoded by duplicated multicopper oxidase genes that have a complex evolutionary history including multiple gene duplication, loss, and ancient horizontal transfer events. The Mn(II) oxidase was most abundant in the extracellular fraction, where it cooccurs with a putative hemolysin-type Ca2+-binding peroxidase. Regulatory elements governing the cellular response to Fe and Mn concentration were identified, and 39 targets of these regulators were detected. The putative Mn(II) oxidase genes were not among the predicted targets, indicating that regulation of Mn(II) oxidation is controlled by other factors yet to be identified. Overall, our results provide novel insights into the physiology and biochemistry of Mn(II) oxidation and reveal a genome specialized for life at the oxic/anoxic interface.

Genomic resources for songbird research and their use in characterizing gene expression during brain development

Vocal learning and neuronal replacement have been studied extensively in songbirds, but until recently, few molecular and genomic tools for songbird research existed. Here we describe new molecular/genomic resources developed in our laboratory. We made cDNA libraries from zebra finch (Taeniopygia guttata) brains at different developmental stages. A total of 11,000 cDNA clones from these libraries, representing 5,866 unique gene transcripts, were randomly picked and sequenced from the 3′ ends. A web-based database was established for clone tracking, sequence analysis, and functional annotations. Our cDNA libraries were not normalized. Sequencing ESTs without normalization produced many developmental stage-specific sequences, yielding insights into patterns of gene expression at different stages of brain development. In particular, the cDNA library made from brains at posthatching day 30–50, corresponding to the period of rapid song system development and song learning, has the most diverse and richest set of genes expressed. We also identified five microRNAs whose sequences are highly conserved between zebra finch and other species. We printed cDNA microarrays and profiled gene expression in the high vocal center of both adult male zebra finches and canaries (Serinus canaria). Genes differentially expressed in the high vocal center were identified from the microarray hybridization results. Selected genes were validated by in situ hybridization. Networks among the regulated genes were also identified. These resources provide songbird biologists with tools for genome annotation, comparative genomics, and microarray gene expression analysis.