Extraintestinal pathogenic Escherichia coli (ExPEC) strains are important pathogens that cause diverse diseases in humans and poultry. Some E. coli isolates from chicken feces contain ExPEC-associated virulence genes, so appear potentially pathogenic; they conceivably could be transmitted to humans through handling and/or consumption of contaminated meat. However, the actual extraintestinal virulence potential of chicken-source fecal E. coli is poorly understood. Here, we assessed whether fecal E. coli isolates from healthy production chickens could cause diseases in a chicken model of avian colibacillosis and three rodent models of ExPEC-associated human infections. From 304 E. coli isolates from chicken fecal samples, 175 E. coli isolates were screened by PCR for virulence genes associated with human-source ExPEC or avian pathogenic E. coli (APEC), an ExPEC subset that causes extraintestinal infections in poultry. Selected isolates genetically identified as ExPEC and non-ExPEC isolates were assessed in vitro for virulence-associated phenotypes, and in vivo for disease-causing ability in animal models of colibacillosis, sepsis, meningitis, and urinary tract infection. Among the study isolates, 13% (40/304) were identified as ExPEC; the majority of these were classified as APEC and uropathogenic E. coli, but none as neonatal meningitis E. coli. Multiple chicken-source fecal ExPEC isolates resembled avian and human clinical ExPEC isolates in causing one or more ExPEC-associated illnesses in experimental animal infection models. Additionally, some isolates that were classified as non-ExPEC were able to cause ExPEC-associated illnesses in animal models, and thus future studies are needed to elucidate their mechanisms of virulence. These findings show that E. coli isolates from chicken feces contain ExPEC-associated genes, exhibit ExPEC-associated in vitro phenotypes, and can cause ExPEC-associated infections in animal models, and thus may pose a health threat to poultry and consumers.

Bismuth drugs, despite being clinically used for decades, surprisingly remain in use and effective for the treatment of Helicobacter pylori infection, even for resistant strains when co-administrated with antibiotics. However, the molecular mechanisms underlying the clinically sustained susceptibility of H. pylori to bismuth drugs remain elusive. Herein, we report that integration of in-house metalloproteomics and quantitative proteomics allows comprehensive uncovering of the bismuth-associated proteomes, including 63 bismuth-binding and 119 bismuth-regulated proteins from Helicobacter pylori, with over 60% being annotated with catalytic functions. Through bioinformatics analysis in combination with bioassays, we demonstrated that bismuth drugs disrupted multiple essential pathways in the pathogen, including ROS defence and pH buffering, by binding and functional perturbation of a number of key enzymes. Moreover, we discovered that HpDnaK may serve as a new target of bismuth drugs to inhibit bacterium-host cell adhesion. The integrative approach we report, herein, provides a novel strategy to unveil the molecular mechanisms of antimicrobial metals against pathogens in general. This study sheds light on the design of new types of antimicrobial agents with multiple targets to tackle the current crisis of antimicrobial resistance.

Competing endogenous RNAs (ceRNAs) are RNA molecules that sequester shared microRNAs (miRNAs) thereby affecting the expression of other targets of the miRNAs. Whether genetic variants in ceRNA can affect its biological function and disease development is still an open question. Here we identified a large number of genetic variants that are associated with ceRNA's function using Geuvaids RNA-seq data for 462 individuals from the 1000 Genomes Project. We call these loci competing endogenous RNA expression quantitative trait loci or ‘cerQTL’, and found that a large number of them were unexplored in conventional eQTL mapping. We identified many cerQTLs that have undergone recent positive selection in different human populations, and showed that single nucleotide polymorphisms in gene 3΄UTRs at the miRNA seed binding regions can simultaneously regulate gene expression changes in both cis and trans by the ceRNA mechanism. We also discovered that cerQTLs are significantly enriched in traits/diseases associated variants reported from genome-wide association studies in the miRNA binding sites, suggesting that disease susceptibilities could be attributed to ceRNA regulation. Further in vitro functional experiments demonstrated that a cerQTL rs11540855 can regulate ceRNA function. These results provide a comprehensive catalog of functional non-coding regulatory variants that may be responsible for ceRNA crosstalk at the post-transcriptional level.

