Menu

Auto Tag: Bioinformatics

September 22, 2019

Constant conflict between Gypsy LTR retrotransposons and CHH methylation within a stress-adapted mangrove genome.

The evolutionary dynamics of the conflict between transposable elements (TEs) and their host genome remain elusive. This conflict will be intense in stress-adapted plants as stress can often reactivate TEs. Mangroves reduce TE load convergently in their adaptation to intertidal environments and thus provide a unique opportunity to address the host-TE conflict and its interaction with stress adaptation. Using the mangrove Rhizophora apiculata as a model, we investigated methylation and short interfering RNA (siRNA) targeting patterns in relation to the abundance and age of long terminal repeat (LTR) retrotransposons. We also examined the distance of LTR retrotransposons to genes, the impact on neighboring gene expression and population frequencies. We found differential accumulation amongst classes of LTR retrotransposons despite high overall methylation levels. This can be attributed to 24-nucleotide siRNA-mediated CHH methylation preferentially targeting Gypsy elements, particularly in their LTR regions. Old Gypsy elements possess unusually abundant siRNAs which show cross-mapping to young copies. Gypsy elements appear to be closer to genes and under stronger purifying selection than other classes. Our results suggest a continuous host-TE battle masked by the TE load reduction in R. apiculata. This conflict may enable mangroves, such as R. apiculata, to maintain genetic diversity and thus evolutionary potential during stress adaptation.© 2018 The Authors. New Phytologist © 2018 New Phytologist Trust.

Read more
September 22, 2019

Comparative genomics of Staphylococcus reveals determinants of speciation and diversification of antimicrobial defense.

The bacterial genus Staphylococcus comprises diverse species with most being described as colonizers of human and animal skin. A relational analysis of features that discriminate its species and contribute to niche adaptation and survival remains to be fully described. In this study, an interspecies, whole-genome comparative analysis of 21 Staphylococcus species was performed based on their orthologues. Three well-defined multi-species groups were identified: group A (including aureus/epidermidis); group B (including saprophyticus/xylosus) and group C (including pseudintermedius/delphini). The machine learning algorithm Random Forest was applied to prioritize orthologs that drive formation of the Staphylococcus species groups A-C. Orthologues driving staphylococcal intrageneric diversity comprised regulatory, metabolic and antimicrobial resistance proteins. Notably, the BraSR (NsaRS) two-component system (TCS) and its associated BraDE transporters that regulate antimicrobial resistance showed limited distribution in the genus and their presence was most closely associated with a subset of Staphylococcus species dominated by those that colonize human skin. Divergence of BraSR and GraSR antimicrobial peptide survival TCS and their associated transporters was observed across the staphylococci, likely reflecting niche specific evolution of these TCS/transporters and their specificities for AMPs. Experimental evolution, with selection for resistance to the lantibiotic nisin, revealed multiple routes to resistance and differences in the selection outcomes of the BraSR-positive species S. hominis and S. aureus. Selection supported a role for GraSR in nisin survival responses of the BraSR-negative species S. saprophyticus. Our study reveals diversification of antimicrobial-sensing TCS across the staphylococci and hints at differential relationships between GraSR and BraSR in those species positive for both TCS.

Read more
September 22, 2019

Whole-Genome Analysis of an Extensively Drug-Resistant Acinetobacter baumannii Strain XDR-BJ83: Insights into the Mechanisms of Resistance of an ST368 Strain from a Tertiary Care Hospital in China.

Acinetobacter baumannii is an important pathogen of nosocomial infections. Nosocomial outbreaks caused by antibiotic-resistant A. baumannii remain a significant challenge. Understanding the antibiotic resistance mechanism of A. baumannii is critical for clinical treatment. The purpose of this study was to determine the whole-genome sequence (WGS) of an extensively drug-resistant (XDR) A. baumannii strain, XDR-BJ83, which was associated with a nosocomial outbreak in a tertiary care hospital of China, and to investigate the antibiotic resistance mechanism of this strain. The WGS of XDR-BJ83 was performed using single-molecule real-time sequencing. The complete genome of XDR-BJ83 consisted of a 4,011,552-bp chromosome and a 69,069-bp plasmid. The sequence type of XDR-BJ83 was ST368, which belongs to clonal complex 92 (CC92). The chromosome of XDR-BJ83 carried multiple antibiotic resistance genes, antibiotic efflux pump genes, and mobile genetic elements, including insertion sequences, transposons, integrons, and resistance islands. The plasmid of XDR-BJ83 (pBJ83) was a conjugative plasmid carrying type IV secretion system. These results indicate that the presence of multiple antibiotic resistance genes, efflux pumps, and mobile genetic elements is likely associated with resistance to various antibiotics in XDR-BJ83.

