ALBERT

Description: Clonal populations accumulate mutations over time, resulting in different haplotypes. Deep sequencing of such a population in principle provides information to reconstruct these haplotypes and the frequency at which the haplotypes occur. However, this reconstruction is technically not trivial, especially not in clonal systems with a relatively low mutation frequency. The low number of segregating sites in those systems adds ambiguity to the haplotype phasing and thus obviates the reconstruction of genome-wide haplotypes based on sequence overlap information. Therefore, we present EVORhA, a haplotype reconstruction method that complements phasing information in the non-empty read overlap with the frequency estimations of inferred local haplotypes. As was shown with simulated data, as soon as read lengths and/or mutation rates become restrictive for state-of-the-art methods, the use of this additional frequency information allows EVORhA to still reliably reconstruct genome-wide haplotypes. On real data, we show the applicability of the method in reconstructing the population composition of evolved bacterial populations and in decomposing mixed bacterial infections from clinical samples.

Description: Data on biological mechanisms of aging are mostly obtained from cross-sectional study designs. An inherent disadvantage of this design is that inter-individual differences can mask small but biologically significant age-dependent changes. A serially sampled design (same individual at different time points) would overcome this problem but is often limited by the relatively small numbers of available paired samples and the statistics being used. To overcome these limitations, we have developed a new vector-based approach, termed three-component analysis, which incorporates temporal distance, signal intensity and variance into one single score for gene ranking and is combined with gene set enrichment analysis. We tested our method on a unique age-based sample set of human skin fibroblasts and combined genome-wide transcription, DNA methylation and histone methylation (H3K4me3 and H3K27me3) data. Importantly, our method can now for the first time demonstrate a clear age-dependent decrease in expression of genes coding for proteins involved in translation and ribosome function. Using analogies with data from lower organisms, we propose a model where age-dependent down-regulation of protein translation-related components contributes to extend human lifespan.

Description: A major challenge in the field of shotgun metagenomics is the accurate identification of organisms present within a microbial community, based on classification of short sequence reads. Though existing microbial community profiling methods have attempted to rapidly classify the millions of reads output from modern sequencers, the combination of incomplete databases, similarity among otherwise divergent genomes, errors and biases in sequencing technologies, and the large volumes of sequencing data required for metagenome sequencing has led to unacceptably high false discovery rates (FDR). Here, we present the application of a novel, gene-independent and signature-based metagenomic taxonomic profiling method with significantly and consistently smaller FDR than any other available method. Our algorithm circumvents false positives using a series of non-redundant signature databases and examines G enomic O rigins T hrough T axonomic CHA llenge (GOTTCHA). GOTTCHA was tested and validated on 20 synthetic and mock datasets ranging in community composition and complexity, was applied successfully to data generated from spiked environmental and clinical samples, and robustly demonstrates superior performance compared with other available tools.

Description: Four antibiotics (pamamycin, oligomycin A, oligomycin B and echinosporin) were isolated and characterized from the fermentation broth of the marine Streptomyces strains B8496 and B8739. Bioassays revealed that each of these compounds impaired motility and caused subsequent lysis of P. viticola zoospores in a dose- and time-dependent manner. Pamamycin displayed the strongest motility inhibitory and lytic activities (IC 50 0.1 μg mL –1 ) followed by oligomycin B (IC 50 0.15 and 0.2 μg mL –1 ) and oligomycin F (IC 50 0.3 and 0.5 μg mL –1 ). Oligomycin A and echinosporin also showed motility inhibitory activities against the zoospores with IC 50 values of 3.0 and 10.0 μg mL –1 , respectively. This is the first report of motility inhibitory and lytic activities of these antibiotics against zoospores of a phytopathogenic peronosporomycete. Structures of all the isolated compounds were determined based on detailed spectroscopic analysis.

Description: In sulfidic environments, microbes oxidize reduced sulfur compounds via several pathways. We used metagenomics to investigate sulfur metabolic pathways from microbial mat communities in two subterranean sulfidic streams in Lower Kane Cave, WY, USA and from Glenwood Hot Springs, CO, USA. Both unassembled and targeted recA gene assembly analyses revealed that these streams were dominated by Epsilonproteobacteria and Gammaproteobacteria , including groups related to Sulfurovum , Sulfurospirillum , Thiothrix and an epsilonproteobacterial group with no close cultured relatives. Genes encoding sulfide:quinone oxidoreductase (SQR) were abundant at all sites, but the specific SQR type and the taxonomic affiliation of each type differed between sites. The abundance of thiosulfate oxidation pathway genes (Sox) was not consistent between sites, although overall they were less abundant than SQR genes. Furthermore, the Sox pathway appeared to be incomplete in all samples. This work reveals both variations in sulfur metabolism within and between taxonomic groups found in these systems, and the presence of novel epsilonproteobacterial groups.

Description: Pseudomonas aeruginosa is an opportunistic pathogen with high resistance to a wide variety of antimicrobials. The multidrug resistance pump MexAB-OprM promotes the efflux of various antibiotics, mostly when mutations accumulate in the transcriptional regulators MexR, NalC and NalD, thereby causing MexAB-OprM overexpression. In this work, a characterization of 50 P. aeruginosa isolates obtained from Brazilian agricultural soils to determine the reasons of their resistance to aztreonam was done. The majority of the isolates showed higher aztreonam resistance than wild-type strain by MIC method. DNA sequence analysis of mexR , nalC and nalD genes from 13 of these isolates showed the amino acid substitution in NalC for all tested isolates, just one mutation was detected in MexR and none in NalD. Furthermore, an increase in the level of mexA expression by real-time RT-PCR analysis in eight isolates harboring mutations in NalC was found. Although there was not a relationship between MIC of aztreonam and the level of mexA expression, on the other hand, the results presented here suggest that novel mutations in NalC, including Arg 97 -Gly and Ala 186 -Thr, are related to MexAB-OprM overexpression causing aztreonam resistance in P. aeruginosa environmental isolates.

Description: Sedge-dominated wetlands on the Qinghai–Tibetan Plateau are methane emission centers. Methanotrophs at these sites play a role in reducing methane emissions, but relatively little is known about the composition of active methanotrophs in these wetlands. Here, we used DNA stable isotope probing to identify the key active aerobic methanotrophs in three sedge-dominated wetlands on the plateau. We found that Methylocystis species were active in two peatlands, Hongyuan and Dangxiong. Methylobacter species were found to be active only in Dangxiong peat. Hongyuan peat had the highest methane oxidation rate, and cross-feeding of carbon from methanotrophs to methylotrophic Hyphomicrobium species was observed. Owing to a low methane oxidation rate during the incubation, the labeling of methanotrophs in Maduo wetland samples was not detected. Our results indicate that there are large differences in the activity of methanotrophs in the wetlands of this region.

Description: Here we present the generation and function of two sets of bacterial plasmids that harbor fluorescent genes encoding either blue, cyan, yellow or red fluorescent proteins. In the first set, protein expression is controlled by the strong and constitutive nptII promoter whereas in the second set, the strong tac promoter was chosen that underlies LacI q regulation. Furthermore, the plasmids are mobilizable, contain Tn 7 transposons and a temperature-sensitive origin of replication. Using Escherichia coli S17-1 as donor strain, the plasmids allow fast and convenient Tn 7 -transposon delivery into many enterobacterial hosts, such as the here-used E. coli O157:H7. This procedure omits the need of preparing competent recipient cells and antibiotic resistances are only transiently conferred to the recipients. As the fluorescence proteins show little to no overlap in fluorescence emission, the constructs are well suited for the study of multicolored synthetic bacterial communities during biofilm production or in host colonization studies, e.g. of plant surfaces. Furthermore, tac promoter-reporter constructs allow the generation of so-called reproductive success reporters, which allow to estimate past doublings of bacterial individuals after introduction into environments, emphasizing the role of individual cells during colonization.

Description: Spa -typing and microarray techniques were used to study epidemiological changes in methicillin-resistant Staphylococcus aureus (MRSA) in South-East Austria. The population structure of 327 MRSA isolated between 2002 and 2012 was investigated. MRSA was assigned to 58 different spa types and 14 different MLST CC (multilocus sequence type clonal complexes); in particular, between 2007 and 2012, an increasing diversity in MRSA clones could be observed. The most abundant clonal complex was CC5. On the respective SCC mec cassettes, the CC5 isolates differed clearly within this decade and CC5/SCC mec I, the South German MRSA, predominant in 2002, was replaced by CC5/SCC mec II, the Rhine-Hesse MRSA in 2012. Whereas in many European countries MLST CC22-MRSA (EMRSA 15, the Barnim epidemic MRSA) is predominant, this clone occurred in Austria nearly 10 years later than in neighbouring countries. CC45, the Berlin EMRSA, epidemic in Germany, was only sporadically found in South-East Austria. The Irish ST8-MRSA-II represented by spa -type t190 was frequently found in 2002 and 2007, but disappeared in 2012. Our results demonstrate clonal replacement of MRSA clones within the last years in Austria. Ongoing surveillance is warranted for detection of changes within the MRSA population.

Description: This study aimed to investigate the effects of dietary fibre sources on the gut microbiota in suckling piglets, and to test the hypothesis that a moderate increase of dietary fibre may affect the gut microbiota during the suckling period. Suckling piglets were fed different fibre-containing diets or a control diet from postnatal day 7 to 22. Digesta samples from cecum, proximal colon and distal colon were used for Pig Intestinal Tract Chip analysis. The data showed that the effects of fibre-containing diet on the gut microbiota differed in the fibre source and gut location. The alfalfa diet increased Clostridium cluster XIVb and Sporobacter termitidis in the cecum compared to the pure cellulose diet. Compared to the control diet, the alfalfa diet also increased Coprococcus eutactus in the distal colon, while the pure cellulose diet decreased Eubacterium pyruvativorans in the cecum. The pure cellulose diet increased Prevotella ruminicola compared to the wheat bran diet. Interestingly, the alfalfa group had the lowest abundance of the potential pathogen Streptococcus suis in the cecum and distal colon. These results indicated that a moderate increase in dietary fibres affected the microbial composition in suckling piglets, and that the alfalfa inclusion produced some beneficial effects on the microbial communities.

Description: One function of the gut microbiota gaining recent attention, especially in herbivorous mammals and insects, is the metabolism of plant secondary metabolites (PSMs). We investigated whether this function exists within the gut communities of a specialist avian herbivore. We sequenced the cecal metagenome of the Greater Sage-Grouse ( Centrocercus urophasianus ), which specializes on chemically defended sagebrush ( Artemisia spp.). We predicted that the cecal metagenome of the sage-grouse would be enriched in genes associated with the metabolism of PSMs when compared to the metagenome of the domestic chicken. We found that representation of microbial genes associated with ‘xenobiotic degradation and metabolism’ was 3-fold higher in the sage-grouse cecal metagenomes when compared to that of the domestic chicken. Further, we identified a complete metabolic pathway for the degradation of phenol to pyruvate, which was not detected in the metagenomes of the domestic chicken, bovine rumen or 14 species of mammalian herbivores. Evidence of monoterpene degradation (a major class of PSMs in sagebrush) was less definitive, although we did detect genes for several enzymes associated with this process. Overall, our results suggest that the gut microbiota of specialist avian herbivores plays a similar role to the microbiota of mammalian and insect herbivores in degrading PSMs.

Description: Intracellular endosymbiotic bacteria are common and can play a crucial role for insect pathology. Therefore, such bacteria could be a potential key to our understanding of major losses of Western honey bees ( Apis mellifera ) colonies. However, the transmission and potential effects of endosymbiotic bacteria in A. mellifera and other Apis spp. are poorly understood. Here, we explore the prevalence and transmission of the genera Arsenophonus , Wolbachia , Spiroplasma and Rickettsia in Apis spp. Colonies of A. mellifera ( N = 33, with 20 eggs from worker brood cells and 100 adult workers each) as well as mated honey bee queens of A. cerana , A. dorsata and A. florea ( N = 12 each) were screened using PCR. While Wolbachia , Spiroplasma and Rickettsia were not detected, Arsenophonus spp. were found in 24.2% of A. mellifera colonies and respective queens as well as in queens of A. dorsata (8.3%) and A. florea (8.3%), but not in A. cerana . The absence of Arsenophonus spp. from reproductive organs of A. mellifera queens and surface-sterilized eggs does not support transovarial vertical transmission. Instead, horizontal transmission is most likely.

