Non-canonical transcriptional start sites in E. coli O157:H7 EDL933 are regulated and appear in surprisingly high numbers

Analysis of genome wide transcription start sites (TSSs) revealed an unexpected complexity since not only canonical TSS of annotated genes are recognized by RNA polymerase. Non-canonical TSS were detected antisense to, or within, annotated genes as well new intergenic (orphan) TSS, not associated with known genes. Previously, it was hypothesized that many such signals represent noise or pervasive transcription, not associated with a biological function. Here, a modified Cappable-seq protocol allows determining the primary transcriptome of the enterohemorrhagic E. coli O157:H7 EDL933 (EHEC). We used four different growth media, both in exponential and stationary growth phase, replicated each thrice. This yielded 19,975 EHEC canonical and non-canonical TSS, which reproducibly occurring in three biological replicates. This questions the hypothesis of experimental noise or pervasive transcription. Accordingly, conserved promoter motifs were found upstream indicating proper TSSs. More than 50% of 5,567 canonical and between 32% and 47% of 10,355 non-canonical TSS were differentially expressed in different media and growth phases, providing evidence for a potential biological function also of non-canonical TSS. Thus, reproducible and environmentally regulated expression suggests that a substantial number of the non-canonical TSSs may be of unknown function rather than being the result of noise or pervasive transcription.

Unified Workflow for the Rapid and In-Depth Characterization of Bacterial Proteomes

Bacteria are the most abundant and diverse organisms among the kingdoms of life. Due to this excessive variance, finding a unified, comprehensive, and safe workflow for quantitative bacterial proteomics is challenging. In this study, we have systematically evaluated and optimized sample preparation, mass spectrometric data acquisition, and data analysis strategies in bacterial proteomics. We investigated workflow performances on six representative species with highly different physiologic properties to mimic bacterial diversity. The best sample preparation strategy was a cell lysis protocol in 100% trifluoroacetic acid followed by an in-solution digest. Peptides were separated on a 30-min linear microflow liquid chromatography gradient and analyzed in data-independent acquisition mode. Data analysis was performed with DIA-NN using a predicted spectral library. Performance was evaluated according to the number of identified proteins, quantitative precision, throughput, costs, and biological safety. With this rapid workflow, over 40% of all encoded genes were detected per bacterial species. We demonstrated the general applicability of our workflow on a set of 23 taxonomically and physiologically diverse bacterial species. We could confidently identify over 45,000 proteins in the combined dataset, of which 30,000 have not been experimentally validated before. Our work thereby provides a valuable resource for the microbial scientific community. Finally, we grew Escherichia coli and Bacillus cereus in replicates under 12 different cultivation conditions to demonstrate the high-throughput suitability of the workflow. The proteomic workflow we present in this manuscript does not require any specialized equipment or commercial software and can be easily applied by other laboratories to support and accelerate the proteomic exploration of the bacterial kingdom.

Three Novel Antisense Overlapping Genes in E. coli O157:H7 EDL933

The abundance of long overlapping genes in prokaryotic genomes is likely to be significantly underestimated. To date, only a few examples of such genes are fully established. Using RNA sequencing and ribosome profiling, we found expression of novel overlapping open reading frames in Escherichia coli O157:H7 EDL933 (EHEC). Indeed, the overlapping candidate genes are equipped with typical structural elements required for transcription and translation, i.e., promoters, transcription start sites, as well as terminators, all of which were experimentally verified. Translationally arrested mutants, unable to produce the overlapping encoded protein, were found to have a growth disadvantage when grown competitively against the wild type. Thus, the phenotypes found imply biological functionality of the genes at the level of proteins produced. The addition of 3 more examples of prokaryotic overlapping genes to the currently limited, yet constantly growing pool of such genes emphasizes the underestimated coding capacity of bacterial genomes. IMPORTANCE The abundance of long overlapping genes in prokaryotic genomes is likely to be significantly underestimated, since such genes are not allowed in genome annotations. However, ribosome profiling catches mRNA in the moment of being template for protein production. Using this technique and subsequent experiments, we verified 3 novel overlapping genes encoded in antisense of known genes. This adds more examples of prokaryotic overlapping genes to the currently limited, yet constantly growing pool of such genes.

Spotlight on alternative frame coding: Two long overlapping genes in Pseudomonas aeruginosa are translated and under purifying selection

The existence of overlapping genes (OLGs) with significant coding overlaps revolutionizes our understanding of genomic complexity. We report two exceptionally long (957 nt and 1536 nt), evolutionarily novel, translated antisense open reading frames (ORFs) embedded within annotated genes in the pathogenic Gram-negative bacterium Pseudomonas aeruginosa. Both OLG pairs show sequence features consistent with being genes and transcriptional signals in RNA sequencing. Translation of both OLGs was confirmed by ribosome profiling and mass spectrometry. Quantitative proteomics of samples taken during different phases of growth revealed regulation of protein abundances, implying biological functionality. Both OLGs are taxonomically restricted, and likely arose by overprinting within the genus. Evidence for purifying selection further supports functionality. The OLGs reported here, designated olg1 and olg2, are the longest yet proposed in prokaryotes and are among the best attested in terms of translation and evolutionary constraint. These results highlight a potentially large unexplored dimension of prokaryotic genomes.

•

Two novel, very long, overlapping genes were found in Pseudomonas aeruginosa

•

Both overlapping genes, olg1 and olg2, are transcribed, translated, and regulated

•

Mass spectrometry verifies translation of the overlapping and their mother genes

•

Both overlapping genes are taxonomically restricted, but under purifying selection

Biological factors in the synthetic construction of overlapping genes

BACKGROUND

Overlapping genes (OLGs) with long protein-coding overlapping sequences are disallowed by standard genome annotation programs, outside of viruses. Recently however they have been discovered in Archaea, diverse Bacteria, and Mammals. The biological factors underlying life's ability to create overlapping genes require more study, and may have important applications in understanding evolution and in biotechnology. A previous study claimed that protein domains from viruses were much better suited to forming overlaps than those from other cellular organisms - in this study we assessed this claim, in order to discover what might underlie taxonomic differences in the creation of gene overlaps.

RESULTS

After overlapping arbitrary Pfam domain pairs and evaluating them with Hidden Markov Models we find OLG construction to be much less constrained than expected. For instance, close to 10% of the constructed sequences cannot be distinguished from typical sequences in their protein family. Most are also indistinguishable from natural protein sequences regarding identity and secondary structure. Surprisingly, contrary to a previous study, virus domains were much less suitable for designing OLGs than bacterial or eukaryotic domains were. In general, the amount of amino acid change required to force a domain to overlap is approximately equal to the variation observed within a typical domain family. The resulting high similarity between natural sequences and those altered so as to overlap is mostly due to the combination of high redundancy in the genetic code and the evolutionary exchangeability of many amino acids.

CONCLUSIONS

Synthetic overlapping genes which closely resemble natural gene sequences, as measured by HMM profiles, are remarkably easy to construct, and most arbitrary domain pairs can be altered so as to overlap while retaining high similarity to the original sequences. Future work however will need to assess important factors not considered such as intragenic interactions which affect protein folding. While the analysis here is not sufficient to guarantee functional folding proteins, further analysis of constructed OLGs will improve our understanding of the origin of these remarkable genetic elements across life and opens up exciting possibilities for synthetic biology.

Facklamia lactis sp. nov., isolated from raw milk

Two strains of a Gram-staining-positive species were isolated from German bulk tank milk. On the basis of their 16S rRNA sequences they were affiliated to the genus Facklamia but could not be assigned to any species with a validly published name. Facklamia miroungae ATCC BAA-466^T (97.3 % 16S rRNA sequence similarity), Facklamia languida CCUG 37842^T (96.9 %), and Facklamia hominis CCUG 36813^T (96.6 %) are the closest relatives. In the 16S rRNA phylogeny and in the core-genome phylogeny strains WS 5301^T and WS 5302 form a well-supported, separate lineage. Pairwise average nucleotide identity calculated using MUMmer (ANIm) between WS 5301^T and type strains of other Facklamia species is well below the species cut-off (95 %) and ranges from 83.4 to 87.7 %. The DNA G+C content of the type strain is 36.4 mol% and the assembly size of the genome is 2.2 Mb. Cells of WS 5301^T are non-motile, non-endospore-forming, oxidase-negative, catalase-negative and facultatively anaerobic cocci. The fastidious species grows at 10-40 °C and with up to 7.0 % (w/v) NaCl in BHI supplemented with 5 g l^-1 yeast extract. Major polar lipids are phosphatidylglycerol, diphosphatidylglycerol and two glycolipids. Predominant fatty acids are C_16 : 1ω9c and C_18 : 1ω9c. On the basis of their genomic, physiological and chemotaxonomic characteristics the strains examined in this study represent the same, hitherto unknown species. We propose the name Facklamia lactis sp. nov. for which WS 5301^T (=DSM 111018^T=LMG 31861^T) is the type strain and WS 5302 (=DSM 111019=LMG 31862) is an additional strain of this novel species.

A Strong Synergy Between the Thiopeptide Bacteriocin Micrococcin P1 and Rifampicin Against MRSA in a Murine Skin Infection Model

