Hal: an automated pipeline for phylogenetic analyses of genomic data

Important innate differences in determining symbiotic responsiveness in host and non-hosts of arbuscular mycorrhiza

Genetic components that regulate arbuscular mycorrhizal (AM) interactions in hosts and non-hosts are not completely known. Comparative transcriptomic analysis was combined with phylogenetic studies to identify the factors that distinguish AM host from non-host. Mycorrhized host, non-mycorrhized host and non-host cultivars of tomato (Solanum lycopersicum) were subjected to RNA seq analysis. The top 10 differentially expressed genes were subjected to extensive in silico phylogenetic analysis along with 10 more candidate genes that have been previously reported for AM-plant interactions. Seven distantly related hosts and four non-hosts were selected to identify structural differences in selected gene/protein candidates. The screened genes/proteins were subjected to MEME, CODEML and DIVERGE analysis to identify evolutionary patterns that differentiate hosts from non-hosts. Based on the results, candidate genes were categorized as highly influenced (SYMRK and CCaMK), moderately influenced and minimally influenced by evolutionary constraints. We propose that the amino acid and nucleotide changes specific to non-hosts are likely to correspond to aberrations in functionality towards AM symbiosis. This study paves way for future research aimed at understanding innate differences in genetic make-up of AM hosts and non-hosts, in addition to the theory of gene losses from the "AM-symbiotic toolkit".

The genome of a novel isolate of Prochlorococcus from the Red Sea contains transcribed genes for compatible solute biosynthesis

Marine microbes possess genomic and physiological adaptations to cope with varying environmental conditions. So far, the effects of high salinity on the most abundant marine photoautotrophic organism, Prochlorococcus, in marine oligotrophic environments, are mostly unknown. Here, we report the isolation of a new Prochlorococcus strain (RSP50) belonging to high-light (HL) clade II from the Red Sea, one of the warmest and most saline bodies of water in the global oceans. A comparative genomic analysis identified a set of 59 genes that were exclusive to RSP50 relative to currently available Prochlorococcus genomes, the majority of which (70%) encode for hypothetical proteins of unknown function. However, three of the unique genes encode for a complete pathway for the biosynthesis of the compatible solute glucosylglycerol, and are homologous to enzymes found in the sister lineage Synechococcus. Metatranscriptomic analyses of this metabolic pathway in the water column of the Red Sea revealed that the corresponding genes were constitutively transcribed, independent of depth and light, suggesting that osmoregulation using glucosylglycerol is a general feature of HL II Prochlorococcus in the Red Sea.

OneTwoTree: An online tool for phylogeny reconstruction

Phylogeny reconstruction is a key instrument in numerous biological analyses, ranging from evolutionary and ecology research, to conservation and systems biology. The increasing accumulation of genomic data makes it possible to reconstruct phylogenies with both high accuracy and at increasingly finer resolution. Yet, taking advantage of the enormous amount of sequence data available requires the use of computational tools for efficient data retrieval and processing, or else the process could quickly become an error-prone endeavour. Here, we present OneTwoTree (http://onetwotree.tau.ac.il/), a Web-based tool for tree reconstruction based on the supermatrix paradigm. Given a list of taxa names of interest as the sole input requirement, OneTwoTree retrieves all available sequence data from NCBI GenBank, clusters these into orthology groups, identifies the most informative set of markers, searches for an appropriate outgroup, and assembles a partitioned sequence matrix that is then used for the final phylogeny reconstruction step. OneTwoTree further allows users to control various steps of the process, such as the merging of sequences from similar clusters, or phylogeny reconstruction based on markers from a specific genome type. By comparing the performance of OneTwoTree to a manually reconstructed phylogeny of the Antirrhineae tribe, we show that the use of OneTwoTree resulted in substantially higher data coverage in terms of both taxon sampling and the number of informative markers assembled. OneTwoTree provides a flexible online tool for species-tree reconstruction, aimed to assist researchers ranging in their level of prior expertise in the task of phylogeny reconstruction.

Optofluidic Single-Cell Genome Amplification of Sub-micron Bacteria in the Ocean Subsurface

Optofluidic single-cell genome amplification was used to obtain genome sequences from sub-micron cells collected from the euphotic and mesopelagic zones of the northwestern Sargasso Sea. Plankton cells were visually selected and manually sorted with an optical trap, yielding 20 partial genome sequences representing seven bacterial phyla. Two organisms, E01-9C-26 (Gammaproteobacteria), represented by four single cell genomes, and Opi.OSU.00C, an uncharacterized Verrucomicrobia, were the first of their types retrieved by single cell genome sequencing and were studied in detail. Metagenomic data showed that E01-9C-26 is found throughout the dark ocean, while Opi.OSU.00C was observed to bloom transiently in the nutrient-depleted euphotic zone of the late spring and early summer. The E01-9C-26 genomes had an estimated size of 4.76-5.05 Mbps, and contained "O" and "W"-type monooxygenase genes related to methane and ammonium monooxygenases that were previously reported from ocean metagenomes. Metabolic reconstruction indicated E01-9C-26 are likely versatile methylotrophs capable of scavenging C1 compounds, methylated compounds, reduced sulfur compounds, and a wide range of amines, including D-amino acids. The genome sequences identified E01-9C-26 as a source of "O" and "W"-type monooxygenase genes related to methane and ammonium monooxygenases that were previously reported from ocean metagenomes, but are of unknown function. In contrast, Opi.OSU.00C genomes encode genes for catabolizing carbohydrate compounds normally associated with eukaryotic phytoplankton. This exploration of optofluidics showed that it was effective for retrieving diverse single-cell bacterioplankton genomes and has potential advantages in microbiology applications that require working with small sample volumes or targeting cells by their morphology.

Genome mining for peptidases in heat-tolerant and mesophilic fungi and putative adaptations for thermostability

BACKGROUND

Peptidases (EC 3.4) consist of a large group of hydrolytic enzymes that catalyze the hydrolysis of proteins accounting for approximately 65% of the total worldwide enzyme production. Peptidases from thermophilic fungi have adaptations to high temperature that makes them adequate for biotechnological application. In the present study, we profiled the genomes of heat-tolerant fungi and phylogenetically related mesophilic species for genes encoding for peptidases and their putative adaptations for thermostability.

RESULTS

We generated an extensive catalogue of these enzymes ranging from 241 to 820 peptidase genes in the genomes of 23 fungi. Thermophilic species presented the smallest number of peptidases encoding genes in relation to mesophilic species, and the peptidases families with a greater number of genes were the most affected. We observed differences in peptidases in thermophilic species in comparison to mesophilic counterparts, at (i) the genome level: a great reduction in the number of peptidases encoding genes that harbored a higher number of copies; (ii) in the primary protein structure: shifts in proportion of single or groups of amino acids; and (iii) in the three-dimensional structure: reduction in the number of internal cavities. Similar results were reported for extremely thermophilic proteins, but here we show for the first time that several changes also occurred on the moderate thermophilic enzymes of fungi. In regards to the amino acids composition, peptidases from thermophilic species in relation to the mesophilic ones, contained a larger proportion of Ala, Glu, Gly, Pro, Arg and Val residues and a lower number of Cys, His, Ile, Lys, Met, Asn, Gln, Ser, Thr and Trp residues (P < 0.05). Moreover, we observed an increase in the proportion of hydrophobic and charged amino acids and a decrease in polar amino acids.

CONCLUSIONS

Although thermophilic fungi present less genes encoding for peptidases, these have adaptations that could play a role in thermal resistance from genome to protein structure level.

SAR202 Genomes from the Dark Ocean Predict Pathways for the Oxidation of Recalcitrant Dissolved Organic Matter

Deep-ocean regions beyond the reach of sunlight contain an estimated 615 Pg of dissolved organic matter (DOM), much of which persists for thousands of years. It is thought that bacteria oxidize DOM until it is too dilute or refractory to support microbial activity. We analyzed five single-amplified genomes (SAGs) from the abundant SAR202 clade of dark-ocean bacterioplankton and found they encode multiple families of paralogous enzymes involved in carbon catabolism, including several families of oxidative enzymes that we hypothesize participate in the degradation of cyclic alkanes. The five partial genomes encoded 152 flavin mononucleotide/F420-dependent monooxygenases (FMNOs), many of which are predicted to be type II Baeyer-Villiger monooxygenases (BVMOs) that catalyze oxygen insertion into semilabile alicyclic alkanes. The large number of oxidative enzymes, as well as other families of enzymes that appear to play complementary roles in catabolic pathways, suggests that SAR202 might catalyze final steps in the biological oxidation of relatively recalcitrant organic compounds to refractory compounds that persist.IMPORTANCE Carbon in the ocean is massively sequestered in a complex mixture of biologically refractory molecules that accumulate as the chemical end member of biological oxidation and diagenetic change. However, few details are known about the biochemical machinery of carbon sequestration in the deep ocean. Reconstruction of the metabolism of a deep-ocean microbial clade, SAR202, led to postulation of new biochemical pathways that may be the penultimate stages of DOM oxidation to refractory forms that persist. These pathways are tied to a proliferation of oxidative enzymes. This research illuminates dark-ocean biochemistry that is broadly consequential for reconstructing the global carbon cycle.

De-novo transcriptome assembly for gene identification, analysis, annotation, and molecular marker discovery in Onobrychis viciifolia

BACKGROUND

Sainfoin (Onobrychis viciifolia) is a highly nutritious tannin-containing forage legume. In the diet of ruminants sainfoin can have anti-parasitic effects and reduce methane emissions under in vitro conditions. Many of these benefits have been attributed to condensed tannins or proanthocyanidins in sainfoin. A combination of increased use of industrially produced nitrogen fertilizer, issues with establishment and productivity in the first year and more reliable alternatives, such as red clover led to a decline in the use of sainfoin since the middle of the last century. In recent years there has been a resurgence of interest in sainfoin due to its potential beneficial nutraceutical and environmental attributes. However, genomic resources are scarce, thus hampering progress in genetic analysis and improvement. To address this we have used next generation RNA sequencing technology to obtain the first transcriptome of sainfoin. We used the library to identify gene-based simple sequence repeats (SSRs) and potential single nucleotide polymorphisms (SNPs).

RESULTS

One genotype from each of five sainfoin accessions was sequenced. Paired-end (PE) sequences were generated from cDNA libraries of RNA extracted from 7 day old seedlings. A combined assembly of 92,772 transcripts was produced de novo using the Trinity programme. About 18,000 transcripts were annotated with at least one GO (gene ontology) term. A total of 63 transcripts were annotated as involved in the tannin biosynthesis pathway. We identified 3786 potential SSRs. SNPs were identified by mapping the reads of the individual assemblies against the combined assembly. After stringent filtering a total of 77,000 putative SNPs were identified. A phylogenetic analysis of single copy number genes showed that sainfoin was most closely related to red clover and Medicago truncatula, while Lotus japonicus, bean and soybean are more distant relatives.

CONCLUSIONS

This work describes the first transcriptome assembly in sainfoin. The 92 K transcripts provide a rich source of SNP and SSR polymorphisms for future use in genetic studies of this crop. Annotation of genes involved in the condensed tannin biosynthesis pathway has provided the basis for further studies of the genetic control of this important trait in sainfoin.

Ectomycorrhizal ecology is imprinted in the genome of the dominant symbiotic fungus Cenococcum geophilum

The most frequently encountered symbiont on tree roots is the ascomycete Cenococcum geophilum, the only mycorrhizal species within the largest fungal class Dothideomycetes, a class known for devastating plant pathogens. Here we show that the symbiotic genomic idiosyncrasies of ectomycorrhizal basidiomycetes are also present in C. geophilum with symbiosis-induced, taxon-specific genes of unknown function and reduced numbers of plant cell wall-degrading enzymes. C. geophilum still holds a significant set of genes in categories known to be involved in pathogenesis and shows an increased genome size due to transposable elements proliferation. Transcript profiling revealed a striking upregulation of membrane transporters, including aquaporin water channels and sugar transporters, and mycorrhiza-induced small secreted proteins (MiSSPs) in ectomycorrhiza compared with free-living mycelium. The frequency with which this symbiont is found on tree roots and its possible role in water and nutrient transport in symbiosis calls for further studies on mechanisms of host and environmental adaptation.

Structural and functional analysis of the finished genome of the recently isolated toxic Anabaena sp. WA102

Background

Very few closed genomes of the cyanobacteria that commonly produce toxic blooms in lakes and reservoirs are available, limiting our understanding of the properties of these organisms. A new anatoxin-a-producing member of the Nostocaceae, Anabaena sp. WA102, was isolated from a freshwater lake in Washington State, USA, in 2013 and maintained in non-axenic culture.

Results

The Anabaena sp. WA102 5.7 Mbp genome assembly has been closed with long-read, single-molecule sequencing and separately a draft genome assembly has been produced with short-read sequencing technology. The closed and draft genome assemblies are compared, showing a correlation between long repeats in the genome and the many gaps in the short-read assembly. Anabaena sp. WA102 encodes anatoxin-a biosynthetic genes, as does its close relative Anabaena sp. AL93 (also introduced in this study). These strains are distinguished by differences in the genes for light-harvesting phycobilins, with Anabaena sp. AL93 possessing a phycoerythrocyanin operon. Biologically relevant structural variants in the Anabaena sp. WA102 genome were detected only by long-read sequencing: a tandem triplication of the anaBCD promoter region in the anatoxin-a synthase gene cluster (not triplicated in Anabaena sp. AL93) and a 5-kbp deletion variant present in two-thirds of the population. The genome has a large number of mobile elements (160). Strikingly, there was no synteny with the genome of its nearest fully assembled relative, Anabaena sp. 90.

Conclusion

Structural and functional genome analyses indicate that Anabaena sp. WA102 has a flexible genome. Genome closure, which can be readily achieved with long-read sequencing, reveals large scale (e.g., gene order) and local structural features that should be considered in understanding genome evolution and function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2738-7) contains supplementary material, which is available to authorized users.