Infection after renal transplantation remains a major cause of morbidity and death, especially infection from the extensively drug-resistant bacteria, A. baumannii. A total of fourteen A. baumannii isolates were isolated from the donors’ preserved fluid from DCD (donation after cardiac death) renal transplantation and four isolates in the recipients’ draining liquid at the Kidney Disease Center, The First Affiliated Hospital, College of Medicine, Zhejiang University, from March 2013 to November 2014. An outbreak of A. baumannii emerging after DCD renal transplantation was tracked to understand the transmission of the pathogen. PFGE displayed similar DNA patterns between isolates from the same hospital. Antimicrobial susceptibility tests against thirteen antimicrobial agents were determined using the K-B diffusion method and eTest. Whole-genome sequencing was applied to investigate the genetic relationship of the isolates. With the clinical data and research results, we concluded that the A. baumannii isolates 3R1 and 3R2 was probably transmitted from the donor who acquired the bacteria during his stay in the ICU, while isolate 4R1 was transmitted from 3R1 and 3R2 via medical manipulation. This study demonstrated the value of integration of clinical profiles with molecular methods in outbreak investigation and their importance in controlling infection and preventing serious complications after DCD transplantation.

Whole genome sequencing (WGS) is a promising strategy to unravel variants or genes responsible for human diseases and traits. However, there is a lack of robust platforms for a comprehensive downstream analysis. In the present study, we first proposed three novel algorithms, sequence gap-filled gene feature annotation, bit-block encoded genotypes and sectional fast access to text lines to address three fundamental problems. The three algorithms then formed the infrastructure of a robust parallel computing framework, KGGSeq, for integrating downstream analysis functions for whole genome sequencing data. KGGSeq has been equipped with a comprehensive set of analysis functions for quality control, filtration, annotation, pathogenic prediction and statistical tests. In the tests with whole genome sequencing data from 1000 Genomes Project, KGGSeq annotated several thousand more reliable non-synonymous variants than other widely used tools (e.g. ANNOVAR and SNPEff). It took only around half an hour on a small server with 10 CPUs to access genotypes of ∼60 million variants of 2504 subjects, while a popular alternative tool required around one day. KGGSeq's bit-block genotype format used 1.5% or less space to flexibly represent phased or unphased genotypes with multiple alleles and achieved a speed of over 1000 times faster to calculate genotypic correlation.

Modeling of transcriptional regulatory networks (TRNs) has been increasingly used to dissect the nature of gene regulation. Inference of regulatory relationships among transcription factors (TFs) and genes, especially among multiple TFs, is still challenging. In this study, we introduced an integrative method, LogicTRN, to decode TF–TF interactions that form TF logics in regulating target genes. By combining cis-regulatory logics and transcriptional kinetics into one single model framework, LogicTRN can naturally integrate dynamic gene expression data and TF-DNA-binding signals in order to identify the TF logics and to reconstruct the underlying TRNs. We evaluated the newly developed methodology using simulation, comparison and application studies, and the results not only show their consistence with existing knowledge, but also demonstrate its ability to accurately reconstruct TRNs in biological complex systems.

Salmonella enterica serovar Typhimurium genome encodes 13 fimbrial operons. Most of the fimbriae encoded by these operons are not produced under laboratory conditions but are likely to be synthesized in vivo. We used an in vivo expression technology (IVET) strategy to identify four fimbrial operons, agf, saf, sti, and stc that are expressed in the spleen. When any three of these operons were deleted, the strain retained wild-type virulence. However, when all four operons were deleted, the resulting strain was completely attenuated, indicating that these four fimbriae play functionally redundant roles critical for virulence. In mice, oral doses of as low as 1 × 10[superscript 5] CFU of the strain with four fimbrial operons deleted provided 100% protection against challenge with 1 × 10[superscript 9] CFU of wild-type S. Typhimurium. We also examined the possible effect of these fimbriae on the ability of a Salmonella vaccine strain to deliver a guest antigen. We modified one of our established attenuated vaccine strains, χ9088, to delete three fimbrial operons while the fourth operon was constitutively expressed. Each derivative was modified to express the Streptococcus pneumoniae antigen PspA. Strains that constitutively expressed saf or stc elicited a strong Th1 response with significantly greater levels of anti-PspA serum IgG and greater protective efficacy than strains carrying saf or stc deletions. The isogenic strain in which all four operons were deleted generated the lowest anti-PspA levels and did not protect against challenge with virulent S. pneumoniae. Our results indicate that these fimbriae play important roles, as yet not understood, in Salmonella virulence and immunogenicity.