Antibiotic-resistant bacterial pathogens have become a serious threat worldwide. One of these pathogens is methicillin-resistant Staphylococcus aureus (MRSA), a major cause of skin and soft tissue infections. In this study we identified a strain of Staphylococcus equorum producing a substance with high antimicrobial activity against many Gram-positive bacteria, including MRSA. By mass spectrometry and whole genome sequencing the antimicrobial substance was identified as the thiopeptide bacteriocin micrococcin P1 (MP1). Based on its properties we developed a one-step purification protocol resulting in high yield (15 mg/L) and high purity (98%) of MP1. For shorter incubation times (5-7 h) MP1 was very potent against MRSA but the inhibitory effect was overshadowed by resistance development during longer incubation time (24h or more). To overcome this problem a synergy study was performed with a number of commercially available antibiotics. Among the antibiotics tested, the combination of MP1 and rifampicin gave the best synergistic effect, with MIC values 25 and 60 times lower than for the individual drugs, respectively. To assess the therapeutic potential of the MP1-rifampicin combination, we used a murine skin infection model based on the use of the multidrug-resistant luciferase-tagged MRSA strain Xen31. As expected, neither of the single antimicrobials (MP1 or rifampicin) could eradicate Xen31 from the wounds. By contrary, the MP1-rifampicin combination was efficient not only to eradicate but also to prevent the recurrence of Xen31 infection. Furthermore, compared to fucidin cream, which is commonly used in skin infection treatments, MP1-rifampicin combination was superior in terms of preventing resistance development. Our results show that combining MP1, and probably other thiopeptides, with antibiotics can be a promising strategy to treat SSTIs caused by MRSA and likely many other Gram-positive bacteria.

Amplicon-sequencing of raw milk microbiota: impact of DNA extraction and library-PCR

Abstract

The highly complex raw milk matrix challenges the sample preparation for amplicon-sequencing due to low bacterial counts and high amounts of eukaryotic DNA originating from the cow. In this study, we optimized the extraction of bacterial DNA from raw milk for microbiome analysis and evaluated the impact of cycle numbers in the library-PCR. The selective lysis of eukaryotic cells by proteinase K and digestion of released DNA before bacterial lysis resulted in a high reduction of mostly eukaryotic DNA and increased the proportion of bacterial DNA. Comparative microbiome analysis showed that a combined enzymatic and mechanical lysis procedure using the DNeasy^® PowerFood^® Microbial Kit with a modified protocol was best suitable to achieve high DNA quantities after library-PCR and broad coverage of detected bacterial biodiversity. Increasing cycle numbers during library-PCR systematically altered results for species and beta-diversity with a tendency to overrepresentation or underrepresentation of particular taxa. To limit PCR bias, high cycle numbers should thus be avoided. An optimized DNA extraction yielding sufficient bacterial DNA and enabling higher PCR efficiency is fundamental for successful library preparation. We suggest that a protocol using ethylenediaminetetraacetic acid (EDTA) to resolve casein micelles, selective lysis of somatic cells, extraction of bacterial DNA with a combination of mechanical and enzymatic lysis, and restriction of PCR cycles for analysis of raw milk microbiomes is optimal even for samples with low bacterial numbers.

Key points

• Sample preparation for high-throughput 16S rRNA gene sequencing of raw milk microbiota.

• Reduction of eukaryotic DNA by enzymatic digestion.

• Shift of detected microbiome caused by high cycle numbers in library-PCR.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-021-11353-4.

Simultaneous quantification of the most common and proteolytic Pseudomonas species in raw milk by multiplex qPCR

Abstract

The heat-stable peptidase AprX, secreted by psychrotolerant Pseudomonas species in raw milk, is a major cause of destabilization and premature spoilage of ultra-high temperature (UHT) milk and milk products. To enable rapid detection and quantification of seven frequent and proteolytic Pseudomonas species (P. proteolytica, P. gessardii, P. lactis, P. fluorescens, P. protegens, P. lundensis, and P. fragi) in raw milk, we developed two triplex qPCR assays taking into account species-dependent differences in AprX activity. Besides five species-specific hydrolysis probes, targeting the aprX gene, a universal rpoB probe was included in the assay to determine the total Pseudomonas counts. For all six probes, linear regression lines between C_q value and target DNA concentration were obtained in singleplex as well as in multiplex approaches, yielding R² values of > 0.975 and amplification efficiencies of 85–97%. Moreover, high specificity was determined using genomic DNA of 75 Pseudomonas strains, assigned to 57 species, and 40 other bacterial species as templates in the qPCR. Quantification of the target species and total Pseudomonas counts resulted in linear detection ranges of approx. 10³–10⁷ cfu/ml, which correspond well to common Pseudomonas counts in raw milk. Application of the assay using 60 raw milk samples from different dairies showed good agreement of total Pseudomonas counts calculated by qPCR with cell counts derived from cultivation. Furthermore, a remarkably high variability regarding the species composition was observed for each milk sample, whereby P. lundensis and P. proteolytica/P. gessardii were the predominant species detected.

Key points

• Multiplex qPCR for quantification of seven proteolytic Pseudomonas species and total Pseudomonas counts in raw milk

• High specificity and sensitivity via hydrolysis probes against aprX and rpoB

• Rapid method to determine Pseudomonas contamination in raw milk and predict spoilage potential

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-021-11109-0.

Pseudomonas cremoris sp. nov., a novel proteolytic species isolated from cream

During a study investigating the microbiota of raw milk and its semi-finished products, strains WS 5106^T and WS 5096 were isolated from cream and skimmed milk concentrate. They could be assigned to the genus Pseudomonas by their 16S rRNA sequences, but not to any validly named species. In this work, a polyphasic approach was used to characterize the novel strains and to investigate their taxonomic status. Examinations based on the topology of core genome phylogenomy as well as average nucleotide identity (ANIm) comparisons suggested a novel Pseudomonas species within the Pseudomonas fluorescens subgroup. With pairwise ANIm values of 90.1 and 89.8 %, WS 5106^T was most closely related to Pseudomonas nabeulensis CECT 9765^T and Pseudomonas kairouanensis CECT 9766^T. The G+C content of strain WS 5106^T was 60.1 mol%. Morphologic analyses revealed Gram-stain-negative, aerobic, catalase and oxidase positive, rod-shaped and motile cells. Proteolysis on skimmed milk agar as well as lipolysis on tributyrin agar occurred at both 28 and 6 °C. Tolerated growth conditions were temperatures between 4 and 34 °C, pH values between 6.0 and 8.0, and salt concentrations of up to 5 %. Fatty acid profiles showed a pattern typical for Pseudomonas, with C_16 : 0 as the dominant component. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol and the dominating quinone was Q-9. Based on these results, it is proposed to classify the strains as a novel species, Pseudomonas cremoris sp. nov., with WS 5106^T (=DSM 111143^T=LMG 31863^T) as type strain and WS 5096 (=DSM 111129=LMG 31864) as an additional strain.

Are Antisense Proteins in Prokaryotes Functional?

Many prokaryotic RNAs are transcribed from loci outside of annotated protein coding genes. Across bacterial species hundreds of short open reading frames antisense to annotated genes show evidence of both transcription and translation, for instance in ribosome profiling data. Determining the functional fraction of these protein products awaits further research, including insights from studies of molecular interactions and detailed evolutionary analysis. There are multiple lines of evidence, however, that many of these newly discovered proteins are of use to the organism. Condition-specific phenotypes have been characterized for a few. These proteins should be added to genome annotations, and the methods for predicting them standardized. Evolutionary analysis of these typically young sequences also may provide important insights into gene evolution. This research should be prioritized for its exciting potential to uncover large numbers of novel proteins with extremely diverse potential practical uses, including applications in synthetic biology and responding to pathogens.

Fundicoccus ignavus gen. nov., sp. nov., a novel genus of the family Aerococcaceae isolated from bulk tank milk

Three strains of a Gram-stain-positive, catalase-negative, facultative anaerobic, and coccoid species were isolated from German bulk tank milk. Phylogenetic analyses based on the 16S rRNA gene sequences indicated that the three strains (WS4937^T, WS4759 and WS5303) constitute an independent phylogenetic lineage within the family Aerococcaceae with Facklamia hominis CCUG 36813^T (93.7-94.1 %) and Eremococcus coleocola M1831/95/2^T (93.5 %) as most closely related type species. The unclassified strains demonstrated variable growth with 6.5 % (w/v) NaCl and tolerated pH 6.5-9.5. Growth was observed from 12 to 39 °C. Their cell-wall peptidoglycan belongs to the A1α type (l-Lys-direct) consisting of alanine, glutamic acid and lysine. The predominant fatty acids were C_16 : 1 ω9c, C_16 : 0 and C_18 : 1 ω9c and in the polar lipids profile three glycolipids, a phospholipid, phosphatidylglycerol, phosphoglycolipid and diphosphatidylglycerol were found. The G+C content of strain WS4937^T was 37.4 mol% with a genome size of ~3.0 Mb. Based on phylogenetic, phylogenomic and biochemical characterizations, the isolates can be demarcated from all other genera of the family Aerococcaceae and, therefore, the novel genus Fundicoccus gen. nov. is proposed. The type species of the novel genus is Fundicoccus ignavus gen. nov., sp. nov. WS4937^T (=DSM 109652^T=LMG 31441^T).