A molecular phylogenetic framework for Bacillus subtilis using genome sequences and its application to Bacillus subtilis subspecies stecoris strain D7XPN1, an isolate from a commercial food-waste degrading bioreactor

A thermophilic, heterotrophic and facultatively anaerobic bacterium designated strain D7XPN1 was isolated from Baku BakuKing™, a commercial food-waste degrading bioreactor (composter). The strain grew optimally at 45 °C (growth range between 24 and 50 °C) and pH 7 (growth pH range between pH 5 and 9) in Luria Broth supplemented with 0.3 % glucose. Strain D7XPN1 tolerated up to 7 % NaCl and showed amylolytic and xylanolytic activities. 16S rRNA gene analysis placed strain D7XPN1 in the cluster represented by Bacillus subtilis and the genome analysis of the 4.1 Mb genome sequence determined using RAST (Rapid Annotation using Subsystem Technology) indicated a total of 5116 genomic features were present of which 2320 features could be grouped into several subsystem categories. Of these, 615 features were related to carbohydrate metabolism which included a range of enzymes with potential in the biodegradation of food wastes, a property consistent with the ecological habitat of the isolate. ANIb (Average Nucleotide Identity based on BLAST) analysis with 49 Bacillus subtilis genomes indicated that it was distantly related to the three currently taxonomically validated B. subtilis subspecies namely B. subtilis subsp. subtilis (95.6 %), B. subtilis subsp. spizizenii (93 %) and B. subtilis subsp. inaquosorum (92 %) and based on our current knowledge warranted that it be included as a separate cluster together with strain JS which it was closely related (98.69 %). The close relationship of strains D7XPN1 and JS is also supported from our results from electronic DNA–DNA Hybridization (e-DDH) studies. Furthermore, our additional in-depth phylogenomic analyses using three different datasets unequivocally supported the creation of a fourth B. subtilis subspecies to include strains D7XPN1 and JS for which we propose strain D7XPN1^T (=KCTC 33554^T, JCM 30051^T) as the type strain, and designate it as B. subtilis subsp. stecoris.

Genome Sequencing and Transposon Mutagenesis of Burkholderia seminalis TC3.4.2R3 Identify Genes Contributing to Suppression of Orchid Necrosis Caused by B. gladioli

From a screen of 36 plant-associated strains of Burkholderia spp., we identified 24 strains that suppressed leaf and pseudobulb necrosis of orchid caused by B. gladioli. To gain insights into the mechanisms of disease suppression, we generated a draft genome sequence from one suppressive strain, TC3.4.2R3. The genome is an estimated 7.67 megabases in size, with three replicons, two chromosomes, and the plasmid pC3. Using a combination of multilocus sequence analysis and phylogenomics, we identified TC3.4.2R3 as B. seminalis, a species within the Burkholderia cepacia complex that includes opportunistic human pathogens and environmental strains. We generated and screened a library of 3,840 transposon mutants of strain TC3.4.2R3 on orchid leaves to identify genes contributing to plant disease suppression. Twelve mutants deficient in suppression of leaf necrosis were selected and the transposon insertions were mapped to eight loci. One gene is in a wcb cluster that is related to synthesis of extracellular polysaccharide, a key determinant in bacterial-host interactions in other systems, and the other seven are highly conserved among Burkholderia spp. The fundamental information developed in this study will serve as a resource for future research aiming to identify mechanisms contributing to biological control.

The Complete Sequence and Comparative Analysis of a Multidrug-Resistance and Virulence Multireplicon IncFII Plasmid pEC302/04 from an Extraintestinal Pathogenic Escherichia coli EC302/04 Indicate Extensive Diversity of IncFII Plasmids

Extraintestinal pathogenic Escherichia coli (ExPEC) that causes extraintestinal infections often harbor plasmids encoding fitness traits such as resistance and virulence determinants that are of clinical importance. We determined the complete nucleotide sequence of plasmid pEC302/04 from a multidrug-resistant E. coli EC302/04 which was isolated from the tracheal aspirate of a patient in Malaysia. In addition, we also performed comparative sequence analyses of 18 related IncFIIA plasmids to determine the phylogenetic relationship and diversity of these plasmids. The 140,232 bp pEC302/04 is a multireplicon plasmid that bears three replication systems (FII, FIA, and FIB) with subtype of F2:A1:B1. The plasmid is self-transmissible with a complete transfer region. pEC302/04 also carries antibiotic resistance genes such as bla_TEM−1 and a class I integron containing sul1, cml and aadA resistance genes, conferring multidrug resistance (MDR) to its host, E. coli EC302/04. Besides, two iron acquisition systems (SitABCD and IutA-IucABCD) which are the conserved virulence determinants of ExPEC-colicin V or B and M (ColV/ColBM)-producing plasmids were identified in pEC302/04. Multiple toxin-antitoxin (TA)-based addiction systems (i.e., PemI/PemK, VagC/VagD, CcdA/CcdB, and Hok/Sok) and a plasmid partitioning system, ParAB, and PsiAB, which are important for plasmid maintenance were also found. Comparative plasmid analysis revealed only one conserved gene, the repA1 as the core genome, showing that there is an extensive diversity among the IncFIIA plasmids. The phylogenetic relationship of 18 IncF plasmids based on the core regions revealed that ColV/ColBM-plasmids and non-ColV/ColBM plasmids were separated into two distinct groups. These plasmids, which carry highly diverse genetic contents, are also mosaic in nature. The atypical combination of genetic materials, i.e., the MDR- and ColV/ColBM-plasmid-virulence encoding regions in a single ExPEC plasmid is rare but of clinical importance. Such phenomenon is bothersome when the plasmids are transmissible, facilitating the spread of virulence and resistance plasmids among pathogenic bacteria. Notably, certain TA systems are more commonly found in particular ExPEC plasmid types, indicating the possible relationships between certain TA systems and ExPEC pathogenesis.

