ProbeMatch: rapid alignment of oligonucleotides to genome allowing both gaps and mismatches

Localization and symbiotic status of probiotics in the coral holobiont

Corals establish symbiotic relationships with microorganisms, especially endosymbiotic photosynthetic algae. Although other microbes have been commonly detected in coral tissues, their identity and beneficial functions for their host are unclear. Here, we confirm the beneficial outcomes of the inoculation of bacteria selected as probiotics and use fluorescence in situ hybridization (FISH) to define their localization in the coral Pocillopora damicornis. Our results show the first evidence of the inherent presence of Halomonas sp. and Cobetia sp. in native coral tissues, even before their inoculation. Furthermore, the relative enrichment of these coral tissue-associated bacteria through their inoculation in corals correlates with health improvements, such as increases in photosynthetic potential, and productivity. Our study suggests the symbiotic status of Halomonas sp. and Cobetia sp. in corals by indicating their localization within coral gastrodermis and epidermis and correlating their increased relative abundance through active inoculation with beneficial outcomes for the holobiont. This knowledge is crucial to facilitate the screening and application of probiotics that may not be transient members of the coral microbiome.

IMPORTANCE

Despite the promising results indicating the beneficial outcomes associated with the application of probiotics in corals and some scarce knowledge regarding the identity of bacterial cells found within the coral tissue, the correlation between these two aspects is still missing. This gap limits our understanding of the actual diversity of coral-associated bacteria and whether these symbionts are beneficial. Some researchers, for example, have been suggesting that probiotic screening should only focus on the very few known tissue-associated bacteria, such as Endozoicomonas sp., assuming that the currently tested probiotics are not tissue-associated. Here, we provide specific FISH probes for Halomonas sp. and Cobetia sp., expand our knowledge of the identity of coral-associated bacteria and confirm the probiotic status of the tested probiotics. The presence of these beneficial microorganisms for corals (BMCs) inside host tissues and gastric cavities also supports the notion that direct interactions with the host may underpin their probiotic role. This is a new breakthrough; these results argue against the possibility that the positive effects of BMCs are due to factors that are not related to a direct symbiotic interaction, for example, that the host simply feeds on inoculated bacteria or that the bacteria change the water quality.

Phylogenetics and biomineralization of a novel magnetotactic Gammaproteobacterium from a freshwater lake in Beijing, China

Magnetotactic bacteria (MTB) have the remarkable capability of producing intracellularly membrane-enveloped magnetic nanocrystals (i.e. magnetosomes) and swimming along geomagnetic field lines. Despite more than 50 years of research, bacterial diversity and magnetosome biomineralization within MTB are relatively less known in the Gammaproteobacteria class than other groups. This is incompatible with the status of Gammaproteobacteria as the most diverse class of gram-negative bacteria with a number of ecologically important bacteria. Here, we identify a novel MTB strain YYHR-1 affiliated with the Gammaproteobacteria class of the Pseudomonadota phylum from a freshwater lake. In YYHR-1, most magnetosome crystals are organized into a long chain aligned along the cell long axis; unusually, a few small superparamagnetic crystals are located at the side of the chain, off the main chain axis. Micromagnetic simulations indicate that magnetostatic interactions among adjacent crystals within a chain reduce the Gibbs energy to enhance chain stability. Genomic analysis suggests that duplication of magnetosome gene clusters may result in off-chain magnetosomes formation. By integrating available genomic data from Gammaproteobacteria, the phylogenetic position of MTB in this class is reassigned here. Our new findings expand knowledge about MTB diversity and magnetosome biomineralization, and deepen understanding of the phylogenetics of the Gammaproteobacteria.

Intracellular silicification by early-branching magnetotactic bacteria

Biosilicification-the formation of biological structures composed of silica-has a wide distribution among eukaryotes; it plays a major role in global biogeochemical cycles, and has driven the decline of dissolved silicon in the oceans through geological time. While it has long been thought that eukaryotes are the only organisms appreciably affecting the biogeochemical cycling of Si, the recent discoveries of silica transporter genes and marked silicon accumulation in bacteria suggest that prokaryotes may play an underappreciated role in the Si cycle, particularly in ancient times. Here, we report a previously unidentified magnetotactic bacterium that forms intracellular, amorphous silica globules. This bacterium, phylogenetically affiliated with the phylum Nitrospirota, belongs to a deep-branching group of magnetotactic bacteria that also forms intracellular magnetite magnetosomes and sulfur inclusions. This contribution reveals intracellularly controlled silicification within prokaryotes and suggests a previously unrecognized influence on the biogeochemical Si cycle that was operational during early Earth history.

A Novel Magnetotactic Alphaproteobacterium Producing Intracellular Magnetite and Calcium-Bearing Minerals