Complementary Use of Cultivation and High-Throughput Amplicon Sequencing Reveals High Biodiversity Within Raw Milk Microbiota

Raw milk microbiota are complex communities with a significant impact on the hygienic, sensory and technological quality of milk products. However, there is a lack of knowledge on factors determining their composition. In the present study, four bulk tank milk samples of two farms at two different time points were analyzed in detail for their microbiota using cultivation and 16S rRNA amplicon sequencing. Diversity in samples from the first time point was assessed via cultivation of 500 aerobic mesophilic bacterial isolates in each sample. A high biodiversity of 70 and 110 species per sample was determined, of which 25–28% corresponded to yet unknown taxa. The isolates were dominated by Gram-positive members of the genera Staphylococcus, Corynebacterium, Streptococcus, or Janibacter, whilst Chryseobacterium and Acinetobacter were most abundant among the Gram-negative taxa. At the second time point, samples of the same farms were analyzed via both cultivation (1,500 individual colonies each) and high-throughput 16S rRNA gene amplicon sequencing. The latter revealed a threefold higher biodiversity at the genus level, as anaerobic or fastidious species were also detected. However, cultivation identified genera not captured by sequencing, indicating that both approaches are complementary. Using amplicon sequencing, the relative abundance of a few genera was distorted, which seems to be an artifact of sample preparation. Therefore, attention needs to be paid to the library preparation procedure with special emphasis on cell lysis and PCR.

Genetic Organization of the aprX-lipA2 Operon Affects the Proteolytic Potential of Pseudomonas Species in Milk

Psychrotolerant Pseudomonas species are a main cause of proteolytic spoilage of ultra-high temperature (UHT) milk products due to the secretion of the heat-resistant metallopeptidase AprX, which is encoded by the first gene of the aprX-lipA2 operon. While the proteolytic property has been characterized for many different Pseudomonas isolates, the underlying aprX-lipA2 gene organization was only described for a few strains so far. In this study, the phylogenomic analysis of 185 Pseudomonas type strains revealed that the presence of aprX is strongly associated to a monophylum composed of 81 species, of which 83% carried the aprX locus. Furthermore, almost all type strains of known milk-relevant species were shown to be members of the three monophyletic groups P. fluorescens, P. gessardii, and P. fragi. In total, 22 different types of aprX-lipA2 genetic organizations were identified in the genus, whereby 31% of the species tested carried the type 1 operon structure consisting of eight genes (aprXIDEF prtAB lipA2). Other genetic structures differed from type 1 mainly in the presence and location of genes coding for two lipases (lipA1 and lipA2) and putative autotransporters (prtA and prtB). The peptidase activity of 129 strains, as determined on skim milk agar and in UHT-milk, correlated largely with different aprX-lipA2 gene compositions. Particularly, isolates harboring the type 1 operon were highly proteolytic, while strains with other operon types, especially ones lacking prtA and prtB, exhibited significantly lower peptidase activities. In conclusion, the phylogenomic position and the aprX-lipA2 gene organization specify the proteolytic potential of Pseudomonas isolates. In addition, however, an interplay of several environmental factors and intrinsic traits influences production and activity of AprX, leading to strain-specific proteolytic phenotypes.

A Novel pH-Regulated, Unusual 603 bp Overlapping Protein Coding Gene pop Is Encoded Antisense to ompA in Escherichia coli O157:H7 (EHEC)

Antisense transcription is well known in bacteria. However, translation of antisense RNAs is typically not considered, as the implied overlapping coding at a DNA locus is assumed to be highly improbable. Therefore, such overlapping genes are systematically excluded in prokaryotic genome annotation. Here we report an exceptional 603 bp long open reading frame completely embedded in antisense to the gene of the outer membrane protein ompA. An active σ70 promoter, transcription start site (TSS), Shine-Dalgarno motif and rho-independent terminator were experimentally validated, providing evidence that this open reading frame has all the structural features of a functional gene. Furthermore, ribosomal profiling revealed translation of the mRNA, the protein was detected in Western blots and a pH-dependent phenotype conferred by the protein was shown in competitive overexpression growth experiments of a translationally arrested mutant versus wild type. We designate this novel gene pop (pH-regulated overlapping protein-coding gene), thus adding another example to the growing list of overlapping, protein coding genes in bacteria.

Proposal of Lysobacter pythonis sp. nov. isolated from royal pythons (Python regius)

The bacterial strains 4284/11T and 812/17 isolated from the respiratory tract of two royal pythons in 2011 and 2017, respectively were subjected to taxonomic characterization. The 16S rRNA gene sequences of the two strains were identical and showed highest sequence similarities to Lysobacter tolerans UM1^T (97.2%) and Luteimonas aestuarii DSM 19680^T (96.7 %). The two strains were identical in the sequences of the 16S-23S rRNA internal transcribed spacer (ITS) and partial groEL gene sequences and almost identical in genomic fingerprints. In the ITS sequence Ly. tolerans DSM 28473^T and in the groEL nucleotide sequence Luteimonas mephitis DSM 12574^T showed the highest similarity. In silico DDH analyses using genome sequence based ANIb and gANI similarity coefficients demonstrated that strain 4284/11^T represents a novel species and revealed Ly. tolerans UM1^T as the next relative (ANIb = 76.2 %, gANI = 78.0 %). Based on the topology of a core gene phylogeny strain 4284/11^T could be assigned to the genus Lysobacter. Chemotaxonomic characteristics including polyamine pattern, quinone system, polar lipid profile and fatty acid profile were in accordance with the characteristics of the genera Lysobacter and Luteimonas. Strains 4284/11^T and 812/17 could be differentiated from the type strains of the most closely related species by several physiological tests. In conclusion we are here proposing the novel species Lysobacter pythonis sp. nov. The type strain is 4284/11^T (= CCM 8829^T = CCUG 72164^T = LMG 30630^T) and strain 812/17 (CCM 8830) is a second strain of this species.

Dynamic Proteome Alteration and Functional Modulation of Human Saliva Induced by Dietary Chemosensory Stimuli

Saliva flow measurements and SDS-PAGE separation of human whole saliva freshly collected after oral stimulation with citric acid (sour), aspartame (sweet), iso-α-acids (bitter), mono sodium l-glutamate (umami), NaCl (salty), 6-gingerol (pungent), hydroxy-α-sanshool (tingling), and hydroxy-β-sanshool (numbing), followed by tryptic digestion, nano-HPLC-MS/MS, and label-free protein quantitation demonstrated a stimulus- and time-dependent influence of the dietary chemosensates on salivation and the salivary proteome composition. Gene ontology enrichment analysis showed evidence for stimulus-induced alterations of the saliva proteome to boot an efficient molecular defense network of the oral cavity, e.g., 6-gingerol increased salivary lactoperoxidase activity, catalyzing the oxidation of thiocyanate to produce the antimicrobial and antifungal hypothiocyanate, from 0.37 ± 0.02 to 0.91 ± 0.05 mU/mL 45 s after stimulation. In comparison, oral citric acid stimulation induced an increase of myeloperoxidase activity, catalyzing the chloride oxidation to generate antimicrobial hypochloride in saliva, from 0.24 ± 0.04 to 0.70 ± 0.1 mU/mL as well as an increase of salivary levels of lysozyme, exhibiting antimicrobial activity on Gram-positive bacteria, from 6.0-10 to 100-150 μg/mL. Finally, microbial growth experiments clearly demonstrated for the first time that the increase of the salivary lysozyme abundance upon oral citric acid stimulation translates into an enhanced biological function, that is an almost complete growth inhibition of the two lysozyme-sensitive Gram-positive bacteria tested.

Oxidative Stress and Cardiovascular Risk in Type 1 Diabetes Mellitus: Insights From the DCCT/EDIC Study

Background

Hyperglycemia leading to increased oxidative stress is implicated in the increased risk for the development of macrovascular and microvascular complications in patients with type 1 diabetes mellitus.

Methods and Results

A random subcohort of 349 participants was selected from the DCCT / EDIC (Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications) cohort. This included 320 controls and 29 cardiovascular disease cases that were augmented with 98 additional known cases to yield a case cohort of 447 participants (320 controls, 127 cases). Biosamples from DCCT baseline, year 1, and closeout of DCCT , and 1 to 2 years post‐ DCCT ( EDIC years 1 and 2) were measured for markers of oxidative stress, including plasma myeloperoxidase, paraoxonase activity, urinary F 2α isoprostanes, and its metabolite, 2,3 dinor‐8 iso prostaglandin F 2α . Following adjustment for glycated hemoblobin and weighting the observations inversely proportional to the sampling selection probabilities, higher paraoxonase activity, reflective of antioxidant activity, and 2,3 dinor‐8 iso prostaglandin F 2α , an oxidative marker, were significantly associated with lower risk of cardiovascular disease (−4.5% risk for 10% higher paraoxonase, P <0.003; −5.3% risk for 10% higher 2,3 dinor‐8 iso prostaglandin F 2α , P =0.0092). In contrast, the oxidative markers myeloperoxidase and F 2α isoprostanes were not significantly associated with cardiovascular disease after adjustment for glycated hemoblobin. There were no significant differences between DCCT intensive and conventional treatment groups in the change in all biomarkers across time segments.

Conclusions