Description: Wood-rotting fungi possess remarkably diverse extracellular oxidation mechanisms, including enzymes, such as laccase and peroxidases, and Fenton chemistry. The ability to biologically drive Fenton chemistry by the redox cycling of quinones has previously been reported to be present in both ecologically diverging main groups of wood-rotting basidiomycetes. Therefore, we investigated whether it is even more widespread among fungal organisms. Screening of a diverse selection of a total of 18 ascomycetes and basidiomycetes for reduction of the model compound 2,6-dimethoxy benzoquinone revealed that all investigated strains were capable of reducing it to its corresponding hydroquinone. In a second step, depolymerization of the synthetic polymer polystyrene sulfonate was used as a proxy for quinone-dependent Fenton-based biodegradation capabilities. A diverse subset of the strains, including environmentally ubiquitous molds, white-rot fungi, as well as peatland and aquatic isolates, caused substantial depolymerization indicative for the effective employment of quinone redox cycling as biodegradation tool. Our results may also open up new paths to utilize diverse fungi for the bioremediation of recalcitrant organic pollutants.

Description: Ice-binding proteins (IBPs), such as antifreeze proteins (AFPs) and ice-nucleating proteins (INPs), have been described in diverse cold-adapted organisms, and their potential applications in biotechnology have been recognized in various fields. Currently, both IBPs are being applied to biotechnological processes, primarily in medicine and the food industry. However, our knowledge regarding the diversity of bacterial IBPs is limited; few studies have purified and characterized AFPs and INPs from bacteria. Phenotypically verified IBPs have been described in members belonging to Gammaproteobacteria, Actinobacteria and Flavobacteriia classes, whereas putative IBPs have been found in Gammaproteobacteria, Alphaproteobacteria and Bacilli classes. Thus, the main goal of this minireview is to summarize the current information on bacterial IBPs and their application in biotechnology, emphasizing the potential application in less explored fields such as agriculture. Investigations have suggested the use of INP-producing bacteria antagonists and AFPs-producing bacteria (or their AFPs) as a very attractive strategy to prevent frost damages in crops. UniProt database analyses of reported IBPs (phenotypically verified) and putative IBPs also show the limited information available on bacterial IBPs and indicate that major studies are required.

Description: Triazophos is a broad-spectrum and highly effective insecticide, and the residues of triazophos have been frequently detected in the environment. A triazophos-degrading bacterium, Burkholderia sp. SZL-1, was isolated from a long-term triazophos-polluted soil. Strain SZL-1 could hydrolyze triazophos to 1-phenyl-3-hydroxy-1,2,4-triazole, which was further utilized as the carbon sources for growth. The triazophos hydrolase gene trhA , cloned from strain SZL-1, was expressed and homogenously purified using Ni-nitrilotriacetic acid affinity chromatography. TrhA is 55 kDa and displays maximum activity at 25°C, pH 8.0. This enzyme still has nearly 60% activity at the range of 15°C–50°C for 30 min. TrhA was mutated by sequential error prone PCR and screened for improved activity for triazophos degradation. One purified variant protein (Val89-Gly89) named TrhA-M1 showed up to 3-fold improvement in specific activity against triazophos, and the specificity constants of K cat and K cat / K m for TrhA-M1 were improved up to 2.3- and 8.28-fold, respectively, compared to the wild-type enzyme. The results in this paper provided potential material for the contaminated soil remediation and hydrolase genetic structure research.

Description: The metal mining industry faces many large challenges in future years, among which is the increasing need to process low-grade ores as accessible higher grade ores become depleted. This is against a backdrop of increasing global demands for base and precious metals, and rare earth elements. Typically about 99% of solid material hauled to, and ground at, the land surface currently ends up as waste (rock dumps and mineral tailings). Exposure of these to air and water frequently leads to the formation of acidic, metal-contaminated run-off waters, referred to as acid mine drainage, which constitutes a severe threat to the environment. Formation of acid drainage is a natural phenomenon involving various species of lithotrophic (literally ‘rock-eating’) bacteria and archaea, which oxidize reduced forms of iron and/or sulfur. However, other microorganisms that reduce inorganic sulfur compounds can essentially reverse this process. These microorganisms can be applied on industrial scale to precipitate metals from industrial mineral leachates and acid mine drainage streams, resulting in a net improvement in metal recovery, while minimizing the amounts of leachable metals to the tailings storage dams. Here, we advocate that more extensive exploitation of microorganisms in metal mining operations could be an important way to green up the industry, reducing environmental risks and improving the efficiency and the economy of metal recovery.

Description: Chromosome-long haplotyping of human genomes is important to identify genetic variants with differing gene expression, in human evolution studies, clinical diagnosis, and other biological and medical fields. Although several methods have realized haplotyping based on sequencing technologies or population statistics, accuracy and cost are factors that prohibit their wide use. Borrowing ideas from group testing theories, we proposed a clone-based haplotyping method by overlapping pool sequencing. The clones from a single individual were pooled combinatorially and then sequenced. According to the distinct pooling pattern for each clone in the overlapping pool sequencing, alleles for the recovered variants could be assigned to their original clones precisely. Subsequently, the clone sequences could be reconstructed by linking these alleles accordingly and assembling them into haplotypes with high accuracy. To verify the utility of our method, we constructed 130 110 clones in silico for the individual NA12878 and simulated the pooling and sequencing process. Ultimately, 99.9% of variants on chromosome 1 that were covered by clones from both parental chromosomes were recovered correctly, and 112 haplotype contigs were assembled with an N50 length of 3.4 Mb and no switch errors. A comparison with current clone-based haplotyping methods indicated our method was more accurate.

Description: Sexual differentiation of malaria parasites into gametocytes in the vertebrate host and subsequent gamete fertilization in mosquitoes is essential for the spreading of the disease. The molecular processes orchestrating these transitions are far from fully understood. Here, we report the first transcriptome analysis of male and female Plasmodium falciparum gametocytes coupled with a comprehensive proteome analysis. In male gametocytes there is an enrichment of proteins involved in the formation of flagellated gametes; proteins involved in DNA replication, chromatin organization and axoneme formation. On the other hand, female gametocytes are enriched in proteins required for zygote formation and functions after fertilization; protein-, lipid- and energy-metabolism. Integration of transcriptome and proteome data revealed 512 highly expressed maternal transcripts without corresponding protein expression indicating large scale translational repression in P. falciparum female gametocytes for the first time. Despite a high degree of conservation between Plasmodium species, 260 of these ‘repressed transcripts’ have not been previously described. Moreover, for some of these genes, protein expression is only reported in oocysts and sporozoites indicating that repressed transcripts can be partitioned into short- and long-term storage. Finally, these data sets provide an essential resource for identification of vaccine/drug targets and for further mechanistic studies.

Description: DNA methylation is a mechanism for long-term transcriptional regulation and is required for normal cellular differentiation. Failure to properly establish or maintain DNA methylation patterns leads to cell dysfunction and diseases such as cancer. Identifying DNA methylation signatures in complex tissues can be challenging owing to inaccurate cell enrichment methods and low DNA yields. We have developed a technique called laser capture microdissection-reduced representation bisulfite sequencing (LCM-RRBS) for the multiplexed interrogation of the DNA methylation status of cytosine–guanine dinucleotide islands and promoters. LCM-RRBS accurately and reproducibly profiles genome-wide methylation of DNA extracted from microdissected fresh frozen or formalin-fixed paraffin-embedded tissue samples. To demonstrate the utility of LCM-RRBS, we characterized changes in DNA methylation associated with gonadectomy-induced adrenocortical neoplasia in the mouse. Compared with adjacent normal tissue, the adrenocortical tumors showed reproducible gains and losses of DNA methylation at genes involved in cell differentiation and organ development. LCM-RRBS is a rapid, cost-effective, and sensitive technique for analyzing DNA methylation in heterogeneous tissues and will facilitate the investigation of DNA methylation in cancer and organ development.

Description: The advent in high-throughput-sequencing (HTS) technologies has revolutionized conventional biodiversity research by enabling parallel capture of DNA sequences possessing species-level diagnosis. However, polymerase chain reaction (PCR)-based implementation is biased by the efficiency of primer binding across lineages of organisms. A PCR-free HTS approach will alleviate this artefact and significantly improve upon the multi-locus method utilizing full mitogenomes. Here we developed a novel multiplex sequencing and assembly pipeline allowing for simultaneous acquisition of full mitogenomes from pooled animals without DNA enrichment or amplification. By concatenating assemblies from three de novo assemblers, we obtained high-quality mitogenomes for all 49 pooled taxa, with 36 species 〉15 kb and the remaining 〉10 kb, including 20 complete mitogenomes and nearly all protein coding genes (99.6%). The assembly quality was carefully validated with Sanger sequences, reference genomes and conservativeness of protein coding genes across taxa. The new method was effective even for closely related taxa, e.g. three Drosophila spp., demonstrating its broad utility for biodiversity research and mito-phylogenomics. Finally, the in silico simulation showed that by recruiting multiple mito-loci, taxon detection was improved at a fixed sequencing depth. Combined, these results demonstrate the plausibility of a multi-locus mito-metagenomics approach as the next phase of the current single-locus metabarcoding method.

Description: Chromatin modifiers and histone modifications are components of a chromatin-signaling network involved in transcription and its regulation. The interactions between chromatin modifiers and histone modifications are often unknown, are based on the analysis of few genes or are studied in vitro . Here, we apply computational methods to recover interactions between chromatin modifiers and histone modifications from genome-wide ChIP-Seq data. These interactions provide a high-confidence backbone of the chromatin-signaling network. Many recovered interactions have literature support; others provide hypotheses about yet unknown interactions. We experimentally verified two of these predicted interactions, leading to a link between H4K20me1 and members of the Polycomb Repressive Complexes 1 and 2. Our results suggest that our computationally derived interactions are likely to lead to novel biological insights required to establish the connectivity of the chromatin-signaling network involved in transcription and its regulation.

Description: With read lengths of currently up to 2 x 300 bp, high throughput and low sequencing costs Illumina's MiSeq is becoming one of the most utilized sequencing platforms worldwide. The platform is manageable and affordable even for smaller labs. This enables quick turnaround on a broad range of applications such as targeted gene sequencing, metagenomics, small genome sequencing and clinical molecular diagnostics. However, Illumina error profiles are still poorly understood and programs are therefore not designed for the idiosyncrasies of Illumina data. A better knowledge of the error patterns is essential for sequence analysis and vital if we are to draw valid conclusions. Studying true genetic variation in a population sample is fundamental for understanding diseases, evolution and origin. We conducted a large study on the error patterns for the MiSeq based on 16S rRNA amplicon sequencing data. We tested state-of-the-art library preparation methods for amplicon sequencing and showed that the library preparation method and the choice of primers are the most significant sources of bias and cause distinct error patterns. Furthermore we tested the efficiency of various error correction strategies and identified quality trimming (Sickle) combined with error correction (BayesHammer) followed by read overlapping (PANDAseq) as the most successful approach, reducing substitution error rates on average by 93%.

Description: RNA-seq is a sensitive and accurate technique to compare steady-state levels of RNA between different cellular states. However, as it does not provide an account of transcriptional activity per se , other technologies are needed to more precisely determine acute transcriptional responses. Here, we have developed an easy, sensitive and accurate novel computational method, iRNA-seq , for genome-wide assessment of transcriptional activity based on analysis of intron coverage from total RNA-seq data. Comparison of the results derived from iRNA-seq analyses with parallel results derived using current methods for genome-wide determination of transcriptional activity, i.e. global run-on (GRO)-seq and RNA polymerase II (RNAPII) ChIP-seq, demonstrate that iRNA-seq provides similar results in terms of number of regulated genes and their fold change. However, unlike the current methods that are all very labor-intensive and demanding in terms of sample material and technologies, iRNA-seq is cheap and easy and requires very little sample material. In conclusion, iRNA-seq offers an attractive novel alternative to current methods for determination of changes in transcriptional activity at a genome-wide level.

Description: Identifying large-scale structural variation in cancer genomes continues to be a challenge to researchers. Current methods rely on genome alignments based on a reference that can be a poor fit to highly variant and complex tumor genomes. To address this challenge we developed a method that uses available breakpoint information to generate models of structural variations. We use these models as references to align previously unmapped and discordant reads from a genome. By using these models to align unmapped reads, we show that our method can help to identify large-scale variations that have been previously missed.

Description: Differential inhibitors are important for measuring the relative contributions of microbial groups, such as ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA), to biogeochemical processes in environmental samples. In particular, 2-phenyl-4,4,5,5-tetramethylimidazoline-1-oxyl 3-oxide (PTIO) represents a nitric oxide scavenger used for the specific inhibition of AOA, implicating nitric oxide as an intermediate of thaumarchaeotal ammonia oxidation. This study investigated four alternative nitric oxide scavengers for their ability to differentially inhibit AOA and AOB in comparison to PTIO. Caffeic acid, curcumin, methylene blue hydrate and trolox were tested on Nitrosopumilus maritimus , two unpublished AOA representatives (AOA-6f and AOA-G6) as well as the AOB representative Nitrosomonas europaea . All four scavengers inhibited ammonia oxidation by AOA at lower concentrations than for AOB. In particular, differential inhibition of AOA and AOB by caffeic acid (100 μM) and methylene blue hydrate (3 μM) was comparable to carboxy-PTIO (100 μM) in pure and enrichment culture incubations. However, when added to aquarium sponge biofilm microcosms, both scavengers were unable to inhibit ammonia oxidation consistently, likely due to degradation of the inhibitors themselves. This study provides evidence that a variety of nitric oxide scavengers result in differential inhibition of ammonia oxidation in AOA and AOB, and provides support to the proposed role of nitric oxide as a key intermediate in the thaumarchaeotal ammonia oxidation pathway.