Magnetotactic bacteria (MTB) are prokaryotes that form intracellular magnetite (Fe₃O₄) or greigite (Fe₃S₄) nanocrystals with tailored sizes, often in chain configurations. Such magnetic particles are each surrounded by a lipid bilayer membrane, called a magnetosome, and provide a model system for studying the formation and function of specialized internal structures in prokaryotes. Using fluorescence-coupled scanning electron microscopy, we identified a novel magnetotactic spirillum, XQGS-1, from freshwater Xingqinggong Lake, Xi'an City, Shaanxi Province, China. Phylogenetic analyses based on 16S rRNA gene sequences indicate that strain XQGS-1 represents a novel genus of the Alphaproteobacteria class in the Proteobacteria phylum. Transmission electron microscopy analyses reveal that strain XQGS-1 forms on average 17 ± 3 magnetite magnetosome particles with an ideal truncated octahedral morphology, with an average length and width of 88.3 ± 11.7 nm and 83.3 ± 11.0 nm, respectively. They are tightly organized into a single chain along the cell long axis close to the concave side of the cell. Intrachain magnetic interactions likely result in these large equidimensional magnetite crystals behaving as magnetically stable single-domain particles that enable bacterial magnetotaxis. Combined structural and chemical analyses demonstrate that XQGS-1 cells also biomineralize intracellular amorphous calcium phosphate (2 to 3 granules per cell; 90.5- ± 19.3-nm average size) and weakly crystalline calcium carbonate (2 to 3 granules per cell; 100.4- ± 21.4-nm average size) in addition to magnetite. Our results expand the taxonomic diversity of MTB and provide evidence for intracellular calcium phosphate biomineralization in MTB. IMPORTANCE Biomineralization is a widespread process in eukaryotes that form shells, teeth, or bones. It also occurs commonly in prokaryotes, resulting in more than 60 known minerals formed by different bacteria under wide-ranging conditions. Among them, magnetotactic bacteria (MTB) are remarkable because they might represent the earliest organisms that biomineralize intracellular magnetic iron minerals (i.e., magnetite [Fe₃O₄] or greigite [Fe₃S₄]). Here, we report a novel magnetotactic spirillum (XQGS-1) that is phylogenetically affiliated with the Alphaproteobacteria class. In addition to magnetite crystals, XQGS-1 cells form intracellular submicrometer calcium carbonate and calcium phosphate granules. This finding supports the view that MTB are also an important microbial group for intracellular calcium carbonate and calcium phosphate biomineralization.

Identification and characterization of magnetotactic Gammaproteobacteria from a salt evaporation pool, Bohai Bay, China

Magnetotactic bacteria (MTB) are phylogenetically diverse prokaryotes that can produce intracellular chain-assembled nanocrystals of magnetite (Fe₃ O₄ ) or greigite (Fe₃ S₄ ). Compared with their wide distribution in the Alpha-, Eta- and Delta-proteobacteria classes, few MTB strains have been identified in the Gammaproteobacteria class, resulting in limited knowledge of bacterial diversity and magnetosome biomineralization within this phylogenetic branch. Here, we identify two magnetotactic Gammaproteobacteria strains (tentatively named FZSR-1 and FZSR-2 respectively) from a salt evaporation pool in Bohai Bay, at the Fuzhou saltern, Dalian City, eastern China. Phylogenetic analysis indicates that strain FZSR-2 is the same species as strains SHHR-1 and SS-5, which were discovered previously from brackish and hypersaline environments respectively. Strain FZSR-1 represents a novel species. Compared with strains FZSR-2, SHHR-1 and SS-5 in which magnetite particles are assembled into a single chain, FZSR-1 cells form relatively narrower magnetite nanoparticles that are often organized into double chains. We find a good relationship between magnetite morphology within strains FZSR-2, SHHR-1 and SS-5 and the salinity of the environment in which they live. This study expands the bacterial diversity of magnetotactic Gammaproteobacteria and provides new insights into magnetosome biomineralization within magnetotactic Gammaproteobacteria.

Characterisation of an atypical manifestation of black band disease on Porites lutea in the Western Indian Ocean