Read more
September 22, 2019

Combining probabilistic alignments with read pair information improves accuracy of split-alignments.

Split-alignments provide base-pair-resolution evidence of genomic rearrangements. In practice, they are found by first computing high-scoring local alignments, parts of which are then combined into a split-alignment. This approach is challenging when aligning a short read to a large and repetitive reference, as it tends to produce many spurious local alignments leading to ambiguities in identifying the correct split-alignment. This problem is further exacerbated by the fact that rearrangements tend to occur in repeat-rich regions.We propose a split-alignment technique that combats the issue of ambiguous alignments by combining information from probabilistic alignment with positional information from paired-end reads. We demonstrate that our method finds accurate split-alignments, and that this translates into improved performance of variant-calling tools that rely on split-alignments.An open-source implementation is freely available at: https://bitbucket.org/splitpairedend/last-split-pe.Supplementary data are available at Bioinformatics online.

Read more
September 22, 2019

Thermosipho spp. immune system differences affect variation in genome size and geographical distributions.

Thermosipho species inhabit thermal environments such as marine hydrothermal vents, petroleum reservoirs, and terrestrial hot springs. A 16S rRNA phylogeny of available Thermosipho spp. sequences suggested habitat specialists adapted to living in hydrothermal vents only, and habitat generalists inhabiting oil reservoirs, hydrothermal vents, and hotsprings. Comparative genomics of 15 Thermosipho genomes separated them into three distinct species with different habitat distributions: The widely distributed T. africanus and the more specialized, T. melanesiensis and T. affectus. Moreover, the species can be differentiated on the basis of genome size (GS), genome content, and immune system composition. For instance, the T. africanus genomes are largest and contained the most carbohydrate metabolism genes, which could explain why these isolates were obtained from ecologically more divergent habitats. Nonetheless, all the Thermosipho genomes, like other Thermotogae genomes, show evidence of genome streamlining. GS differences between the species could further be correlated to differences in defense capacities against foreign DNA, which influence recombination via HGT. The smallest genomes are found in T. affectus that contain both CRISPR-cas Type I and III systems, but no RM system genes. We suggest that this has caused these genomes to be almost devoid of mobile elements, contrasting the two other species genomes that contain a higher abundance of mobile elements combined with different immune system configurations. Taken together, the comparative genomic analyses of Thermosipho spp. revealed genetic variation allowing habitat differentiation within the genus as well as differentiation with respect to invading mobile DNA.

Read more
September 22, 2019

Prevalence, antimicrobial resistance and phylogenetic characterization of Yersinia enterocolitica in retail poultry meat and swine feces in parts of China

Yersinia enterocolitica is an enteropathogen transmitted by contaminated food. In this study, a total of 500 retail poultry meat samples from 4 provinces and 145 swine feces samples from 12 provinces in China was tested for Y. enterocolitica and 26 isolates were obtained for further bio-serotyping, testing with antimicrobial susceptibility testing to a panel of antimicrobial compounds, and genetically characterization based on the whole genome sequencing. Higher prevalence (4.8%) of Y. enterocolitica contamination in retail poultry meat than that in swine feces (2.76%) was observed. No difference in bio-serotypes, multilocus sequence typing (MLST) and virulence genes distribution between swine and poultry origin were found. All isolates were resistant to ampicillin, amoxicillin/clavulanic acid, and cefazolin and were multi-drug resistant (MDR). The most predominant drug-resistance profile was AMP-CFZ-AMC-FOX (42.31%). A pathogenic isolate with bio-serotype 3/O:3 and ST135 was cultured from retail fresh chicken meat for the first time in China. Based on the whole-genome single nucleotide polymorphisms (SNPs) tree analysis, pathogenic isolates clustered closely, while nonpathogenic isolates exhibited high genetic heterogeneity. These indicated that pathogenic isolates were conserved on genetic level. The whole-genome SNP tree also revealed that Y. enterocolitica of swine, chicken and duck origin may share a common ancestor. The findings highlight the emergence of drug-resistant pathogenic Y. entrocoliticas in retailed poultry meats in China.