Comprehensive Genomic Analyses of the OM43 Clade, Including a Novel Species from the Red Sea, Indicate Ecotype Differentiation among Marine Methylotrophs

The OM43 clade within the family Methylophilaceae of Betaproteobacteria represents a group of methylotrophs that play important roles in the metabolism of C1 compounds in marine environments and other aquatic environments around the globe. Using dilution-to-extinction cultivation techniques, we successfully isolated a novel species of this clade (here designated MBRS-H7) from the ultraoligotrophic open ocean waters of the central Red Sea. Phylogenomic analyses indicate that MBRS-H7 is a novel species that forms a distinct cluster together with isolate KB13 from Hawaii (Hawaii-Red Sea [H-RS] cluster) that is separate from the cluster represented by strain HTCC2181 (from the Oregon coast). Phylogenetic analyses using the robust 16S-23S internal transcribed spacer revealed a potential ecotype separation of the marine OM43 clade members, which was further confirmed by metagenomic fragment recruitment analyses that showed trends of higher abundance in low-chlorophyll and/or high-temperature provinces for the H-RS cluster but a preference for colder, highly productive waters for the HTCC2181 cluster. This potential environmentally driven niche differentiation is also reflected in the metabolic gene inventories, which in the case of the H-RS cluster include those conferring resistance to high levels of UV irradiation, temperature, and salinity. Interestingly, we also found different energy conservation modules between these OM43 subclades, namely, the existence of the NADH:quinone oxidoreductase complex I (NUO) system in the H-RS cluster and the nonhomologous NADH:quinone oxidoreductase (NQR) system in the HTCC2181 cluster, which might have implications for their overall energetic yields.

Red clover (Trifolium pratense L.) draft genome provides a platform for trait improvement

Red clover (Trifolium pratense L.) is a globally significant forage legume in pastoral livestock farming systems. It is an attractive component of grassland farming, because of its high yield and protein content, nutritional value and ability to fix atmospheric nitrogen. Enhancing its role further in sustainable agriculture requires genetic improvement of persistency, disease resistance, and tolerance to grazing. To help address these challenges, we have assembled a chromosome-scale reference genome for red clover. We observed large blocks of conserved synteny with Medicago truncatula and estimated that the two species diverged ~23 million years ago. Among the 40,868 annotated genes, we identified gene clusters involved in biochemical pathways of importance for forage quality and livestock nutrition. Genotyping by sequencing of a synthetic population of 86 genotypes show that the number of markers required for genomics-based breeding approaches is tractable, making red clover a suitable candidate for association studies and genomic selection.

Whole genome annotation and comparative genomic analyses of bio-control fungus Purpureocillium lilacinum

Background

The fungus Purpureocillium lilacinum is widely known as a biological control agent against plant parasitic nematodes. This research article consists of genomic annotation of the first draft of whole genome sequence of P. lilacinum. The study aims to decipher the putative genetic components of the fungus involved in nematode pathogenesis by performing comparative genomic analysis with nine closely related fungal species in Hypocreales.

Results

de novo genomic assembly was done and a total of 301 scaffolds were constructed for P. lilacinum genomic DNA. By employing structural genome prediction models, 13, 266 genes coding for proteins were predicted in the genome. Approximately 73 % of the predicted genes were functionally annotated using Blastp, InterProScan and Gene Ontology. A 14.7 % fraction of the predicted genes shared significant homology with genes in the Pathogen Host Interactions (PHI) database. The phylogenomic analysis carried out using maximum likelihood RAxML algorithm provided insight into the evolutionary relationship of P. lilacinum. In congruence with other closely related species in the Hypocreales namely, Metarhizium spp., Pochonia chlamydosporia, Cordyceps militaris, Trichoderma reesei and Fusarium spp., P. lilacinum has large gene sets coding for G-protein coupled receptors (GPCRs), proteases, glycoside hydrolases and carbohydrate esterases that are required for degradation of nematode-egg shell components. Screening of the genome by Antibiotics & Secondary Metabolite Analysis Shell (AntiSMASH) pipeline indicated that the genome potentially codes for a variety of secondary metabolites, possibly required for adaptation to heterogeneous lifestyles reported for P. lilacinum. Significant up-regulation of subtilisin-like serine protease genes in presence of nematode eggs in quantitative real-time analyses suggested potential role of serine proteases in nematode pathogenesis.

Conclusions

The data offer a better understanding of Purpureocillium lilacinum genome and will enhance our understanding on the molecular mechanism involved in nematophagy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2229-2) contains supplementary material, which is available to authorized users.

Transcriptome Analysis of the Emerald Ash Borer (EAB), Agrilus planipennis: De Novo Assembly, Functional Annotation and Comparative Analysis

Background

The Emerald ash borer (EAB), Agrilus planipennis, is an invasive phloem-feeding insect pest of ash trees. Since its initial discovery near the Detroit, US- Windsor, Canada area in 2002, the spread of EAB has had strong negative economic, social and environmental impacts in both countries. Several transcriptomes from specific tissues including midgut, fat body and antenna have recently been generated. However, the relatively low sequence depth, gene coverage and completeness limited the usefulness of these EAB databases.

Methodology and Principal Findings

High-throughput deep RNA-Sequencing (RNA-Seq) was used to obtain 473.9 million pairs of 100 bp length paired-end reads from various life stages and tissues. These reads were assembled into 88,907 contigs using the Trinity strategy and integrated into 38,160 unigenes after redundant sequences were removed. We annotated 11,229 unigenes by searching against the public nr, Swiss-Prot and COG. The EAB transcriptome assembly was compared with 13 other sequenced insect species, resulting in the prediction of 536 unigenes that are Coleoptera-specific. Differential gene expression revealed that 290 unigenes are expressed during larval molting and 3,911 unigenes during metamorphosis from larvae to pupae, respectively (FDR< 0.01 and log2 FC>2). In addition, 1,167 differentially expressed unigenes were identified from larval and adult midguts, 435 unigenes were up-regulated in larval midgut and 732 unigenes were up-regulated in adult midgut. Most of the genes involved in RNA interference (RNAi) pathways were identified, which implies the existence of a system RNAi in EAB.

Conclusions and Significance

This study provides one of the most fundamental and comprehensive transcriptome resources available for EAB to date. Identification of the tissue- stage- or species- specific unigenes will benefit the further study of gene functions during growth and metamorphosis processes in EAB and other pest insects.

The Genome of Winter Moth (Operophtera brumata) Provides a Genomic Perspective on Sexual Dimorphism and Phenology

The winter moth (Operophtera brumata) belongs to one of the most species-rich families in Lepidoptera, the Geometridae (approximately 23,000 species). This family is of great economic importance as most species are herbivorous and capable of defoliating trees. Genome assembly of the winter moth allows the study of genes and gene families, such as the cytochrome P450 gene family, which is known to be vital in plant secondary metabolite detoxification and host-plant selection. It also enables exploration of the genomic basis for female brachyptery (wing reduction), a feature of sexual dimorphism in winter moth, and for seasonal timing, a trait extensively studied in this species. Here we present a reference genome for the winter moth, the first geometrid and largest sequenced Lepidopteran genome to date (638 Mb) including a set of 16,912 predicted protein-coding genes. This allowed us to assess the dynamics of evolution on a genome-wide scale using the P450 gene family. We also identified an expanded gene family potentially linked to female brachyptery, and annotated the genes involved in the circadian clock mechanism as main candidates for involvement in seasonal timing. The genome will contribute to Lepidopteran genomic resources and comparative genomics. In addition, the genome enhances our ability to understand the genetic and molecular basis of insect seasonal timing and thereby provides a reference for future evolutionary and population studies on the winter moth.

The genome of the truffle-parasite Tolypocladium ophioglossoides and the evolution of antifungal peptaibiotics

Background

Two major mycoparasitic lineages, the family Hypocreaceae and the genus Tolypocladium, exist within the fungal order, Hypocreales. Peptaibiotics are a group of secondary metabolites almost exclusively described from Trichoderma species of Hypocreaceae. Peptaibiotics are produced by nonribosomal peptide synthetases (NRPSs) and have antibiotic and antifungal activities. Tolypocladium species are mainly truffle parasites, but a few species are insect pathogens.

Results