Description: Any given human individual carries multiple genetic variants that disrupt protein-coding genes, through structural variation, as well as nucleotide variants and indels. Predicting the phenotypic consequences of a gene disruption remains a significant challenge. Current approaches employ information from a range of biological networks to predict which human genes are haploinsufficient (meaning two copies are required for normal function) or essential (meaning at least one copy is required for viability). Using recently available study gene sets, we show that these approaches are strongly biased towards providing accurate predictions for well-studied genes. By contrast, we derive a haploinsufficiency score from a combination of unbiased large-scale high-throughput datasets, including gene co-expression and genetic variation in over 6000 human exomes. Our approach provides a haploinsufficiency prediction for over twice as many genes currently unassociated with papers listed in Pubmed as three commonly-used approaches, and outperforms these approaches for predicting haploinsufficiency for less-studied genes. We also show that fine-tuning the predictor on a set of well-studied ‘gold standard’ haploinsufficient genes does not improve the prediction for less-studied genes. This new score can readily be used to prioritize gene disruptions resulting from any genetic variant, including copy number variants, indels and single-nucleotide variants.

Description: Meta-analysis of gene expression has enabled numerous insights into biological systems, but current methods have several limitations. We developed a method to perform a meta-analysis using the elastic net, a powerful and versatile approach for classification and regression. To demonstrate the utility of our method, we conducted a meta-analysis of lung cancer gene expression based on publicly available data. Using 629 samples from five data sets, we trained a multinomial classifier to distinguish between four lung cancer subtypes. Our meta-analysis-derived classifier included 58 genes and achieved 91% accuracy on leave-one-study-out cross-validation and on three independent data sets. Our method makes meta-analysis of gene expression more systematic and expands the range of questions that a meta-analysis can be used to address. As the amount of publicly available gene expression data continues to grow, our method will be an effective tool to help distill these data into knowledge.

Description: Peatlands of all latitudes play an integral role in global climate change by serving as a carbon sink and a primary source of atmospheric methane; however, the microbial ecology of mid-latitude peatlands is vastly understudied. Herein, next generation Illumina amplicon sequencing of small subunit rRNA genes was utilized to elucidate the microbial communities in three southern Appalachian peatlands. In contrast to northern peatlands, Proteobacteria dominated over Acidobacteria in all three sites. An average of 11 bacterial phyla was detected at relative abundance values 〉1%, with three candidate divisions (OP3, WS3 and NC10) represented, indicating high phylogenetic diversity. Physiological traits of isolates within the candidate alphaproteobacterial order, Ellin 329, obtained here and in previous studies indicate that bacteria of this order may be involved in hydrolysis of poly-, di- and monosaccharides. Community analyses indicate that Ellin 329 is the third most abundant order and is most abundant near the surface layers where plant litter decomposition should be primarily occurring. In sum, members of Ellin 329 likely play important roles in organic matter decomposition, in southern Appalachian peatlands and should be investigated further in other peatlands and ecosystem types.

Description: Marine viruses are the most abundant biological entity in the oceans, the majority of which infect bacteria and are known as bacteriophages. Yet, the bulk of bacteriophages form part of the vast uncultured dark matter of the microbial biosphere. In spite of the paucity of cultured marine bacteriophages, it is known that marine bacteriophages have major impacts on microbial population structure and the biogeochemical cycling of key elements. Despite the ecological relevance of marine bacteriophages, there are relatively few isolates with complete genome sequences. This minireview focuses on knowledge gathered from these genomes put in the context of viral metagenomic data and highlights key advances in the field, particularly focusing on genome structure and auxiliary metabolic genes.

Description: The fynbos biome in South Africa is globally recognised as a plant biodiversity hotspot. However, very little is known about the bacterial communities associated with fynbos plants, despite interactions between primary producers and bacteria having an impact on the physiology of both partners and shaping ecosystem diversity. This study reports on the structure, phylogenetic composition and potential roles of the endophytic bacterial communities located in the stems of three fynbos plants ( Erepsia anceps , Phaenocoma prolifera and Leucadendron laureolum ). Using Illumina MiSeq 16S rRNA sequencing we found that different subpopulations of Deinococcus-Thermus, Alphaproteobacteria, Acidobacteria and Firmicutes dominated the endophytic bacterial communities. Alphaproteobacteria and Actinobacteria were prevalent in P. prolifera , whereas Deinococcus-Thermus dominated in L. laureolum , revealing species-specific host–bacteria associations. Although a high degree of variability in the endophytic bacterial communities within hosts was observed, we also detected a core microbiome across the stems of the three plant species, which accounted for 72% of the sequences. Altogether, it seems that both deterministic and stochastic processes shaped microbial communities. Endophytic bacterial communities harboured putative plant growth-promoting bacteria, thus having the potential to influence host health and growth.

Description: The functioning of many natural and engineered environments is dependent on long distance electron transfer mediated through electrical currents. These currents have been observed in exoelectrogenic biofilms and it has been proposed that microbial biofilms can mediate electron transfer via electrical currents on the centimeter scale. However, direct evidence to confirm this hypothesis has not been demonstrated and the longest known electrical transfer distance for single species exoelectrogenic biofilms is limited to 100 μm. In the present study, biofilms were developed on electrodes with electrically non-conductive gaps from 50 μm to 1 mm and the in situ conductance of biofilms was evaluated over time. Results demonstrated that the exoelectrogenic mixed species biofilms in the present study possess the ability to transfer electrons through electrical currents over a distance of up to 1 mm, 10 times further than previously observed. Results indicate the possibility of interspecies interactions playing an important role in the spatial development of exoelectrogenic biofilms, suggesting that these biological networks might remain conductive even at longer distance. These findings have significant implications in regards to future optimization of microbial electrochemical systems.

Description: Advanced sequencing technologies have generated a plethora of data for many chromatin marks in multiple tissues and cell types, yet there is lack of a generalized tool for optimal utility of those data. A major challenge is to quantitatively model the epigenetic dynamics across both the genome and many cell types for understanding their impacts on differential gene regulation and disease. We introduce IDEAS, an i ntegrative and d iscriminative e pigenome a nnotation s ystem, for jointly characterizing epigenetic landscapes in many cell types and detecting differential regulatory regions. A key distinction between our method and existing state-of-the-art algorithms is that IDEAS integrates epigenomes of many cell types simultaneously in a way that preserves the position-dependent and cell type-specific information at fine scales, thereby greatly improving segmentation accuracy and producing comparable annotations across cell types.

Description: Bacteriophages are increasingly being used as water quality indicators. Two groups of phages infecting Escherichia coli , somatic and F-specific coliphages, are being considered as indicators of fecal and viral contamination for several types of water around the world. However, some uncertainties remain regarding which coliphages to assess. Recently, E. coli strain CB390 has been reported to be suitable for simultaneous detection of both groups, which seems to be more informative than determining only one of the groups. Here, a significant number of samples from different settings, mostly those where F-specific phages have been reported to outnumber somatic coliphages, are analyzed for somatic coliphages, F-specific RNA phages by standardized methods and coliphages detected by host strain CB390. The results presented here confirm that the numbers of phages counted using CB390 are equivalent to the sum of the somatic and F-specific coliphages counted independently in all settings. Hence the usefulness of this strain for simultaneous detection of somatic and F-specific coliphages is confirmed. Also, sets of data on the presence of coliphages in reclaimed and groundwater are reported.

Description: It is well known that Methylosinus trichosporium OB3b has two forms of methane monooxygenase (MMO) responsible for the initial conversion of methane to methanol, a cytoplasmic (soluble) methane monooxygenase and a membrane-associated (particulate) methane monooxygenase, and that copper strongly regulates expression of these alternative forms of MMO. More recently, it has been discovered that M. trichosporium OB3b has multiple types of the methanol dehydrogenase (MeDH), i.e. the Mxa-type MeDH (Mxa-MeDH) and Xox-type MeDH (Xox-MeDH), and the expression of these two forms is regulated by the availability of the rare earth element (REE), cerium. Here, we extend these studies and show that lanthanum, praseodymium, neodymium and samarium also regulate expression of alternative forms of MeDH. The effect of these REEs on MeDH expression, however, was only observed in the absence of copper. Further, a mutant of M. trichosporium OB3b, where the Mxa-MeDH was knocked out, was able to grow in the presence of lanthanum, praseodymium and neodymium, but was not able to grow in the presence of samarium. Collectively, these data suggest that multiple levels of gene regulation by metals exist in M. trichosporium OB3b, but that copper overrides the effect of other metals by an as yet unknown mechanism.

Description: A majority of large-scale bacterial genome rearrangements involve mobile genetic elements such as insertion sequence (IS) elements. Here we report novel insertions and excisions of IS elements and recombination between homologous IS elements identified in a large collection of Escherichia coli mutation accumulation lines by analysis of whole genome shotgun sequencing data. Based on 857 identified events (758 IS insertions, 98 recombinations and 1 excision), we estimate that the rate of IS insertion is 3.5 x 10 –4 insertions per genome per generation and the rate of IS homologous recombination is 4.5 x 10 –5 recombinations per genome per generation. These events are mostly contributed by the IS elements IS 1 , IS 2 , IS 5 and IS 186 . Spatial analysis of new insertions suggest that transposition is biased to proximal insertions, and the length spectrum of IS-caused deletions is largely explained by local hopping. For any of the ISs studied there is no region of the circular genome that is favored or disfavored for new insertions but there are notable hotspots for deletions. Some elements have preferences for non-coding sequence or for the beginning and end of coding regions, largely explained by target site motifs. Interestingly, transposition and deletion rates remain constant across the wild-type and 12 mutant E. coli lines, each deficient in a distinct DNA repair pathway. Finally, we characterized the target sites of four IS families, confirming previous results and characterizing a highly specific pattern at IS 186 target-sites, 5'-GGGG(N6/N7)CCCC-3'. We also detected 48 long deletions not involving IS elements.

Description: Polysulfides (S x 2– ) are sulfide oxidation intermediates that are important for a variety of environmentally relevant processes including pyrite formation, organic matter sulfidization, isotope exchange among reduced sulfur species, and metal chelation. In addition to their chemical reactivity, laboratory experiments with microbial cultures and enzymes indicate both indirect and direct roles for microorganisms in affecting polysulfide chemistry in natural environments through production and consumption. As polysulfides have been detected in a wide array of natural systems ranging from microbial mats to hydrothermal vents, constraining their biogeochemical cycling has broad impacts. However, many questions remain regarding the processes responsible for polysulfide dynamics in these environments and the precise role that microorganisms play in these processes. This review provides a summary of laboratory experiments investigating the role of polysulfides in microbial metabolism, and observations of polysulfides in the environment in order to provide further insight into and highlight open questions about this significant component of the sulfur cycle.

Description: Extensive and multi-dimensional data sets generated from recent cancer omics profiling projects have presented new challenges and opportunities for unraveling the complexity of cancer genome landscapes. In particular, distinguishing the unique complement of genes that drive tumorigenesis in each patient from a sea of passenger mutations is necessary for translating the full benefit of cancer genome sequencing into the clinic. We address this need by presenting a data integration framework (OncoIMPACT) to nominate patient-specific driver genes based on their phenotypic impact. Extensive in silico and in vitro validation helped establish OncoIMPACT's robustness, improved precision over competing approaches and verifiable patient and cell line specific predictions (2/2 and 6/7 true positives and negatives, respectively). In particular, we computationally predicted and experimentally validated the gene TRIM24 as a putative novel amplified driver in a melanoma patient. Applying OncoIMPACT to more than 1000 tumor samples, we generated patient-specific driver gene lists in five different cancer types to identify modes of synergistic action. We also provide the first demonstration that computationally derived driver mutation signatures can be overall superior to single gene and gene expression based signatures in enabling patient stratification and prognostication. Source code and executables for OncoIMPACT are freely available from http://sourceforge.net/projects/oncoimpact .