Heightened antioxidant activity (rather than diminished oxidative stress markers) is associated with lower cardiovascular disease risk in type 1 diabetes mellitus, but these biomarkers did not change over time with intensification of glycemic control.

Clinical Trial Registration

URL : https://www.clinicaltrials.gov . Unique identifiers: NCT 00360815 and NCT 00360893.

The Novel Anaerobiosis-Responsive Overlapping Gene ano Is Overlapping Antisense to the Annotated Gene ECs2385 of Escherichia coli O157:H7 Sakai

Current notion presumes that only one protein is encoded at a given bacterial genetic locus. However, transcription and translation of an overlapping open reading frame (ORF) of 186 bp length were discovered by RNAseq and RIBOseq experiments. This ORF is almost completely embedded in the annotated L,D-transpeptidase gene ECs2385 of Escherichia coli O157:H7 Sakai in the antisense reading frame -3. The ORF is transcribed as part of a bicistronic mRNA, which includes the annotated upstream gene ECs2384, encoding a murein lipoprotein. The transcriptional start site of the operon resides 38 bp upstream of the ECs2384 start codon and is driven by a predicted σ⁷⁰ promoter, which is constitutively active under different growth conditions. The bicistronic operon contains a ρ-independent terminator just upstream of the novel gene, significantly decreasing its transcription. The novel gene can be stably expressed as an EGFP-fusion protein and a translationally arrested mutant of ano, unable to produce the protein, shows a growth advantage in competitive growth experiments compared to the wild type under anaerobiosis. Therefore, the novel antisense overlapping gene is named ano (anaerobiosis responsive overlapping gene). A phylostratigraphic analysis indicates that ano originated very recently de novo by overprinting after the Escherichia/Shigella clade separated from other enterobacteria. Therefore, ano is one of the very rare cases of overlapping genes known in the genus Escherichia.

A novel short L-arginine responsive protein-coding gene (laoB) antiparallel overlapping to a CadC-like transcriptional regulator in Escherichia coli O157:H7 Sakai originated by overprinting

Background

Due to the DNA triplet code, it is possible that the sequences of two or more protein-coding genes overlap to a large degree. However, such non-trivial overlaps are usually excluded by genome annotation pipelines and, thus, only a few overlapping gene pairs have been described in bacteria. In contrast, transcriptome and translatome sequencing reveals many signals originated from the antisense strand of annotated genes, of which we analyzed an example gene pair in more detail.

Results

A small open reading frame of Escherichia coli O157:H7 strain Sakai (EHEC), designated laoB (L-arginine responsive overlapping gene), is embedded in reading frame −2 in the antisense strand of ECs5115, encoding a CadC-like transcriptional regulator. This overlapping gene shows evidence of transcription and translation in Luria-Bertani (LB) and brain-heart infusion (BHI) medium based on RNA sequencing (RNAseq) and ribosomal-footprint sequencing (RIBOseq). The transcriptional start site is 289 base pairs (bp) upstream of the start codon and transcription termination is 155 bp downstream of the stop codon. Overexpression of LaoB fused to an enhanced green fluorescent protein (EGFP) reporter was possible. The sequence upstream of the transcriptional start site displayed strong promoter activity under different conditions, whereas promoter activity was significantly decreased in the presence of L-arginine. A strand-specific translationally arrested mutant of laoB provided a significant growth advantage in competitive growth experiments in the presence of L-arginine compared to the wild type, which returned to wild type level after complementation of laoB in trans. A phylostratigraphic analysis indicated that the novel gene is restricted to the Escherichia/Shigella clade and might have originated recently by overprinting leading to the expression of part of the antisense strand of ECs5115.

Conclusions

Here, we present evidence of a novel small protein-coding gene laoB encoded in the antisense frame −2 of the annotated gene ECs5115. Clearly, laoB is evolutionarily young and it originated in the Escherichia/Shigella clade by overprinting, a process which may cause the de novo evolution of bacterial genes like laoB.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1134-0) contains supplementary material, which is available to authorized users.

Discovery of numerous novel small genes in the intergenic regions of the Escherichia coli O157:H7 Sakai genome

In the past, short protein-coding genes were often disregarded by genome annotation pipelines. Transcriptome sequencing (RNAseq) signals outside of annotated genes have usually been interpreted to indicate either ncRNA or pervasive transcription. Therefore, in addition to the transcriptome, the translatome (RIBOseq) of the enteric pathogen Escherichia coli O157:H7 strain Sakai was determined at two optimal growth conditions and a severe stress condition combining low temperature and high osmotic pressure. All intergenic open reading frames potentially encoding a protein of ≥ 30 amino acids were investigated with regard to coverage by transcription and translation signals and their translatability expressed by the ribosomal coverage value. This led to discovery of 465 unique, putative novel genes not yet annotated in this E. coli strain, which are evenly distributed over both DNA strands of the genome. For 255 of the novel genes, annotated homologs in other bacteria were found, and a machine-learning algorithm, trained on small protein-coding E. coli genes, predicted that 89% of these translated open reading frames represent bona fide genes. The remaining 210 putative novel genes without annotated homologs were compared to the 255 novel genes with homologs and to 250 short annotated genes of this E. coli strain. All three groups turned out to be similar with respect to their translatability distribution, fractions of differentially regulated genes, secondary structure composition, and the distribution of evolutionary constraint, suggesting that both novel groups represent legitimate genes. However, the machine-learning algorithm only recognized a small fraction of the 210 genes without annotated homologs. It is possible that these genes represent a novel group of genes, which have unusual features dissimilar to the genes of the machine-learning algorithm training set.

Complete Circular Genome Sequence and Temperature Independent Adaptation to Anaerobiosis of Listeria weihenstephanensis DSM 24698

The aim of this study was to analyze the adaptation of the environmental Listeria weihenstephanensis DSM 24698 to anaerobiosis. The complete circular genome sequence of this species is reported and the adaptation of L. weihenstephanensis DSM 24698 to oxygen availability was investigated by global transcriptional analyses via RNAseq at 18 and 34°C. A list of operons was created based on the transcriptional data. Forty-two genes were upregulated anaerobically and 62 genes were downregulated anaerobically. The oxygen dependent gene expression of selected genes was further validated via qPCR. Many of the differentially regulated genes encode metabolic enzymes indicating broad metabolic adaptations with respect to oxygen availability. Genes showing the strongest oxygen-dependent adaption encoded nitrate (narGHJI) and nitrite (nirBD) reductases. Together with the observation that nitrate supported anaerobic growth, these data indicate that L. weihenstephanensis DSM 24698 performs anaerobic nitrate respiration. The wide overlap between the oxygen-dependent transcriptional regulation at 18 and 34°C suggest that temperature does not play a key role in the oxygen-dependent transcriptional regulation of L. weihenstephanensis DSM 24698.

Pseudomonas lactis sp. nov. and Pseudomonas paralactis sp. nov., isolated from bovine raw milk

Five strains, designated WS 4672T, WS 4998, WS 4992T, WS 4997 and WS 5000, isolated from bovine raw milk formed two individual groups in a phylogenetic analysis. The most similar species on the basis of 16S rRNA gene sequences were Pseudomonas azotoformans IAM 1603T, Pseudomonas gessardii CIP 105469T and Pseudomonas libanensis CIP 105460T showing 99.7-99.6 % similarity. Using rpoD gene sequences Pseudomonas veronii LMG 17761T (93.3 %) was most closely related to strain WS 4672T and Pseudomonas libanensis CIP 105460T to strain WS 4992T (93.3 %). The five strains could be differentiated from their closest relatives and from each other by phenotypic and chemotaxonomic characterization and ANIb values calculated from draft genome assemblies. ANIb values of strains WS 4992T and WS4671T to the closest relatives are lower than 90 %. The major cellular polar lipids of both strains are phosphatidylethanolamine, phosphatidylglycerol, a phospholipid and diphosphatidylglycerol, and their major quinone is Q-9. The DNA G+C content of strains WS 4992T and WS 4672T were 60.0  and 59.7  mol%, respectively. Based on these genotypic and phenotypic traits two novel species of the genus Pseudomonas are proposed: Pseudomonas lactis sp. nov. [with type strain WS 4992T (=DSM 29167T=LMG 28435T) and the additional strains WS 4997 and WS 5000], and Pseudomonasparalactis sp. nov. [with type strain WS 4672T (=DSM 29164T=LMG 28439T) and additional strain WS 4998].

Intranasal oxytocin increases facial expressivity, but not ratings of trustworthiness, in patients with schizophrenia and healthy controls

BACKGROUND

Blunted facial affect is a common negative symptom of schizophrenia. Additionally, assessing the trustworthiness of faces is a social cognitive ability that is impaired in schizophrenia. Currently available pharmacological agents are ineffective at improving either of these symptoms, despite their clinical significance. The hypothalamic neuropeptide oxytocin has multiple prosocial effects when administered intranasally to healthy individuals and shows promise in decreasing negative symptoms and enhancing social cognition in schizophrenia. Although two small studies have investigated oxytocin's effects on ratings of facial trustworthiness in schizophrenia, its effects on facial expressivity have not been investigated in any population.

METHOD

