Menu

Auto Tag: Bioinformatics

September 22, 2019

Adaptation of Pseudomonas aeruginosa to phage PaP1 predation via O-antigen polymerase mutation.

Adaptation of bacteria to phage predation poses a major obstacle for phage therapy. Bacteria adopt multiple mechanisms, such as inhibition of phage adsorption and CRISPR/Cas systems, to resist phage infection. Here, a phage-resistant mutant of Pseudomonas aeruginosa strain PA1 under the infection of lytic phage PaP1 was selected for further study. The PaP1-resistant variant, termed PA1RG, showed decreased adsorption to PaP1 and was devoid of long chain O-antigen on its cell envelope. Whole genome sequencing and comparative analysis revealed a single nucleotide mutation in the gene PA1S_08510, which encodes the O-antigen polymerase Wzy that is involved in lipopolysaccharide (LPS) biosynthesis. PA1_Wzy was classified into the O6 serotype based on sequence homology analysis and adopts a transmembrane topology similar to that seem with P. aeruginosa strain PAO1. Complementation of gene wzy in trans enabled the mutant PA1RG to produce the normal LPS pattern with long chain O-antigen and restored the susceptibility of PA1RG to phage PaP1 infection. While wzy mutation did not affect bacterial growth, mutant PA1RG exhibited decreased biofilm production, suggesting a fitness cost of PA1 associated with resistance of phage PaP1 predation. This study uncovered the mechanism responsible for PA1RG resistance to phage PaP1 via wzy mutation and revealed the role of phages in regulating bacterial behavior.

Read more
September 22, 2019

Assembly of chromosome-scale contigs by efficiently resolving repetitive sequences with long reads

Due to the large number of repetitive sequences in complex eukaryotic genomes, fragmented and incompletely assembled genomes lose value as reference sequences, often due to short contigs that cannot be anchored or mispositioned onto chromosomes. Here we report a novel method Highly Efficient Repeat Assembly (HERA), which includes a new concept called a connection graph as well as algorithms for constructing the graph. HERA resolves repeats at high efficiency with single-molecule sequencing data, and enables the assembly of chromosome-scale contigs by further integrating genome maps and Hi-C data. We tested HERA with the genomes of rice R498, maize B73, human HX1 and Tartary buckwheat Pinku1. HERA can correctly assemble most of the tandemly repetitive sequences in rice using single-molecule sequencing data only. Using the same maize and human sequencing data published by Jiao et al. (2017) and Shi et al. (2016), respectively, we dramatically improved on the sequence contiguity compared with the published assemblies, increasing the contig N50 from 1.3 Mb to 61.2 Mb in maize B73 assembly and from 8.3 Mb to 54.4 Mb in human HX1 assembly with HERA. We provided a high-quality maize reference genome with 96.9% of the gaps filled (only 76 gaps left) and several incorrectly positioned sequences fixed compared with the B73 RefGen_v4 assembly. Comparisons between the HERA assembly of HX1 and the human GRCh38 reference genome showed that many gaps in GRCh38 could be filled, and that GRCh38 contained some potential errors that could be fixed. We assembled the Pinku1 genome into 12 scaffolds with a contig N50 size of 27.85 Mb. HERA serves as a new genome assembly/phasing method to generate high quality sequences for complex genomes and as a curation tool to improve the contiguity and completeness of existing reference genomes, including the correction of assembly errors in repetitive regions.

Read more
September 22, 2019

C-di-GMP turnover influences motility and biofilm formation in Bacillus amyloliquefaciens PG12.