Description: Analysis of rewired upstream subnetworks impacting downstream differential gene expression aids the delineation of evolving molecular mechanisms. Cumulative statistics based on conventional differential correlation are limited for subnetwork rewiring analysis since rewiring is not necessarily equivalent to change in correlation coefficients. Here we present a computational method ChiNet to quantify subnetwork rewiring by statistical heterogeneity that enables detection of potential genotype changes causing altered transcription regulation in evolving organisms. Given a differentially expressed downstream gene set, ChiNet backtracks a rewired upstream subnetwork from a super-network including gene interactions known to occur under various molecular contexts. We benchmarked ChiNet for its high accuracy in distinguishing rewired artificial subnetworks, in silico yeast transcription-metabolic subnetworks, and rewired transcription subnetworks for Candida albicans versus Saccharomyces cerevisiae , against two differential-correlation based subnetwork rewiring approaches. Then, using transcriptome data from tolerant S. cerevisiae strain NRRL Y-50049 and a wild-type intolerant strain, ChiNet identified 44 metabolic pathways affected by rewired transcription subnetworks anchored to major adaptively activated transcription factor genes YAP1 , RPN4 , SFP1 and ROX1 , in response to toxic chemical challenges involved in lignocellulose-to-biofuels conversion. These findings support the use of ChiNet in rewiring analysis of subnetworks where differential interaction patterns resulting from divergent nonlinear dynamics abound.

Description: Ribosome biogenesis, a central and essential cellular process, occurs through sequential association and mutual co-folding of protein–RNA constituents in a well-defined assembly pathway. Here, we construct a network of co-evolving nucleotide/amino acid residues within the ribosome and demonstrate that assembly constraints are strong predictors of co-evolutionary patterns. Predictors of co-evolution include a wide spectrum of structural reconstitution events, such as cooperativity phenomenon, protein-induced rRNA reconstitutions, molecular packing of different rRNA domains, protein–rRNA recognition, etc. A correlation between folding rate of small globular proteins and their topological features is known. We have introduced an analogous topological characteristic for co-evolutionary network of ribosome, which allows us to differentiate between rRNA regions subjected to rapid reconstitutions from those hindered by kinetic traps. Furthermore, co-evolutionary patterns provide a biological basis for deleterious mutation sites and further allow prediction of potential antibiotic targeting sites. Understanding assembly pathways of multicomponent macromolecules remains a key challenge in biophysics. Our study provides a ‘proof of concept’ that directly relates co-evolution to biophysical interactions during multicomponent assembly and suggests predictive power to identify candidates for critical functional interactions as well as for assembly-blocking antibiotic target sites.

Description: The emergence of new sequencing technologies has facilitated the use of bacterial whole genome alignments for evolutionary studies and outbreak analyses. These datasets, of increasing size, often include examples of multiple different mechanisms of horizontal sequence transfer resulting in substantial alterations to prokaryotic chromosomes. The impact of these processes demands rapid and flexible approaches able to account for recombination when reconstructing isolates’ recent diversification. Gubbins is an iterative algorithm that uses spatial scanning statistics to identify loci containing elevated densities of base substitutions suggestive of horizontal sequence transfer while concurrently constructing a maximum likelihood phylogeny based on the putative point mutations outside these regions of high sequence diversity. Simulations demonstrate the algorithm generates highly accurate reconstructions under realistically parameterized models of bacterial evolution, and achieves convergence in only a few hours on alignments of hundreds of bacterial genome sequences. Gubbins is appropriate for reconstructing the recent evolutionary history of a variety of haploid genotype alignments, as it makes no assumptions about the underlying mechanism of recombination. The software is freely available for download at github.com/sanger-pathogens/Gubbins , implemented in Python and C and supported on Linux and Mac OS X.

Description: Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for ‘targeted’ resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a ‘kmer’ strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.

Description: A total of 65 spore-forming mercury-resistant bacteria were isolated from natural environments worldwide in order to understand the acquisition of additional genes by and dissemination of mercury resistance transposons across related Bacilli genera by horizontal gene movement. PCR amplification using a single primer complementary to the inverted repeat sequence of Tn MERI1 -like transposons showed that 12 of 65 isolates had a transposon-like structure. There were four types of amplified fragments: Tn 5084 , Tn 5085 , Tn d MER3 (a newly identified deleted transposon-like fragment) and Tn 6294 (a newly identified transposon). Tn d MER3 is a 3.5-kb sequence that carries a merRETPA operon with no merB or transposase genes. It is related to the mer operon of Bacillus licheniformis strain FA6-12 from Russia. DNA homology analysis shows that Tn 6294 is an 8.5-kb sequence that is possibly derived from Tn d MER3 by integration of a Tn MERI1 -type transposase and resolvase genes and in addition the merR2 and merB1 genes. Bacteria harboring Tn 6294 exhibited broad-spectrum mercury resistance to organomercurial compounds, although Tn 6294 had only merB1 and did not have the merB2 and merB3 sequences for organomercurial lyases found in Tn 5084 of B. cereus strain RC607. Strains with Tn 6294 encode mercuric reductase (MerA) of less than 600 amino acids in length with a single N-terminal mercury-binding domain, whereas MerA encoded by strains MB1 and RC607 has two tandem domains. Thus, Tn d MER3 and Tn 6294 are shorter prototypes for Tn MERI1 -like transposons. Identification of Tn 6294 in Bacillus sp. from Taiwan and in Paenibacillus sp. from Antarctica indicates the wide horizontal dissemination of Tn MERI1 -like transposons across bacterial species and geographical barriers.

Description: Fungi may play an important role in the production of the greenhouse gas nitrous oxide (N 2 O). Bipolaris sorokiniana is a ubiquitous saprobe found in soils worldwide, yet denitrification by this fungal strain has not previously been reported. We aimed to test if B. sorokiniana would produce N 2 O and CO 2 in the presence of organic and inorganic forms of nitrogen (N) under microaerobic and anaerobic conditions. Nitrogen source (organic-N, inorganic-N, no-N control) significantly affected N 2 O and CO 2 production both in the presence and absence of oxygen, which contrasts with bacterial denitrification. Inorganic N addition increased denitrification of N 2 O (from 0 to 0.3 μg N 2 0-N h –1  g –1 biomass) and reduced respiration of CO 2 (from 0.1 to 0.02 mg CO 2 h –1  g –1 biomass). Isotope analyses indicated that nitrite, rather than ammonium or glutamine, was transformed to N 2 O. Results suggest the source of N may play a larger role in fungal N 2 O production than oxygen status.

Description: Genetic variants in or near miRNA genes can have profound effects on miRNA expression and targeting. As user-friendly software for the impact prediction of miRNA variants on a large scale is still lacking, we created a tool called miRVaS. miRVaS automates this prediction by annotating the location of the variant relative to functional regions within the miRNA hairpin (seed, mature, loop, hairpin arm, flanks) and by annotating all predicted structural changes within the miRNA due to the variant. In addition, the tool defines the most important region that is predicted to have structural changes and calculates a conservation score that is indicative of the reliability of the structure prediction. The output is presented in a tab-separated file, which enables fast screening, and in an html file, which allows visual comparison between wild-type and variant structures. All separate images are provided for downstream use. Finally, we tested two different approaches on a small test set of published functionally validated genetic variants for their capacity to predict the impact of variants on miRNA expression.

Description: Analysis of RNA-seq data often detects numerous ‘non-co-linear’ (NCL) transcripts, which comprised sequence segments that are topologically inconsistent with their corresponding DNA sequences in the reference genome. However, detection of NCL transcripts involves two major challenges: removal of false positives arising from alignment artifacts and discrimination between different types of NCL transcripts ( trans -spliced, circular or fusion transcripts). Here, we developed a new NCL-transcript-detecting method (‘NCLscan’), which utilized a stepwise alignment strategy to almost completely eliminate false calls (〉98% precision) without sacrificing true positives, enabling NCLscan outperform 18 other publicly-available tools (including fusion- and circular-RNA-detecting tools) in terms of sensitivity and precision, regardless of the generation strategy of simulated dataset, type of intragenic or intergenic NCL event, read depth of coverage, read length or expression level of NCL transcript. With the high accuracy, NCLscan was applied to distinguishing between trans -spliced, circular and fusion transcripts on the basis of poly(A)- and nonpoly(A)-selected RNA-seq data. We showed that circular RNAs were expressed more ubiquitously, more abundantly and less cell type-specifically than trans -spliced and fusion transcripts. Our study thus describes a robust pipeline for the discovery of NCL transcripts, and sheds light on the fundamental biology of these non-canonical RNA events in human transcriptome.

Description: Legionella pneumophila is a pathogenic bacterium commonly found in water and responsible for severe pneumonia. Free-living amoebae are protozoa also found in water, which feed on bacteria by phagocytosis. Under favorable conditions, some L. pneumophila are able to resist phagocytic digestion and even multiply within amoebae. However, it is not clear whether L. pneumophila could infect at a same rate a large range of amoebae or if there is some selectivity towards specific amoebal genera or strains. Also, most studies have been performed using collection strains and not with freshly isolated strains. In our study, we assess the permissiveness of freshly isolated environmental strains of amoebae, belonging to three common genera (i.e. Acanthamoeba, Naegleria and Vermamoeba ), for growth of L. pneumophila at three different temperatures. Our results indicated that all the tested strains of amoebae were permissive to L. pneumophila Lens and that there was no significant difference between the strains. Intracellular proliferation was more efficient at a temperature of 40°C. In conclusion, our work suggests that, under favorable conditions, virulent strains of L. pneumophila could equally infect a large number of isolates of common freshwater amoeba genera.

Description: Detecting allelic biases from high-throughput sequencing data requires an approach that maximises sensitivity while minimizing false positives. Here, we present Allelome.PRO, an automated user-friendly bioinformatics pipeline, which uses high-throughput sequencing data from reciprocal crosses of two genetically distinct mouse strains to detect allele-specific expression and chromatin modifications. Allelome.PRO extends approaches used in previous studies that exclusively analyzed imprinted expression to give a complete picture of the ‘allelome’ by automatically categorising the allelic expression of all genes in a given cell type into imprinted, strain-biased, biallelic or non-informative. Allelome.PRO offers increased sensitivity to analyze lowly expressed transcripts, together with a robust false discovery rate empirically calculated from variation in the sequencing data. We used RNA-seq data from mouse embryonic fibroblasts from F1 reciprocal crosses to determine a biologically relevant allelic ratio cutoff, and define for the first time an entire allelome. Furthermore, we show that Allelome.PRO detects differential enrichment of H3K4me3 over promoters from ChIP-seq data validating the RNA-seq results. This approach can be easily extended to analyze histone marks of active enhancers, or transcription factor binding sites and therefore provides a powerful tool to identify candidate cis regulatory elements genome wide.

Description: Identity by descent (IBD) can be reliably detected for long shared DNA segments, which are found in related individuals. However, many studies contain cohorts of unrelated individuals that share only short IBD segments. New sequencing technologies facilitate identification of short IBD segments through rare variants, which convey more information on IBD than common variants. Current IBD detection methods, however, are not designed to use rare variants for the detection of short IBD segments. Short IBD segments reveal genetic structures at high resolution. Therefore, they can help to improve imputation and phasing, to increase genotyping accuracy for low-coverage sequencing and to increase the power of association studies. Since short IBD segments are further assumed to be old, they can shed light on the evolutionary history of humans. We propose HapFABIA, a computational method that applies biclustering to identify very short IBD segments characterized by rare variants. HapFABIA is designed to detect short IBD segments in genotype data that were obtained from next-generation sequencing, but can also be applied to DNA microarray data. Especially in next-generation sequencing data, HapFABIA exploits rare variants for IBD detection. HapFABIA significantly outperformed competing algorithms at detecting short IBD segments on artificial and simulated data with rare variants. HapFABIA identified 160 588 different short IBD segments characterized by rare variants with a median length of 23 kb (mean 24 kb) in data for chromosome 1 of the 1000 Genomes Project. These short IBD segments contain 752 000 single nucleotide variants (SNVs), which account for 39% of the rare variants and 23.5% of all variants. The vast majority—152 000 IBD segments—are shared by Africans, while only 19 000 and 11 000 are shared by Europeans and Asians, respectively. IBD segments that match the Denisova or the Neandertal genome are found significantly more often in Asians and Europeans but also, in some cases exclusively, in Africans. The lengths of IBD segments and their sharing between continental populations indicate that many short IBD segments from chromosome 1 existed before humans migrated out of Africa. Thus, rare variants that tag these short IBD segments predate human migration from Africa. The software package HapFABIA is available from Bioconductor. All data sets, result files and programs for data simulation, preprocessing and evaluation are supplied at http://www.bioinf.jku.at/research/short-IBD .

Description: Glioma is the most common and fatal primary brain tumour with poor prognosis; however, the functional roles of miRNAs in glioma malignant progression are insufficiently understood. Here, we used an integrated approach to identify miRNA functional targets during glioma malignant progression by combining the paired expression profiles of miRNAs and mRNAs across 160 Chinese glioma patients, and further constructed the functional miRNA–mRNA regulatory network. As a result, most tumour-suppressive miRNAs in glioma progression were newly discovered, whose functions were widely involved in gliomagenesis. Moreover, three miRNA signatures, with different combinations of hub miRNAs (regulations≥30) were constructed, which could independently predict the survival of patients with all gliomas, high-grade glioma and glioblastoma. Our network-based method increased the ability to identify the prognostic biomarkers, when compared with the traditional method and random conditions. Hsa-miR-524-5p and hsa-miR-628-5p, shared by these three signatures, acted as protective factors and their expression decreased gradually during glioma progression. Functional analysis of these miRNA signatures highlighted their critical roles in cell cycle and cell proliferation in glioblastoma malignant progression, especially hsa-miR-524-5p and hsa-miR-628-5p exhibited dominant regulatory activities. Therefore, network-based biomarkers are expected to be more effective and provide deep insights into the molecular mechanism of glioma malignant progression.