The draft genome sequence of the truffle parasite Tolypocladium ophioglossoides was generated and numerous secondary metabolite clusters were discovered, many of which have no known putative product. However, three large peptaibiotic gene clusters were identified using phylogenetic analyses. Peptaibiotic genes are absent from the predominantly plant and insect pathogenic lineages of Hypocreales, and are therefore exclusive to the largely mycoparasitic lineages. Using NRPS adenylation domain phylogenies and reconciliation of the domain tree with the organismal phylogeny, it is demonstrated that the distribution of these domains is likely not the product of horizontal gene transfer between mycoparasitic lineages, but represents independent losses in insect pathogenic lineages. Peptaibiotic genes are less conserved between species of Tolypocladium and are the product of complex patterns of lineage sorting and module duplication. In contrast, these genes are more conserved within the genus Trichoderma and consistent with diversification through speciation.

Conclusions

Peptaibiotic NRPS genes are restricted to mycoparasitic lineages of Hypocreales, based on current sampling. Phylogenomics and comparative genomics can provide insights into the evolution of secondary metabolite genes, their distribution across a broader range of taxa, and their possible function related to host specificity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1777-9) contains supplementary material, which is available to authorized users.

SAR11 lipid renovation in response to phosphate starvation

Phytoplankton inhabiting oligotrophic ocean gyres actively reduce their phosphorus demand by replacing polar membrane phospholipids with those lacking phosphorus. Although the synthesis of nonphosphorus lipids is well documented in some heterotrophic bacterial lineages, phosphorus-free lipid synthesis in oligotrophic marine chemoheterotrophs has not been directly demonstrated, implying they are disadvantaged in phosphate-deplete ecosystems, relative to phytoplankton. Here, we show the SAR11 clade chemoheterotroph Pelagibacter sp. str. HTCC7211 renovates membrane lipids when phosphate starved by replacing a portion of its phospholipids with monoglucosyl- and glucuronosyl-diacylglycerols and by synthesizing new ornithine lipids. Lipid profiles of cells grown with excess phosphate consisted entirely of phospholipids. Conversely, up to 40% of the total lipids were converted to nonphosphorus lipids when cells were starved for phosphate, or when growing on methylphosphonate. Cells sequentially limited by phosphate and methylphosphonate transformed >75% of their lipids to phosphorus-free analogs. During phosphate starvation, a four-gene cluster was significantly up-regulated that likely encodes the enzymes responsible for lipid renovation. These genes were found in Pelagibacterales strains isolated from a phosphate-deficient ocean gyre, but not in other strains from coastal environments, suggesting alternate lipid synthesis is a specific adaptation to phosphate scarcity. Similar gene clusters are found in the genomes of other marine α-proteobacteria, implying lipid renovation is a common strategy used by heterotrophic cells to reduce their requirement for phosphorus in oligotrophic habitats.

Phylogenomic Analyses Indicate that Early Fungi Evolved Digesting Cell Walls of Algal Ancestors of Land Plants

As decomposers, fungi are key players in recycling plant material in global carbon cycles. We hypothesized that genomes of early diverging fungi may have inherited pectinases from an ancestral species that had been able to extract nutrients from pectin-containing land plants and their algal allies (Streptophytes). We aimed to infer, based on pectinase gene expansions and on the organismal phylogeny, the geological timing of the plant-fungus association. We analyzed 40 fungal genomes, three of which, including Gonapodya prolifera, were sequenced for this study. In the organismal phylogeny from 136 housekeeping loci, Rozella diverged first from all other fungi. Gonapodya prolifera was included among the flagellated, predominantly aquatic fungal species in Chytridiomycota. Sister to Chytridiomycota were the predominantly terrestrial fungi including zygomycota I and zygomycota II, along with the ascomycetes and basidiomycetes that comprise Dikarya. The Gonapodya genome has 27 genes representing five of the seven classes of pectin-specific enzymes known from fungi. Most of these share a common ancestry with pectinases from Dikarya. Indicating functional and sequence similarity, Gonapodya, like many Dikarya, can use pectin as a carbon source for growth in pure culture. Shared pectinases of Dikarya and Gonapodya provide evidence that even ancient aquatic fungi had adapted to extract nutrients from the plants in the green lineage. This implies that 750 million years, the estimated maximum age of origin of the pectin-containing streptophytes represents a maximum age for the divergence of Chytridiomycota from the lineage including Dikarya.

BIR Pipeline for Preparation of Phylogenomic Data

SUMMARY

We present a pipeline named BIR (Blast, Identify and Realign) developed for phylogenomic analyses. BIR is intended for the identification of gene sequences applicable for phylogenomic inference. The pipeline allows users to apply their own manually curated sequence alignments (seed) in search for homologous genes in sequence databases and available genomes. BIR automatically adds the identified sequences from these databases to the seed alignments and reconstruct a phylogenetic tree from each. The BIR pipeline is an efficient tool for the identification of orthologous gene copies because it expands user-defined sequence alignments and conducts massive parallel phylogenetic reconstruction. The application is also particularly useful for large-scale sequencing projects that require management of a large number of single-gene alignments for gene comparison, functional annotation, and evolutionary analyses.

AVAILABILITY

The BIR user manual is available at http://www.bioportal.no/ and can be accessed through Lifeportal at https://lifeportal.uio.no. Access is free but requires a user account registration using the link "Register for BIR access" from the Lifeportal homepage.

Metagenome sequence of Elaphomyces granulatus from sporocarp tissue reveals Ascomycota ectomycorrhizal fingerprints of genome expansion and a Proteobacteria-rich microbiome

Many obligate symbiotic fungi are difficult to maintain in culture, and there is a growing need for alternative approaches to obtaining tissue and subsequent genomic assemblies from such species. In this study, the genome of Elaphomyces granulatus was sequenced from sporocarp tissue. The genome assembly remains on many contigs, but gene space is estimated to be mostly complete. Phylogenetic analyses revealed that the Elaphomyces lineage is most closely related to Talaromyces and Trichocomaceae s.s. The genome of E. granulatus is reduced in carbohydrate-active enzymes, despite a large expansion in genome size, both of which are consistent with what is seen in Tuber melanosporum, the other sequenced ectomycorrhizal ascomycete. A large number of transposable elements are predicted in the E. granulatus genome, especially Gypsy-like long terminal repeats, and there has also been an expansion in helicases. The metagenome is a complex community dominated by bacteria in Bradyrhizobiaceae, and there is evidence to suggest that the community may be reduced in functional capacity as estimated by KEGG pathways. Through the sequencing of sporocarp tissue, this study has provided insights into Elaphomyces phylogenetics, genomics, metagenomics and the evolution of the ectomycorrhizal association.

Synechococcus: 3 billion years of global dominance

Cyanobacteria are among the most important primary producers on the Earth. However, the evolutionary forces driving cyanobacterial species diversity remain largely enigmatic due to both their distinction from macro-organisms and an undersampling of sequenced genomes. Thus, we present a new genome of a Synechococcus-like cyanobacterium from a novel evolutionary lineage. Further, we analyse all existing 16S rRNA sequences and genomes of Synechococcus-like cyanobacteria. Chronograms showed extremely polyphyletic relationships in Synechococcus, which has not been observed in any other cyanobacteria. Moreover, most Synechococcus lineages bifurcated after the Great Oxidation Event, including the most abundant marine picoplankton lineage. Quantification of horizontal gene transfer among 70 cyanobacterial genomes revealed significant differences among studied genomes. Horizontal gene transfer levels were not correlated with ecology, genome size or phenotype, but were correlated with the age of divergence. All findings were synthetized into a novel model of cyanobacterial evolution, characterized by serial convergence of the features, that is multicellularity and ecology.

Use of whole genome sequences to develop a molecular phylogenetic framework for Rhodococcus fascians and the Rhodococcus genus

The accurate diagnosis of diseases caused by pathogenic bacteria requires a stable species classification. Rhodococcus fascians is the only documented member of its ill-defined genus that is capable of causing disease on a wide range of agriculturally important plants. Comparisons of genome sequences generated from isolates of Rhodococcus associated with diseased plants revealed a level of genetic diversity consistent with them representing multiple species. To test this, we generated a tree based on more than 1700 homologous sequences from plant-associated isolates of Rhodococcus, and obtained support from additional approaches that measure and cluster based on genome similarities. Results were consistent in supporting the definition of new Rhodococcus species within clades containing phytopathogenic members. We also used the genome sequences, along with other rhodococcal genome sequences to construct a molecular phylogenetic tree as a framework for resolving the Rhodococcus genus. Results indicated that Rhodococcus has the potential for having 20 species and also confirmed a need to revisit the taxonomic groupings within Rhodococcus.

Genome sequencing of four Aureobasidium pullulans varieties: biotechnological potential, stress tolerance, and description of new species

BACKGROUND

Aureobasidium pullulans is a black-yeast-like fungus used for production of the polysaccharide pullulan and the antimycotic aureobasidin A, and as a biocontrol agent in agriculture. It can cause opportunistic human infections, and it inhabits various extreme environments. To promote the understanding of these traits, we performed de-novo genome sequencing of the four varieties of A. pullulans.

RESULTS

The 25.43-29.62 Mb genomes of these four varieties of A. pullulans encode between 10266 and 11866 predicted proteins. Their genomes encode most of the enzyme families involved in degradation of plant material and many sugar transporters, and they have genes possibly associated with degradation of plastic and aromatic compounds. Proteins believed to be involved in the synthesis of pullulan and siderophores, but not of aureobasidin A, are predicted. Putative stress-tolerance genes include several aquaporins and aquaglyceroporins, large numbers of alkali-metal cation transporters, genes for the synthesis of compatible solutes and melanin, all of the components of the high-osmolarity glycerol pathway, and bacteriorhodopsin-like proteins. All of these genomes contain a homothallic mating-type locus.

CONCLUSIONS