Recent surveys conducted on Reunion Island coral reefs revealed an atypical manifestation of black band disease on the main framework building coral, Porites lutea. This BBD manifestation (PorBBD) presented a thick lighter-colored band, which preceded the typical BBD lesion. Whilst BBD aetiology has been intensively described worldwide, it remains unclear if corals with apparently similar lesions across coral reefs are affected by the same pathogens. Therefore, a multidisciplinary approach involving field surveys, gross lesion monitoring, histopathology and 454-pyrosequencing was employed to provide the first comprehensive characterization of this particular manifestation. Surveys conducted within two geomorphological zones over two consecutive summers and winters showed spatial and seasonal patterns consistent with those found for typical BBD. Genetic analyses suggested an uncharacteristically high level of Vibrio spp. bacterial infection within PorBBD. However, microscopic analysis revealed high densities of cyanobacteria, penetrating the compromised tissue as well as the presence of basophilic bodies resembling bacterial aggregates in the living tissue, adjacent to the bacterial mat. Additionally, classical BBD-associated cyanobacterial strains, genetically related to Pseudoscillatoria coralii and Roseofilum reptotaenium were identified and isolated and the presence of sulfate-reducers or sulfide-oxidizers such as Desulfovibrio and Arcobacter, previously shown to be associated with anoxic microenvironment within typical BBD was also observed, confirming that PorBBD is a manifestation of classical BBD.

MOSAIK: a hash-based algorithm for accurate next-generation sequencing short-read mapping