Bis-(3′?5′) cyclic dimeric guanosine monophosphate (c-di-GMP) is defined as a highly versatile secondary messenger in bacteria, coordinating diverse aspects of bacterial growth and behavior, including motility and biofilm formation. Bacillus amyloliquefaciens PG12 is an effective biocontrol agent against apple ring rot caused by Botryosphaeria dothidea. In this study, we characterized the core regulators of c-di-GMP turnover in B. amyloliquefaciens PG12. Using bioinformatic analysis, heterologous expression and biochemical characterization of knockout and overexpression derivatives, we identified and characterized two active diguanylate cyclases (which catalyze c-di-GMP biosynthesis), YhcK and YtrP and one active c-di-GMP phosphodiesterase (which degrades c-di-GMP), YuxH. Furthermore, we showed that elevating c-di-GMP levels up to a certain threshold inhibited the swimming motility of B. amyloliquefaciens PG12. Although yhcK, ytrP and yuxH knockout mutants did not display defects in biofilm formation, significant increases in c-di-GMP levels induced by YtrP or YuxH overexpression stimulated biofilm formation in B. amyloliquefaciens PG12. Our results indicate that B. amyloliquefaciens possesses a functional c-di-GMP signaling system that influences the bacterium’s motility and ability to form biofilms. Since motility and biofilm formation influence the efficacy of biological control agent, our work provides a basis for engineering a more effective strain of B. amyloliquefaciens PG12. Copyright © 2018 Institut Pasteur. Published by Elsevier Masson SAS. All rights reserved.

Read more
September 22, 2019

Saliniramus fredricksonii gen. nov., sp. nov., a heterotrophic halophile isolated from Hot Lake, Washington, a member of a novel lineage (Salinarimonadaceae fam. nov.) within the order Rhizobiales, and reclassification of the genus Salinarimonas Liu et al. 2010 into Salinarimonadaceae.

There was an error in the proposed genus name in the published article, in that the genus ‘Salinivirga’ was effectively published while this article was in review. Therefore, the genus ‘Salinivirga’ should be replaced with ‘Saliniramus’. For the convenience of future readers, we have included the complete corrected article below, in which all occurrences of the incorrect genus name have been amended: A halophilic bacterial strain, HL-109T, was isolated from the unicyanobacterial consortium UCC-O, which was obtained from the photosynthetic mat of Hot Lake (Washington, USA). A polyphasic approach using phenotypic, genotypic and chemotaxonomic data was used to classify the strain within the order Rhizobiales. The organism stained Gram-negative and was a moderate thermophile with a growth optimum of 45?°C. It was obligately aerobic, heterotrophic and halophilic, growing in both NaCl and MgSO4 brines. The novel isolate had a polymorphic cellular morphology of short rods with occasional branching, and cells were monotrichous. The major fatty acids detected were C18?:?1, C18?:?0, C16?:?0 and C18?:?cyc. Phylogenetic analysis of the 16S rRNA gene placed the strain in the order Rhizobiales and it shared 94?% identity with the type strain of its nearest relative, Salinarimonas ramus. Morphological, chemotaxonomic and phylogenetic results did not affiliate the novel organism with any of the families in the Rhizobiales; therefore, HL-109T is representative of a new lineage, for which the name Saliniramus fredricksonii gen. nov., sp. nov. is proposed, with the type strain HL-109T (=JCM 31876T=DSM 102886T). In addition, examination of the phylogenetics of strain HL-109T and its nearest relatives, Salinarimonas ramus and Salinarimonasrosea, demonstrates that these halophiles form a clade distinct from the described families of the Rhizobiales. We further propose the establishment of a new family, Salinarimonadaceae fam. nov., to accommodate the genera Saliniramus and Salinarimonas (the type genus of the family).

Read more
September 22, 2019

Recurrent loss, horizontal transfer, and the obscure origins of mitochondrial introns in diatoms (Bacillariophyta).

We sequenced mitochondrial genomes from five diverse diatoms (Toxarium undulatum, Psammoneis japonica, Eunotia naegelii, Cylindrotheca closterium, and Nitzschia sp.), chosen to fill important phylogenetic gaps and help us characterize broadscale patterns of mitochondrial genome evolution in diatoms. Although gene content was strongly conserved, intron content varied widely across species. The vast majority of introns were of group II type and were located in the cox1 or rnl genes. Although recurrent intron loss appears to be the principal underlying cause of the sporadic distributions of mitochondrial introns across diatoms, phylogenetic analyses showed that intron distributions superficially consistent with a recurrent-loss model were sometimes more complicated, implicating horizontal transfer as a likely mechanism of intron acquisition as well. It was not clear, however, whether diatoms were the donors or recipients of horizontally transferred introns, highlighting a general challenge in resolving the evolutionary histories of many diatom mitochondrial introns. Although some of these histories may become clearer as more genomes are sampled, high rates of intron loss suggest that the origins of many diatom mitochondrial introns are likely to remain unclear.