Antibiotic-resistant bacterial infections are difficult to treat, producing a burden on healthcare and the economy. Extraintestinal pathogenic Escherichia coli (ExPEC) strains frequently carry antibiotic resistance genes, cause infections outside of the intestine, and are causative agents of hospital-acquired infections. Developing a prevention strategy against this pathogen is challenging due to its antibiotic resistance and antigenic diversity. E. coli common pilus (ECP) is frequently found in ExPEC strains and may serve as a common antigen to induce protection against several ExPEC serotypes. In addition, live recombinant attenuated Salmonella vaccine (RASV) strains have been used to prevent Salmonella infection and can also be modified to deliver foreign antigens. Thus, the objective of this study was to design a RASV to produce ECP on its surface and assess its ability to provide protection against ExPEC infections. To constitutively display ECP in a RASV strain, we genetically engineered a vector (pYA4428) containing aspartate-β-semialdehyde dehydrogenase and E. coli ecp genes and introduced it into RASV χ9558. RASV χ9558 containing an empty vector (pYA3337) was used as a control to assess protection conferred by the RASV strain without ECP.

We assessed vaccine efficacy in in vitro bacterial inhibition assays and mouse models of ExPEC-associated human infections. We found that RASV χ9558(pYA4428) synthesized the major pilin (EcpA) and tip pilus adhesin (EcpD) on the bacterial surface. Mice orally vaccinated with RASV χ9558(pYA3337) without ECP or χ9558(pYA4428) with ECP, produced anti-Salmonella LPS and anti-E. coli EcpA and EcpD IgG and IgA antibodies. RASV strains showed protective potential against some E. coli and Salmonella strains as assessed using in vitro assays. In mouse sepsis and urinary tract infection challenge models, both vaccines had significant protection in some internal organs. Overall, this work showed that RASVs can elicit an immune response to E. coli and Salmonella antigens in some mice, provide significant protection in some internal organs during ExPEC challenge, and thus this study is a promising initial step toward developing a vaccine for prevention of ExPEC infections. Future studies should optimize the ExPEC antigens displayed by the RASV strain for a more robust immune response and enhanced protection against ExPEC infection.

Linnorm is a novel normalization and transformation method for the analysis of single cell RNA sequencing (scRNA-seq) data. Linnorm is developed to remove technical noises and simultaneously preserve biological variations in scRNA-seq data, such that existing statistical methods can be improved. Using real scRNA-seq data, we compared Linnorm with existing normalization methods, including NODES, SAMstrt, SCnorm, scran, DESeq and TMM. Linnorm shows advantages in speed, technical noise removal and preservation of cell heterogeneity, which can improve existing methods in the discovery of novel subtypes, pseudo-temporal ordering of cells, clustering analysis, etc. Linnorm also performs better than existing DEG analysis methods, including BASiCS, NODES, SAMstrt, Seurat and DESeq2, in false positive rate control and accuracy.

It remains challenging to predict regulatory variants in particular tissues or cell types due to highly context-specific gene regulation. By connecting large-scale epigenomic profiles to expression quantitative trait loci (eQTLs) in a wide range of human tissues/cell types, we identify critical chromatin features that predict variant regulatory potential. We present cepip, a joint likelihood framework, for estimating a variant’s regulatory probability in a context-dependent manner. Our method exhibits significant GWAS signal enrichment and is superior to existing cell type-specific methods. Furthermore, using phenotypically relevant epigenomes to weight the GWAS single-nucleotide polymorphisms, we improve the statistical power of the gene-based association test.