We investigated the effects of oxytocin on facial emotional expressivity while participants performed a facial trustworthiness rating task in 33 individuals with schizophrenia and 35 age-matched healthy controls using a double-blind, placebo-controlled, cross-over design. Participants rated the trustworthiness of presented faces interspersed with emotionally evocative photographs while being video-recorded. Participants' facial expressivity in these videos was quantified by blind raters using a well-validated manualized approach (i.e. the Facial Expression Coding System; FACES).

RESULTS

While oxytocin administration did not affect ratings of facial trustworthiness, it significantly increased facial expressivity in individuals with schizophrenia (Z = -2.33, p = 0.02) and at trend level in healthy controls (Z = -1.87, p = 0.06).

CONCLUSIONS

These results demonstrate that oxytocin administration can increase facial expressivity in response to emotional stimuli and suggest that oxytocin may have the potential to serve as a treatment for blunted facial affect in schizophrenia.

Simulating Intestinal Growth Conditions Enhances Toxin Production of Enteropathogenic Bacillus cereus

Bacillus cereus is a ubiquitous bacterial pathogen increasingly reported to be the causative agent of foodborne infections and intoxications. Since the enterotoxins linked to the diarrheal form of food poising are foremost produced in the human intestine, the toxic potential of enteropathogenic B. cereus strains is difficult to predict from studies carried out under routine cultivation procedures. In this study, toxigenic properties of a panel of strains (n = 19) of diverse origin were compared using cell culture medium pre-incubated with CaCo-2 cells to mimic intestinal growth conditions. Shortly after contact of the bacteria with the simulated host environment, enterotoxin gene expression was activated and total protein secretion of all strains was accelerated. Although the signal stimulating enterotoxin production still needs to be elucidated, it could be shown that it originated from the CaCo-2 cells. Overall, our study demonstrates that the currently used methods in B. cereus diagnostics, based on standard culture medium, are not allowing a conclusive prediction of the potential health risk related to a certain strain. Thus, these methods should be complemented by cultivation procedures that are simulating intestinal host conditions.

Differentiation of ncRNAs from small mRNAs in Escherichia coli O157:H7 EDL933 (EHEC) by combined RNAseq and RIBOseq - ryhB encodes the regulatory RNA RyhB and a peptide, RyhP

Background

While NGS allows rapid global detection of transcripts, it remains difficult to distinguish ncRNAs from short mRNAs. To detect potentially translated RNAs, we developed an improved protocol for bacterial ribosomal footprinting (RIBOseq). This allowed distinguishing ncRNA from mRNA in EHEC. A high ratio of ribosomal footprints per transcript (ribosomal coverage value, RCV) is expected to indicate a translated RNA, while a low RCV should point to a non-translated RNA.

Results

Based on their low RCV, 150 novel non-translated EHEC transcripts were identified as putative ncRNAs, representing both antisense and intergenic transcripts, 74 of which had expressed homologs in E. coli MG1655. Bioinformatics analysis predicted statistically significant target regulons for 15 of the intergenic transcripts; experimental analysis revealed 4-fold or higher differential expression of 46 novel ncRNA in different growth media. Out of 329 annotated EHEC ncRNAs, 52 showed an RCV similar to protein-coding genes, of those, 16 had RIBOseq patterns matching annotated genes in other enterobacteriaceae, and 11 seem to possess a Shine-Dalgarno sequence, suggesting that such ncRNAs may encode small proteins instead of being solely non-coding. To support that the RIBOseq signals are reflecting translation, we tested the ribosomal-footprint covered ORF of ryhB and found a phenotype for the encoded peptide in iron-limiting condition.

Conclusion

Determination of the RCV is a useful approach for a rapid first-step differentiation between bacterial ncRNAs and small mRNAs. Further, many known ncRNAs may encode proteins as well.

Electronic supplementary material

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

Spoilage of Microfiltered and Pasteurized Extended Shelf Life Milk Is Mainly Induced by Psychrotolerant Spore-Forming Bacteria that often Originate from Recontamination

Premature spoilage and varying product quality due to microbial contamination still constitute major problems in the production of microfiltered and pasteurized extended shelf life (ESL) milk. Spoilage-associated bacteria may enter the product either as part of the raw milk microbiota or as recontaminants in the dairy plant. To identify spoilage-inducing bacteria and their routes of entry, we analyzed end products for their predominant microbiota as well as the prevalence and biodiversity of psychrotolerant spores in bulk tank milk. Process analyses were performed to determine the removal of psychrotolerant spores at each production step. To detect transmission and recontamination events, strain typing was conducted with isolates obtained from all process stages. Microbial counts in 287 ESL milk packages at the end of shelf life were highly diverse ranging from <1 to 7.9 log cfu/mL. In total, 15% of samples were spoiled. High G+C Gram-positive bacteria were the most abundant taxonomic group, but were responsible for only 31% of spoilage. In contrast, psychrotolerant spores were isolated from 55% of spoiled packages. In 90% of samples with pure cultures of Bacillus cereus sensu lato and Paenibacillus spp., counts exceeded 6 log cfu/mL. In bulk tank milk, the concentration of psychrotolerant spores was low, accounting for merely 0.5 ± 0.8 MPN/mL. Paenibacillus amylolyticus/xylanexedens was by far the most dominant species in bulk tank milk (48% of all isolates), but was never detected in ESL milk, pointing to efficient removal during manufacturing. Six large-scale process analyses confirmed a high removal rate for psychrotolerant spores (reduction by nearly 4 log-units). B. cereus sensu lato, on the contrary, was frequently found in spoiled end products, but was rarely detected in bulk tank milk. Due to low counts in bulk tank samples and efficient spore removal during production, we suggest that shelf life is influenced only to a minor extent by raw-milk-associated factors. In contrast, recontamination with spores, particularly from the B. cereus complex, seems to occur. To enhance milk quality throughout the entire shelf life, improved plant sanitation and disinfection that target the elimination of spores are necessary.

Transcriptional and translational regulation by RNA thermometers, riboswitches and the sRNA DsrA in Escherichia coli O157:H7 Sakai under combined cold and osmotic stress adaptation

The enteric pathogen Escherichia coli O157:H7 Sakai (EHEC) is able to grow at lower temperatures compared to commensal E. coli Growth at environmental conditions displays complex challenges different to those in a host. EHEC was grown at 37°C and at 14°C with 4% NaCl, a combination of cold and osmotic stress as present in the food chain. Comparison of RNAseq and RIBOseq data provided a snap shot of ongoing transcription and translation, differentiating transcriptional and post-transcriptional gene regulation, respectively. Indeed, cold and osmotic stress related genes are simultaneously regulated at both levels, but translational regulation clearly dominates. Special emphasis was given to genes regulated by RNA secondary structures in their 5^'UTRs, such as RNA thermometers and riboswitches, or genes controlled by small RNAs encoded in trans The results reveal large differences in gene expression between short-time shock compared to adaptation in combined cold and osmotic stress. Whereas the majority of cold shock proteins, such as CspA, are translationally downregulated after adaptation, many osmotic stress genes are still significantly upregulated mainly translationally, but several also transcriptionally.

Permanent colonization of creek sediments, creek water and limnic water plants by four Listeria species in low population densities

During a 1-year longitudinal study, water, sediment and water plants from two creeks and one pond were sampled monthly and analyzed for the presence of Listeria species. A total of 90 % of 30 sediment samples, 84 % of 31 water plant samples and 67 % of 36 water samples were tested positive. Generally, most probable number counts ranged between 1 and 40 g-1, only occasionally >110 cfu g-1 were detected. Species differentiation based on FT-IR spectroscopy and multiplex PCR of a total of 1220 isolates revealed L. innocua (46 %), L. seeligeri (27 %), L. monocytogenes (25 %) and L. ivanovii (2 %). Titers and species compositions were similar during all seasons. While the species distributions in sediments and associated Ranunculus fluitans plants appeared to be similar in both creeks, RAPD typing did not provide conclusive evidence that the populations of these environments were connected. It is concluded that (i) the fresh-water sediments and water plants are year-round populated by Listeria, (ii) no clear preference for growth in habitats as different as sediments and water plants was found and (iii) the RAPD-based intraspecific biodiversity is high compared to the low population density.

Human genome meeting 2016 : Houston, TX, USA. 28 February - 2 March 2016