The differences between these four varieties of A. pullulans are large enough to justify their redefinition as separate species: A. pullulans, A. melanogenum, A. subglaciale and A. namibiae. The redundancy observed in several gene families can be linked to the nutritional versatility of these species and their particular stress tolerance. The availability of the genome sequences of the four Aureobasidium species should improve their biotechnological exploitation and promote our understanding of their stress-tolerance mechanisms, diverse lifestyles, and pathogenic potential.

Genomic differentiation among two strains of the PS1 clade isolated from geographically separated marine habitats

Using dilution-to-extinction cultivation, we isolated a strain affiliated with the PS1 clade from surface waters of the Red Sea. Strain RS24 represents the second isolate of this group of marine Alphaproteobacteria after IMCC14465 that was isolated from the East (Japan) Sea. The PS1 clade is a sister group to the OCS116 clade, together forming a putatively novel order closely related to Rhizobiales. While most genomic features and most of the genetic content are conserved between RS24 and IMCC14465, their average nucleotide identity (ANI) is < 81%, suggesting two distinct species of the PS1 clade. Next to encoding two different variants of proteorhodopsin genes, they also harbor several unique genomic islands that contain genes related to degradation of aromatic compounds in IMCC14465 and in polymer degradation in RS24, possibly reflecting the physicochemical differences in the environment they were isolated from. No clear differences in abundance of the genomic content of either strain could be found in fragment recruitment analyses using different metagenomic datasets, in which both genomes were detectable albeit as minor part of the communities. The comparative genomic analysis of both isolates of the PS1 clade and the fragment recruitment analysis provide first insights into the ecology of this group.

Mitochondrial genome of Phlebia radiata is the second largest (156 kbp) among fungi and features signs of genome flexibility and recent recombination events

Mitochondria are eukaryotic organelles supporting individual life-style via generation of proton motive force and cellular energy, and indispensable metabolic pathways. As part of genome sequencing of the white rot Basidiomycota species Phlebia radiata, we first assembled its mitochondrial genome (mtDNA). So far, the 156 348 bp mtDNA is the second largest described for fungi, and of considerable size among eukaryotes. The P. radiata mtDNA assembled as single circular dsDNA molecule containing genes for the large and small ribosomal RNAs, 28 transfer RNAs, and over 100 open reading frames encoding the 14 fungal conserved protein subunits of the mitochondrial complexes I, III, IV, and V. Two genes (atp6 and tRNA-IleGAU) were duplicated within 6.1 kbp inverted region, which is a unique feature of the genome. The large mtDNA size, however, is explained by the dominance of intronic and intergenic regions (sum 80% of mtDNA sequence). The intergenic DNA stretches harness short (≤ 200 nt) repetitive, dispersed and overlapping sequence elements in abundance. Long self-splicing introns of types I and II interrupt eleven of the conserved genes (cox1,2,3; cob; nad1,2,4,4L,5; rnl; rns). The introns embrace a total of 57 homing endonucleases with LAGLIDADGD and GYI-YIG core motifs, which makes P. radiata mtDNA to one of the largest known reservoirs of intron-homing endonucleases. The inverted duplication, intergenic stretches, and intronic features are indications of dynamics and genetic flexibility of the mtDNA, not fully recognized to this extent in fungal mitochondrial genomes previously, thus giving new insights for the evolution of organelle genomes in eukaryotes.

Discovery of a SAR11 growth requirement for thiamin's pyrimidine precursor and its distribution in the Sargasso Sea

Vitamin traffic, the production of organic growth factors by some microbial community members and their use by other taxa, is being scrutinized as a potential explanation for the variation and highly connected behavior observed in ocean plankton by community network analysis. Thiamin (vitamin B1), a cofactor in many essential biochemical reactions that modify carbon-carbon bonds of organic compounds, is distributed in complex patterns at subpicomolar concentrations in the marine surface layer (0-300 m). Sequenced genomes from organisms belonging to the abundant and ubiquitous SAR11 clade of marine chemoheterotrophic bacteria contain genes coding for a complete thiamin biosynthetic pathway, except for thiC, encoding the 4-amino-5-hydroxymethyl-2-methylpyrimidine (HMP) synthase, which is required for de novo synthesis of thiamin's pyrimidine moiety. Here we demonstrate that the SAR11 isolate 'Candidatus Pelagibacter ubique', strain HTCC1062, is auxotrophic for the thiamin precursor HMP, and cannot use exogenous thiamin for growth. In culture, strain HTCC1062 required 0.7 zeptomoles per cell (ca. 400 HMP molecules per cell). Measurements of dissolved HMP in the Sargasso Sea surface layer showed that HMP ranged from undetectable (detection limit: 2.4 pM) to 35.7 pM, with maximum concentrations coincident with the deep chlorophyll maximum. In culture, some marine cyanobacteria, microalgae and bacteria exuded HMP, and in the Western Sargasso Sea, HMP profiles changed between the morning and evening, suggesting a dynamic biological flux from producers to consumers.

Building a phylogenomic pipeline for the eukaryotic tree of life - addressing deep phylogenies with genome-scale data

Background Understanding the evolutionary relationships of all eukaryotes on Earth remains a paramount goal of modern biology, yet analyzing homologous sequences across 1.8 billion years of eukaryotic evolution is challenging. Many existing tools for identifying gene orthologs are inadequate when working with heterogeneous rates of evolution and endosymbiotic/lateral gene transfer. Moreover, genomic-scale sequencing, which was once the domain of large sequencing centers, has advanced to the point where small laboratories can now generate the data needed for phylogenomic studies. This has opened the door for increased taxonomic sampling as individual research groups have the ability to conduct genome-scale projects on their favorite non-model organism. Results Here we present some of the tools developed, and insights gained, as we created a pipeline that combines data-mining from public databases and our own transcriptome data to study the eukaryotic tree of life. The first steps of a phylogenomic pipeline involve choosing taxa and loci, and making decisions about how to handle alleles, paralogs and non-overlapping sequences. Next, orthologs are aligned for analyses including gene tree reconstruction and concatenation for supermatrix approaches. To build our pipeline, we created scripts written in Python that integrate third-party tools with custom methods. As a test case, we present the placement of five amoebae on the eukaryotic tree of life based on analyses of transcriptome data. Our scripts available on GitHUb and may be used as-is for automated analyses of large scale phylogenomics, or adapted for use in other types of studies. Conclusion Analyses on the scale of all eukaryotes present challenges not necessarily found in studies of more closely related organisms. Our approach will be of relevance to others for whom existing third-party tools fail to fully answer desired phylogenetic questions.

Genome sequencing provides insight into the reproductive biology, nutritional mode and ploidy of the fern pathogen Mixia osmundae

Mixia osmundae (Basidiomycota, Pucciniomycotina) represents a monotypic class containing an unusual fern pathogen with incompletely understood biology. We sequenced and analyzed the genome of M. osmundae, focusing on genes that may provide some insight into its mode of pathogenicity and reproductive biology. Mixia osmundae has the smallest plant pathogenic basidiomycete genome sequenced to date, at 13.6 Mb, with very few repeats, high gene density, and relatively few significant gene family gains. The genome shows that the yeast state of M. osmundae is haploid and the lack of segregation of mating genes suggests that the spores produced on Osmunda spp. fronds are probably asexual. However, our finding of a complete complement of mating and meiosis genes suggests the capacity to undergo sexual reproduction. Analyses of carbohydrate active enzymes suggest that this fungus is a biotroph with the ability to break down several plant cell wall components. Analyses of publicly available sequence data show that other Mixia members may exist on other plant hosts and with a broader distribution than previously known.

Reconciliation of gene and species trees

The first part of the paper briefly overviews the problem of gene and species trees reconciliation with the focus on defining and algorithmic construction of the evolutionary scenario. Basic ideas are discussed for the aspects of mapping definitions, costs of the mapping and evolutionary scenario, imposing time scales on a scenario, incorporating horizontal gene transfers, binarization and reconciliation of polytomous trees, and construction of species trees and scenarios. The review does not intend to cover the vast diversity of literature published on these subjects. Instead, the authors strived to overview the problem of the evolutionary scenario as a central concept in many areas of evolutionary research. The second part provides detailed mathematical proofs for the solutions of two problems: (i) inferring a gene evolution along a species tree accounting for various types of evolutionary events and (ii) trees reconciliation into a single species tree when only gene duplications and losses are allowed. All proposed algorithms have a cubic time complexity and are mathematically proved to find exact solutions. Solving algorithms for problem (ii) can be naturally extended to incorporate horizontal transfers, other evolutionary events, and time scales on the species tree.

Single-cell enabled comparative genomics of a deep ocean SAR11 bathytype

Bacterioplankton of the SAR11 clade are the most abundant microorganisms in marine systems, usually representing 25% or more of the total bacterial cells in seawater worldwide. SAR11 is divided into subclades with distinct spatiotemporal distributions (ecotypes), some of which appear to be specific to deep water. Here we examine the genomic basis for deep ocean distribution of one SAR11 bathytype (depth-specific ecotype), subclade Ic. Four single-cell Ic genomes, with estimated completeness of 55%-86%, were isolated from 770 m at station ALOHA and compared with eight SAR11 surface genomes and metagenomic datasets. Subclade Ic genomes dominated metagenomic fragment recruitment below the euphotic zone. They had similar COG distributions, high local synteny and shared a large number (69%) of orthologous clusters with SAR11 surface genomes, yet were distinct at the 16S rRNA gene and amino-acid level, and formed a separate, monophyletic group in phylogenetic trees. Subclade Ic genomes were enriched in genes associated with membrane/cell wall/envelope biosynthesis and showed evidence of unique phage defenses. The majority of subclade Ic-specfic genes were hypothetical, and some were highly abundant in deep ocean metagenomic data, potentially masking mechanisms for niche differentiation. However, the evidence suggests these organisms have a similar metabolism to their surface counterparts, and that subclade Ic adaptations to the deep ocean do not involve large variations in gene content, but rather more subtle differences previously observed deep ocean genomic data, like preferential amino-acid substitutions, larger coding regions among SAR11 clade orthologs, larger intergenic regions and larger estimated average genome size.