Read more
September 22, 2019

Isolation, characterization, genomic sequencing, and GFP-marked insertional mutagenesis of a high-performance nitrogen-fixing bacterium, Kosakonia radicincitans GXGL-4A and visualization of bacterial colonization on cucumber roots.

A gram-negative bacterium GXGL-4A was originally isolated from maize roots. It displayed nitrogen-fixing (NF) ability under nitrogen-free culture condition, and had a significant promotion effect on cucumber growth in the pot inoculation test. The preliminary physiological and biochemical traits of GXGL-4A were characterized. Furthermore, a phylogenetic tree was constructed based on 16S ribosomal DNA (rDNA) sequences of genetically related species. To determine the taxonomic status of GXGL-4A and further utilize its nitrogen-fixing potential, genome sequence was obtained using PacBio RS II technology. The analyses of average nucleotide identity based on BLAST+ (ANIb) and correlation indexes of tetra-nucleotide signatures (Tetra) showed that the NF isolate GXGL-4A is closely related to the Kosakonia radicincitans type strain DSM 16656. Therefore, the isolate GXGL-4A was eventually classified into the species of Kosakonia radicincitans and designated K. radicincitans GXGL-4A. A high consistency in composition and gene arrangement of nitrogen-fixing gene cluster I (nif cluster I) was found between K. radicincitans GXGL-4A and other Kosakonia NF strains. The mutants tagged with green fluorescence protein (GFP) were obtained by transposon Tn5 mutagenesis, and then, the colonization of gfp-marked K. radicincitans GXGL-4A cells on cucumber seedling root were observed under fluorescence microscopy. The preferential sites of the labeled GXGL-4A cell population were the lateral root junctions, the differentiation zone, and the elongation zone. All these results should benefit for the deep exploration of nitrogen fixation mechanism of K. radicincitans GXGL-4A and will definitely facilitate the genetic modification process of this NF bacterium in sustainable agriculture.

Read more
September 22, 2019

The unique evolution of the pig LRC, a single KIR but expansion of LILR and a novel Ig receptor family.

The leukocyte receptor complex (LRC) encodes numerous immunoglobulin (Ig)-like receptors involved in innate immunity. These include the killer-cell Ig-like receptors (KIR) and the leukocyte Ig-like receptors (LILR) which can be polymorphic and vary greatly in number between species. Using the recent long-read genome assembly, Sscrofa11.1, we have characterized the porcine LRC on chromosome 6. We identified a ~?197-kb region containing numerous LILR genes that were missing in previous assemblies. Out of 17 such LILR genes and fragments, six encode functional proteins, of which three are inhibitory and three are activating, while the majority of pseudogenes had the potential to encode activating receptors. Elsewhere in the LRC, between FCAR and GP6, we identified a novel gene that encodes two Ig-like domains and a long inhibitory intracellular tail. Comparison with two other porcine assemblies revealed a second, nearly identical, non-functional gene encoding a short intracellular tail with ambiguous function. These novel genes were found in a diverse range of mammalian species, including a pseudogene in humans, and typically consist of a single long-tailed receptor and a variable number of short-tailed receptors. Using porcine transcriptome data, both the novel inhibitory gene and the LILR were highly expressed in peripheral blood, while the single KIR gene, KIR2DL1, was either very poorly expressed or not at all. These observations are a prerequisite for improved understanding of immune cell functions in the pig and other species.

Read more
September 22, 2019

How complete are “complete” genome assemblies?-An avian perspective.

The genomics revolution has led to the sequencing of a large variety of nonmodel organisms often referred to as “whole” or “complete” genome assemblies. But how complete are these, really? Here, we use birds as an example for nonmodel vertebrates and find that, although suitable in principle for genomic studies, the current standard of short-read assemblies misses a significant proportion of the expected genome size (7% to 42%; mean 20 ± 9%). In particular, regions with strongly deviating nucleotide composition (e.g., guanine-cytosine-[GC]-rich) and regions highly enriched in repetitive DNA (e.g., transposable elements and satellite DNA) are usually underrepresented in assemblies. However, long-read sequencing technologies successfully characterize many of these underrepresented GC-rich or repeat-rich regions in several bird genomes. For instance, only ~2% of the expected total base pairs are missing in the last chicken reference (galGal5). These assemblies still contain thousands of gaps (i.e., fragmented sequences) because some chromosomal structures (e.g., centromeres) likely contain arrays of repetitive DNA that are too long to bridge with currently available technologies. We discuss how to minimize the number of assembly gaps by combining the latest available technologies with complementary strengths. At last, we emphasize the importance of knowing the location, size and potential content of assembly gaps when making population genetic inferences about adjacent genomic regions.© 2018 The Authors. Molecular Ecology Resources Published by John Wiley & Sons Ltd.