O1 The metabolomics approach to autism: identification of biomarkers for early detection of autism spectrum disorder A. K. Srivastava, Y. Wang, R. Huang, C. Skinner, T. Thompson, L. Pollard, T. Wood, F. Luo, R. Stevenson O2 Phenome-wide association study for smoking- and drinking-associated genes in 26,394 American women with African, Asian, European, and Hispanic descents R. Polimanti, J. Gelernter O3 Effects of prenatal environment, genotype and DNA methylation on birth weight and subsequent postnatal outcomes: findings from GUSTO, an Asian birth cohort X. Lin, I. Y. Lim, Y. Wu, A. L. Teh, L. Chen, I. M. Aris, S. E. Soh, M. T. Tint, J. L. MacIsaac, F. Yap, K. Kwek, S. M. Saw, M. S. Kobor, M. J. Meaney, K. M. Godfrey, Y. S. Chong, J. D. Holbrook, Y. S. Lee, P. D. Gluckman, N. Karnani, GUSTO study group O4 High-throughput identification of specific qt interval modulating enhancers at the SCN5A locus A. Kapoor, D. Lee, A. Chakravarti O5 Identification of extracellular matrix components inducing cancer cell migration in the supernatant of cultivated mesenchymal stem cells C. Maercker, F. Graf, M. Boutros O6 Single cell allele specific expression (ASE) IN T21 and common trisomies: a novel approach to understand DOWN syndrome and other aneuploidies G. Stamoulis, F. Santoni, P. Makrythanasis, A. Letourneau, M. Guipponi, N. Panousis, M. Garieri, P. Ribaux, E. Falconnet, C. Borel, S. E. Antonarakis O7 Role of microRNA in LCL to IPSC reprogramming S. Kumar, J. Curran, J. Blangero O8 Multiple enhancer variants disrupt gene regulatory network in Hirschsprung disease S. Chatterjee, A. Kapoor, J. Akiyama, D. Auer, C. Berrios, L. Pennacchio, A. Chakravarti O9 Metabolomic profiling for the diagnosis of neurometabolic disorders T. R. Donti, G. Cappuccio, M. Miller, P. Atwal, A. Kennedy, A. Cardon, C. Bacino, L. Emrick, J. Hertecant, F. Baumer, B. Porter, M. Bainbridge, P. Bonnen, B. Graham, R. Sutton, Q. Sun, S. Elsea O10 A novel causal methylation network approach to Alzheimer’s disease Z. Hu, P. Wang, Y. Zhu, J. Zhao, M. Xiong, David A Bennett O11 A microRNA signature identifies subtypes of triple-negative breast cancer and reveals MIR-342-3P as regulator of a lactate metabolic pathway A. Hidalgo-Miranda, S. Romero-Cordoba, S. Rodriguez-Cuevas, R. Rebollar-Vega, E. Tagliabue, M. Iorio, E. D’Ippolito, S. Baroni O12 Transcriptome analysis identifies genes, enhancer RNAs and repetitive elements that are recurrently deregulated across multiple cancer types B. Kaczkowski, Y. Tanaka, H. Kawaji, A. Sandelin, R. Andersson, M. Itoh, T. Lassmann, the FANTOM5 consortium, Y. Hayashizaki, P. Carninci, A. R. R. Forrest O13 Elevated mutation and widespread loss of constraint at regulatory and architectural binding sites across 11 tumour types C. A. Semple O14 Exome sequencing provides evidence of pathogenicity for genes implicated in colorectal cancer E. A. Rosenthal, B. Shirts, L. Amendola, C. Gallego, M. Horike-Pyne, A. Burt, P. Robertson, P. Beyers, C. Nefcy, D. Veenstra, F. Hisama, R. Bennett, M. Dorschner, D. Nickerson, J. Smith, K. Patterson, D. Crosslin, R. Nassir, N. Zubair, T. Harrison, U. Peters, G. Jarvik, NHLBI GO Exome Sequencing Project O15 The tandem duplicator phenotype as a distinct genomic configuration in cancer F. Menghi, K. Inaki, X. Woo, P. Kumar, K. Grzeda, A. Malhotra, H. Kim, D. Ucar, P. Shreckengast, K. Karuturi, J. Keck, J. Chuang, E. T. Liu O16 Modeling genetic interactions associated with molecular subtypes of breast cancer B. Ji, A. Tyler, G. Ananda, G. Carter O17 Recurrent somatic mutation in the MYC associated factor X in brain tumors H. Nikbakht, M. Montagne, M. Zeinieh, A. Harutyunyan, M. Mcconechy, N. Jabado, P. Lavigne, J. Majewski O18 Predictive biomarkers to metastatic pancreatic cancer treatment J. B. Goldstein, M. Overman, G. Varadhachary, R. Shroff, R. Wolff, M. Javle, A. Futreal, D. Fogelman O19 DDIT4 gene expression as a prognostic marker in several malignant tumors L. Bravo, W. Fajardo, H. Gomez, C. Castaneda, C. Rolfo, J. A. Pinto O20 Spatial organization of the genome and genomic alterations in human cancers K. C. Akdemir, L. Chin, A. Futreal, ICGC PCAWG Structural Alterations Group O21 Landscape of targeted therapies in solid tumors S. Patterson, C. Statz, S. Mockus O22 Genomic analysis reveals novel drivers and progression pathways in skin basal cell carcinoma S. N. Nikolaev, X. I. Bonilla, L. Parmentier, B. King, F. Bezrukov, G. Kaya, V. Zoete, V. Seplyarskiy, H. Sharpe, T. McKee, A. Letourneau, P. Ribaux, K. Popadin, N. Basset-Seguin, R. Ben Chaabene, F. Santoni, M. Andrianova, M. Guipponi, M. Garieri, C. Verdan, K. Grosdemange, O. Sumara, M. Eilers, I. Aifantis, O. Michielin, F. de Sauvage, S. Antonarakis O23 Identification of differential biomarkers of hepatocellular carcinoma and cholangiocarcinoma via transcriptome microarray meta-analysis S. Likhitrattanapisal O24 Clinical validity and actionability of multigene tests for hereditary cancers in a large multi-center study S. Lincoln, A. Kurian, A. Desmond, S. Yang, Y. Kobayashi, J. Ford, L. Ellisen O25 Correlation with tumor ploidy status is essential for correct determination of genome-wide copy number changes by SNP array T. L. Peters, K. R. Alvarez, E. F. Hollingsworth, D. H. Lopez-Terrada O26 Nanochannel based next-generation mapping for interrogation of clinically relevant structural variation A. Hastie, Z. Dzakula, A. W. Pang, E. T. Lam, T. Anantharaman, M. Saghbini, H. Cao, BioNano Genomics O27 Mutation spectrum in a pulmonary arterial hypertension (PAH) cohort and identification of associated truncating mutations in TBX4 C. Gonzaga-Jauregui, L. Ma, A. King, E. Berman Rosenzweig, U. Krishnan, J. G. Reid, J. D. Overton, F. Dewey, W. K. Chung O28 NORTH CAROLINA macular dystrophy (MCDR1): mutations found affecting PRDM13 K. Small, A. DeLuca, F. Cremers, R. A. Lewis, V. Puech, B. Bakall, R. Silva-Garcia, K. Rohrschneider, M. Leys, F. S. Shaya, E. Stone O29 PhenoDB and genematcher, solving unsolved whole exome sequencing data N. L. Sobreira, F. Schiettecatte, H. Ling, E. Pugh, D. Witmer, K. Hetrick, P. Zhang, K. Doheny, D. Valle, A. Hamosh O30 Baylor-Johns Hopkins Center for Mendelian genomics: a four year review S. N. Jhangiani, Z. Coban Akdemir, M. N. Bainbridge, W. Charng, W. Wiszniewski, T. Gambin, E. Karaca, Y. Bayram, M. K. Eldomery, J. Posey, H. Doddapaneni, J. Hu, V. R. Sutton, D. M. Muzny, E. A. Boerwinkle, D. Valle, J. R. Lupski, R. A. Gibbs O31 Using read overlap assembly to accurately identify structural genetic differences in an ashkenazi jewish trio S. Shekar, W. Salerno, A. English, A. Mangubat, J. Bruestle O32 Legal interoperability: a sine qua non for international data sharing A. Thorogood, B. M. Knoppers, Global Alliance for Genomics and Health - Regulatory and Ethics Working Group O33 High throughput screening platform of competent sineups: that can enhance translation activities of therapeutic target H. Takahashi, K. R. Nitta, A. Kozhuharova, A. M. Suzuki, H. Sharma, D. Cotella, C. Santoro, S. Zucchelli, S. Gustincich, P. Carninci O34 The undiagnosed diseases network international (UDNI): clinical and laboratory research to meet patient needs J. J. Mulvihill, G. Baynam, W. Gahl, S. C. Groft, K. Kosaki, P. Lasko, B. Melegh, D. Taruscio O36 Performance of computational algorithms in pathogenicity predictions for activating variants in oncogenes versus loss of function mutations in tumor suppressor genes R. Ghosh, S. Plon O37 Identification and electronic health record incorporation of clinically actionable pharmacogenomic variants using prospective targeted sequencing S. Scherer, X. Qin, R. Sanghvi, K. Walker, T. Chiang, D. Muzny, L. Wang, J. Black, E. Boerwinkle, R. Weinshilboum, R. Gibbs O38 Melanoma reprogramming state correlates with response to CTLA-4 blockade in metastatic melanoma T. Karpinets, T. Calderone, K. Wani, X. Yu, C. Creasy, C. Haymaker, M. Forget, V. Nanda, J. Roszik, J. Wargo, L. Haydu, X. Song, A. Lazar, J. Gershenwald, M. Davies, C. Bernatchez, J. Zhang, A. Futreal, S. Woodman O39 Data-driven refinement of complex disease classification from integration of heterogeneous functional genomics data in GeneWeaver E. J. Chesler, T. Reynolds, J. A. Bubier, C. Phillips, M. A. Langston, E. J. Baker O40 A general statistic framework for genome-based disease risk prediction M. Xiong, L. Ma, N. Lin, C. Amos O41 Integrative large-scale causal network analysis of imaging and genomic data and its application in schizophrenia studies N. Lin, P. Wang, Y. Zhu, J. Zhao, V. Calhoun, M. Xiong O42 Big data and NGS data analysis: the cloud to the rescue O. Dobretsberger, M. Egger, F. Leimgruber O43 Cpipe: a convergent clinical exome pipeline specialised for targeted sequencing S. Sadedin, A. Oshlack, Melbourne Genomics Health Alliance O44 A Bayesian classification of biomedical images using feature extraction from deep neural networks implemented on lung cancer data V. A. A. Antonio, N. Ono, Clark Kendrick C. Go O45 MAV-SEQ: an interactive platform for the Management, Analysis, and Visualization of sequence data Z. Ahmed, M. Bolisetty, S. Zeeshan, E. Anguiano, D. Ucar O47 Allele specific enhancer in EPAS1 intronic regions may contribute to high altitude adaptation of Tibetans C. Zeng, J. Shao O48 Nanochannel based next-generation mapping for structural variation detection and comparison in trios and populations H. Cao, A. Hastie, A. W. Pang, E. T. Lam, T. Liang, K. Pham, M. Saghbini, Z. Dzakula O49 Archaic introgression in indigenous populations of Malaysia revealed by whole genome sequencing Y. Chee-Wei, L. Dongsheng, W. Lai-Ping, D. Lian, R. O. Twee Hee, Y. Yunus, F. Aghakhanian, S. S. Mokhtar, C. V. Lok-Yung, J. Bhak, M. Phipps, X. Shuhua, T. Yik-Ying, V. Kumar, H. Boon-Peng O50 Breast and ovarian cancer prevention: is it time for population-based mutation screening of high risk genes? I. Campbell, M.-A. Young, P. James, Lifepool O53 Comprehensive coverage from low DNA input using novel NGS library preparation methods for WGS and WGBS C. Schumacher, S. Sandhu, T. Harkins, V. Makarov O54 Methods for large scale construction of robust PCR-free libraries for sequencing on Illumina HiSeqX platform H. DoddapaneniR. Glenn, Z. Momin, B. Dilrukshi, H. Chao, Q. Meng, B. Gudenkauf, R. Kshitij, J. Jayaseelan, C. Nessner, S. Lee, K. Blankenberg, L. Lewis, J. Hu, Y. Han, H. Dinh, S. Jireh, K. Walker, E. Boerwinkle, D. Muzny, R. Gibbs O55 Rapid capture methods for clinical sequencing J. Hu, K. Walker, C. Buhay, X. Liu, Q. Wang, R. Sanghvi, H. Doddapaneni, Y. Ding, N. Veeraraghavan, Y. Yang, E. Boerwinkle, A. L. Beaudet, C. M. Eng, D. M. Muzny, R. A. Gibbs O56 A diploid personal human genome model for better genomes from diverse sequence data K. C. C. Worley, Y. Liu, D. S. T. Hughes, S. C. Murali, R. A. Harris, A. C. English, X. Qin, O. A. Hampton, P. Larsen, C. Beck, Y. Han, M. Wang, H. Doddapaneni, C. L. Kovar, W. J. Salerno, A. Yoder, S. Richards, J. Rogers, J. R. Lupski, D. M. Muzny, R. A. Gibbs O57 Development of PacBio long range capture for detection of pathogenic structural variants Q. Meng, M. Bainbridge, M. Wang, H. Doddapaneni, Y. Han, D. Muzny, R. Gibbs O58 Rhesus macaques exhibit more non-synonymous variation but greater impact of purifying selection than humans R. A. Harris, M. Raveenedran, C. Xue, M. Dahdouli, L. Cox, G. Fan, B. Ferguson, J. Hovarth, Z. Johnson, S. Kanthaswamy, M. Kubisch, M. Platt, D. Smith, E. Vallender, R. Wiseman, X. Liu, J. Below, D. Muzny, R. Gibbs, F. Yu, J. Rogers O59 Assessing RNA structure disruption induced by single-nucleotide variation J. Lin, Y. Zhang, Z. Ouyang P1 A meta-analysis of genome-wide association studies of mitochondrial dna copy number A. Moore, Z. Wang, J. Hofmann, M. Purdue, R. Stolzenberg-Solomon, S. Weinstein, D. Albanes, C.-S. Liu, W.-L. Cheng, T.-T. Lin, Q. Lan, N. Rothman, S. Berndt P2 Missense polymorphic genetic combinations underlying down syndrome susceptibility E. S. Chen P4 The evaluation of alteration of ELAM-1 expression in the endometriosis patients H. Bahrami, A. Khoshzaban, S. Heidari Keshal P5 Obesity and the incidence of apolipoprotein E polymorphisms in an assorted population from Saudi Arabia population K. K. R. Alharbi P6 Genome-associated personalized antithrombotical therapy for patients with high risk of thrombosis and bleeding M. Zhalbinova, A. Akilzhanova, S. Rakhimova, M. Bekbosynova, S. Myrzakhmetova P7 Frequency of Xmn1 polymorphism among sickle cell carrier cases in UAE population M. Matar P8 Differentiating inflammatory bowel diseases by using genomic data: dimension of the problem and network organization N. Mili, R. Molinari, Y. Ma, S. Guerrier P9 Vulnerability of genetic variants to the risk of autism among Saudi children N. Elhawary, M. Tayeb, N. Bogari, N. Qotb P10 Chromatin profiles from ex vivo purified dopaminergic neurons establish a promising model to support studies of neurological function and dysfunction S. A. McClymont, P. W. Hook, L. A. Goff, A. McCallion P11 Utilization of a sensitized chemical mutagenesis screen to identify genetic modifiers of retinal dysplasia in homozygous Nr2e3^rd7 mice Y. Kong, J. R. Charette, W. L. Hicks, J. K. Naggert, L. Zhao, P. M. Nishina P12 Ion torrent next generation sequencing of recessive polycystic kidney disease in Saudi patients B. M. Edrees, M. Athar, F. A. Al-Allaf, M. M. Taher, W. Khan, A. Bouazzaoui, N. A. Harbi, R. Safar, H. Al-Edressi, A. Anazi, N. Altayeb, M. A. Ahmed, K. Alansary, Z. Abduljaleel P13 Digital expression profiling of Purkinje neurons and dendrites in different subcellular compartments A. Kratz, P. Beguin, S. Poulain, M. Kaneko, C. Takahiko, A. Matsunaga, S. Kato, A. M. Suzuki, N. Bertin, T. Lassmann, R. Vigot, P. Carninci, C. Plessy, T. Launey P14 The evolution of imperfection and imperfection of evolution: the functional and functionless fractions of the human genome D. Graur P16 Species-independent identification of known and novel recurrent genomic entities in multiple cancer patients J. Friis-Nielsen, J. M. Izarzugaza, S. Brunak P18 Discovery of active gene modules which are densely conserved across multiple cancer types reveal their prognostic power and mutually exclusive mutation patterns B. S. Soibam P19 Whole exome sequencing of dysplastic leukoplakia tissue indicates sequential accumulation of somatic mutations from oral precancer to cancer D. Das, N. Biswas, S. Das, S. Sarkar, A. Maitra, C. Panda, P. Majumder P21 Epigenetic mechanisms of carcinogensis by hereditary breast cancer genes J. J. Gruber, N. Jaeger, M. Snyder P22 RNA direct: a novel RNA enrichment strategy applied to transcripts associated with solid tumors K. Patel, S. Bowman, T. Davis, D. Kraushaar, A. Emerman, S. Russello, N. Henig, C. Hendrickson P23 RNA sequencing identifies gene mutations for neuroblastoma K. Zhang P24 Participation of SFRP1 in the modulation of TMPRSS2-ERG fusion gene in prostate cancer cell lines M. Rodriguez-Dorantes, C. D. Cruz-Hernandez, C. D. P. Garcia-Tobilla, S. Solorzano-Rosales P25 Targeted Methylation Sequencing of Prostate Cancer N. Jäger, J. Chen, R. Haile, M. Hitchins, J. D. Brooks, M. Snyder P26 Mutant TPMT alleles in children with acute lymphoblastic leukemia from México City and Yucatán, Mexico S. Jiménez-Morales, M. Ramírez, J. Nuñez, V. Bekker, Y. Leal, E. Jiménez, A. Medina, A. Hidalgo, J. Mejía P28 Genetic modifiers of Alström syndrome J. Naggert, G. B. Collin, K. DeMauro, R. Hanusek, P. M. Nishina P31 Association of genomic variants with the occurrence of angiotensin-converting-enzyme inhibitor (ACEI)-induced coughing among Filipinos E. M. Cutiongco De La Paz, R. Sy, J. Nevado, P. Reganit, L. Santos, J. D. Magno, F. E. Punzalan , D. Ona , E. Llanes, R. L. Santos-Cortes , R. Tiongco, J. Aherrera, L. Abrahan, P. Pagauitan-Alan; Philippine Cardiogenomics Study Group P32 The use of “humanized” mouse models to validate disease association of a de novo GARS variant and to test a novel gene therapy strategy for Charcot-Marie-Tooth disease type 2D K. H. Morelli, J. S. Domire, N. Pyne, S. Harper, R. Burgess P34 Molecular regulation of chondrogenic human induced pluripotent stem cells M. A. Gari, A. Dallol, H. Alsehli, A. Gari, M. Gari, A. Abuzenadah P35 Molecular profiling of hematologic malignancies: implementation of a variant assessment algorithm for next generation sequencing data analysis and clinical reporting M. Thomas, M. Sukhai, S. Garg, M. Misyura, T. Zhang, A. Schuh, T. Stockley, S. Kamel-Reid P36 Accessing genomic evidence for clinical variants at NCBI S. Sherry, C. Xiao, D. Slotta, K. Rodarmer, M. Feolo, M. Kimelman, G. Godynskiy, C. O’Sullivan, E. Yaschenko P37 NGS-SWIFT: a cloud-based variant analysis framework using control-accessed sequencing data from DBGAP/SRA C. Xiao, E. Yaschenko, S. Sherry P38 Computational assessment of drug induced hepatotoxicity through gene expression profiling C. Rangel-Escareño, H. Rueda-Zarate P40 Flowr: robust and efficient pipelines using a simple language-agnostic approach;ultraseq; fast modular pipeline for somatic variation calling using flowr S. Seth, S. Amin, X. Song, X. Mao, H. Sun, R. G. Verhaak, A. Futreal, J. Zhang P41 Applying “Big data” technologies to the rapid analysis of heterogenous large cohort data S. J. Whiite, T. Chiang, A. English, J. Farek, Z. Kahn, W. Salerno, N. Veeraraghavan, E. Boerwinkle, R. Gibbs P42 FANTOM5 web resource for the large-scale genome-wide transcription start site activity profiles of wide-range of mammalian cells T. Kasukawa, M. Lizio, J. Harshbarger, S. Hisashi, J. Severin, A. Imad, S. Sahin, T. C. Freeman, K. Baillie, A. Sandelin, P. Carninci, A. R. R. Forrest, H. Kawaji, The FANTOM Consortium P43 Rapid and scalable typing of structural variants for disease cohorts W. Salerno, A. English, S. N. Shekar, A. Mangubat, J. Bruestle, E. Boerwinkle, R. A. Gibbs P44 Polymorphism of glutathione S-transferases and sulphotransferases genes in an Arab population A. H. Salem, M. Ali, A. Ibrahim, M. Ibrahim P46 Genetic divergence of CYP3A5*3 pharmacogenomic marker for native and admixed Mexican populations J. C. Fernandez-Lopez, V. Bonifaz-Peña, C. Rangel-Escareño, A. Hidalgo-Miranda, A. V. Contreras P47 Whole exome sequence meta-analysis of 13 white blood cell, red blood cell, and platelet traits L. Polfus, CHARGE and NHLBI Exome Sequence Project Working Groups P48 Association of adipoq gene with type 2 diabetes and related phenotypes in african american men and women: The jackson heart study S. Davis, R. Xu, S. Gebeab, P Riestra, A Gaye, R. Khan, J. Wilson, A. Bidulescu P49 Common variants in casr gene are associated with serum calcium levels in koreans S. H. Jung, N. Vinayagamoorthy, S. H. Yim, Y. J. Chung P50 Inference of multiple-wave population admixture by modeling decay of linkage disequilibrium with multiple exponential functions Y. Zhou, S. Xu P51 A Bayesian framework for generalized linear mixed models in genome-wide association studies X. Wang, V. Philip, G. Carter P52 Targeted sequencing approach for the identification of the genetic causes of hereditary hearing impairment A. A. Abuzenadah, M. Gari, R. Turki, A. Dallol P53 Identification of enhancer sequences by ATAC-seq open chromatin profiling A. Uyar, A. Kaygun, S. Zaman, E. Marquez, J. George, D. Ucar P54 Direct enrichment for the rapid preparation of targeted NGS libraries C. L. Hendrickson, A. Emerman, D. Kraushaar, S. Bowman, N. Henig, T. Davis, S. Russello, K. Patel P56 Performance of the Agilent D5000 and High Sensitivity D5000 ScreenTape assays for the Agilent 4200 Tapestation System R. Nitsche, L. Prieto-Lafuente P57 ClinVar: a multi-source archive for variant interpretation M. Landrum, J. Lee, W. Rubinstein, D. Maglott P59 Association of functional variants and protein physical interactions of human MUTY homolog linked with familial adenomatous polyposis and colorectal cancer syndrome Z. Abduljaleel, W. Khan, F. A. Al-Allaf, M. Athar , M. M. Taher, N. Shahzad P60 Modification of the microbiom constitution in the gut using chicken IgY antibodies resulted in a reduction of acute graft-versus-host disease after experimental bone marrow transplantation A. Bouazzaoui, E. Huber, A. Dan, F. A. Al-Allaf, W. Herr, G. Sprotte, J. Köstler, A. Hiergeist, A. Gessner, R. Andreesen, E. Holler P61 Compound heterozygous mutation in the LDLR gene in Saudi patients suffering severe hypercholesterolemia F. Al-Allaf, A. Alashwal, Z. Abduljaleel, M. Taher, A. Bouazzaoui, H. Abalkhail, A. Al-Allaf, R. Bamardadh, M. Athar