Phylogenomic test of the hypotheses for the evolutionary origin of eukaryotes

The evolutionary origin of eukaryotes is a question of great interest for which many different hypotheses have been proposed. These hypotheses predict distinct patterns of evolutionary relationships for individual genes of the ancestral eukaryotic genome. The availability of numerous completely sequenced genomes covering the three domains of life makes it possible to contrast these predictions with empirical data. We performed a systematic analysis of the phylogenetic relationships of ancestral eukaryotic genes with archaeal and bacterial genes. In contrast with previous studies, we emphasize the critical importance of methods accounting for statistical support, horizontal gene transfer, and gene loss, and we disentangle the processes underlying the phylogenomic pattern we observe. We first recover a clear signal indicating that a fraction of the bacteria-like eukaryotic genes are of alphaproteobacterial origin. Then, we show that the majority of bacteria-related eukaryotic genes actually do not point to a relationship with a specific bacterial taxonomic group. We also provide evidence that eukaryotes branch close to the last archaeal common ancestor. Our results demonstrate that there is no phylogenetic support for hypotheses involving a fusion with a bacterium other than the ancestor of mitochondria. Overall, they leave only two possible interpretations, respectively, based on the early-mitochondria hypotheses, which suppose an early endosymbiosis of an alphaproteobacterium in an archaeal host and on the slow-drip autogenous hypothesis, in which early eukaryotic ancestors were particularly prone to horizontal gene transfers.

Proteomic and transcriptomic analyses of "Candidatus Pelagibacter ubique" describe the first PII-independent response to nitrogen limitation in a free-living Alphaproteobacterium

Nitrogen is one of the major nutrients limiting microbial productivity in the ocean, and as a result, most marine microorganisms have evolved systems for responding to nitrogen stress. The highly abundant alphaproteobacterium "Candidatus Pelagibacter ubique," a cultured member of the order Pelagibacterales (SAR11), lacks the canonical GlnB, GlnD, GlnK, and NtrB/NtrC genes for regulating nitrogen assimilation, raising questions about how these organisms respond to nitrogen limitation. A survey of 266 Alphaproteobacteria genomes found these five regulatory genes nearly universally conserved, absent only in intracellular parasites and members of the order Pelagibacterales, including "Ca. Pelagibacter ubique." Global differences in mRNA and protein expression between nitrogen-limited and nitrogen-replete cultures were measured to identify nitrogen stress responses in "Ca. Pelagibacter ubique" strain HTCC1062. Transporters for ammonium (AmtB), taurine (TauA), amino acids (YhdW), and opines (OccT) were all elevated in nitrogen-limited cells, indicating that they devote increased resources to the assimilation of nitrogenous organic compounds. Enzymes for assimilating amine into glutamine (GlnA), glutamate (GltBD), and glycine (AspC) were similarly upregulated. Differential regulation of the transcriptional regulator NtrX in the two-component signaling system NtrY/NtrX was also observed, implicating it in control of the nitrogen starvation response. Comparisons of the transcriptome and proteome supported previous observations of uncoupling between transcription and translation in nutrient-deprived "Ca. Pelagibacter ubique" cells. Overall, these data reveal a streamlined, PII-independent response to nitrogen stress in "Ca. Pelagibacter ubique," and likely other Pelagibacterales, and show that they respond to nitrogen stress by allocating more resources to the assimilation of nitrogen-rich organic compounds.

IMPORTANCE

Pelagibacterales are extraordinarily abundant and play a pivotal role in marine geochemical cycles, as one of the major recyclers of labile dissolved organic matter. They are also models for understanding how streamlining selection can reshape chemoheterotroph metabolism. Streamlining and its broad importance to environmental microbiology are emerging slowly from studies that reveal the complete genomes of uncultured organisms. Here, we report another remarkable example of streamlined metabolism in Pelagibacterales, this time in systems that control nitrogen assimilation. Pelagibacterales are major contributors to metatranscriptomes and metaproteomes from ocean systems, where patterns of gene expression are used to gain insight into ocean conditions and geochemical cycles. The data presented here supply background that is essential to interpreting data from field studies.

Comparative and phylogenomic evidence that the alphaproteobacterium HIMB59 is not a member of the oceanic SAR11 clade

SAR11 is a globally abundant group of Alphaproteobacteria in the oceans that is taxonomically not well defined. It has been suggested SAR11 should be classified into the novel order Pelagibacterales. Features such as conservation of gene content and synteny have been taken as evidence that also the divergent member HIMB59 should be included in the order. However, this proposition is controversial since phylogenetic analyses have questioned the monophyly of this grouping. Here, we performed phylogenetic analyses and reinvestigated the genomic similarity of SAR11 and HIMB59. Our phylogenetic analysis confirmed that HIMB59 is not a sister group to the other SAR11 strains. By placing the comparison in the context of the evolution of the Alphaproteobacteria, we found that none of the measures of genomic similarity supports a clustering of HIMB59 and SAR11 to the exclusion of other Alphaproteobacteria. First, pairwise sequence similarity measures for the SAR11 and HIMB59 genomes were within the range observed for unrelated pairs of Alphaproteobacteria. Second, pairwise comparisons of gene contents revealed a higher similarity of SAR11 to several other alphaproteobacterial genomes than to HIMB59. Third, the SAR11 genomes are not more similar in gene order to the HIMB59 genome than what they are to several other alphaproteobacterial genomes. Finally, in contrast to earlier reports, we observed no sequence similarity between the hypervariable region HVR2 in the SAR11 genomes and the region located at the corresponding position in the HIMB59 genome. Based on these observations, we conclude that the alphaproteobacterium HIMB59 is not monophyletic with the SAR11 strains and that genome streamlining has evolved multiple times independently in Alphaproteobacteria adapted to the upper surface waters of the oceans.

Comparative genomic analysis indicates that niche adaptation of terrestrial Flavobacteria is strongly linked to plant glycan metabolism

Flavobacteria are important members of aquatic and terrestrial bacterial communities, displaying extreme variations in lifestyle, geographical distribution and genome size. They are ubiquitous in soil, but are often strongly enriched in the rhizosphere and phyllosphere of plants. In this study, we compared the genome of a root-associated Flavobacterium that we recently isolated, physiologically characterized and sequenced, to 14 additional Flavobacterium genomes, in order to pinpoint characteristics associated with its high abundance in the rhizosphere. Interestingly, flavobacterial genomes vary in size by approximately two-fold, with terrestrial isolates having predominantly larger genomes than those from aquatic environments. Comparative functional gene analysis revealed that terrestrial and aquatic Flavobacteria generally segregated into two distinct clades. Members of the aquatic clade had a higher ratio of peptide and protein utilization genes, whereas members of the terrestrial clade were characterized by a significantly higher abundance and diversity of genes involved in metabolism of carbohydrates such as xylose, arabinose and pectin. Interestingly, genes encoding glycoside hydrolase (GH) families GH78 and GH106, responsible for rhamnogalacturonan utilization (exclusively associated with terrestrial plant hemicelluloses), were only present in terrestrial clade genomes, suggesting adaptation of the terrestrial strains to plant-related carbohydrate metabolism. The Peptidase/GH ratio of aquatic clade Flavobacteria was significantly higher than that of terrestrial strains (1.7±0.7 and 9.7±4.7, respectively), supporting the concept that this relation can be used to infer Flavobacterium lifestyles. Collectively, our research suggests that terrestrial Flavobacteria are highly adapted to plant carbohydrate metabolism, which appears to be a key to their profusion in plant environments.

Genome and transcriptome sequencing of the halophilic fungus Wallemia ichthyophaga: haloadaptations present and absent

Background

The basidomycete Wallemia ichthyophaga from the phylogenetically distinct class Wallemiomycetes is the most halophilic fungus known to date. It requires at least 10% NaCl and thrives in saturated salt solution. To investigate the genomic basis of this exceptional phenotype, we obtained a de-novo genome sequence of the species type-strain and analysed its transcriptomic response to conditions close to the limits of its lower and upper salinity range.

Results

The unusually compact genome is 9.6 Mb large and contains 1.67% repetitive sequences. Only 4884 predicted protein coding genes cover almost three quarters of the sequence. Of 639 differentially expressed genes, two thirds are more expressed at lower salinity. Phylogenomic analysis based on the largest dataset used to date (whole proteomes) positions Wallemiomycetes as a 250-million-year-old sister group of Agaricomycotina. Contrary to the closely related species Wallemia sebi, W. ichthyophaga appears to have lost the ability for sexual reproduction. Several protein families are significantly expanded or contracted in the genome. Among these, there are the P-type ATPase cation transporters, but not the sodium/ hydrogen exchanger family. Transcription of all but three cation transporters is not salt dependent. The analysis also reveals a significant enrichment in hydrophobins, which are cell-wall proteins with multiple cellular functions. Half of these are differentially expressed, and most contain an unusually large number of acidic amino acids. This discovery is of particular interest due to the numerous applications of hydrophobines from other fungi in industry, pharmaceutics and medicine.

