With the advent of high-throughput sequencing technologies, the rapid generation and accumulation of large amounts of sequencing data pose an insurmountable demand for efficient algorithms for constructing whole-genome phylogenies. The existing phylogenomic methods all use assembled sequences, which are often not available owing to the difficulty of assembling short-reads; this obstructs phylogenetic investigations on species without a reference genome. In this report, we present co-phylog, an assembly-free phylogenomic approach that creates a ‘micro-alignment’ at each ‘object’ in the sequence using the ‘context’ of the object and calculates pairwise distances before reconstructing the phylogenetic tree based on those distances. We explored the parameters’ usages and the optimal working range of co-phylog, assessed co-phylog using the simulated next-generation sequencing (NGS) data and the real NGS raw data. We also compared co-phylog method with traditional alignment and alignment-free methods and illustrated the advantages and limitations of co-phylog method. In conclusion, we demonstrated that co-phylog is efficient algorithm and that it delivers high resolution and accurate phylogenies using whole-genome unassembled sequencing data, especially in the case of closely related organisms, thereby significantly alleviating the computational burden in the genomic era.

No abstract

Chronic kidney disease (CKD) is a complex gene-environmental disease affecting close to 10% of the US population. Genome-wide association studies (GWASs) have identified sequence variants, localized to non-coding genomic regions, associated with kidney function. Despite these robust observations, the mechanism by which variants lead to CKD remains a critical unanswered question. Expression quantitative trait loci (eQTL) analysis is a method to identify genetic variation associated with gene expression changes in specific tissue types. We hypothesized that an integrative analysis combining CKD GWAS and kidney eQTL results can identify candidate genes for CKD. We performed eQTL analysis by correlating genotype with RNA-seq-based gene expression levels in 96 human kidney samples. Applying stringent statistical criteria, we detected 1,886 genes whose expression differs with the sequence variants. Using direct overlap and Bayesian methods, we identified new potential target genes for CKD. With respect to one of the target genes, lysosomal beta A mannosidase (MANBA), we observed that genetic variants associated with MANBA expression in the kidney showed statistically significant colocalization with variants identified in CKD GWASs, indicating that MANBA is a potential target gene for CKD. The expression of MANBA was significantly lower in kidneys of subjects with risk alleles. Suppressing manba expression in zebrafish resulted in renal tubule defects and pericardial edema, phenotypes typically induced by kidney dysfunction. Our analysis shows that gene-expression changes driven by genetic variation in the kidney can highlight potential new target genes for CKD development.

BackgroundThe development of multidrug resistance is a major problem in the treatment of pathogenic microorganisms by distinct antimicrobial agents. Characterizing the genetic variation among plasmids from different bacterial species or strains is a key step towards understanding the mechanism of virulence and their evolution.ResultsWe applied a deep sequencing approach to 206 clinical strains of Klebsiella pneumoniae collected from 2002 to 2008 to understand the genetic variation of multidrug resistance plasmids, and to reveal the dynamic change of drug resistance over time. First, we sequenced three plasmids (70 Kb, 94 Kb, and 147 Kb) from a clonal strain of K. pneumoniae using Sanger sequencing. Using the Illumina sequencing technology, we obtained more than 17 million of short reads from two pooled plasmid samples. We mapped these short reads to the three reference plasmid sequences, and identified a large number of single nucleotide polymorphisms (SNPs) in these pooled plasmids. Many of these SNPs are present in drug-resistance genes. We also found that a significant fraction of short reads could not be mapped to the reference sequences, indicating a high degree of genetic variation among the collection of K. pneumoniae isolates. Moreover, we identified that plasmid conjugative transfer genes and antibiotic resistance genes are more likely to suffer from positive selection, as indicated by the elevated rates of nonsynonymous substitution.ConclusionThese data represent the first large-scale study of genetic variation in multidrug resistance plasmids and provide insight into the mechanisms of plasmid diversification and the genetic basis of antibiotic resistance.

BackgroundFlorfenicol is a derivative of chloramphenicol that is used only for the treatment of animal diseases. A key resistance gene for florfenicol, floR, can spread among bacteria of the same and different species or genera through horizontal gene transfer. To analyze the potential transmission of resistance genes between animal and human pathogens, we investigated floR in Klebsiella pneumoniae isolates from patient samples. floR in human pathogens may originate from animal pathogens and would reflect the risk to human health of using antimicrobial agents in animals.MethodsPCR was used to identify floR-positive strains. The floR genes were cloned, and the minimum inhibitory concentrations (MICs) were determined to assess the relative resistance levels of the genes and strains. Sequencing and comparative genomics methods were used to analyze floR gene-related sequence structure as well as the molecular mechanism of resistance dissemination.ResultsOf the strains evaluated, 20.42% (67/328) were resistant to florfenicol, and 86.96% (20/23) of the floR-positive strains demonstrated high resistance to florfenicol with MICs ≥512 μg/mL. Conjugation experiments showed that transferrable plasmids carried the floR gene in three isolates. Sequencing analysis of a plasmid approximately 125 kb in size (pKP18–125) indicated that the floR gene was flanked by multiple copies of mobile genetic elements. Comparative genomics analysis of a 9-kb transposon-like fragment of pKP18–125 showed that an approximately 2-kb sequence encoding lysR-floR-virD2 was conserved in the majority (79.01%, 83/105) of floR sequences collected from NCBI nucleotide database. Interestingly, the most similar sequence was a 7-kb fragment of plasmid pEC012 from an Escherichia coli strain isolated from a chicken.ConclusionsIdentified on a transferable plasmid in the human pathogen K. pneumoniae, the floR gene may be disseminated through horizontal gene transfer from animal pathogens. Studies on the molecular mechanism of resistance gene dissemination in different bacterial species of animal origin could provide useful information for preventing or controlling the spread of resistance between animal and human pathogens.Electronic supplementary materialThe online version of this article (10.1186/s13756-018-0415-0) contains supplementary material, which is available to authorized users.