Comparative Bioinformatics and Experimental Analysis of the Intergenic Regulatory Regions of Bacillus cereus hbl and nhe Enterotoxin Operons and the Impact of CodY on Virulence Heterogeneity

Bacillus cereus is a food contaminant with greatly varying enteropathogenic potential. Almost all known strains harbor the genes for at least one of the three enterotoxins Nhe, Hbl, and CytK. While some strains show no cytotoxicity, others have caused outbreaks, in rare cases even with lethal outcome. The reason for these differences in cytotoxicity is unknown. To gain insight into the origin of enterotoxin expression heterogeneity in different strains, the architecture and role of 5′ intergenic regions (5′ IGRs) upstream of the nhe and hbl operons was investigated. In silico comparison of 142 strains of all seven phylogenetic groups of B. cereus sensu lato proved the presence of long 5′ IGRs upstream of the nheABC and hblCDAB operons, which harbor recognition sites for several transcriptional regulators, including the virulence regulator PlcR, redox regulators ResD and Fnr, the nutrient-sensitive regulator CodY as well as the master regulator for biofilm formation SinR. By determining transcription start sites, unusually long 5′ untranslated regions (5′ UTRs) upstream of the nhe and hbl start codons were identified, which are not present upstream of cytK-1 and cytK-2. Promoter fusions lacking various parts of the nhe and hbl 5′ UTR in B. cereus INRA C3 showed that the entire 331 bp 5′ UTR of nhe is necessary for full promoter activity, while the presence of the complete 606 bp hbl 5′ UTR lowers promoter activity. Repression was caused by a 268 bp sequence directly upstream of the hbl transcription start. Luciferase activity of reporter strains containing nhe and hbl 5′ IGR lux fusions provided evidence that toxin gene transcription is upregulated by the depletion of free amino acids. Electrophoretic mobility shift assays showed that the branched-chain amino acid sensing regulator CodY binds to both nhe and hbl 5′ UTR downstream of the promoter, potentially acting as a nutrient-responsive roadblock repressor of toxin gene transcription. PlcR binding sites are highly conserved among all B. cereus sensu lato strains, indicating that this regulator does not significantly contribute to the heterogeneity in virulence potentials. The CodY recognition sites are far less conserved, perhaps conferring varying strengths of CodY binding, which might modulate toxin synthesis in a strain-specific manner.