Conclusions

W. ichthyophaga is an extremophilic specialist that shows only low levels of adaptability and genetic recombination. This is reflected in the characteristics of its genome and its transcriptomic response to salt. No unusual traits were observed in common salt-tolerance mechanisms, such as transport of inorganic ions or synthesis of compatible solutes. Instead, various data indicate a role of the cell wall of W. ichthyophaga in its response to salt. Availability of the genomic sequence is expected to facilitate further research into this unique species, and shed more light on adaptations that allow it to thrive in conditions lethal to most other eukaryotes.

Workflow management systems for gene sequence analysis and evolutionary studies - A Review

Post 'omic' era has resulted in the development of many primary, secondary and derived databases. Many analytical and visualization bioinformatics tools have been developed to manage and analyze the data available through large sequencing projects. Availability of heterogeneous databases and tools make it difficult for researchers to access information from varied sources and run different bioinformatics tools to get desired analysis done. Building integrated bioinformatics platforms is one of the most challenging tasks that bioinformatics community is facing. Integration of various databases, tools and algorithm is a challenging problem to deal with. This article describes the bioinformatics analysis workflow management systems that are developed in the area of gene sequence analysis and phylogeny. This article will be useful for biotechnologists, molecular biologists, computer scientists and statisticians engaged in computational biology and bioinformatics research.

Comparative genomic analysis of six bacteria belonging to the genus Novosphingobium: insights into marine adaptation, cell-cell signaling and bioremediation

Background

Bacteria belonging to the genus Novosphingobium are known to be metabolically versatile and occupy different ecological niches. In the absence of genomic data and/or analysis, knowledge of the bacteria that belong to this genus is currently limited to biochemical characteristics. In this study, we analyzed the whole genome sequencing data of six bacteria in the Novosphingobium genus and provide evidence to show the presence of genes that are associated with salt tolerance, cell-cell signaling and aromatic compound biodegradation phenotypes. Additionally, we show the taxonomic relationship between the sequenced bacteria based on phylogenomic analysis, average amino acid identity (AAI) and genomic signatures.

Results

The taxonomic clustering of Novosphingobium strains is generally influenced by their isolation source. AAI and genomic signature provide strong support the classification of Novosphingobium sp. PP1Y as Novosphingobium pentaromaticivorans PP1Y. The identification and subsequent functional annotation of the unique core genome in the marine Novosphingobium bacteria show that ectoine synthesis may be the main contributing factor in salt water adaptation. Genes coding for the synthesis and receptor of the cell-cell signaling molecules, of the N-acyl-homoserine lactones (AHL) class are identified. Notably, a solo luxR homolog was found in strain PP1Y that may have been recently acquired via horizontal gene transfer as evident by the presence of multiple mobile elements upstream of the gene. Additionally, phylogenetic tree analysis and sequence comparison with functionally validated aromatic ring hydroxylating dioxygenases (ARDO) revealed the presence of several ARDOs (oxygenase) in Novosphingobium bacteria with the majority of them belonging to the Groups II and III of the enzyme.

Conclusions

The combination of prior knowledge on the distinctive phenotypes of Novosphingobium strains and meta-analysis of their whole genomes enables the identification of several genes that are relevant in industrial applications and bioremediation. The results from such targeted but comprehensive comparative genomics analysis have the potential to contribute to the understanding of adaptation, cell-cell communication and bioremediation properties of bacteria belonging to the genus Novosphingobium.

The genome of tolypocladium inflatum: evolution, organization, and expression of the cyclosporin biosynthetic gene cluster

The ascomycete fungus Tolypocladium inflatum, a pathogen of beetle larvae, is best known as the producer of the immunosuppressant drug cyclosporin. The draft genome of T. inflatum strain NRRL 8044 (ATCC 34921), the isolate from which cyclosporin was first isolated, is presented along with comparative analyses of the biosynthesis of cyclosporin and other secondary metabolites in T. inflatum and related taxa. Phylogenomic analyses reveal previously undetected and complex patterns of homology between the nonribosomal peptide synthetase (NRPS) that encodes for cyclosporin synthetase (simA) and those of other secondary metabolites with activities against insects (e.g., beauvericin, destruxins, etc.), and demonstrate the roles of module duplication and gene fusion in diversification of NRPSs. The secondary metabolite gene cluster responsible for cyclosporin biosynthesis is described. In addition to genes necessary for cyclosporin biosynthesis, it harbors a gene for a cyclophilin, which is a member of a family of immunophilins known to bind cyclosporin. Comparative analyses support a lineage specific origin of the cyclosporin gene cluster rather than horizontal gene transfer from bacteria or other fungi. RNA-Seq transcriptome analyses in a cyclosporin-inducing medium delineate the boundaries of the cyclosporin cluster and reveal high levels of expression of the gene cluster cyclophilin. In medium containing insect hemolymph, weaker but significant upregulation of several genes within the cyclosporin cluster, including the highly expressed cyclophilin gene, was observed. T. inflatum also represents the first reference draft genome of Ophiocordycipitaceae, a third family of insect pathogenic fungi within the fungal order Hypocreales, and supports parallel and qualitatively distinct radiations of insect pathogens. The T. inflatum genome provides additional insight into the evolution and biosynthesis of cyclosporin and lays a foundation for further investigations of the role of secondary metabolite gene clusters and their metabolites in fungal biology.

Mutualistic co-evolution of type III effector genes in Sinorhizobium fredii and Bradyrhizobium japonicum

Two diametric paradigms have been proposed to model the molecular co-evolution of microbial mutualists and their eukaryotic hosts. In one, mutualist and host exhibit an antagonistic arms race and each partner evolves rapidly to maximize their own fitness from the interaction at potential expense of the other. In the opposing model, conflicts between mutualist and host are largely resolved and the interaction is characterized by evolutionary stasis. We tested these opposing frameworks in two lineages of mutualistic rhizobia, Sinorhizobium fredii and Bradyrhizobium japonicum. To examine genes demonstrably important for host-interactions we coupled the mining of genome sequences to a comprehensive functional screen for type III effector genes, which are necessary for many Gram-negative pathogens to infect their hosts. We demonstrate that the rhizobial type III effector genes exhibit a surprisingly high degree of conservation in content and sequence that is in contrast to those of a well characterized plant pathogenic species. This type III effector gene conservation is particularly striking in the context of the relatively high genome-wide diversity of rhizobia. The evolution of rhizobial type III effectors is inconsistent with the molecular arms race paradigm. Instead, our results reveal that these loci are relatively static in rhizobial lineages and suggest that fitness conflicts between rhizobia mutualists and their host plants have been largely resolved.

De novo transcriptome sequencing and analysis for Venturia inaequalis, the devastating apple scab pathogen

Venturia inaequalis is the causal agent of apple scab, one of the most devastating diseases of apple. Due to several distinct features, it has emerged as a model fungal pathogen to study various aspects of hemibiotrophic plant pathogen interactions. The present study reports de novo assembling, annotation and characterization of the transcriptome of V. inaequalis. Venturia transcripts expressed during its growth on laboratory medium and that expressed during its biotrophic stage of infection on apple were sequenced using Illumina RNAseq technology. A total of 94,350,055 reads (50 bp read length) specific to Venturia were obtained after filtering. The reads were assembled into 62,061 contigs representing 24,571 unique genes. GO analysis suggested prevalence of genes associated with biological process categories like metabolism, transport and response to stimulus. Genes associated with molecular function like binding, catalytic activities and transferase activities were found in majority. EC and KEGG pathway analyses suggested prevalence of genes encoding kinases, proteases, glycoside hydrolases, cutinases, cytochrome P450 and transcription factors. The study has identified several putative pathogenicity determinants and candidate effectors in V. inaequalis. A large number of transcripts encoding membrane transporters were identified and comparative analysis revealed that the number of transporters encoded by Venturia is significantly more as compared to that encoded by several other important plant fungal pathogens. Phylogenomics analysis indicated that V. inaequalis is closely related to Pyrenophora tritici-repentis (the causal organism of tan spot of wheat). In conclusion, the findings from this study provide a better understanding of the biology of the apple scab pathogen and have identified candidate genes/functions required for its pathogenesis. This work lays the foundation for facilitating further research towards understanding this host-pathogen interaction.

Comparative genomics of a plant-pathogenic fungus, Pyrenophora tritici-repentis, reveals transduplication and the impact of repeat elements on pathogenicity and population divergence

Pyrenophora tritici-repentis is a necrotrophic fungus causal to the disease tan spot of wheat, whose contribution to crop loss has increased significantly during the last few decades. Pathogenicity by this fungus is attributed to the production of host-selective toxins (HST), which are recognized by their host in a genotype-specific manner. To better understand the mechanisms that have led to the increase in disease incidence related to this pathogen, we sequenced the genomes of three P. tritici-repentis isolates. A pathogenic isolate that produces two known HSTs was used to assemble a reference nuclear genome of approximately 40 Mb composed of 11 chromosomes that encode 12,141 predicted genes. Comparison of the reference genome with those of a pathogenic isolate that produces a third HST, and a nonpathogenic isolate, showed the nonpathogen genome to be more diverged than those of the two pathogens. Examination of gene-coding regions has provided candidate pathogen-specific proteins and revealed gene families that may play a role in a necrotrophic lifestyle. Analysis of transposable elements suggests that their presence in the genome of pathogenic isolates contributes to the creation of novel genes, effector diversification, possible horizontal gene transfer events, identified copy number variation, and the first example of transduplication by DNA transposable elements in fungi. Overall, comparative analysis of these genomes provides evidence that pathogenicity in this species arose through an influx of transposable elements, which created a genetically flexible landscape that can easily respond to environmental changes.