Read more
September 22, 2019

Whole genome sequencing of Klebsiella pneumoniae strain unravels a new model for the development of extensive drug resistance in Enterobacteriaceae.

Increased incidence of carbapenem-resistant Enterobacteriaceae (CRE) has been reported worldwide. The WHO warns about the imminent risk to global health if the spread of resistant bacteria is not contained.Here, single molecule real time sequencing was used to analyse the whole genome and resistome of SKGH01, a strain of Klebsiella pneumoniae .The data showed that SKGH01 was resistant to all commercially available antibiotics. A complete account of extensively drug-resistant (XDR) CRE at a genomic level and the entire location map of all antibiotic resistance components are here presented. Additionally, this work proposes a model of XDR acquisition in Enterobacteriaceae.

Read more
September 22, 2019

High-Resolution Full-Length HLA Typing Method Using Third Generation (Pac-Bio SMRT) Sequencing Technology.

The human HLA genes are among the most polymorphic genes in the human genome. Therefore, it is very difficult to find two unrelated individuals with identical HLA molecules. As a result, HLA Class I and Class II genes are routinely sequenced or serotyped for organ transplantation, autoimmune disease-association studies, drug hypersensitivity research, and other applications. However, these methods were able to give two or four digit data, which was not sufficient enough to understand the completeness of haplotypes of HLA genes. To overcome these limitations, we here described end-to-end workflow for sequencing of HLA class I and class II genes using third generation sequencing, SMRT technology. This method produces fully-phased, unambiguous, allele-level information on the PacBio System.

Read more
September 22, 2019

Genome Assembly.

Genome assembly uses sequence similarity to go from sequencing reads to longer contiguous sequences (contigs). Scaffolds are contigs linked together by gaps where the order and orientation of the contigs is known but the exact sequence connecting two contigs is unknown, represented by Ns which estimate the gap length. Here we describe recommendations for genome assembly for different sequencing technologies, describe organelle assembly, and review how to perform assembly quality control.

Read more
September 22, 2019

Fungal Epigenomics: Detection and Analysis.

Across Eukaryota, DNA modifications play an important role in regulation of gene expression. While 5-methylcytosine (5mC) has been explored in depth, other modifications such as 6-methyladenine (6 mA) have historically been overlooked, in part due to technical difficulties in collecting/analyzing these data. However, recent technological advances have enabled exploration of these marks with much greater detail and on a larger scale. In this chapter, we discuss multiple methods for identifying and analyzing both 5mC and 6 mA across fungi.

Read more
September 22, 2019

Computational Modeling of Multidrug-Resistant Bacteria

Understanding how complex phenotypes arise from individual molecules and their interactions is a primary challenge in biology, and computational approaches have been increasingly employed to tackle this task. In this chapter, we describe current efforts by FIOCRUZ and partners to develop integrated computational models of multidrug-resistant bacteria. The bacterium chosen as the main focus of this effort is Pseudomonas aeruginosa, an opportunistic pathogen associated with a broad spectrum of infections in humans. Nowadays, P. aeruginosa is one of the main problems of healthcare-associated infections (HAI) in the world, because of its great capacity of survival in hospital environments and its intrinsic resistance to many antibiotics. Our overall research objective is to use integrated computational models to accurately predict a wide range of observable cellular behaviors of multidrug-resistant P. aeruginosa CCBH4851, which is a strain belonging to the clone ST277, endemic in Brazil. In this chapter, after a brief introduction to P. aeruginosa biology, we discuss the construction of metabolic and gene regulatory networks of P. aeruginosa CCBH 4851 from its genome. We also illustrate how these networks can be integrated into a single model, and we discuss methods for identifying potential therapeutic targets through integrated models.

Read more
September 22, 2019

Long-read whole genome sequencing and comparative analysis of six strains of the human pathogen Orientia tsutsugamushi.