MOSAIK is a stable, sensitive and open-source program for mapping second and third-generation sequencing reads to a reference genome. Uniquely among current mapping tools, MOSAIK can align reads generated by all the major sequencing technologies, including Illumina, Applied Biosystems SOLiD, Roche 454, Ion Torrent and Pacific BioSciences SMRT. Indeed, MOSAIK was the only aligner to provide consistent mappings for all the generated data (sequencing technologies, low-coverage and exome) in the 1000 Genomes Project. To provide highly accurate alignments, MOSAIK employs a hash clustering strategy coupled with the Smith-Waterman algorithm. This method is well-suited to capture mismatches as well as short insertions and deletions. To support the growing interest in larger structural variant (SV) discovery, MOSAIK provides explicit support for handling known-sequence SVs, e.g. mobile element insertions (MEIs) as well as generating outputs tailored to aid in SV discovery. All variant discovery benefits from an accurate description of the read placement confidence. To this end, MOSAIK uses a neural-network based training scheme to provide well-calibrated mapping quality scores, demonstrated by a correlation coefficient between MOSAIK assigned and actual mapping qualities greater than 0.98. In order to ensure that studies of any genome are supported, a training pipeline is provided to ensure optimal mapping quality scores for the genome under investigation. MOSAIK is multi-threaded, open source, and incorporated into our command and pipeline launcher system GKNO (http://gkno.me).

MOSAIK: a hash-based algorithm for accurate next-generation sequencing short-read mapping

MOSAIK is a stable, sensitive and open-source program for mapping second and third-generation sequencing reads to a reference genome. Uniquely among current mapping tools, MOSAIK can align reads generated by all the major sequencing technologies, including Illumina, Applied Biosystems SOLiD, Roche 454, Ion Torrent and Pacific BioSciences SMRT. Indeed, MOSAIK was the only aligner to provide consistent mappings for all the generated data (sequencing technologies, low-coverage and exome) in the 1000 Genomes Project. To provide highly accurate alignments, MOSAIK employs a hash clustering strategy coupled with the Smith-Waterman algorithm. This method is well-suited to capture mismatches as well as short insertions and deletions. To support the growing interest in larger structural variant (SV) discovery, MOSAIK provides explicit support for handling known-sequence SVs, e.g. mobile element insertions (MEIs) as well as generating outputs tailored to aid in SV discovery. All variant discovery benefits from an accurate description of the read placement confidence. To this end, MOSAIK uses a neural-network based training scheme to provide well-calibrated mapping quality scores, demonstrated by a correlation coefficient between MOSAIK assigned and actual mapping qualities greater than 0.98. In order to ensure that studies of any genome are supported, a training pipeline is provided to ensure optimal mapping quality scores for the genome under investigation. MOSAIK is multi-threaded, open source, and incorporated into our command and pipeline launcher system GKNO (http://gkno.me).

WITHDRAWN: Evaluation of next-generation sequencing software in mapping and assembly

Next-generation high-throughput DNA sequencing technologies have advanced progressively in sequence-based genomic research and novel biological applications with the promise of sequencing DNA at unprecedented speed. These new non-Sanger-based technologies feature several advantages, when compared with traditional sequencing methods in terms of higher sequencing speed, lower per run cost and higher accuracy. However, reads from next-generation sequencing (NGS) platforms, such as 454/Roche, ABI/SOLiD and Illumina/Solexa, are usually short, thereby restricting the applications of NGS platforms in genome assembly and annotation. We presented an overview of the challenges that these novel technologies meet and particularly illustrated various bioinformatics attempts on mapping and assembly for problem solving. We then compared the performance of several programs in these two fields and further provided advices on selecting suitable tools for specific biological applications.Journal of Human Genetics advance online publication, 16 June 2011; doi:10.1038/jhg.2011.62.

Evaluation of next-generation sequencing software in mapping and assembly

Next-generation high-throughput DNA sequencing technologies have advanced progressively in sequence-based genomic research and novel biological applications with the promise of sequencing DNA at unprecedented speed. These new non-Sanger-based technologies feature several advantages when compared with traditional sequencing methods in terms of higher sequencing speed, lower per run cost and higher accuracy. However, reads from next-generation sequencing (NGS) platforms, such as 454/Roche, ABI/SOLiD and Illumina/Solexa, are usually short, thereby restricting the applications of NGS platforms in genome assembly and annotation. We presented an overview of the challenges that these novel technologies meet and particularly illustrated various bioinformatics attempts on mapping and assembly for problem solving. We then compared the performance of several programs in these two fields, and further provided advices on selecting suitable tools for specific biological applications.

HiTEC: accurate error correction in high-throughput sequencing data

MOTIVATION

High-throughput sequencing technologies produce very large amounts of data and sequencing errors constitute one of the major problems in analyzing such data. Current algorithms for correcting these errors are not very accurate and do not automatically adapt to the given data.

RESULTS

We present HiTEC, an algorithm that provides a highly accurate, robust and fully automated method to correct reads produced by high-throughput sequencing methods. Our approach provides significantly higher accuracy than previous methods. It is time and space efficient and works very well for all read lengths, genome sizes and coverage levels.

AVAILABILITY

The source code of HiTEC is freely available at www.csd.uwo.ca/~ilie/HiTEC/.

Methylation analysis by DNA immunoprecipitation

DNA methylation regulates gene expression primarily through modification of chromatin structure. Global methylation studies have revealed biologically relevant patterns of DNA methylation in the human genome affecting sequences such as gene promoters, gene bodies, and repetitive elements. Disruption of normal methylation patterns and subsequent gene expression changes have been observed in several diseases especially in human cancers. Immunoprecipitation (IP)-based methods to evaluate methylation status of DNA have been instrumental in such genome-wide methylation studies. This review describes techniques commonly used to identify and quantify methylated DNA with emphasis on IP based platforms. In an effort to consolidate the wealth of information and highlight critical aspects of methylated DNA analysis, sample considerations, experimental and bioinformatic approaches for analyzing genome-wide methylation profiles, and the benefit of integrating DNA methylation data with complementary dimensions of genomic data are discussed.

PerM: efficient mapping of short sequencing reads with periodic full sensitive spaced seeds

MOTIVATION

The explosion of next-generation sequencing data has spawned the design of new algorithms and software tools to provide efficient mapping for different read lengths and sequencing technologies. In particular, ABI's sequencer (SOLiD system) poses a big computational challenge with its capacity to produce very large amounts of data, and its unique strategy of encoding sequence data into color signals.

RESULTS

We present the mapping software, named PerM (Periodic Seed Mapping) that uses periodic spaced seeds to significantly improve mapping efficiency for large reference genomes when compared with state-of-the-art programs. The data structure in PerM requires only 4.5 bytes per base to index the human genome, allowing entire genomes to be loaded to memory, while multiple processors simultaneously map reads to the reference. Weight maximized periodic seeds offer full sensitivity for up to three mismatches and high sensitivity for four and five mismatches while minimizing the number random hits per query, significantly speeding up the running time. Such sensitivity makes PerM a valuable mapping tool for SOLiD and Solexa reads.

AVAILABILITY

http://code.google.com/p/perm/

Fast and accurate short read alignment with Burrows-Wheeler transform

MOTIVATION

The enormous amount of short reads generated by the new DNA sequencing technologies call for the development of fast and accurate read alignment programs. A first generation of hash table-based methods has been developed, including MAQ, which is accurate, feature rich and fast enough to align short reads from a single individual. However, MAQ does not support gapped alignment for single-end reads, which makes it unsuitable for alignment of longer reads where indels may occur frequently. The speed of MAQ is also a concern when the alignment is scaled up to the resequencing of hundreds of individuals.

RESULTS

We implemented Burrows-Wheeler Alignment tool (BWA), a new read alignment package that is based on backward search with Burrows-Wheeler Transform (BWT), to efficiently align short sequencing reads against a large reference sequence such as the human genome, allowing mismatches and gaps. BWA supports both base space reads, e.g. from Illumina sequencing machines, and color space reads from AB SOLiD machines. Evaluations on both simulated and real data suggest that BWA is approximately 10-20x faster than MAQ, while achieving similar accuracy. In addition, BWA outputs alignment in the new standard SAM (Sequence Alignment/Map) format. Variant calling and other downstream analyses after the alignment can be achieved with the open source SAMtools software package.

AVAILABILITY

http://maq.sourceforge.net.