Optimized Illumina PCR-free library preparation for bacterial whole genome sequencing and analysis of factors influencing de novo assembly

Background

Next-generation sequencing (NGS) technology has paved the way for rapid and cost-efficient de novo sequencing of bacterial genomes. In particular, the introduction of PCR-free library preparation procedures (LPPs) lead to major improvements as PCR bias is largely reduced. However, in order to facilitate the assembly of Illumina paired-end sequence data and to enhance assembly performance, an increase of insert sizes to facilitate the repeat bridging and resolution capabilities of current state of the art assembly tools is needed. In addition, information concerning the relationships between genomic GC content, library insert size and sequencing quality as well as the influence of library insert size, read length and sequencing depth on assembly performance would be helpful to specifically target sequencing projects.

Results

Optimized DNA fragmentation settings and fine-tuned resuspension buffer to bead buffer ratios during fragment size selection were integrated in the Illumina TruSeq^® DNA PCR-free LPP in order to produce sequencing libraries varying in average insert size for bacterial genomes within a range of 35.4–73.0 % GC content. The modified protocol consumes only half of the reagents per sample, thus doubling the number of preparations possible with a kit. Examination of different libraries revealed that sequencing quality decreases with increased genomic GC content and with larger insert sizes. The estimation of assembly performance using assembly metrics like corrected NG50 and NGA50 showed that libraries with larger insert sizes can result in substantial assembly improvements as long as appropriate assembly tools are chosen. However, such improvements seem to be limited to genomes with a low to medium GC content. A positive trend between read length and assembly performance was observed while sequencing depth is less important, provided a minimum coverage is reached.

Conclusions

Based on the optimized protocol developed, sequencing libraries with flexible insert sizes and lower reagent costs can be generated. Furthermore, increased knowledge about the interplay of sequencing quality, insert size, genomic GC content, read length, sequencing depth and the assembler used will help molecular biologists to set up an optimal experimental and analytical framework with respect to Illumina next-generation sequencing of bacterial genomes.

Electronic supplementary material

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