Orientia tsutsugamushi is a clinically important but neglected obligate intracellular bacterial pathogen of the Rickettsiaceae family that causes the potentially life-threatening human disease scrub typhus. In contrast to the genome reduction seen in many obligate intracellular bacteria, early genetic studies of Orientia have revealed one of the most repetitive bacterial genomes sequenced to date. The dramatic expansion of mobile elements has hampered efforts to generate complete genome sequences using short read sequencing methodologies, and consequently there have been few studies of the comparative genomics of this neglected species.We report new high-quality genomes of O. tsutsugamushi, generated using PacBio single molecule long read sequencing, for six strains: Karp, Kato, Gilliam, TA686, UT76 and UT176. In comparative genomics analyses of these strains together with existing reference genomes from Ikeda and Boryong strains, we identify a relatively small core genome of 657 genes, grouped into core gene islands and separated by repeat regions, and use the core genes to infer the first whole-genome phylogeny of Orientia.Complete assemblies of multiple Orientia genomes verify initial suggestions that these are remarkable organisms. They have larger genomes compared with most other Rickettsiaceae, with widespread amplification of repeat elements and massive chromosomal rearrangements between strains. At the gene level, Orientia has a relatively small set of universally conserved genes, similar to other obligate intracellular bacteria, and the relative expansion in genome size can be accounted for by gene duplication and repeat amplification. Our study demonstrates the utility of long read sequencing to investigate complex bacterial genomes and characterise genomic variation.

Read more
September 22, 2019

The mutation rate and the age of the sex chromosomes in Silene latifolia.

Many aspects of sex chromosome evolution are common to both plants and animals [1], but the process of Y chromosome degeneration, where genes on the Y become non-functional over time, may be much slower in plants due to purifying selection against deleterious mutations in the haploid gametophyte [2, 3]. Testing for differences in Y degeneration between the kingdoms has been hindered by the absence of accurate age estimates for plant sex chromosomes. Here, we used genome resequencing to estimate the spontaneous mutation rate and the age of the sex chromosomes in white campion (Silene latifolia). Screening of single nucleotide polymorphisms (SNPs) in parents and 10 F1 progeny identified 39 de novo mutations and yielded a rate of 7.31 × 10-9 (95% confidence interval: 5.20 × 10-9 – 8.00 × 10-9) mutations per site per haploid genome per generation. Applying this mutation rate to the synonymous divergence between homologous X- and Y-linked genes (gametologs) gave age estimates of 11.00 and 6.32 million years for the old and young strata, respectively. Based on SNP segregation patterns, we inferred which genes were Y-linked and found that at least 47% are already dysfunctional. Applying our new estimates for the age of the sex chromosomes indicates that the rate of Y degeneration in S. latifolia is nearly 2-fold slower when compared to animal sex chromosomes of a similar age. Our revised estimates support Y degeneration taking place more slowly in plants, a discrepancy that may be explained by differences in the life cycles of animals and plants. Copyright © 2018 Elsevier Ltd. All rights reserved.

Read more
September 22, 2019

Homogenization of sub-genome secretome gene expression patterns in the allodiploid fungus Verticillium longisporum

Allopolyploidization, genome duplication through interspecific hybridization, is an important evolutionary mechanism that can enable organisms to adapt to environmental changes or stresses. The increased adaptive potential of allopolyploids can be particularly relevant for plant pathogens in their ongoing quest for host immune response evasion. To this end, plant pathogens secrete a plethora of molecules that enable host colonization. Allodiploidization has resulted in the new plant pathogen Verticillium longisporum that infects different hosts than haploid Verticillium species. To reveal the impact of allodiploidization on plant pathogen evolution, we studied the genome and transcriptome dynamics of V. longisporum using next-generation sequencing. V. longisporum genome evolution is characterized by extensive chromosomal rearrangements, between as well as within parental chromosome sets, leading to a mosaic genome structure. In comparison to haploid Verticillium species, V. longisporum genes display stronger signs of positive selection. The expression patterns of the two sub-genomes show remarkable resemblance, suggesting that the parental gene expression patterns homogenized upon hybridization. Moreover, whereas V. longisporum genes encoding secreted proteins frequently display differential expression between the parental sub-genomes in culture medium, expression patterns homogenize upon plant colonization. Collectively, our results illustrate of the adaptive potential of allodiploidy mediated by the interplay of two sub-genomes. Author summary Hybridization followed by whole-genome duplication, so-called allopolyploidization, provides genomic flexibility that is beneficial for survival under stressful conditions or invasiveness into new habitats. Allopolyploidization has mainly been studied in plants, but also occurs in other organisms, including fungi. Verticillium longisporum, an emerging fungal pathogen on brassicaceous plants, arose by allodiploidization between two Verticillium spp. We used comparative genomics to reveal the plastic nature of the V. longisporum genomes, showing that parental chromosome sets recombined extensively, resulting in a mosaic genome pattern. Furthermore, we show that non-synonymous substitutions frequently occurred in V. longisporum. Moreover, we reveal that expression patterns of genes encoding secreted proteins homogenized between the V. longisporum sub-genomes upon plant colonization. In conclusion, our results illustrate the large adaptive potential upon genome hybridization for fungi mediated by genomic plasticity and interaction between sub-genomes.