Read more
September 22, 2019

Computational tools to unmask transposable elements.

A substantial proportion of the genome of many species is derived from transposable elements (TEs). Moreover, through various self-copying mechanisms, TEs continue to proliferate in the genomes of most species. TEs have contributed numerous regulatory, transcript and protein innovations and have also been linked to disease. However, notwithstanding their demonstrated impact, many genomic studies still exclude them because their repetitive nature results in various analytical complexities. Fortunately, a growing array of methods and software tools are being developed to cater for them. This Review presents a summary of computational resources for TEs and highlights some of the challenges and remaining gaps to perform comprehensive genomic analyses that do not simply ‘mask’ repeats.

Read more
September 22, 2019

Genomic characterization reveals significant divergence within Chlorella sorokiniana (Chlorellales, Trebouxiophyceae)

Selection of highly productive algal strains is crucial for establishing economically viable biomass and biopro- duct cultivation systems. Characterization of algal genomes, including understanding strain-specific differences in genome content and architecture is a critical step in this process. Using genomic analyses, we demonstrate significant differences between three strains of Chlorella sorokiniana (strain 1228, UTEX 1230, and DOE1412). We found that unique, strain-specific genes comprise a substantial proportion of each genome, and genomic regions with> 80% local nucleotide identity constitute <15% of each genome among the strains, indicating substantial strain specific evolution. Furthermore, cataloging of meiosis and other sex-related genes in C. sor- okiniana strains suggests strategic breeding could be utilized to improve biomass and bioproduct yields if a sexual cycle can be characterized. Finally, preliminary investigation of epigenetic machinery suggests the pre- sence of potentially unique transcriptional regulation in each strain. Our data demonstrate that these three C. sorokiniana strains represent significantly different genomic content. Based on these findings, we propose in- dividualized assessment of each strain for potential performance in cultivation systems.

Read more
September 22, 2019

TranSurVeyor: an improved database-free algorithm for finding non-reference transpositions in high-throughput sequencing data.

Transpositions transfer DNA segments between different loci within a genome; in particular, when a transposition is found in a sample but not in a reference genome, it is called a non-reference transposition. They are important structural variations that have clinical impact. Transpositions can be called by analyzing second generation high-throughput sequencing datasets. Current methods follow either a database-based or a database-free approach. Database-based methods require a database of transposable elements. Some of them have good specificity; however this approach cannot detect novel transpositions, and it requires a good database of transposable elements, which is not yet available for many species. Database-free methods perform de novo calling of transpositions, but their accuracy is low. We observe that this is due to the misalignment of the reads; since reads are short and the human genome has many repeats, false alignments create false positive predictions while missing alignments reduce the true positive rate. This paper proposes new techniques to improve database-free non-reference transposition calling: first, we propose a realignment strategy called one-end remapping that corrects the alignments of reads in interspersed repeats; second, we propose a SNV-aware filter that removes some incorrectly aligned reads. By combining these two techniques and other techniques like clustering and positive-to-negative ratio filter, our proposed transposition caller TranSurVeyor shows at least 3.1-fold improvement in terms of F1-score over existing database-free methods. More importantly, even though TranSurVeyor does not use databases of prior information, its performance is at least as good as existing database-based methods such as MELT, Mobster and Retroseq. We also illustrate that TranSurVeyor can discover transpositions that are not known in the current database.

Read more
September 22, 2019

FRI-4 carbapenemase-producing Enterobacter cloacae complex isolated in Tokyo, Japan.