Description: A major challenge in cancer genomics is uncovering genes with an active role in tumorigenesis from a potentially large pool of mutated genes across patient samples. Here we focus on the interactions that proteins make with nucleic acids, small molecules, ions and peptides, and show that residues within proteins that are involved in these interactions are more frequently affected by mutations observed in large-scale cancer genomic data than are other residues. We leverage this observation to predict genes that play a functionally important role in cancers by introducing a computational pipeline ( http://canbind.princeton.edu ) for mapping large-scale cancer exome data across patients onto protein structures, and automatically extracting proteins with an enriched number of mutations affecting their nucleic acid, small molecule, ion or peptide binding sites. Using this computational approach, we show that many previously known genes implicated in cancers are enriched in mutations within the binding sites of their encoded proteins. By focusing on functionally relevant portions of proteins—specifically those known to be involved in molecular interactions—our approach is particularly well suited to detect infrequent mutations that may nonetheless be important in cancer, and should aid in expanding our functional understanding of the genomic landscape of cancer.

Description: Genome duplication with hybridization, or allopolyploidization, occurs commonly in plants, and is considered to be a strong force for generating new species. However, genome-wide quantification of homeolog expression ratios was technically hindered because of the high homology between homeologous gene pairs. To quantify the homeolog expression ratio using RNA-seq obtained from polyploids, a new method named HomeoRoq was developed, in which the genomic origin of sequencing reads was estimated using mismatches between the read and each parental genome. To verify this method, we first assembled the two diploid parental genomes of Arabidopsis halleri subsp. gemmifera and Arabidopsis lyrata subsp. petraea ( Arabidopsis petraea subsp. umbrosa ), then generated a synthetic allotetraploid, mimicking the natural allopolyploid Arabidopsis kamchatica . The quantified ratios corresponded well to those obtained by Pyrosequencing. We found that the ratios of homeologs before and after cold stress treatment were highly correlated ( r = 0.870). This highlights the presence of nonstochastic polyploid gene regulation despite previous research identifying stochastic variation in expression. Moreover, our new statistical test incorporating overdispersion identified 226 homeologs (1.11% of 20 369 expressed homeologs) with significant ratio changes, many of which were related to stress responses. HomeoRoq would contribute to the study of the genes responsible for polyploid-specific environmental responses.

Description: Methods to interpret personal genome sequences are increasingly required. Here, we report a novel framework (EvoTol) to identify disease-causing genes using patient sequence data from within protein coding-regions. EvoTol quantifies a gene's intolerance to mutation using evolutionary conservation of protein sequences and can incorporate tissue-specific gene expression data. We apply this framework to the analysis of whole-exome sequence data in epilepsy and congenital heart disease, and demonstrate EvoTol's ability to identify known disease-causing genes is unmatched by competing methods. Application of EvoTol to the human interactome revealed networks enriched for genes intolerant to protein sequence variation, informing novel polygenic contributions to human disease.

Description: Detecting in vivo transcription factor (TF) binding is important for understanding gene regulatory circuitries. ChIP-seq is a powerful technique to empirically define TF binding in vivo . However, the multitude of distinct TFs makes genome-wide profiling for them all labor-intensive and costly. Algorithms for in silico prediction of TF binding have been developed, based mostly on histone modification or DNase I hypersensitivity data in conjunction with DNA motif and other genomic features. However, technical limitations of these methods prevent them from being applied broadly, especially in clinical settings. We conducted a comprehensive survey involving multiple cell lines, TFs, and methylation types and found that there are intimate relationships between TF binding and methylation level changes around the binding sites. Exploiting the connection between DNA methylation and TF binding, we proposed a novel supervised learning approach to predict TF–DNA interaction using data from base-resolution whole-genome methylation sequencing experiments. We devised beta-binomial models to characterize methylation data around TF binding sites and the background. Along with other static genomic features, we adopted a random forest framework to predict TF–DNA interaction. After conducting comprehensive tests, we saw that the proposed method accurately predicts TF binding and performs favorably versus competing methods.

Description: Enriching target sequences in sequencing libraries via capture hybridization to bait/probes is an efficient means of leveraging the capabilities of next-generation sequencing for obtaining sequence data from target regions of interest. However, homologous sequences from non-target regions may also be enriched by such methods. Here we investigate the fidelity of capture enrichment for complete mitochondrial DNA (mtDNA) genome sequencing by analyzing sequence data for nuclear copies of mtDNA (NUMTs). Using capture-enriched sequencing data from a mitochondria-free cell line and the parental cell line, and from samples previously sequenced from long-range PCR products, we demonstrate that NUMT alleles are indeed present in capture-enriched sequence data, but at low enough levels to not influence calling the authentic mtDNA genome sequence. However, distinguishing NUMT alleles from true low-level mutations (e.g. heteroplasmy) is more challenging. We develop here a computational method to distinguish NUMT alleles from heteroplasmies, using sequence data from artificial mixtures to optimize the method.

Description: Single nucleotide polymorphisms (SNPs) are increasingly used to tag genetic loci associated with phenotypes such as risk of complex diseases. Technically, this is done genome-wide without prior restriction or knowledge of biological feasibility in scans referred to as genome-wide association studies (GWAS). Depending on the linkage disequilibrium (LD) structure at a particular locus, such tagSNPs may be surrogates for many thousands of other SNPs, and it is difficult to distinguish those that may play a functional role in the phenotype from those simply genetically linked. Because a large proportion of tagSNPs have been identified within non-coding regions of the genome, distinguishing functional from non-functional SNPs has been an even greater challenge. A strategy was recently proposed that prioritizes surrogate SNPs based on non-coding chromatin and epigenomic mapping techniques that have become feasible with the advent of massively parallel sequencing. Here, we introduce an R/Bioconductor software package that enables the identification of candidate functional SNPs by integrating information from tagSNP locations, lists of linked SNPs from the 1000 genomes project and locations of chromatin features which may have functional significance. Availability: FunciSNP is available from Bioconductor (bioconductor.org).

Description: Reconstructing the evolutionary relationships of species is a major goal in biology. Despite the increasing number of completely sequenced genomes, a large number of phylogenetic projects rely on targeted sequencing and analysis of a relatively small sample of marker genes. The selection of these phylogenetic markers should ideally be based on accurate predictions of their combined, rather than individual, potential to accurately resolve the phylogeny of interest. Here we present and validate a new phylogenomics strategy to efficiently select a minimal set of stable markers able to reconstruct the underlying species phylogeny. In contrast to previous approaches, our methodology does not only rely on the ability of individual genes to reconstruct a known phylogeny, but it also explores the combined power of sets of concatenated genes to accurately infer phylogenetic relationships of species not previously analyzed. We applied our approach to two broad sets of cyanobacterial and ascomycetous fungal species, and provide two minimal sets of six and four genes, respectively, necessary to fully resolve the target phylogenies. This approach paves the way for the informed selection of phylogenetic markers in the effort of reconstructing the tree of life.

Description: Visualization of chromosomal dynamics is important for understanding many fundamental intra-nuclear processes. Efficient and reliable live-cell multicolor labeling of chromosomal loci can realize this goal. However, the current methods are constrained mainly by insufficient labeling throughput, efficiency, flexibility as well as photostability. Here we have developed a new approach to realize dual-color chromosomal loci imaging based on a modified single-guide RNA (sgRNA) of the CRISPR/Cas9 system. The modification of sgRNA was optimized by structure-guided engineering of the original sgRNA, consisting of RNA aptamer insertions that bind fluorescent protein-tagged effectors. By labeling and tracking telomeres, centromeres and genomic loci, we demonstrate that the new approach is easy to implement and enables robust dual-color imaging of genomic elements. Importantly, our data also indicate that the fast exchange rate of RNA aptamer binding effectors makes our sgRNA-based labeling method much more tolerant to photobleaching than the Cas9-based labeling method. This is crucial for continuous, long-term tracking of chromosomal dynamics. Lastly, as our method is complementary to other live-cell genomic labeling systems, it is therefore possible to combine them into a plentiful palette for the study of native chromatin organization and genome ultrastructure dynamics in living cells.

Description: Growth media have been developed to facilitate the enrichment and isolation of acidophilic and acid-tolerant sulfate-reducing bacteria (aSRB) from environmental and industrial samples, and to allow their cultivation in vitro . The main features of the ‘standard’ solid and liquid devised media are as follows: (i) use of glycerol rather than an aliphatic acid as electron donor; (ii) inclusion of stoichiometric concentrations of zinc ions to both buffer pH and to convert potentially harmful hydrogen sulphide produced by the aSRB to insoluble zinc sulphide; (iii) inclusion of Acidocella aromatica (an heterotrophic acidophile that does not metabolize glycerol or yeast extract) in the gel underlayer of double layered (overlay) solid media, to remove acetic acid produced by aSRB that incompletely oxidize glycerol and also aliphatic acids (mostly pyruvic) released by acid hydrolysis of the gelling agent used (agarose). Colonies of aSRB are readily distinguished from those of other anaerobes due to their deposition and accumulation of metal sulphide precipitates. Data presented illustrate the effectiveness of the overlay solid media described for isolating aSRB from acidic anaerobic sediments and low pH sulfidogenic bioreactors.

Description: A common dye of prussian blue (PB) as an indicator was used to develop a colorimetric method for detecting the efficacy of the antibiotics in vitro. Considering the electronic production capacity of microbial respiration, ferricyanide was employed in transferring electrons from target microorganism of Escherichia coli ( E. coli ) to produce ferrocyanide. Subsequently, ferrocyanide reacted with ferric ions to form PB. In view of relationship between the PB yield and the bacterial activity, the efficacy of the antibiotics on E. coli was directly detected at 700 nm of PB absorption. When the 5% activity of antibiotics on 20 isolates of E. coli was quantified as 5% efficacy, the applied concentrations of eight antibiotics, such as cefepime, ceftriaxone sodium, cefoperazone sodium, piperacillin sodium, amoxicillin, gentamicin, amikacin and levofloxacin were 2, 2, 4, 4, 10, 4, 8 and 8 μg mL –1 , respectively. To compare with minimum inhibitory concentration results obtained by Clinical and Laboratory Standards Institute broth macrodilution method, the results of PB methods showed good agreements except with gentamicin. Paired t- test result ( P ) also showed that difference between two methods was statistically significant ( P = 0.006).

Description: Metarhizium acridum is an entomopathogenic fungus commonly used as a bioinsecticide. The conidium is the fungal stage normally employed as field inoculum in biological control programs and must survive under field conditions such as high ultraviolet-B (UV-B) exposure. Light, which is an important stimulus for many fungi, has been shown to induce the production of M. robertsii conidia with increased stress tolerance. Here we show that a two-hour exposure to white or blue/UV-A light of fast-growing mycelium induces tolerance to subsequent UV-B irradiation. Red light, however, does not have the same effect. In addition, we established that this induction can take place with as little as 1 min of white-light exposure. This brief illumination scheme could be relevant in future studies of M. acridum photobiology and for the production of UV-B resistant mycelium used in mycelium-based formulations for biological control.

Description: Photorhabdus (Enterobacteriaceae) bacteria are pathogenic to insects and mutualistic with entomopathogenic Heterorhabditis nematodes . Photorhabdus luminescens subsp. akhurstii LN2, associated with Heterorhabditis indica LN2, shows nematicidal activity against H. bacteriophora H06 infective juveniles (IJs). In the present study, an rpoS mutant of P. luminescens LN2 was generated through allelic exchange to examine the effects of rpoS deletion on the nematicidal activity and nematode development. The results showed that P. luminescens LN2 required rpoS for nematicidal activity against H06 nematodes, normal IJ recovery and development of H. indica LN2, however, not for the bacterial colonization in LN2 and H06 IJs. This provides cues for further understanding the role of rpoS in the mutualistic association between entomopathogenic nematodes and their symbionts.

Description: Catechol 2, 3-dioxygenase (C23O) is the key enzyme for aerobic aromatic degradation. Based on clone libraries and quantitative real-time polymerase chain reaction, we characterized diversity and distribution patterns of C23O genes in surface sediments of the Bohai Sea. The results showed that sediments of the Bohai Sea were dominated by genes related to C23O subfamily I.2.A. The samples from wastewater discharge area (DG) and aquaculture farm (KL) showed distinct composition of C23O genes when compared to the samples from Bohai Bay (BH), and total organic carbon was a crucial determinant accounted for the composition variation. C6BH12-38 and C2BH2-35 displayed the highest gene copies and highest ratios to the 16S rRNA genes in KL, and they might prefer biologically labile aromatic hydrocarbons via aquaculture inputs. Meanwhile, C7BH3-48 showed the highest gene copies and highest ratios to the 16S rRNA genes in DG, and this could be selective effect of organic loadings from wastewater discharge. An evident increase in C6BH12-38 and C7BH3-48 gene copies and reduction in diversity of C23O genes in DG and KL indicated composition perturbations of C23O genes and potential loss in functional redundancy. We suggest that ecological habitat and trophic specificity could shape the distribution of C23O genes in the Bohai Sea sediments.

Description: LysR-type transcriptional regulators (LTTRs) regulate various cellular processes in bacteria. pnpR is an LTTR-encoding gene involved in the regulation of hydroquinone (HQ) degradation, and its effects on the cellular processes of Pseudomonas putida DLL-E4 were investigated at the physiological, biochemical and molecular levels. Reverse transcription polymerase chain reaction revealed that pnpR positively regulated its own expression and that of the pnpC1C2DECX1X2 operon; additionally, pnpR partially regulated the expression of pnpA when P. putida was grown on para -nitrophenol (PNP) or HQ. Strains DLL-E4 and DLL- pnpR exhibited similar cellular morphologies and growth rates. Transcriptome analysis revealed that pnpR regulated the expression of genes in addition to those involved in PNP degradation. A total of 20 genes were upregulated and 19 genes were downregulated by at least 2-fold in strain DLL- pnpR relative to strain DLL-E4. Bioinformatic analysis revealed putative PnpR-binding sites located in the upstream regions of genes involved in PNP degradation, carbon catabolite repression and other cellular processes. The utilization of L-aspartic acid, L-histidine, L-pyroglutamic acid, L-serine, -aminobutyric acid, D,L-lactic acid, D-saccharic acid, succinic acid and L-alaninamide was increased at least 1.3-fold in strain DLL- pnpR as shown by BIOLOG assays, indicating that pnpR plays a potential negative regulation role in the utilization of carbon sources.

Description: Landfills are significant global sources of atmospheric methane, but little is known about the ecology and community structure of methanogens in these sites. Here, we investigated the methanogen community based on methyl coenzyme M reductase A gene amplicons in the vertical profiles of three different sites at a municipal landfill complex in China. Links between methanogen communities and refuse properties were explored using multivariate analysis. Clone library results showed that most clones (92%) were related to the hydrogenotrophic methanogens, Methanomicrobiales. Almost all of the Methanomicrobiales clones retrieved in this study are members of the genus Methanoculleus . Eight clones were affiliated with the genus Methanofollis . The remaining clones were clustered within the genus Methanosarcina . Terminal restriction fragment length polymorphism profiles showed that the landfill was predominated by 22 taxa, making up 69%–96% of the community. Of these, a single taxon comprised 36%–65% of the communities across all sites and depths. Principal components analysis separated the methanogen community into three groups, irrespective of site or depth. Redundancy analysis suggested that total phosphorus and pH play roles in structuring methanogen communities in landfills.

Description: Here we report a newly identified ‘Chalky back’ phenomenon in banana prawns ( Fenneropenaeus merguiensis ) farmed in North Queensland, Australia. This was characterized by localized white discoloured segmentation of the cervical groove, moreover, after cooking the prawns exploded, making them unfit for commercial sale. Histological examination revealed breakdown of gut and abdominal muscle tissue in some moribund specimens. We selectively isolated Vibrio spp., which are known prawn pathogens, from healthy and Chalky back specimens. Isolated bacteria were identified, typed and tested for the presence of eight virulence genes (VGs), biofilm formation, adherence and cytotoxicity to fish cells. In all, 32 isolates were recovered and identified as Vibrio harveyi , V. owensii , V. sinaloensis -like, V. campbellii , V. shilonii , Vibrio sp. and Photobacterium damselae using 16S rRNA gene sequencing. All V. harveyi carried VGs coding for haemolysin, tox R and flagella; formed biofilm; and adhered to both cell lines. This was similar to the V. sinaloensis -like strains that were only isolated from Chalky back specimens. Our data suggest that Vibrio spp. may play a role in the pathogenesis of Chalky back. This study is the first report of Chalky back phenomenon in farmed banana prawns that needs to be closely monitored by the industry.

Description: Environmental biosurveillance and microbial ecology studies use PCR-based assays to detect and quantify microbial taxa and gene sequences within a complex background of microorganisms. However, the fragmentary nature and growing quantity of DNA-sequence data make group-specific assay design challenging. We solved this problem by developing a software platform that enables PCR-assay design at an unprecedented scale. As a demonstration, we developed quantitative PCR assays for a globally widespread, ecologically important bacterial group in soil, Acidobacteria Group 1. A total of 33 684 Acidobacteria 16S rRNA gene sequences were used for assay design. Following 1 week of computation on a 376-core cluster, 83 assays were obtained. We validated the specificity of the top three assays, collectively predicted to detect 42% of the Acidobacteria Group 1 sequences, by PCR amplification and sequencing of DNA from soil. Based on previous analyses of 16S rRNA gene sequencing, Acidobacteria Group 1 species were expected to decrease in response to elevated atmospheric CO 2 . Quantitative PCR results, using the Acidobacteria Group 1-specific PCR assays, confirmed the expected decrease and provided higher statistical confidence than the 16S rRNA gene-sequencing data. These results demonstrate a powerful capacity to address previously intractable assay design challenges.

Description: Due to dramatic advances in DNA technology, quantitative measures of annotation data can now be obtained in continuous coordinates across the entire genome, allowing various heterogeneous ‘genomic landscapes’ to emerge. Although much effort has been devoted to comparing DNA sequences, not much attention has been given to comparing these large quantities of data comprehensively. In this article, we introduce a method for rapidly detecting local regions that show high correlations between genomic landscapes. We overcame the size problem for genome-wide data by converting the data into series of symbols and then carrying out sequence alignment. We also decomposed the oscillation of the landscape data into different frequency bands before analysis, since the real genomic landscape is a mixture of embedded and confounded biological processes working at different scales in the cell nucleus. To verify the usefulness and generality of our method, we applied our approach to well investigated landscapes from the human genome, including several histone modifications. Furthermore, by applying our method to over 20 genomic landscapes in human and 12 in mouse, we found that DNA replication timing and the density of Alu insertions are highly correlated genome-wide in both species, even though the Alu elements have amplified independently in the two genomes. To our knowledge, this is the first method to align genomic landscapes at multiple scales according to their shape.

Description: The numerous discovered cases of domesticated transposable element (TE) proteins led to the recognition that TEs are a significant source of evolutionary innovation. However, much less is known about the reverse process, whether and to what degree the evolution of TEs is influenced by the genome of their hosts. We addressed this issue by searching for cases of incorporation of host genes into the sequence of TEs and examined the systems-level properties of these genes using the Saccharomyces cerevisiae and Drosophila melanogaster genomes. We identified 51 cases where the evolutionary scenario was the incorporation of a host gene fragment into a TE consensus sequence, and we show that both the yeast and fly homologues of the incorporated protein sequences have central positions in the cellular networks. An analysis of selective pressure (Ka/Ks ratio) detected significant selection in 37% of the cases. Recent research on retrovirus-host interactions shows that virus proteins preferentially target hubs of the host interaction networks enabling them to take over the host cell using only a few proteins. We propose that TEs face a similar evolutionary pressure to evolve proteins with high interacting capacities and take some of the necessary protein domains directly from their hosts.

Description: Plant-growth-promoting bacteria belonging to Azospirillum and Pseudomonas genera are major inhabitants of the rhizosphere. Both are increasingly commercialized as crops inoculants. Interspecific interaction in the rhizosphere is critical for inoculants aptness. The objective of this work was to evaluate Azospirillum and Pseudomonas interaction in mixed biofilms by co-cultivation of the model strains Azospirillum brasilense Sp245 and Pseudomonas protegens CHA0. The results revealed enhanced growth of both strains when co-cultured in static conditions. Moreover, Sp245 biofilm formed in plastic surfaces was increased 2-fold in the presence of CHA0. Confocal microscopy revealed highly structured mixed biofilms showing Sp245 mainly on the bottom and CHA0 towards the biofilm surface. In addition, A. brasilense biofilm was thicker and denser when co-cultured with P. protegens. In a colony–colony interaction assay, Sp245 changed nearby CHA0 producing small colony phenotype, which accounts for a diffusible metabolite mediator; though CHA0 spent medium did not affect Sp245 colony phenotype. Altogether, these results point to a cooperative interaction between A. brasilense Sp245 and P. protegens CHA0 in which both strains increase their static growth and produce structured mixed biofilms with a strain-specific distribution.

Description: Presence of excess unaltered, wild-type (WT) DNA providing no information of biological or clinical value often masks rare alterations containing diagnostic or therapeutic clues in cancer, prenatal diagnosis, infectious diseases or organ transplantation. With the surge of high-throughput technologies there is a growing demand for removing unaltered DNA over large pools-of-sequences. Here we present nuclease-assisted minor-allele enrichment with probe-overlap (NaME-PrO), a single-step approach with broad genome coverage that can remove WT-DNA from numerous sequences simultaneously, prior to genomic analysis. NaME-PrO employs a double-strand-DNA-specific nuclease and overlapping oligonucleotide-probes interrogating WT-DNA targets and guiding nuclease digestion to these sites. Mutation-containing DNA creates probe-DNA mismatches that inhibit digestion, thus subsequent DNA-amplification magnifies DNA-alterations at all selected targets. We demonstrate several-hundred-fold mutation enrichment in diverse human samples on multiple clinically relevant targets including tumor samples and circulating DNA in 50-plex reactions. Enrichment enables routine mutation detection at 0.01% abundance while by adjusting conditions it is possible to sequence mutations down to 0.00003% abundance, or to scan tumor-suppressor genes for rare mutations. NaME-PrO introduces a simple and highly parallel process to remove un-informative DNA sequences and unmask clinically and biologically useful alterations.

Description: Knowledge about the factors shaping the rumen microbiota in wild animals is limited. Therefore, the aim of this study was to compare the microbiota from the three cervid species moose ( Alces alces , n = 5), red deer ( Cervus elaphus , n = 4) and roe deer ( Capreolus capreolus , n = 12), sharing the same habitat. Using deep 16S rRNA gene sequencing, we found that the largest species moose had the highest number of unique operational taxonomic units. Furthermore, red deer and moose shared more of the microbiota, compared with the smallest species, roe deer, with Firmicutes and Euryarchaeota being significantly overrepresented for the shared microbiota. These differences could not be explained by diet or range. The animals largely shared the same range, and there are no systematic differences in diet. We therefore believe rumen physiology can be one of the main contributing factors to the observed distribution of the rumen microbiota in cervid species.

Description: Mesorhizobium loti MAFF303099 has a functional Type III secretion system (T3SS) that is involved in the determination of competitiveness for legume nodulation. Here we demonstrate that the transcriptional factor TtsI, which positively regulates T3SS genes expression, is involved in a negative regulation of M. loti swimming motility in soft-agar. Conditions that induce T3SS expression affect flagella production. The same conditions also affect promoter activity of M. loti visN gene, a homolog to the positive regulator of flagellar genes that has been described in other rhizobia. Defects in T3SS complex assembly at membranes limited the negative regulation of motility by the expression of TtsI.

Description: Frequent burning is commonly undertaken to maintain diversity in temperate grasslands of southern Australia. How burning affects below-ground fungal community diversity remains unknown. We show, using a fungal rDNA metabarcoding approach (Illumina MiSeq), that the fungal community composition was influenced by fire regime (frequency) but not time-since-fire. Fungal community composition was resilient to direct fire effects, most likely because grassland fires transfer little heat to the soil. Differences in the fungal community composition due to fire regime was likely due to associated changes that occur in vegetation with recurrent fire, via the break up of obligate symbiotic relationships. However, fire history only partially explains the observed dissimilarity in composition among the soil samples, suggesting a distinctiveness in composition in each grassland site. The importance of considering changes in soil microbe communities when managing vegetation with fire is highlighted.

Description: Soil is thought to be important both as a source and a sink of carbonyl sulfide (COS) in the troposphere, but the mechanism affecting COS uptake, especially for fungi, remains uncertain. Fungal isolates that were collected randomly from forest soil showed COS-degrading ability at high frequencies: 38 out of 43 isolates grown on potato dextrose agar showed degradation of 30 ppmv COS within 24 h. Of these isolates, eight degraded 30 ppmv of COS to below the detection limit within 2 h. These isolates also showed an ability to degrade COS included in ambient air (around 500 pptv) and highly concentrated (12 500 ppmv) level, even though the latter is higher than the lethal level for mammals. COS-degrading activity was estimated by using ergosterol as a biomass index for fungi. Trichoderma sp. THIF08 had the highest COS-degrading activity of all the isolates. Interestingly, Umbelopsis/Mortierella spp. THIF09 and THIF13 were unable to degrade 30 ppmv COS within 24 h, and actually emitted COS during the cultivation in ambient air. These results indicate a fungal contribution to the flux of COS between the terrestrial and atmospheric environments.

Description: A composite transposon is a mobile genetic element consisting of two insertion sequences (ISs) flanking a segment of cargo DNA often containing antibiotic resistance (AR) genes. Composite transposons can move as a discreet unit. There have been recently several reports on a novel mechanism of movement of an IS 26 -based composite transposon through the formation of a translocatable unit (TU), carrying the internal DNA segment of a composite transposon and one copy of a flanking IS. In this study, we determined the presence of composite transposons and TUs in human oral metagenomic DNA using PCR primers from common IS elements. Analysis of resulting amplicons showed four different IS 1216 composite transposons and one IS 257 composite transposon in our metagenomic sample. As our PCR strategy would also detect TUs, PCR was carried out to detect circular TUs predicted to originate from these composite transposons. We confirmed the presence of two novel TUs, one containing an experimentally proven antiseptic resistance gene and another containing a putative universal stress response protein (UspA) encoding gene. This is the first report of a PCR strategy to amplify the DNA segment on composite transposons and TUs in metagenomic DNA. This can be used to identify AR genes associated with a variety of mobile genetic elements from metagenomes.

Description: Carbonyl sulfide (COS) is an atmospheric trace gas and one of the sources of stratospheric aerosol contributing to climate change. Although one of the major sinks of COS is soil, the distribution of COS degradation ability among bacteria remains unclear. Seventeen out of 20 named bacteria belonging to Actinomycetales had COS degradation activity at mole fractions of 30 parts per million by volume (ppmv) COS. Dietzia maris NBRC 15801 T and Mycobacterium sp. THI405 had the activity comparable to a chemolithoautotroph Thiobacillus thioparus THI115 that degrade COS by COS hydrolase for energy production. Among 12 bacteria manifesting rapid degradation at 30 ppmv COS, D. maris NBRC 15801 T and Streptomyces ambofaciens NBRC 12836 T degraded ambient COS (~500 parts per trillion by volume). Geodermatophilus obscurus NBRC 13315 T and Amycolatopsis orientalis NBRC 12806 T increased COS concentrations. Moreover, six of eight COS-degrading bacteria isolated from soils had partial nucleotide sequences similar to that of the gene encoding clade D of β-class carbonic anhydrase, which included COS hydrolase. These results indicate the potential importance of Actinomycetes in the role of soils as sinks of atmospheric COS.

Description: Neonicotinoids are neurotoxic systemic insecticides used in plant protection worldwide. Unfortunately, application of neonicotinoids affects both beneficial and target insects indiscriminately. Being water soluble and persistent, these pesticides are capable of disrupting both food chains and biogeochemical cycles. This review focuses on the biodegradation of neonicotinoids in soil and water systems by the bacterial community. Several bacterial strains have been isolated and identified as capable of transforming neonicotinoids in the presence of an additional carbon source. Environmental parameters have been established for accelerated transformation in some of these strains. Studies have also indicated that enhanced biotransformation of these pesticides can be accomplished by mixed microbial populations under optimised environmental conditions. Substantial research into the identification of neonicotinoid-mineralising bacterial strains and identification of the genes and enzymes responsible for neonicotinoid degradation is still required to complete the understanding of microbial biodegradation pathways, and advance bioremediation efforts.

Description: Many toxic insecticides used worldwide as well as some chemical warfare agents are phosphotriester derivatives. Therefore, detoxification of organophosphorus compounds has become the subject of many studies and in particular bioremediation, based on the phosphotriesterase catalysed hydrolysis of these compounds, has shown to be an effective and ecological methodology. In order to identify new bacterial phosphotriesterases, a simple and sensitive fluorimetric screening method on solid media was employed that allowed the selection of six strains with phosphotriesterase activity. Since pH and temperature are important parameters for bioremediation of contaminated soils and waters, the influence of these variables on the rate of the enzymatic hydrolysis was assessed. This study afforded notable results, being the most remarkable one the increased activity exhibited by Nocardia asteroides and Streptomyces setonii strains at 50°C, 7 and 30 times higher than at 30°C, respectively. Compared with the results obtained with Brevundimonas diminuta , whose activity is usually considered as reference, an increase of 26 and 75 times is observed, respectively.

Description: Phosphorus (P) is a critical, non-renewable nutrient; yet excess discharges can lead to eutrophication and deterioration of water quality. Thus, P removal from water must be coupled with P recovery to achieve sustainable P management. P-specific proteins provide a novel, promising approach to recover P from water. Bacterial phosphate-binding proteins (PBPs) are able to effectively remove phosphate, achieving extremely low levels in water (i.e. 0.015 mg-P L –1 ). A prerequisite of using PBP for P recovery, however, is not only removal, but also controlled P release, which has not yet been reported. Phosphate release using recombinant PBP-expressing Escherichia coli was explored in this study. Escherichia coli was genetically modified to overexpress PBP in the periplasmic space. The impacts of ionic strength, temperature and pH on phosphate release were assessed. PBP-expressed E. coli demonstrated consistently superior ability to adsorb more phosphate from liquid and release more phosphate under controlled conditions relative to negative controls (unexpressed PBP E. coli and E. coli K12). Lower pH (3.8), higher temperature (35ºC) and higher ionic strength (100 mM KCl) facilitated increased phosphate release, providing a maximum of 2.1% P recovery within 3 h. This study provides proof of concept of the feasibility of using PBP to recover P.

Description: A number of studies have shown that transcriptome analysis in terms of chromosomal location can reveal regions of non-random transcriptional activity within the genome. Genomic clusters of differentially expressed genes can identify genomic patterns of structural organization, underlying copy number variations or long-range epigenetic regulation such as X-chromosome inactivation. Here we apply an integrative bioinformatics analysis to a collection of 315 freely available mouse pluripotent stem cell samples to discover transcriptional clusters in the genome. We show that over half of the analysed samples (56.83%) carry whole or partial-chromosome spanning clusters which recur in genomic regions previously implicated in chromosomal imbalances. Strikingly, we found that the presence of such large-clusters is linked to the differential expression of a limited number of genes, common to all samples carrying clusters irrespectively of the chromosome where the cluster is found. We have used these genes to train and test classification models that can predict samples that carry large-scale clusters on any chromosome with over 90% accuracy. Our findings suggest that there is a common downstream activation in these cells that affects a limited number of nodes. We propose that this effect is linked to selective advantage and identify potential driver genes.

Description: The rapidly growing amount of genomic sequence data being generated and made publicly available necessitate the development of new data storage and archiving methods. The vast amount of data being shared and manipulated also create new challenges for network resources. Thus, developing advanced data compression techniques is becoming an integral part of data production and analysis. The HapMap project is one of the largest public resources of human single-nucleotide polymorphisms (SNPs), characterizing over 3 million SNPs genotyped in over 1000 individuals. The standard format and biological properties of HapMap data suggest that a dedicated genetic compression method can outperform generic compression tools. We propose a compression methodology for genetic data by introducing H ap Z ipper , a lossless compression tool tailored to compress HapMap data beyond benchmarks defined by generic tools such as gzip , bzip2 and lzma . We demonstrate the usefulness of H ap Z ipper by compressing HapMap 3 populations to 〈5% of their original sizes. H ap Z ipper is freely downloadable from https://bitbucket.org/pchanda/hapzipper/downloads/HapZipper.tar.bz2 .

Description: Spliced alignment plays a central role in the precise identification of eukaryotic gene structures. Even though many spliced alignment programs have been developed, recent rapid progress in DNA sequencing technologies demands further improvements in software tools. Benchmarking algorithms under various conditions is an indispensable task for the development of better software; however, there is a dire lack of appropriate datasets usable for benchmarking spliced alignment programs. In this study, we have constructed two types of datasets: simulated sequence datasets and actual cross-species datasets. The datasets are designed to correspond to various real situations, i.e. divergent eukaryotic species, different types of reference sequences, and the wide divergence between query and target sequences. In addition, we have developed an extended version of our program Spaln , which incorporates two additional features to the scoring scheme of the original version, and examined this extended version, Spaln2, together with the original Spaln and other representative aligners based on our benchmark datasets. Although the effects of the modifications are not individually striking, Spaln2 is consistently most accurate and reasonably fast in most practical cases, especially for plants and fungi and for increasingly divergent pairs of target and query sequences.

Description: Various ‘omics’ technologies, including microarrays and gas chromatography mass spectrometry, can be used to identify hundreds of interesting genes, proteins and metabolites, such as differential genes, proteins and metabolites associated with diseases. Identifying metabolic pathways has become an invaluable aid to understanding the genes and metabolites associated with studying conditions. However, the classical methods used to identify pathways fail to accurately consider joint power of interesting gene/metabolite and the key regions impacted by them within metabolic pathways. In this study, we propose a powerful analytical method referred to as Subpathway-GM for the identification of metabolic subpathways. This provides a more accurate level of pathway analysis by integrating information from genes and metabolites, and their positions and cascade regions within the given pathway. We analyzed two colorectal cancer and one metastatic prostate cancer data sets and demonstrated that Subpathway-GM was able to identify disease-relevant subpathways whose corresponding entire pathways might be ignored using classical entire pathway identification methods. Further analysis indicated that the power of a joint genes/metabolites and subpathway strategy based on their topologies may play a key role in reliably recalling disease-relevant subpathways and finding novel subpathways.

Description: Multiplex analytical systems that allow detection of multiple nucleic acid targets in one assay can provide rapid characterization of a sample while still saving cost and resources. However, few systems have proven to offer a solution for mid-plex (e.g. 10- to 50-plex) analysis that is high throughput and cost effective. Here we describe the combined use of fluorescence color and melting temperature (T m ) as a virtual 2D label that enables homogenous detection of one order of magnitude more targets than current strategies on real-time polymerase chain reaction platform. The target was first hybridized with a pair of ligation oligonucleotides, one of which harbored an artificial sequence that had a unique T m when hybridized with a reporter fluorogenic probe. The ligated products were then amplified by a universal primer pair and denatured by a melting curve analysis procedure. The targets were identified by their respective T m values in the corresponding fluorescence detection channels. The proof-of-principle of this approach was validated by genotyping 15 high-risk human papillomaviruses and 48 human single-nucleotide polymorphisms. The robustness of this method was demonstrated by analyzing a large number of clinical samples in both cases. The combined merits of multiplexity, flexibility and simplicity should make this approach suitable for a variety of applications.

Description: The chromatin structure of eukaryotic telomeres plays an essential role in telomere functions. However, their study might be impaired by the presence of interstitial telomeric sequences (ITSs), which have a widespread distribution in different model systems. We have developed a simple approach to study the chromatin structure of Arabidopsis telomeres independently of ITSs by analyzing ChIP-seq data. This approach could be used to study the chromatin structure of telomeres in some other eukaryotes. The analysis of ChIP-seq experiments revealed that Arabidopsis telomeres have higher density of histone H3 than centromeres, which might reflects their short nucleosomal organization. These experiments also revealed that Arabidopsis telomeres have lower levels of heterochromatic marks than centromeres (H3K9 Me2 and H3K27 Me ), higher levels of some euchromatic marks (H3K4 Me2 and H3K9Ac) and similar or lower levels of other euchromatic marks (H3K4 Me3 , H3K36 Me2 , H3K36 Me3 and H3K18Ac). Interestingly, the ChIP-seq experiments also revealed that Arabidopsis telomeres exhibit high levels of H3K27 Me3 , a repressive mark that associates with many euchromatic genes. The epigenetic profile of Arabidopsis telomeres is closely related to the previously defined chromatin state 2. This chromatin state is found in 23% of Arabidopsis genes, many of which are repressed or lowly expressed. At least, in part, this scenario is similar in rice.

Description: Copy-number variants (CNVs) are a major form of genetic variation and a risk factor for various human diseases, so it is crucial to accurately detect and characterize them. It is conceivable that allele-specific reads from high-throughput sequencing data could be leveraged to both enhance CNV detection and produce allele-specific copy number (ASCN) calls. Although statistical methods have been developed to detect CNVs using whole-genome sequence (WGS) and/or whole-exome sequence (WES) data, information from allele-specific read counts has not yet been adequately exploited. In this paper, we develop an integrated method, called AS-GENSENG, which incorporates allele-specific read counts in CNV detection and estimates ASCN using either WGS or WES data. To evaluate the performance of AS-GENSENG, we conducted extensive simulations, generated empirical data using existing WGS and WES data sets and validated predicted CNVs using an independent methodology. We conclude that AS-GENSENG not only predicts accurate ASCN calls but also improves the accuracy of total copy number calls, owing to its unique ability to exploit information from both total and allele-specific read counts while accounting for various experimental biases in sequence data. Our novel, user-friendly and computationally efficient method and a complete analytic protocol is freely available at https://sourceforge.net/projects/asgenseng/ .

Description: Increased sequencing of microbial genomes has revealed that prevailing prokaryotic species assignments can be inconsistent with whole genome information for a significant number of species. The long-standing need for a systematic and scalable species assignment technique can be met by the genome-wide Average Nucleotide Identity (gANI) metric, which is widely acknowledged as a robust measure of genomic relatedness. In this work, we demonstrate that the combination of gANI and the alignment fraction (AF) between two genomes accurately reflects their genomic relatedness. We introduce an efficient implementation of AF,gANI and discuss its successful application to 86.5M genome pairs between 13,151 prokaryotic genomes assigned to 3032 species. Subsequently, by comparing the genome clusters obtained from complete linkage clustering of these pairs to existing taxonomy, we observed that nearly 18% of all prokaryotic species suffer from anomalies in species definition. Our results can be used to explore central questions such as whether microorganisms form a continuum of genetic diversity or distinct species represented by distinct genetic signatures. We propose that this precise and objective AF,gANI-based species definition: the MiSI (Microbial Species Identifier) method, be used to address previous inconsistencies in species classification and as the primary guide for new taxonomic species assignment, supplemented by the traditional polyphasic approach, as required.

Description: Classically or alternatively activated macrophages (M1 and M2, respectively) play distinct and important roles for microbiocidal activity, regulation of inflammation and tissue homeostasis. Despite this, their transcriptional regulatory dynamics are poorly understood. Using promoter-level expression profiling by non-biased deepCAGE we have studied the transcriptional dynamics of classically and alternatively activated macrophages. Transcription factor (TF) binding motif activity analysis revealed four motifs, NFKB1_REL_RELA, IRF1,2, IRF7 and TBP that are commonly activated but have distinct activity dynamics in M1 and M2 activation. We observe matching changes in the expression profiles of the corresponding TFs and show that only a restricted set of TFs change expression. There is an overall drastic and transient up-regulation in M1 and a weaker and more sustainable up-regulation in M2. Novel TFs, such as Thap6, Maff , (M1) and Hivep1, Nfil3, Prdm1 , (M2) among others, were suggested to be involved in the activation processes. Additionally, 52 (M1) and 67 (M2) novel differentially expressed genes and, for the first time, several differentially expressed long non-coding RNA (lncRNA) transcriptome markers were identified. In conclusion, the finding of novel motifs, TFs and protein-coding and lncRNA genes is an important step forward to fully understand the transcriptional machinery of macrophage activation.

Description: Clostridium difficile is both a hospital and community-acquired pathogen. The current study determined if C. difficile could be cultured from clinical laundry facility surfaces. A total of 240 surface samples were collected from dirty areas ( n = 120), which handle soiled clinical linens, and from clean areas ( n = 120), which process and fold the clean linens, within the University of Washington Consolidated Laundry facility in 2015. Sampling was done four times over the course of 1 year. The dirty area was significantly more contaminated than the clean area (21% vs 2%, P 〈 0.001). Clostridium difficile isolates were genetically characterized using multilocus sequence typing and PCR for the detection of genes encoding toxin A and toxin B. The MLST types 1, 2, 3, 15, 26, 34, 35, 39, 42, 43, 44, 53, 63 and 284 were identified and have previously been found in both clinical and community settings. Toxin positive isolates were identified in both the dirty ( n = 16/25) and clean areas ( n = 2/2). Seasonal variation was observed with 40% of the 27 isolates cultured in April 2015. The study suggests that soiled clinical linens may be a source of C. difficile surface contamination.

Description: Polyploidy is a well-described trait in some prokaryotic organisms; however, it is unusual in marine microbes from oligotrophic environments, which typically display a tendency towards genome streamlining. The biogeochemically significant diazotrophic cyanobacterium Trichodesmium is a potential exception. With a relatively large genome and a comparatively high proportion of non-protein-coding DNA, Trichodesmium appears to allocate relatively more resources to genetic material than closely related organisms and microbes within the same environment. Through simultaneous analysis of gene abundance and direct cell counts, we show for the first time that Trichodesmium spp. can also be highly polyploid, containing as many as 100 genome copies per cell in field-collected samples and 〉600 copies per cell in laboratory cultures. These findings have implications for the widespread use of the abundance of the nifH gene (encoding a subunit of the N 2 -fixing enzyme nitrogenase) as an approach for quantifying the abundance and distribution of marine diazotrophs. Moreover, polyploidy may combine with the unusual genomic characteristics of this genus both in reflecting evolutionary dynamics and influencing phenotypic plasticity and ecological resilience.

Description: The RAG1/RAG2 endonuclease initiates V(D)J recombination at antigen receptor loci but also binds to thousands of places outside of these loci. RAG2 localizes directly to lysine 4 trimethylated histone 3 (H3K4me3) through a plant homeodomain (PHD) finger. The relative contribution of RAG2-dependent and RAG1-intrinsic mechanisms in determining RAG1 binding patterns is not known. Through analysis of deep RAG1 ChIP-seq data, we provide a quantitative description of the forces underlying genome-wide targeting of RAG1. Surprisingly, sequence-specific DNA binding contributes minimally to RAG1 targeting outside of antigen receptor loci. Instead, RAG1 binding is driven by two distinct modes of interaction with chromatin: the first is driven by H3K4me3, promoter-focused and dependent on the RAG2 PHD, and the second is defined by H3K27Ac, enhancer-focused and dependent on ‘non-core’ portions of RAG1. Based on this and additional chromatin and genomic features, we formulated a predictive model of RAG1 targeting to the genome. RAG1 binding sites predicted by our model correlate well with observed patterns of RAG1-mediated breaks in human pro-B acute lymphoblastic leukemia. Overall, this study provides an integrative model for RAG1 genome-wide binding and off-target activity and reveals a novel role for the RAG1 non-core region in RAG1 targeting.

Description: This study aimed to evaluate the survival and gene expression of Vibrio harveyi under starvation conditions. The microcosms V. harveyi were incubated in sterilized seawater for 4 weeks at room temperature. Overall, the cell numeration declined rapidly about 10 3 CFU/ml during starvation, with a tiny rebound at day 21. Scanning electron microscopy revealed that rod-shaped cells became sphere with a rippled cell surface. By polymerase chain reaction (PCR) assay, nine genes, named lux R, tox R, vhh B, fla A, top A, fur , rpo S, mre B and fts Z, were detected in the non-starved cells. In the starved cells, the expression levels of the detected genes declined substantially ranging from 0.005-fold to 0.028-fold compared to the non-starved cells performed by reverse transcription quantitative real-time PCR with 16S rRNA as the internal control. In the recovering cells, the expression levels of the detected genes, except lux R and mre B, were upregulated dramatically compared to the wild, especially top A (23.720-fold), fur (39.400-fold) and tox R (9.837-fold), validating that the expressions of both the metabolism and virulence genes were important for growth and survival of V. harveyi. The results may shed a new light on understanding of stress adaptation in bacteria.

Description: Type 1 fimbriae (T1F) are well characterised cell surface organelles expressed by Escherichia coli and required for adherence to mannosylated host tissue. They satisfy molecular Koch's postulates as a virulence determinant and a host-adapted role has been reinforced by reports that T1F expression is repressed at submammalian temperatures. Analysis of a group of 136 environmental and animal E. coli isolates that express T1F at 37°C showed that 28% are also capable of expression at 20°C, in a phase variable manner. The heterogeneous proportions varied widely, and although growth temperature impacted the total proportion expressing T1F, there was no direct correlation between growth at 37°C and 20°C, indicative of differences in thermoregulation of the genetic switch ( fimS ) that controls phase variation. Specificities of the adhesin (FimH) also varied between the isolates: most bound to α-(1-3) mannan and yeast extracts as expected, but some recognised β-(1-4)-mannans and N -linked glycoproteins from plants, and T1F from two of the isolates mediated binding to plant roots. The results expand our view of a well-described adherence factor to show alternative expression profiles and adhesin specificities, which in turn may confer an advantage for certain isolates in alternative hosts and habitats.

Description: If the in situ growth rate of filamentous bacteria in activated sludge can be quantified, researchers can more accurately assess the effect of operating conditions on the growth of filaments and improve the mathematical modeling of filamentous bulking. We developed a method to quantify the in situ specific growth rate of Sphaerotilus natans (a model filament) in activated sludge using the species-specific 16S rRNA:rDNA ratio. Primers targeting the 16S rRNA of S. natans were designed, and real-time PCR and RT-PCR were used to quantify DNA and RNA levels of S. natans , respectively. A positive linear relationship was found between the rRNA:rDNA ratio (from 440 to 4500) and the specific growth rate of S. natans (from 0.036 to 0.172 h –1 ) using chemostat experiments. The in situ growth rates of S. natans in activated sludge samples from three water reclamation facilities were quantified, illustrating how the approach can be applied in a complex environment such as activated sludge.

Description: Animal-associated microbiotas form complex communities, which play crucial functions for their host, including susceptibility to infections. Despite increasing attention to bats as reservoirs of zoonotic pathogens, their microbiota is poorly documented, especially for samples potentially implicated in pathogen transmission such as urine and saliva. Here, using low-biomass individual samples, we examined the composition and structure of bacterial communities excreted by insectivorous bats, focusing on three body habitats (saliva, urine and faeces). We show that niche specialisation occurs as bacterial community composition was distinct across body habitats with the majority of phylotypes being body habitat specific. Our results suggest that urine harbours more diverse bacterial communities than saliva and faeces and reveal potentially zoonotic bacteria such as Leptospira , Rickettsia , Bartonella and Coxiella in all body habitats. Our study emphasised that, in addition to the traditional use of gut-associated samples such as faeces, both urine and saliva are also of interest because of their diverse microbiota and the potential transmission of pathogenic bacteria. Our results represent a critical baseline for future studies investigating the interactions between microbiota and infection dynamics in bats.

Description: R-type bacteriocins are contractile phage tail-like structures that are bactericidal towards related bacterial species. The C-terminal region of the phage tail fiber protein determines target-binding specificity. The mutualistic bacteria Xenorhabdus nematophila and X. bovienii produce R-type bacteriocins (xenorhabdicins) that are selectively active against different Xenorhabdus species. We analyzed the P2-type remnant prophage clusters in draft sequences of nine strains of X. bovienii . The C-terminal tail fiber region in each of the respective strains was unique and consisted of mosaics of modular units. The region between the main tail fiber gene ( xbpH1 ) and the sheath gene ( xbpS1 ) contained a variable number of modules encoding tail fiber fragments. DNA inversion and module exchange between strains was involved in generating tail fiber diversity. Xenorhabdicin-enriched fractions from three different X. bovienii strains isolated from the same nematode species displayed distinct activities against each other. In one set of strains, the strain that produced highly active xenorhabdicin was able to eliminate a sensitive strain. In contrast, xenorhabdicin activity was not a determining factor in the competitive fitness of a second set of strains. These findings suggest that related strains of X. bovienii use xenorhabdicin and additional antagonistic molecules to compete against each other.

Description: Anabaena PCC7120 has two annotated toxin–antitoxin systems: MazEF and HicAB. Overexpression of either of the toxins severely inhibited the growth of Escherichia coli BL21(p lysS )(DE3). Of the two Anabaena toxins, MazF exhibited higher toxicity than HicA as evidenced by (i) 100-fold lower viability upon overexpression of MazF compared to HicA; (ii) complete loss of cell viability within 1 h of induction of MazF expression, as against 〉10 3 colony forming units mL –1 in case of HicA; (iii) inability to maintain the MazF overexpressing plasmid in E. coli cells; and (iv) neutralisation of the toxin was effective at the molar ratio of 1:1.9 for MazF:MazE and 13:1 for HicA:HicB, indicating higher antitoxin requirement for neutralisation of MazF. The growth inhibitory effect of MazF was found to be higher in lag phase cultures compared to mid-logarithmic phase cultures of E. coli , while the reverse was true for HicA. The results suggest possible distinct roles for MazEF and HicAB systems of Anabaena .

Description: Members of subdivision 1 of the phylum Acidobacteria were grown at different pH values in a new medium formulation named PSYL 5, which includes sucrose as a carbon source and other compounds (such as KH 2 PO 4 and MgSO 4 .7H 2 O). Growth rate was nearly constant at pH 5.0 and declined at pH 3–4 and 6–7. However, it was found that effects involving good carbon/nitrogen ratios and pH on the growth of the members of Acidobacteria subdivision 1 were significant, and the strongest effect of these conditions was at pH 5.0. In addition, incubation time of 48, 72, 96 and 120 h was shorter than that described previously for members of Acidobacteria subdivision 1 on solid laboratory media.