Streamlining and core genome conservation among highly divergent members of the SAR11 clade

SAR11 is an ancient and diverse clade of heterotrophic bacteria that are abundant throughout the world’s oceans, where they play a major role in the ocean carbon cycle. Correlations between the phylogenetic branching order and spatiotemporal patterns in cell distributions from planktonic ocean environments indicate that SAR11 has evolved into perhaps a dozen or more specialized ecotypes that span evolutionary distances equivalent to a bacterial order. We isolated and sequenced genomes from diverse SAR11 cultures that represent three major lineages and encompass the full breadth of the clade. The new data expand observations about genome evolution and gene content that previously had been restricted to the SAR11 Ia subclade, providing a much broader perspective on the clade’s origins, evolution, and ecology. We found small genomes throughout the clade and a very high proportion of core genome genes (48 to 56%), indicating that small genome size is probably an ancestral characteristic. In their level of core genome conservation, the members of SAR11 are outliers, the most conserved free-living bacteria known. Shared features of the clade include low GC content, high gene synteny, a large hypervariable region bounded by rRNA genes, and low numbers of paralogs. Variation among the genomes included genes for phosphorus metabolism, glycolysis, and C1 metabolism, suggesting that adaptive specialization in nutrient resource utilization is important to niche partitioning and ecotype divergence within the clade. These data provide support for the conclusion that streamlining selection for efficient cell replication in the planktonic habitat has occurred throughout the evolution and diversification of this clade.

The SAR11 clade is the most abundant group of marine microorganisms worldwide, making them key players in the global carbon cycle. Growing knowledge about their biochemistry and metabolism is leading to a more mechanistic understanding of organic carbon oxidation and sequestration in the oceans. The discovery of small genomes in SAR11 provided crucial support for the theory that streamlining selection can drive genome reduction in low-nutrient environments. Study of isolates in culture revealed atypical organic nutrient requirements that can be attributed to genome reduction, such as conditional auxotrophy for glycine and its precursors, a requirement for reduced sulfur compounds, and evidence for widespread cycling of C1 compounds in marine environments. However, understanding the genetic variation and distribution of such pathways and characteristics like streamlining throughout the group has required the isolation and genome sequencing of diverse SAR11 representatives, an analysis of which we provide here.

Comparative genomics of plant-associated Pseudomonas spp.: insights into diversity and inheritance of traits involved in multitrophic interactions

We provide here a comparative genome analysis of ten strains within the Pseudomonas fluorescens group including seven new genomic sequences. These strains exhibit a diverse spectrum of traits involved in biological control and other multitrophic interactions with plants, microbes, and insects. Multilocus sequence analysis placed the strains in three sub-clades, which was reinforced by high levels of synteny, size of core genomes, and relatedness of orthologous genes between strains within a sub-clade. The heterogeneity of the P. fluorescens group was reflected in the large size of its pan-genome, which makes up approximately 54% of the pan-genome of the genus as a whole, and a core genome representing only 45–52% of the genome of any individual strain. We discovered genes for traits that were not known previously in the strains, including genes for the biosynthesis of the siderophores achromobactin and pseudomonine and the antibiotic 2-hexyl-5-propyl-alkylresorcinol; novel bacteriocins; type II, III, and VI secretion systems; and insect toxins. Certain gene clusters, such as those for two type III secretion systems, are present only in specific sub-clades, suggesting vertical inheritance. Almost all of the genes associated with multitrophic interactions map to genomic regions present in only a subset of the strains or unique to a specific strain. To explore the evolutionary origin of these genes, we mapped their distributions relative to the locations of mobile genetic elements and repetitive extragenic palindromic (REP) elements in each genome. The mobile genetic elements and many strain-specific genes fall into regions devoid of REP elements (i.e., REP deserts) and regions displaying atypical tri-nucleotide composition, possibly indicating relatively recent acquisition of these loci. Collectively, the results of this study highlight the enormous heterogeneity of the P. fluorescens group and the importance of the variable genome in tailoring individual strains to their specific lifestyles and functional repertoire.

We sequenced the genomes of seven strains of the Pseudomonas fluorescens group that colonize plant surfaces and function as biological control agents, protecting plants from disease. In this study, we demonstrated the genomic diversity of the group by comparing these strains to each other and to three other strains that were sequenced previously. Only about half of the genes in each strain are present in all of the other strains, and each strain has hundreds of unique genes that are not present in the other genomes. We mapped the genes that contribute to biological control in each genome and found that most of the biological control genes are in the variable regions of the genome, which are not shared by all of the other strains. This finding is consistent with our knowledge of the distinctive biology of each strain. Finally, we looked for new genes that are likely to confer antimicrobial traits needed to suppress plant pathogens, but have not been identified previously. In each genome, we discovered many of these new genes, which provide avenues for future discovery of new traits with the potential to manage plant diseases in agriculture or natural ecosystems.

The genome of the xerotolerant mold Wallemia sebi reveals adaptations to osmotic stress and suggests cryptic sexual reproduction

Wallemia (Wallemiales, Wallemiomycetes) is a genus of xerophilic Fungi of uncertain phylogenetic position within Basidiomycota. Most commonly found as food contaminants, species of Wallemia have also been isolated from hypersaline environments. The ability to tolerate environments with reduced water activity is rare in Basidiomycota. We sequenced the genome of W. sebi in order to understand its adaptations for surviving in osmotically challenging environments, and we performed phylogenomic and ultrastructural analyses to address its systematic placement and reproductive biology. W. sebi has a compact genome (9.8 Mb), with few repeats and the largest fraction of genes with functional domains compared with other Basidiomycota. We applied several approaches to searching for osmotic stress-related proteins. In silico analyses identified 93 putative osmotic stress proteins; homology searches showed the HOG (High Osmolarity Glycerol) pathway to be mostly conserved. Despite the seemingly reduced genome, several gene family expansions and a high number of transporters (549) were found that also provide clues to the ability of W. sebi to colonize harsh environments. Phylogenetic analyses of a 71-protein dataset support the position of Wallemia as the earliest diverging lineage of Agaricomycotina, which is confirmed by septal pore ultrastructure that shows the septal pore apparatus as a variant of the Tremella-type. Mating type gene homologs were identified although we found no evidence of meiosis during conidiogenesis, suggesting there may be aspects of the life cycle of W. sebi that remain cryptic.

Advances in insect phylogeny at the dawn of the postgenomic era

Most species on Earth are insects and thus, understanding their evolutionary relationships is key to understanding the evolution of life. Insect relationships are increasingly well supported, due largely to technological advances in molecular sequencing and phylogenetic computational analysis. In this postgenomic era, insect systematics will be furthered best by integrative methods aimed at hypothesis corroboration from molecular, morphological, and paleontological evidence. This review of the current consensus of insect relationships provides a foundation for comparative study and offers a framework to evaluate incoming genomic evidence. Notable recent phylogenetic successes include the resolution of Holometabola, including the identification of the enigmatic Strepsiptera as a beetle relative and the early divergence of Hymenoptera; the recognition of hexapods as a crustacean lineage within Pancrustacea; and the elucidation of Dictyoptera orders, with termites placed as social cockroaches. Regions of the tree that require further investigation include the earliest winged insects (Palaeoptera) and Polyneoptera (orthopteroid lineages).

One carbon metabolism in SAR11 pelagic marine bacteria

The SAR11 Alphaproteobacteria are the most abundant heterotrophs in the oceans and are believed to play a major role in mineralizing marine dissolved organic carbon. Their genomes are among the smallest known for free-living heterotrophic cells, raising questions about how they successfully utilize complex organic matter with a limited metabolic repertoire. Here we show that conserved genes in SAR11 subgroup Ia (Candidatus Pelagibacter ubique) genomes encode pathways for the oxidation of a variety of one-carbon compounds and methyl functional groups from methylated compounds. These pathways were predicted to produce energy by tetrahydrofolate (THF)-mediated oxidation, but not to support the net assimilation of biomass from C1 compounds. Measurements of cellular ATP content and the oxidation of ¹⁴C-labeled compounds to ¹⁴CO₂ indicated that methanol, formaldehyde, methylamine, and methyl groups from glycine betaine (GBT), trimethylamine (TMA), trimethylamine N-oxide (TMAO), and dimethylsulfoniopropionate (DMSP) were oxidized by axenic cultures of the SAR11 strain Ca. P. ubique HTCC1062. Analyses of metagenomic data showed that genes for C1 metabolism occur at a high frequency in natural SAR11 populations. In short term incubations, natural communities of Sargasso Sea microbial plankton expressed a potential for the oxidation of ¹⁴C-labeled formate, formaldehyde, methanol and TMAO that was similar to cultured SAR11 cells and, like cultured SAR11 cells, incorporated a much larger percentage of pyruvate and glucose (27–35%) than of C1 compounds (2–6%) into biomass. Collectively, these genomic, cellular and environmental data show a surprising capacity for demethylation and C1 oxidation in SAR11 cultures and in natural microbial communities dominated by SAR11, and support the conclusion that C1 oxidation might be a significant conduit by which dissolved organic carbon is recycled to CO₂ in the upper ocean.