A carbapenem-resistant Enterobacter cloacae complex isolated in Tokyo, Japan, produced a carbapenemase that was detected by a Carba NP test and a modified carbapenem inactivation method, but none of the ‘Big Five’ carbapenemase genes was detected by PCR. This study aimed to identify the carbapenemase.Carbapenemase genes were screened by WGS. Next, we generated a recombinant plasmid in which the carbapenemase gene was inserted. We also extracted the carbapenemase gene-carrying plasmid from the E. cloacae complex. The effects of both plasmids on the antibiotic susceptibility of Escherichia coli were then tested. The carbapenemase gene-carrying plasmid in the E. cloacae complex was completely sequenced.A novel carbapenemase gene, blaFRI-4, encoded an amino acid sequence that was 93.2% identical to French imipenemase (FRI-1). E. coli transformed with blaFRI-4 showed reduced carbapenem susceptibility. A complete sequence of the blaFRI-4-carrying 98?508?bp IncFII/IncR plasmid (pTMTA61661) showed that blaFRI-4 and the surrounding region (18.7?kb) were duplicated.The FRI-4-producing E. cloacae complex was isolated in Japan, whereas all other FRI variants have been found in Europe, suggesting that the spread of FRI carbapenemases is global.

Read more
September 22, 2019

Plasmid and chromosomal integration of four novel blaIMP-carrying transposons from Pseudomonas aeruginosa, Klebsiella pneumoniae and an Enterobacter sp.

To provide detailed genetic characterization of four novel blaIMP-carrying transposons from Pseudomonas aeruginosa, Klebsiella pneumoniae and an Enterobacter sp.P. aeruginosa 60512, K. pneumoniae 447, P. aeruginosa 12939 and Enterobacter sp. A1137 were subjected to genome sequencing. The complete nucleotide sequences of two plasmids (p60512-IMP from the 60512 isolate and p447-IMP from the 447 isolate) and two chromosomes (the 12939 and A1137 isolates) were determined, then a genomic comparison of p60512-IMP, p447-IMP and four novel blaIMP-carrying transposons (Tn6394, Tn6375, Tn6411 and Tn6397) with related sequences was performed. Transferability of the blaIMP gene and bacterial antimicrobial susceptibility were tested.Tn6394 and Tn6375 were located in p60512-IMP and p447-IMP, respectively, while Tn6411 and Tn6397 were integrated into the 12939 and A1137 chromosomes, respectively. Tn6394 was an ISPa17-based transposition unit that harboured the integron In992 (carrying blaIMP-1). In73 (carrying blaIMP-8), In73 and In992, together with the ISEcp1:IS1R-blaCTX-M-14-IS903D unit, the macAB-tolC region and the truncated aacC2-tmrB region, respectively, were integrated into the prototype transposons Tn1722, Tn1696 and Tn7, respectively, generating the Tn3-family unit transposons, Tn6375 and Tn6378, and the Tn7-family unit transposon Tn6411, respectively. Tn6397 was a large integrative and conjugative element carrying Tn6378.Complex events of transposition and homologous recombination have occurred during the original formation and further plasmid and chromosomal integration of these four transposons, promoting accumulation and spread of antimicrobial resistance genes.

Read more
September 22, 2019

Diversity of DHA-1-encoding plasmids in Klebsiella pneumoniae isolates from 16 French hospitals.

To provide new insights into the spread of plasmidic cephalosporinase DHA-1, 16 strains of Klebsiella pneumoniae and a strain of Klebsiella variicola producing DHA-1 were isolated between January 2012 and December 2013 in six regions of France and two French overseas departments and territories.Disc diffusion assays, isoelectric focusing and PCRs were used to characterize the plasmidic DHA-1 ß-lactamase. Plasmid analysis was performed by the method of Kado and Liu and WGS. Virulence of the strains was studied by biofilm formation and the survival of Drosophila.The strains were of low virulence and had one to three plasmids including one of various sizes (~40 to 319?kb) mediating DHA-1. Nine strains belonged to ST11 and possessed a pKPS30-type DHA-1 plasmid of the IncR (incompatibility) group. A strain of ST307 possessed pENVA, a DHA-1 plasmid of the IncH-type group. The seven remaining plasmids were unknown. Three belonged to the IncL/M group. They were closely related and their sequences were determined. One of the four remaining strains was chosen for further investigation. This strain of ST16 had two plasmids, a pUUH239.2-related plasmid and a new DHA-1 plasmid of ~319?kb of IncHI2 type.These findings demonstrate the major role of the pKPS30-type plasmid in the spread of DHA-1 cephalosporinase in France and provide evidence of two new emerging plasmids carrying this enzyme.

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.