Read more
September 22, 2019

Improved de novo genome assembly and analysis of the Chinese cucurbit Siraitia grosvenorii, also known as monk fruit or luo-han-guo.

Luo-han-guo (Siraitia grosvenorii), also called monk fruit, is a member of the Cucurbitaceae family. Monk fruit has become an important area for research because of the pharmacological and economic potential of its noncaloric, extremely sweet components (mogrosides). It is also commonly used in traditional Chinese medicine for the treatment of lung congestion, sore throat, and constipation. Recently, a single reference genome became available for monk fruit, assembled from 36.9x genome coverage reads via Illumina sequencing platforms. This genome assembly has a relatively short (34.2 kb) contig N50 length and lacks integrated annotations. These drawbacks make it difficult to use as a reference in assembling transcriptomes and discovering novel functional genes.Here, we offer a new high-quality draft of the S. grosvenorii genome assembled using 31 Gb (~73.8x) long single molecule real time sequencing reads and polished with ~50 Gb Illumina paired-end reads. The final genome assembly is approximately 469.5 Mb, with a contig N50 length of 432,384 bp, representing a 12.6-fold improvement. We further annotated 237.3 Mb of repetitive sequence and 30,565 consensus protein coding genes with combined evidence. Phylogenetic analysis showed that S. grosvenorii diverged from members of the Cucurbitaceae family approximately 40.9 million years ago. With comprehensive transcriptomic analysis and differential expression testing, we identified 4,606 up-regulated genes in the early fruit compared to the leaf, a number of which were linked to metabolic pathways regulating fruit development and ripening.The availability of this new monk fruit genome assembly, as well as the annotations, will facilitate the discovery of new functional genes and the genetic improvement of monk fruit.

Read more
September 22, 2019

Transcriptional regulation of cysteine and methionine metabolism in Lactobacillus paracasei FAM18149.

Lactobacillus paracasei is common in the non-starter lactic acid bacteria (LAB) community of raw milk cheeses. This species can significantly contribute to flavor formation through amino acid metabolism. In this study, the DNA and RNA of L. paracasei FAM18149 were sequenced using next-generation sequencing technologies to reconstruct the metabolism of the sulfur-containing amino acids cysteine and methionine. Twenty-three genes were found to be involved in cysteine biosynthesis, the conversion of cysteine to methionine and vice versa, the S-adenosylmethionine recycling pathway, and the transport of sulfur-containing amino acids. Additionally, six methionine-specific T-boxes and one cysteine-specific T-box were found. Five of these were located upstream of genes encoding transporter functions. RNA-seq analysis and reverse-transcription quantitative polymerase reaction assays showed that expression of genes located downstream of these T-boxes was affected by the absence of either cysteine or methionine. Remarkably, the cysK2-ctl1-cysE2 operon, which is associated with te methionine-to-cysteine conversion and is upregulated in the absence of cysteine, showed high read coverage in the 5′-untranslated region and an antisense-RNA in the 3′-untranslated region. This indicates that this operon is regulated by the combination of cis- and antisense-mediated regulation mechanisms. The results of this study may help in the selection of L. paracasei strains to control sulfuric flavor formation in cheese.

Read more

Subscribe for blog updates:

Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.