ALBERT

Description: Clonal populations accumulate mutations over time, resulting in different haplotypes. Deep sequencing of such a population in principle provides information to reconstruct these haplotypes and the frequency at which the haplotypes occur. However, this reconstruction is technically not trivial, especially not in clonal systems with a relatively low mutation frequency. The low number of segregating sites in those systems adds ambiguity to the haplotype phasing and thus obviates the reconstruction of genome-wide haplotypes based on sequence overlap information. Therefore, we present EVORhA, a haplotype reconstruction method that complements phasing information in the non-empty read overlap with the frequency estimations of inferred local haplotypes. As was shown with simulated data, as soon as read lengths and/or mutation rates become restrictive for state-of-the-art methods, the use of this additional frequency information allows EVORhA to still reliably reconstruct genome-wide haplotypes. On real data, we show the applicability of the method in reconstructing the population composition of evolved bacterial populations and in decomposing mixed bacterial infections from clinical samples.

Description: Data on biological mechanisms of aging are mostly obtained from cross-sectional study designs. An inherent disadvantage of this design is that inter-individual differences can mask small but biologically significant age-dependent changes. A serially sampled design (same individual at different time points) would overcome this problem but is often limited by the relatively small numbers of available paired samples and the statistics being used. To overcome these limitations, we have developed a new vector-based approach, termed three-component analysis, which incorporates temporal distance, signal intensity and variance into one single score for gene ranking and is combined with gene set enrichment analysis. We tested our method on a unique age-based sample set of human skin fibroblasts and combined genome-wide transcription, DNA methylation and histone methylation (H3K4me3 and H3K27me3) data. Importantly, our method can now for the first time demonstrate a clear age-dependent decrease in expression of genes coding for proteins involved in translation and ribosome function. Using analogies with data from lower organisms, we propose a model where age-dependent down-regulation of protein translation-related components contributes to extend human lifespan.

Description: A major challenge in the field of shotgun metagenomics is the accurate identification of organisms present within a microbial community, based on classification of short sequence reads. Though existing microbial community profiling methods have attempted to rapidly classify the millions of reads output from modern sequencers, the combination of incomplete databases, similarity among otherwise divergent genomes, errors and biases in sequencing technologies, and the large volumes of sequencing data required for metagenome sequencing has led to unacceptably high false discovery rates (FDR). Here, we present the application of a novel, gene-independent and signature-based metagenomic taxonomic profiling method with significantly and consistently smaller FDR than any other available method. Our algorithm circumvents false positives using a series of non-redundant signature databases and examines G enomic O rigins T hrough T axonomic CHA llenge (GOTTCHA). GOTTCHA was tested and validated on 20 synthetic and mock datasets ranging in community composition and complexity, was applied successfully to data generated from spiked environmental and clinical samples, and robustly demonstrates superior performance compared with other available tools.

Description: Invertebrate animal species that can withstand temperatures as high as 37°C, the human body temperature, are limited. In the present study, we utilized the two-spotted cricket, Gryllus bimaculatus , which lives in tropical and subtropical regions, as an animal model of human pathogenic bacterial infection. Injection of Pseudomonas aeruginosa or Staphylococcus aureus into the hemolymph killed crickets. Injected P. aeruginosa or S. aureus proliferated in the hemolymph until the cricket died. The ability of these pathogenic bacteria to kill the crickets was blocked by the administration of antibiotics. S. aureus gene-knockout mutants of virulence factors, including cvfA, agr and srtA , exhibited decreased killing ability compared with the parent strain. The dose at which 50% of crickets were killed by P. aeruginosa or S. aureus was not decreased at 37°C compared with that at 27°C. Injection of Listeria monocytogenes , which upregulates toxin expression at 37°C, killed crickets, and the dose at which 50% of crickets were killed was decreased at 37°C compared with that at 27°C. These findings suggest that the two-spotted cricket is a useful model animal for evaluating the virulence properties of various human pathogenic bacteria at variable temperature including 37°C.

Description: The OmpA-like protein domain has been associated with peptidoglycan-binding proteins, and is often found in virulence factors of bacterial pathogens. The intracellular pathogen Legionella pneumophila encodes for six proteins that contain the OmpA-like domain, among them the highly conserved uncharacterized protein we named CmpA. Here we set out to characterize the CmpA protein and determine its contribution to intracellular survival of L. pneumophila . Secondary structure analysis suggests that CmpA is an inner membrane protein with a peptidoglycan-binding domain at the C-teminus. A cmpA mutant was able to replicate normally in broth, but failed to compete with an isogenic wild-type strain in an intracellular growth competition assay. The cmpA mutant also displayed significant intracellular growth defects in both the protozoan host Acanthamoeba castellanii and in primary bone marrow-derived macrophages, where uptake into the cells was also impaired. The cmpA phenotypes were completely restored upon expression of CmpA in trans . The data presented here establish CmpA as a novel virulence factor of L. pneumophila that is required for efficient intracellular replication in both mammalian and protozoan hosts.

Description: Shiga toxin-encoding Escherichia coli (STEC) regroup strains that carry genes encoding Shiga toxin (Stx). Among intestinal pathogenic E. coli , enterohaemorrhagic E. coli (EHEC) constitute the major subgroup of virulent STEC. EHEC cause serious human disease such as haemorrhagic colitis and haemolytic-uremic syndrome. While EHEC have evolved from enteropathogenic E. coli , hybrids with enteroaggregative E. coli have recently emerged. Of note, some enteroinvasive E. coli also belong to the STEC group. While the LEE (locus of enterocyte effacement) is a key and prominent molecular determinant in the pathogenicity, neither all EHEC nor STEC contain the LEE, suggesting that they possess additional virulence and colonisation factors. Currently, nine protein secretion systems have been described in diderm-lipopolysaccharide bacteria (archetypal Gram-negative) and can be involved in the secretion of extracellular effectors, cell-surface proteins or assembly of cell-surface organelles, such as flagella or pili. In this review, we focus on the secretome of STEC and related enteropathotypes, which are relevant to the colonisation of biotic and abiotic surfaces. Considering the wealth of potential protein trafficking mechanisms, the different combinations of colonisation factors and modulation of their expression is further emphasised with regard to the ecophysiology of STEC.

Description: The activity of levofloxacin against planktonic and biofilm Stenotrophomonas maltophilia cells and the role played by the multidrug efflux pump SmeDEF were evaluated under conditions relevant to the cystic fibrosis (CF) lung. MIC, MBC and MBEC of levofloxacin were assessed, against five CF strains, under ‘standard’ (CLSI-recommended) and ‘CF-like’ (pH 6.8, 5% CO 2 , in a synthetic CF sputum) conditions. Levofloxacin was tested against biofilms at concentrations (10, 50 and 100 μg mL –1 ) corresponding to achievable serum levels and sputum levels by aerosolisation. smeD expression was evaluated, under both conditions, in planktonic and biofilm cells by RT-PCR. The bactericidal effect of levofloxacin was decreased, in three out of five strains tested, under ‘CF-like’ conditions (MBC: 2–4 vs 8–16 μg mL –1 , under ‘standard’ and ‘CF-like’ conditions, respectively). Biofilm was intrinsically resistant to levofloxacin, regardless of conditions tested (MBECs ≥ 100 μg mL –1 for all strains). Only under ‘CF-like’ conditions, smeD expression increased during planktonic-to-biofilm transition, and in biofilm cells compared to stationary planktonic cells. Our findings confirmed that S. maltophilia biofilm is intrinsically resistant to therapeutic concentrations of levofloxacin. Under conditions relevant to CF, smeD overexpression could contribute to levofloxacin resistance. Further studies are warranted to define the clinical relevance of our findings .

Description: Auranofin is an FDA-approved gold-containing compound used for the treatment of rheumatoid arthritis. Recent reports of antimicrobial activity against protozoa and bacteria indicate that auranofin targets the reductive enzyme thioredoxin reductase (TrxR). We evaluated auranofin as well as five auranofin analogs containing N- heterocyclic carbenes (instead of the triethylphosphane present in auranofin) and five gold-carbene controls for their ability to inhibit or kill Helicobacter pylori in vitro . Auranofin completely inhibited bacterial growth at 1.2 μM. Purified H. pylori TrxR was inhibited by auranofin in a cell-free assay (IC 50 ~88 nM). The most active gold(I)- N- heterocyclic carbene compounds exhibited MICs comparable to auranofin against H. pylori (2 μM), while also exhibiting lower toxicities for human embryonic kidney cells (HEK-293T cells). Median toxic concentrations (TC 50 ) were 13–20-fold higher compared to auranofin indicating that they were less cytotoxic. The N- heterocyclic carbene analogs maybe well tolerated, but further evaluation is needed in vivo . Finally, auranofin was synergistic with the antibiotic amoxicillin, suggesting that targeting both the reductive enzyme TrxR and cell wall synthesis may be effective against H. pylori infections.

Description: Legionella feeleii is a Gram-negative pathogenic bacterium that causes Pontiac fever and pneumonia in humans. When L. feeleii serogroup 1 (ATCC 35072) was cultured on BCYE agar plates, two types of colonies were observed and exhibited differences in color, opacity and morphology. Since the two colony types are white rugose and brown translucent, they were termed as white rugose L. feeleii (WRLf) and brown translucent L. feeleii (BTLf), respectively. They exhibited different growth capacities in BYE broth in vitro , and it was found that WRLf could transform to BTLf. Under the electron microscope, it was observed that WRLf secreted materials which could be stained with ruthenium red, which was absent in BTLf. When U937 macrophages and HeLa cells were infected with the bacteria, WRLf manifested stronger internalization ability than BTLf. Intracellular growth in murine macrophages and Acanthamoeba cells was affected by the level of initial phagocytosis. WRLf was more resistant to human serum bactericidal action than BTLf. After being inoculated to guinea pigs, both organisms caused fever in the animals. These results suggest that ruthenium red-stained materials secreted in the surroundings may play a crucial role in determining L. feeleii colony morphology and virulence traits.

Description: Piscirickettsia salmonis is a fastidious intracellular pathogen responsible for high mortality rates in farmed salmonids, with serious economic consequences for the Chilean aquaculture industry. Oxytetracycline and florfenicol are the most frequently used antibiotics against P. salmonis , but routine use could contribute to drug resistance. This study identified differentiated florfenicol susceptibilities in two P. salmonis strains, LF-89 and AUSTRAL-005. The less susceptible isolate, AUSTRAL-005, also showed a high ethidium bromide efflux rate, indicating a higher activity of general efflux pump genes than LF-89. The P. salmonis genome presented resistance nodulation division (RND) family members, a family containing typical multidrug resistance-related efflux pumps in Gram-negative bacteria. Additionally, efflux pump acrAB genes were overexpressed in AUSTRAL-005 following exposure to the tolerated maximal concentration of florfenicol, in contrast to LF-89. These results indicate that tolerated maximum concentrations of florfenicol can modulate RND gene expression and increase efflux pump activity. We propose that the acrAB efflux pump is essential for P. salmonis survival at critical florfenicol concentrations and for the generation of antibiotic-resistant bacterial strains.

Description: Pseudomonas aeruginosa is an opportunistic pathogen, known to develop robust biofilms. Its biofilm development increases when antibiotics are presented at subminimal inhibitory concentrations (MICs) for reasons that remain unclear. In order to identify genes that affect biofilm development under such a sublethal antibiotic stress condition, we screened a transposon (Tn) mutant library of PAO1, a prototype P. aeruginosa strain. Among ~5000 mutants, a fiuA gene mutant was verified to form very defective biofilms in the presence of sub-MIC carbenicillin. The fiuA gene encodes ferrichrome receptor A, involved in the iron acquisition process. Of note, biofilm formation was not decreased in the pchpvd mutant defective in the production of pyochelin and pyoverdine, two well-characterized P. aeruginosa siderophore molecules. Moreover, fiuA , a non-polar fiuA deletion mutant, produced a significantly decreased level of elastase, a major virulence determinant. Mouse airway infection experiments revealed that the mutant expressed significantly less pathogenicity. Our results suggest that the fiuA gene has pleiotropic functions that affect P. aeruginosa biofilm development and virulence. The targeting of FiuA could enable the attenuation of P. aeruginosa virulence and may be suitable for the development of a drug that specifically controls the virulence of this important pathogen.

Description: Chromosome-long haplotyping of human genomes is important to identify genetic variants with differing gene expression, in human evolution studies, clinical diagnosis, and other biological and medical fields. Although several methods have realized haplotyping based on sequencing technologies or population statistics, accuracy and cost are factors that prohibit their wide use. Borrowing ideas from group testing theories, we proposed a clone-based haplotyping method by overlapping pool sequencing. The clones from a single individual were pooled combinatorially and then sequenced. According to the distinct pooling pattern for each clone in the overlapping pool sequencing, alleles for the recovered variants could be assigned to their original clones precisely. Subsequently, the clone sequences could be reconstructed by linking these alleles accordingly and assembling them into haplotypes with high accuracy. To verify the utility of our method, we constructed 130 110 clones in silico for the individual NA12878 and simulated the pooling and sequencing process. Ultimately, 99.9% of variants on chromosome 1 that were covered by clones from both parental chromosomes were recovered correctly, and 112 haplotype contigs were assembled with an N50 length of 3.4 Mb and no switch errors. A comparison with current clone-based haplotyping methods indicated our method was more accurate.

Description: Sexual differentiation of malaria parasites into gametocytes in the vertebrate host and subsequent gamete fertilization in mosquitoes is essential for the spreading of the disease. The molecular processes orchestrating these transitions are far from fully understood. Here, we report the first transcriptome analysis of male and female Plasmodium falciparum gametocytes coupled with a comprehensive proteome analysis. In male gametocytes there is an enrichment of proteins involved in the formation of flagellated gametes; proteins involved in DNA replication, chromatin organization and axoneme formation. On the other hand, female gametocytes are enriched in proteins required for zygote formation and functions after fertilization; protein-, lipid- and energy-metabolism. Integration of transcriptome and proteome data revealed 512 highly expressed maternal transcripts without corresponding protein expression indicating large scale translational repression in P. falciparum female gametocytes for the first time. Despite a high degree of conservation between Plasmodium species, 260 of these ‘repressed transcripts’ have not been previously described. Moreover, for some of these genes, protein expression is only reported in oocysts and sporozoites indicating that repressed transcripts can be partitioned into short- and long-term storage. Finally, these data sets provide an essential resource for identification of vaccine/drug targets and for further mechanistic studies.

Description: DNA methylation is a mechanism for long-term transcriptional regulation and is required for normal cellular differentiation. Failure to properly establish or maintain DNA methylation patterns leads to cell dysfunction and diseases such as cancer. Identifying DNA methylation signatures in complex tissues can be challenging owing to inaccurate cell enrichment methods and low DNA yields. We have developed a technique called laser capture microdissection-reduced representation bisulfite sequencing (LCM-RRBS) for the multiplexed interrogation of the DNA methylation status of cytosine–guanine dinucleotide islands and promoters. LCM-RRBS accurately and reproducibly profiles genome-wide methylation of DNA extracted from microdissected fresh frozen or formalin-fixed paraffin-embedded tissue samples. To demonstrate the utility of LCM-RRBS, we characterized changes in DNA methylation associated with gonadectomy-induced adrenocortical neoplasia in the mouse. Compared with adjacent normal tissue, the adrenocortical tumors showed reproducible gains and losses of DNA methylation at genes involved in cell differentiation and organ development. LCM-RRBS is a rapid, cost-effective, and sensitive technique for analyzing DNA methylation in heterogeneous tissues and will facilitate the investigation of DNA methylation in cancer and organ development.

Description: The advent in high-throughput-sequencing (HTS) technologies has revolutionized conventional biodiversity research by enabling parallel capture of DNA sequences possessing species-level diagnosis. However, polymerase chain reaction (PCR)-based implementation is biased by the efficiency of primer binding across lineages of organisms. A PCR-free HTS approach will alleviate this artefact and significantly improve upon the multi-locus method utilizing full mitogenomes. Here we developed a novel multiplex sequencing and assembly pipeline allowing for simultaneous acquisition of full mitogenomes from pooled animals without DNA enrichment or amplification. By concatenating assemblies from three de novo assemblers, we obtained high-quality mitogenomes for all 49 pooled taxa, with 36 species 〉15 kb and the remaining 〉10 kb, including 20 complete mitogenomes and nearly all protein coding genes (99.6%). The assembly quality was carefully validated with Sanger sequences, reference genomes and conservativeness of protein coding genes across taxa. The new method was effective even for closely related taxa, e.g. three Drosophila spp., demonstrating its broad utility for biodiversity research and mito-phylogenomics. Finally, the in silico simulation showed that by recruiting multiple mito-loci, taxon detection was improved at a fixed sequencing depth. Combined, these results demonstrate the plausibility of a multi-locus mito-metagenomics approach as the next phase of the current single-locus metabarcoding method.

Description: Chromatin modifiers and histone modifications are components of a chromatin-signaling network involved in transcription and its regulation. The interactions between chromatin modifiers and histone modifications are often unknown, are based on the analysis of few genes or are studied in vitro . Here, we apply computational methods to recover interactions between chromatin modifiers and histone modifications from genome-wide ChIP-Seq data. These interactions provide a high-confidence backbone of the chromatin-signaling network. Many recovered interactions have literature support; others provide hypotheses about yet unknown interactions. We experimentally verified two of these predicted interactions, leading to a link between H4K20me1 and members of the Polycomb Repressive Complexes 1 and 2. Our results suggest that our computationally derived interactions are likely to lead to novel biological insights required to establish the connectivity of the chromatin-signaling network involved in transcription and its regulation.

Description: With read lengths of currently up to 2 x 300 bp, high throughput and low sequencing costs Illumina's MiSeq is becoming one of the most utilized sequencing platforms worldwide. The platform is manageable and affordable even for smaller labs. This enables quick turnaround on a broad range of applications such as targeted gene sequencing, metagenomics, small genome sequencing and clinical molecular diagnostics. However, Illumina error profiles are still poorly understood and programs are therefore not designed for the idiosyncrasies of Illumina data. A better knowledge of the error patterns is essential for sequence analysis and vital if we are to draw valid conclusions. Studying true genetic variation in a population sample is fundamental for understanding diseases, evolution and origin. We conducted a large study on the error patterns for the MiSeq based on 16S rRNA amplicon sequencing data. We tested state-of-the-art library preparation methods for amplicon sequencing and showed that the library preparation method and the choice of primers are the most significant sources of bias and cause distinct error patterns. Furthermore we tested the efficiency of various error correction strategies and identified quality trimming (Sickle) combined with error correction (BayesHammer) followed by read overlapping (PANDAseq) as the most successful approach, reducing substitution error rates on average by 93%.

Description: RNA-seq is a sensitive and accurate technique to compare steady-state levels of RNA between different cellular states. However, as it does not provide an account of transcriptional activity per se , other technologies are needed to more precisely determine acute transcriptional responses. Here, we have developed an easy, sensitive and accurate novel computational method, iRNA-seq , for genome-wide assessment of transcriptional activity based on analysis of intron coverage from total RNA-seq data. Comparison of the results derived from iRNA-seq analyses with parallel results derived using current methods for genome-wide determination of transcriptional activity, i.e. global run-on (GRO)-seq and RNA polymerase II (RNAPII) ChIP-seq, demonstrate that iRNA-seq provides similar results in terms of number of regulated genes and their fold change. However, unlike the current methods that are all very labor-intensive and demanding in terms of sample material and technologies, iRNA-seq is cheap and easy and requires very little sample material. In conclusion, iRNA-seq offers an attractive novel alternative to current methods for determination of changes in transcriptional activity at a genome-wide level.

Description: Identifying large-scale structural variation in cancer genomes continues to be a challenge to researchers. Current methods rely on genome alignments based on a reference that can be a poor fit to highly variant and complex tumor genomes. To address this challenge we developed a method that uses available breakpoint information to generate models of structural variations. We use these models as references to align previously unmapped and discordant reads from a genome. By using these models to align unmapped reads, we show that our method can help to identify large-scale variations that have been previously missed.

Description: Activating transcription factor 3 (ATF3) is a stress-induced transcriptional regulator in eukaryote. The role of ATF3 in cancer has been well defined, but how ATF3 functions in bacterial infection is not well understood. Pneumococcal infection has been shown to induce ATF3 expression, which subsequently enhances cytokine production and provides protection from lethal Streptococcus pneumoniae infection, but the role of ATF3 in other Gram-positive (G + ) infections remains unclear. Here, we report that infection with other G + bacteria ( Staphylococcus aureus and Listeria monocytogenes ) and with G – bacteria (uropathogenic Escherichia coli ) also significantly induced ATF3 expression. Moreover, the production of cytokines (tumor necrosis factor alpha [TNF]-α, interleukin [IL]-1β, IL-6 and interferon [IFN]-) was enhanced by ATF3 in S. aureus and L. monocytogenes infection, but decreased in uropathogenic E. coli (UPEC) infection. In addition, in S. aureus and L. monocytogenes infections, ATF3 WT mice cleared bacteria more efficiently and had higher survival rates than ATF3 knockout mice. However, in UPEC infection, no significant difference was found in survival rate. Taken together, these data suggest that ATF3 provides protection from S. aureus and L. monocytogenes infections; however, the role of ATF3 in UPEC infection is more complicated and should be further elucidated.

Description: Any given human individual carries multiple genetic variants that disrupt protein-coding genes, through structural variation, as well as nucleotide variants and indels. Predicting the phenotypic consequences of a gene disruption remains a significant challenge. Current approaches employ information from a range of biological networks to predict which human genes are haploinsufficient (meaning two copies are required for normal function) or essential (meaning at least one copy is required for viability). Using recently available study gene sets, we show that these approaches are strongly biased towards providing accurate predictions for well-studied genes. By contrast, we derive a haploinsufficiency score from a combination of unbiased large-scale high-throughput datasets, including gene co-expression and genetic variation in over 6000 human exomes. Our approach provides a haploinsufficiency prediction for over twice as many genes currently unassociated with papers listed in Pubmed as three commonly-used approaches, and outperforms these approaches for predicting haploinsufficiency for less-studied genes. We also show that fine-tuning the predictor on a set of well-studied ‘gold standard’ haploinsufficient genes does not improve the prediction for less-studied genes. This new score can readily be used to prioritize gene disruptions resulting from any genetic variant, including copy number variants, indels and single-nucleotide variants.

Description: Meta-analysis of gene expression has enabled numerous insights into biological systems, but current methods have several limitations. We developed a method to perform a meta-analysis using the elastic net, a powerful and versatile approach for classification and regression. To demonstrate the utility of our method, we conducted a meta-analysis of lung cancer gene expression based on publicly available data. Using 629 samples from five data sets, we trained a multinomial classifier to distinguish between four lung cancer subtypes. Our meta-analysis-derived classifier included 58 genes and achieved 91% accuracy on leave-one-study-out cross-validation and on three independent data sets. Our method makes meta-analysis of gene expression more systematic and expands the range of questions that a meta-analysis can be used to address. As the amount of publicly available gene expression data continues to grow, our method will be an effective tool to help distill these data into knowledge.

Description: The presence of carbapenemase gene bla KPC-2 in a wide variety of plasmids, especially conjugative plasmids, is key to the rapid, worldwide spread of carbapenemase enzymes. Thirty-eight, non-duplicated, carbapenem-resistant, clinical Klebsiella pneumoniae isolates were collected, all carrying bla KPC-2 -bearing plasmids. Relaxase analysis was used to classify these plasmids; 8 and 30 plasmids belonged to the MOB P3 and MOB F12 subfamilies, respectively. Phylogenetic analysis revealed two genetic subclades in the MOB F12 subfamily and suggested that these subclades might not have originated from the same ancestor. Crossing PCR, used to sequence fully the type IV secretion system (T4SS, essential structures for conjugative plasmids) of the MOB F12 plasmids, found that T4SSs were distinctively different in certain functional genes, e.g. traS and traG. In conclusion, this study delineated the evolution of bla KPC-2 -bearing plasmids at Huashan Hospital, Shanghai, China. The plasmids bearing bla KPC-2 were diverse and the MOB F12 plasmids were dominant in clinical K. pneumoniae isolates.

Description: Advanced sequencing technologies have generated a plethora of data for many chromatin marks in multiple tissues and cell types, yet there is lack of a generalized tool for optimal utility of those data. A major challenge is to quantitatively model the epigenetic dynamics across both the genome and many cell types for understanding their impacts on differential gene regulation and disease. We introduce IDEAS, an i ntegrative and d iscriminative e pigenome a nnotation s ystem, for jointly characterizing epigenetic landscapes in many cell types and detecting differential regulatory regions. A key distinction between our method and existing state-of-the-art algorithms is that IDEAS integrates epigenomes of many cell types simultaneously in a way that preserves the position-dependent and cell type-specific information at fine scales, thereby greatly improving segmentation accuracy and producing comparable annotations across cell types.

Description: Enteroaggregative Escherichia coli (EAEC) is an important diarrhoeal pathogen causing diseases in multiple epidemiological and clinical settings. In developing countries like India, diarrhoeal diseases are one of the major killers among paediatric population and oddly, few studies are available from Indian paediatric population on the variability of EAEC virulence genes. In this study, we examined the distribution of plasmid and chromosomal-encoded virulence determinants in EAEC isolates, and analysed cytokines response generated against EAEC with specific aggregative adherence fimbriae (AAF) type in duodenal biopsies using in vitro organ culture (IVOC) mimicking in vivo conditions. Different virulence marker combinations among strains were reflected as a function of specific adhesins signifying EAEC heterogeneity. fis gene emerged as an important genetic marker apart from aggA and aap . Further, EAEC infection in IVOC showed upregulation of IL-8, IL-1β, IL-6, TNF-α and TLR-5 expression. EAEC with AAFII induced significant TLR-5 and IL-8 response, conceivably owing to more pathogenicity markers. This study sheds light on the pattern of EAEC pathotypes prevalent in North Indian paediatric population and highlights the presence of unique virulence combinations in pathogenic strains. Thus, evident diversity in EAEC virulence and multifaceted bacteria-host crosstalk can provide useful insights for the strategic management of diarrhoeal diseases in India, where diarrhoeal outbreaks are more frequent.

Description: Coagulase-negative staphylococci are thought to act as reservoirs of antibiotic resistance genes that can be transferred to Staphylococcus aureus , thus hindering the combat of this bacterium. In this work, we analyzed the presence of plasmids conferring resistance to the antibiotic mupirocin—widely used to treat and prevent S. aureus infections in hospital environments—in nosocomial S. haemolyticus strains. About 12% of the 75 strains tested were resistant to mupirocin, and this phenotype was correlated with the presence of plasmids. These plasmids were shown to be diverse, being either conjugative or mobilizable, and capable of transferring mupirocin resistance to S. aureus . Our findings reinforce that S. haemolyticus , historically and mistakenly considered as a less important pathogen, is a reservoir of resistance genes which can be transferred to other bacteria, such as S. aureus , emphasizing the necessity of more effective strategies to detect and combat this emergent opportunistic pathogen.

Description: Emergence of races in Fusarium oxysporum f. sp. lycopersici ( Fol ) is caused by loss or mutation of at least one avirulence ( AVR ) gene. The product of AVR1 is a small protein (Avr1) secreted by Fol in tomato xylem sap during infection. This protein triggers Fol race 1 specific resistance (I) in tomato, indicating that AVR1 is an AVR gene. Deletion of AVR1 in race 1 resulted in the emergence of race 2, and an additional mutation in AVR2 generated race 3. Previously, we reported a new biotype of race 3, KoChi-1, in which AVR1 was truncated by a transposon Hormin , which suggested a new route to evolution of races in Fol . However, to date no race 2 isolate carrying Hormin -truncated AVR1 has been reported. In this report, we describe such isolates, represented by Chiba-5, in which Hormin insertion occurred in AVR1 at a position different from that in KoChi-1. AVR1 truncation in both isolates resulted in production of defective Avr1 proteins. Chiba-5 and KoChi-1 belong to different phylogenetic clades, A1 and A2, respectively, suggesting that insertion of Hormin in AVR1 in Chiba-5 and KoChi-1 occurred as independent evolutionary events.

Description: The LBIT-1200 strain of Bacillus thuringiensis was recently isolated from soil, and showed a 6.4 and 9.5 increase in toxicity, against Manduca sexta and Trichoplusia ni , respectively, compared to HD-73. However, LBIT-1200 was still highly similar to HD-73, including the production of bipyramidal crystals containing only one protein of ~130 000 kDa, its flagellin gene sequence related to the kurstaki serotype, plasmid and RepPCR patterns similar to HD-73, no production of β-exotoxin and no presence of VIP genes. Sequencing of its cry gene showed the presence of a cry1Ac -type gene with four amino acid differences, including two amino acid replacements in domain III, compared to Cry1Ac1, which may explain its higher toxicity. In conclusion, the LBIT-1200 strain is a variant of the HD-73 strain but shows a much higher toxicity, which makes this new strain an important candidate to be developed as a bioinsecticide, once it passes other tests, throughout its biotechnological development.

Description: A majority of large-scale bacterial genome rearrangements involve mobile genetic elements such as insertion sequence (IS) elements. Here we report novel insertions and excisions of IS elements and recombination between homologous IS elements identified in a large collection of Escherichia coli mutation accumulation lines by analysis of whole genome shotgun sequencing data. Based on 857 identified events (758 IS insertions, 98 recombinations and 1 excision), we estimate that the rate of IS insertion is 3.5 x 10 –4 insertions per genome per generation and the rate of IS homologous recombination is 4.5 x 10 –5 recombinations per genome per generation. These events are mostly contributed by the IS elements IS 1 , IS 2 , IS 5 and IS 186 . Spatial analysis of new insertions suggest that transposition is biased to proximal insertions, and the length spectrum of IS-caused deletions is largely explained by local hopping. For any of the ISs studied there is no region of the circular genome that is favored or disfavored for new insertions but there are notable hotspots for deletions. Some elements have preferences for non-coding sequence or for the beginning and end of coding regions, largely explained by target site motifs. Interestingly, transposition and deletion rates remain constant across the wild-type and 12 mutant E. coli lines, each deficient in a distinct DNA repair pathway. Finally, we characterized the target sites of four IS families, confirming previous results and characterizing a highly specific pattern at IS 186 target-sites, 5'-GGGG(N6/N7)CCCC-3'. We also detected 48 long deletions not involving IS elements.

Description: Gastrointestinal (GI) leakage in Clostridium difficile -associated diarrhea (CDAD) is well known but is not routinely assessed in clinical practice. Serum (1-〉3)-β-D-glucan (BG), a fungal cell wall component used as a biomarker for invasive fungal disease, was tested in a CDAD mouse model with and without probiotics. Higher serum fluorescein isothiocyanate-dextran (FITC-dextran) and spontaneous gram-negative bacteremia, GI leakage indicators, were frequently found in CDAD mice, which died compared with those which survived. BG, serum macrophage inflammatory protein-2 and FITC-dextran but not quantitative blood bacterial count differentiated the clinical severity. Interestingly, a specific dose of Lactobacillus rhamnosus L34 attenuated CDAD and decreased serum BG and FITC-dextran, but not other parameters. BG also showed a higher area under the receiver operating characteristic curve for 7-day mortality than FITC-dextran. Fifty-five percent of CDAD mice with BG ≥ 60 pg/ml (the human negative cut-off value for invasive fungal disease) at 1 day after C. difficile gavage died within 7 days. In conclusion, s erum BG was elevated in mice with severe CDAD, an established model of GI leakage with a strong association with mortality rate. BG monitoring in patients with CDAD is of interest as both a potential prognostic tool and a therapeutic efficacy indicator.

Description: Chlamydial species are common intracellular parasites that cause various diseases, mainly characterized by persistent infection, which lead to inflammatory responses modulated by pattern recognition receptors (PRRs). The best understood PRRs are the extracellular Toll-like receptors, but recent significant advances have focused on two important proteins, NOD1 and NOD2, which are members of the intracellular nucleotide-binding oligomerization domain receptor family and are capable of triggering the host innate immune signaling pathways. This results in the production of pro-inflammatory cytokines, which is vital for an adequate host defense against intracellular chlamydial infection. NOD1/2 ligands are known to derive from peptidoglycan, and the latest research has resolved the paradox of whether chlamydial species possess this bacterial cell wall component; this finding is likely to promote in-depth investigations into the interaction between the NOD proteins and chlamydial pathogens. In this review, we summarize the basic characteristics and signal transduction functions of NOD1 and NOD2 and highlight the new research on the roles of NOD1 and NOD2 in the host defense against chlamydial infection.

Description: Two strains of Aeromonas salmonicida , YK and BG, were isolated from largemouth bronze gudgeon and northern whitefish in China, and identified as A. salmonicida subsp. salmonicida based on phylogenetic analysis of vapA and 16S rRNA gene sequences. YK and BG originated from freshwater fish, one of which belonged to the cyprinid family, and the strains showed a difference in virulence. Subsequently, we performed whole genome sequencing of the strains, and comparison of their genomic sequences to the genome of the A449 reference strain revealed various genomic rearrangements, including a new variant of the genomic island AsaGEI in BG, designated as AsaGEI2c . This is the first report on a GEI of A. salmonicida strain from China. Furthermore, both YK and BG strains contained a Tn7 transposon inserted at the same position in the chromosome. Finally, IS-dependent rearrangements on pAsa5 are deemed likely to have occurred, with omission of the resD gene in both strains as well as omission of genes related to the IncF conjugal transfer system in the YK isolate. This study demonstrates that A. salmonicida subsp. salmonicida can infect non-salmonids (cyprinids) in addition to salmonids, and that AsaGEI2c might be useful as a geographical indicator of Chinese A. salmonicida subsp. salmonicida isolates.

Description: Burkholderia pseudomallei is an intracellular Gram-negative bacterial pathogen and the causative agent of melioidosis, a widespread disease in Southeast Asia. Reactive nitrogen, in an intermediate form of nitric oxide (NO), is one of the first lines of defense used by host cells to eliminate intracellular pathogens, through the stimulation of inducible nitric oxide synthase (iNOS). Studies in phagocytotic cells have shown that the iNOS response is muted in B. pseudomallei infection, and implicated the rpoS sigma factor as a key regulatory factor mediating suppression. The liver is a main visceral organ affected by B. pseudomallei , and there is little knowledge about the interaction of liver cells and B. pseudomallei . This study investigated the induction of iNOS, as well as autophagic flux and light-chain 3 (LC3) localization in human liver (HC04) cells in response to infection with B. pseudomallei and its rpoS deficient mutant. Results showed that the rpoS mutant was unable to suppress iNOS induction and that the mutant showed less induction of autophagy and lower co-localization with LC3, and this was coupled with a lower intracellular growth rate. Combining these results suggest that B. pseudomallei rpoS is an important factor in establishing infection in liver cells.

Description: Extensive and multi-dimensional data sets generated from recent cancer omics profiling projects have presented new challenges and opportunities for unraveling the complexity of cancer genome landscapes. In particular, distinguishing the unique complement of genes that drive tumorigenesis in each patient from a sea of passenger mutations is necessary for translating the full benefit of cancer genome sequencing into the clinic. We address this need by presenting a data integration framework (OncoIMPACT) to nominate patient-specific driver genes based on their phenotypic impact. Extensive in silico and in vitro validation helped establish OncoIMPACT's robustness, improved precision over competing approaches and verifiable patient and cell line specific predictions (2/2 and 6/7 true positives and negatives, respectively). In particular, we computationally predicted and experimentally validated the gene TRIM24 as a putative novel amplified driver in a melanoma patient. Applying OncoIMPACT to more than 1000 tumor samples, we generated patient-specific driver gene lists in five different cancer types to identify modes of synergistic action. We also provide the first demonstration that computationally derived driver mutation signatures can be overall superior to single gene and gene expression based signatures in enabling patient stratification and prognostication. Source code and executables for OncoIMPACT are freely available from http://sourceforge.net/projects/oncoimpact .

Description: Analysis of rewired upstream subnetworks impacting downstream differential gene expression aids the delineation of evolving molecular mechanisms. Cumulative statistics based on conventional differential correlation are limited for subnetwork rewiring analysis since rewiring is not necessarily equivalent to change in correlation coefficients. Here we present a computational method ChiNet to quantify subnetwork rewiring by statistical heterogeneity that enables detection of potential genotype changes causing altered transcription regulation in evolving organisms. Given a differentially expressed downstream gene set, ChiNet backtracks a rewired upstream subnetwork from a super-network including gene interactions known to occur under various molecular contexts. We benchmarked ChiNet for its high accuracy in distinguishing rewired artificial subnetworks, in silico yeast transcription-metabolic subnetworks, and rewired transcription subnetworks for Candida albicans versus Saccharomyces cerevisiae , against two differential-correlation based subnetwork rewiring approaches. Then, using transcriptome data from tolerant S. cerevisiae strain NRRL Y-50049 and a wild-type intolerant strain, ChiNet identified 44 metabolic pathways affected by rewired transcription subnetworks anchored to major adaptively activated transcription factor genes YAP1 , RPN4 , SFP1 and ROX1 , in response to toxic chemical challenges involved in lignocellulose-to-biofuels conversion. These findings support the use of ChiNet in rewiring analysis of subnetworks where differential interaction patterns resulting from divergent nonlinear dynamics abound.

Description: Ribosome biogenesis, a central and essential cellular process, occurs through sequential association and mutual co-folding of protein–RNA constituents in a well-defined assembly pathway. Here, we construct a network of co-evolving nucleotide/amino acid residues within the ribosome and demonstrate that assembly constraints are strong predictors of co-evolutionary patterns. Predictors of co-evolution include a wide spectrum of structural reconstitution events, such as cooperativity phenomenon, protein-induced rRNA reconstitutions, molecular packing of different rRNA domains, protein–rRNA recognition, etc. A correlation between folding rate of small globular proteins and their topological features is known. We have introduced an analogous topological characteristic for co-evolutionary network of ribosome, which allows us to differentiate between rRNA regions subjected to rapid reconstitutions from those hindered by kinetic traps. Furthermore, co-evolutionary patterns provide a biological basis for deleterious mutation sites and further allow prediction of potential antibiotic targeting sites. Understanding assembly pathways of multicomponent macromolecules remains a key challenge in biophysics. Our study provides a ‘proof of concept’ that directly relates co-evolution to biophysical interactions during multicomponent assembly and suggests predictive power to identify candidates for critical functional interactions as well as for assembly-blocking antibiotic target sites.

Description: The emergence of new sequencing technologies has facilitated the use of bacterial whole genome alignments for evolutionary studies and outbreak analyses. These datasets, of increasing size, often include examples of multiple different mechanisms of horizontal sequence transfer resulting in substantial alterations to prokaryotic chromosomes. The impact of these processes demands rapid and flexible approaches able to account for recombination when reconstructing isolates’ recent diversification. Gubbins is an iterative algorithm that uses spatial scanning statistics to identify loci containing elevated densities of base substitutions suggestive of horizontal sequence transfer while concurrently constructing a maximum likelihood phylogeny based on the putative point mutations outside these regions of high sequence diversity. Simulations demonstrate the algorithm generates highly accurate reconstructions under realistically parameterized models of bacterial evolution, and achieves convergence in only a few hours on alignments of hundreds of bacterial genome sequences. Gubbins is appropriate for reconstructing the recent evolutionary history of a variety of haploid genotype alignments, as it makes no assumptions about the underlying mechanism of recombination. The software is freely available for download at github.com/sanger-pathogens/Gubbins , implemented in Python and C and supported on Linux and Mac OS X.

Description: Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for ‘targeted’ resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a ‘kmer’ strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.

Description: Genetic variants in or near miRNA genes can have profound effects on miRNA expression and targeting. As user-friendly software for the impact prediction of miRNA variants on a large scale is still lacking, we created a tool called miRVaS. miRVaS automates this prediction by annotating the location of the variant relative to functional regions within the miRNA hairpin (seed, mature, loop, hairpin arm, flanks) and by annotating all predicted structural changes within the miRNA due to the variant. In addition, the tool defines the most important region that is predicted to have structural changes and calculates a conservation score that is indicative of the reliability of the structure prediction. The output is presented in a tab-separated file, which enables fast screening, and in an html file, which allows visual comparison between wild-type and variant structures. All separate images are provided for downstream use. Finally, we tested two different approaches on a small test set of published functionally validated genetic variants for their capacity to predict the impact of variants on miRNA expression.

Description: Analysis of RNA-seq data often detects numerous ‘non-co-linear’ (NCL) transcripts, which comprised sequence segments that are topologically inconsistent with their corresponding DNA sequences in the reference genome. However, detection of NCL transcripts involves two major challenges: removal of false positives arising from alignment artifacts and discrimination between different types of NCL transcripts ( trans -spliced, circular or fusion transcripts). Here, we developed a new NCL-transcript-detecting method (‘NCLscan’), which utilized a stepwise alignment strategy to almost completely eliminate false calls (〉98% precision) without sacrificing true positives, enabling NCLscan outperform 18 other publicly-available tools (including fusion- and circular-RNA-detecting tools) in terms of sensitivity and precision, regardless of the generation strategy of simulated dataset, type of intragenic or intergenic NCL event, read depth of coverage, read length or expression level of NCL transcript. With the high accuracy, NCLscan was applied to distinguishing between trans -spliced, circular and fusion transcripts on the basis of poly(A)- and nonpoly(A)-selected RNA-seq data. We showed that circular RNAs were expressed more ubiquitously, more abundantly and less cell type-specifically than trans -spliced and fusion transcripts. Our study thus describes a robust pipeline for the discovery of NCL transcripts, and sheds light on the fundamental biology of these non-canonical RNA events in human transcriptome.

Description: Bacteriophages are viruses that infect bacteria. There are an estimated 10 31 phage on the planet, making them the most abundant form of life. We are rapidly approaching the centenary of their identification, and yet still have only a limited understanding of their role in the ecology and evolution of bacterial populations. Temperate prophage carriage is often associated with increased bacterial virulence. The rise in use of technologies, such as genome sequencing and transcriptomics, has highlighted more subtle ways in which prophages contribute to pathogenicity. This review discusses the current knowledge of the multifaceted effects that phage can exert on their hosts and how this may contribute to bacterial adaptation during infection.

Description: Phthiocerol dimycocerosates (PDIMs) and phenolic glycolipids (PGLs) contribute to the pathogenicity of several mycobacteria. Biosynthesis of these virulence factors requires polyketide synthases and other enzymes that represent potential targets for the development of adjuvant antivirulence drugs. We used six isogenic Mycobacterium marinum mutants, each with a different gene knockout in the PDIM/PGL biosynthetic pathway, to probe the pleiotropy of mutations leading to PDIM – PGL – , PDIM + PGL – or PDIM – PGL + phenotypes. We evaluated the M. marinum mutants for changes in antibiotic susceptibility, cell envelope permeability, biofilm formation, surface properties, sliding motility and virulence in an amoeba model. The analysis also permitted us to begin exploring the hypothesis that different gene knockouts rendering the same PDIM and/or PGL deficiency phenotypes lead to M. marinum mutants with equivalent pleiotropic profiles. Overall, the results of our study revealed a complex picture of pleiotropic patterns emerging from different gene knockouts, uncovered unexpected phenotypic inequalities between mutants, and provided new insight into the phenotypic consequences of gene knockouts in the PDIM/PGL biosynthetic pathway.

Description: Culture medium from an isolate of the fungus Aspergillus candidus was extracted, fractionated and examined to discover compounds antagonistic to plant-parasitic nematodes that are important pathogens of agricultural crops. Column, thin layer and preparative chromatographies and spectral and elemental analyses, were used to isolate and identify two major constituents of an active fraction (Fraction F) obtained from the medium. Compound 1 was identified as 2-hydroxypropane-1, 2, 3-tricarboxylic acid (citric acid). Compound 2 was identified as 3-hydroxy-5-methoxy-3-(methoxycarbonyl)-5-oxopentanoic acid, an isomer of 1, 2-dimethyl citrate. Compound 1 and a citric acid standard, each tested at 50 mg mL –1 in water, decreased hatch from eggs of the plant-parasitic nematode Meloidogyne incognita by more than 94%, and completely immobilized second-stage juveniles after 4–6 days exposure. Fraction F and Compounds 1 and 2 decreased the mobility of adults of the plant-parasitic nematode Ditylenchus destructor in vitro . Fraction F (25 mg mL –1 ) inhibited mobility 〉99% at 72 hrs. Compounds 1 and 2 (50 mg mL –1 ) each inhibited mobility more than 25% at 24 hr and more than 50% at 72 hr. This is the first assignment of nematode-antagonistic properties to specifically identified A. candidus metabolites.

Description: An RNAi system based on T7 RNA polymerase (TRNAP) was designed and examined in Aspergillus fumigatus . This system consists of two elements; an inducible T7RNAP expressing cassette and an AMA1-based episomal RNAi plasmid. These constructs were transformed into the A. fumigatus protoplasts and the efficiency of this system was tested in downregulation of alb1 gene. Upon the induction of T7RNAP expression, the recombinant T7RNAP was able to recognize T7 promoters, which were located on the episomal plasmid and in opposite direction. As a result, the bidirectional transcription of alb1 fragment led to the silencing of the target gene. However, our results demonstrated that this silencing system is unstable and may not be applicable in preparation of RNAi libraries.

Description: Detecting allelic biases from high-throughput sequencing data requires an approach that maximises sensitivity while minimizing false positives. Here, we present Allelome.PRO, an automated user-friendly bioinformatics pipeline, which uses high-throughput sequencing data from reciprocal crosses of two genetically distinct mouse strains to detect allele-specific expression and chromatin modifications. Allelome.PRO extends approaches used in previous studies that exclusively analyzed imprinted expression to give a complete picture of the ‘allelome’ by automatically categorising the allelic expression of all genes in a given cell type into imprinted, strain-biased, biallelic or non-informative. Allelome.PRO offers increased sensitivity to analyze lowly expressed transcripts, together with a robust false discovery rate empirically calculated from variation in the sequencing data. We used RNA-seq data from mouse embryonic fibroblasts from F1 reciprocal crosses to determine a biologically relevant allelic ratio cutoff, and define for the first time an entire allelome. Furthermore, we show that Allelome.PRO detects differential enrichment of H3K4me3 over promoters from ChIP-seq data validating the RNA-seq results. This approach can be easily extended to analyze histone marks of active enhancers, or transcription factor binding sites and therefore provides a powerful tool to identify candidate cis regulatory elements genome wide.

Description: Identity by descent (IBD) can be reliably detected for long shared DNA segments, which are found in related individuals. However, many studies contain cohorts of unrelated individuals that share only short IBD segments. New sequencing technologies facilitate identification of short IBD segments through rare variants, which convey more information on IBD than common variants. Current IBD detection methods, however, are not designed to use rare variants for the detection of short IBD segments. Short IBD segments reveal genetic structures at high resolution. Therefore, they can help to improve imputation and phasing, to increase genotyping accuracy for low-coverage sequencing and to increase the power of association studies. Since short IBD segments are further assumed to be old, they can shed light on the evolutionary history of humans. We propose HapFABIA, a computational method that applies biclustering to identify very short IBD segments characterized by rare variants. HapFABIA is designed to detect short IBD segments in genotype data that were obtained from next-generation sequencing, but can also be applied to DNA microarray data. Especially in next-generation sequencing data, HapFABIA exploits rare variants for IBD detection. HapFABIA significantly outperformed competing algorithms at detecting short IBD segments on artificial and simulated data with rare variants. HapFABIA identified 160 588 different short IBD segments characterized by rare variants with a median length of 23 kb (mean 24 kb) in data for chromosome 1 of the 1000 Genomes Project. These short IBD segments contain 752 000 single nucleotide variants (SNVs), which account for 39% of the rare variants and 23.5% of all variants. The vast majority—152 000 IBD segments—are shared by Africans, while only 19 000 and 11 000 are shared by Europeans and Asians, respectively. IBD segments that match the Denisova or the Neandertal genome are found significantly more often in Asians and Europeans but also, in some cases exclusively, in Africans. The lengths of IBD segments and their sharing between continental populations indicate that many short IBD segments from chromosome 1 existed before humans migrated out of Africa. Thus, rare variants that tag these short IBD segments predate human migration from Africa. The software package HapFABIA is available from Bioconductor. All data sets, result files and programs for data simulation, preprocessing and evaluation are supplied at http://www.bioinf.jku.at/research/short-IBD .

Description: Glioma is the most common and fatal primary brain tumour with poor prognosis; however, the functional roles of miRNAs in glioma malignant progression are insufficiently understood. Here, we used an integrated approach to identify miRNA functional targets during glioma malignant progression by combining the paired expression profiles of miRNAs and mRNAs across 160 Chinese glioma patients, and further constructed the functional miRNA–mRNA regulatory network. As a result, most tumour-suppressive miRNAs in glioma progression were newly discovered, whose functions were widely involved in gliomagenesis. Moreover, three miRNA signatures, with different combinations of hub miRNAs (regulations≥30) were constructed, which could independently predict the survival of patients with all gliomas, high-grade glioma and glioblastoma. Our network-based method increased the ability to identify the prognostic biomarkers, when compared with the traditional method and random conditions. Hsa-miR-524-5p and hsa-miR-628-5p, shared by these three signatures, acted as protective factors and their expression decreased gradually during glioma progression. Functional analysis of these miRNA signatures highlighted their critical roles in cell cycle and cell proliferation in glioblastoma malignant progression, especially hsa-miR-524-5p and hsa-miR-628-5p exhibited dominant regulatory activities. Therefore, network-based biomarkers are expected to be more effective and provide deep insights into the molecular mechanism of glioma malignant progression.

Description: A major challenge in cancer genomics is uncovering genes with an active role in tumorigenesis from a potentially large pool of mutated genes across patient samples. Here we focus on the interactions that proteins make with nucleic acids, small molecules, ions and peptides, and show that residues within proteins that are involved in these interactions are more frequently affected by mutations observed in large-scale cancer genomic data than are other residues. We leverage this observation to predict genes that play a functionally important role in cancers by introducing a computational pipeline ( http://canbind.princeton.edu ) for mapping large-scale cancer exome data across patients onto protein structures, and automatically extracting proteins with an enriched number of mutations affecting their nucleic acid, small molecule, ion or peptide binding sites. Using this computational approach, we show that many previously known genes implicated in cancers are enriched in mutations within the binding sites of their encoded proteins. By focusing on functionally relevant portions of proteins—specifically those known to be involved in molecular interactions—our approach is particularly well suited to detect infrequent mutations that may nonetheless be important in cancer, and should aid in expanding our functional understanding of the genomic landscape of cancer.

Description: Genome duplication with hybridization, or allopolyploidization, occurs commonly in plants, and is considered to be a strong force for generating new species. However, genome-wide quantification of homeolog expression ratios was technically hindered because of the high homology between homeologous gene pairs. To quantify the homeolog expression ratio using RNA-seq obtained from polyploids, a new method named HomeoRoq was developed, in which the genomic origin of sequencing reads was estimated using mismatches between the read and each parental genome. To verify this method, we first assembled the two diploid parental genomes of Arabidopsis halleri subsp. gemmifera and Arabidopsis lyrata subsp. petraea ( Arabidopsis petraea subsp. umbrosa ), then generated a synthetic allotetraploid, mimicking the natural allopolyploid Arabidopsis kamchatica . The quantified ratios corresponded well to those obtained by Pyrosequencing. We found that the ratios of homeologs before and after cold stress treatment were highly correlated ( r = 0.870). This highlights the presence of nonstochastic polyploid gene regulation despite previous research identifying stochastic variation in expression. Moreover, our new statistical test incorporating overdispersion identified 226 homeologs (1.11% of 20 369 expressed homeologs) with significant ratio changes, many of which were related to stress responses. HomeoRoq would contribute to the study of the genes responsible for polyploid-specific environmental responses.

Description: Methods to interpret personal genome sequences are increasingly required. Here, we report a novel framework (EvoTol) to identify disease-causing genes using patient sequence data from within protein coding-regions. EvoTol quantifies a gene's intolerance to mutation using evolutionary conservation of protein sequences and can incorporate tissue-specific gene expression data. We apply this framework to the analysis of whole-exome sequence data in epilepsy and congenital heart disease, and demonstrate EvoTol's ability to identify known disease-causing genes is unmatched by competing methods. Application of EvoTol to the human interactome revealed networks enriched for genes intolerant to protein sequence variation, informing novel polygenic contributions to human disease.

Description: Detecting in vivo transcription factor (TF) binding is important for understanding gene regulatory circuitries. ChIP-seq is a powerful technique to empirically define TF binding in vivo . However, the multitude of distinct TFs makes genome-wide profiling for them all labor-intensive and costly. Algorithms for in silico prediction of TF binding have been developed, based mostly on histone modification or DNase I hypersensitivity data in conjunction with DNA motif and other genomic features. However, technical limitations of these methods prevent them from being applied broadly, especially in clinical settings. We conducted a comprehensive survey involving multiple cell lines, TFs, and methylation types and found that there are intimate relationships between TF binding and methylation level changes around the binding sites. Exploiting the connection between DNA methylation and TF binding, we proposed a novel supervised learning approach to predict TF–DNA interaction using data from base-resolution whole-genome methylation sequencing experiments. We devised beta-binomial models to characterize methylation data around TF binding sites and the background. Along with other static genomic features, we adopted a random forest framework to predict TF–DNA interaction. After conducting comprehensive tests, we saw that the proposed method accurately predicts TF binding and performs favorably versus competing methods.

Description: Enriching target sequences in sequencing libraries via capture hybridization to bait/probes is an efficient means of leveraging the capabilities of next-generation sequencing for obtaining sequence data from target regions of interest. However, homologous sequences from non-target regions may also be enriched by such methods. Here we investigate the fidelity of capture enrichment for complete mitochondrial DNA (mtDNA) genome sequencing by analyzing sequence data for nuclear copies of mtDNA (NUMTs). Using capture-enriched sequencing data from a mitochondria-free cell line and the parental cell line, and from samples previously sequenced from long-range PCR products, we demonstrate that NUMT alleles are indeed present in capture-enriched sequence data, but at low enough levels to not influence calling the authentic mtDNA genome sequence. However, distinguishing NUMT alleles from true low-level mutations (e.g. heteroplasmy) is more challenging. We develop here a computational method to distinguish NUMT alleles from heteroplasmies, using sequence data from artificial mixtures to optimize the method.

Description: Single nucleotide polymorphisms (SNPs) are increasingly used to tag genetic loci associated with phenotypes such as risk of complex diseases. Technically, this is done genome-wide without prior restriction or knowledge of biological feasibility in scans referred to as genome-wide association studies (GWAS). Depending on the linkage disequilibrium (LD) structure at a particular locus, such tagSNPs may be surrogates for many thousands of other SNPs, and it is difficult to distinguish those that may play a functional role in the phenotype from those simply genetically linked. Because a large proportion of tagSNPs have been identified within non-coding regions of the genome, distinguishing functional from non-functional SNPs has been an even greater challenge. A strategy was recently proposed that prioritizes surrogate SNPs based on non-coding chromatin and epigenomic mapping techniques that have become feasible with the advent of massively parallel sequencing. Here, we introduce an R/Bioconductor software package that enables the identification of candidate functional SNPs by integrating information from tagSNP locations, lists of linked SNPs from the 1000 genomes project and locations of chromatin features which may have functional significance. Availability: FunciSNP is available from Bioconductor (bioconductor.org).

Description: Reconstructing the evolutionary relationships of species is a major goal in biology. Despite the increasing number of completely sequenced genomes, a large number of phylogenetic projects rely on targeted sequencing and analysis of a relatively small sample of marker genes. The selection of these phylogenetic markers should ideally be based on accurate predictions of their combined, rather than individual, potential to accurately resolve the phylogeny of interest. Here we present and validate a new phylogenomics strategy to efficiently select a minimal set of stable markers able to reconstruct the underlying species phylogeny. In contrast to previous approaches, our methodology does not only rely on the ability of individual genes to reconstruct a known phylogeny, but it also explores the combined power of sets of concatenated genes to accurately infer phylogenetic relationships of species not previously analyzed. We applied our approach to two broad sets of cyanobacterial and ascomycetous fungal species, and provide two minimal sets of six and four genes, respectively, necessary to fully resolve the target phylogenies. This approach paves the way for the informed selection of phylogenetic markers in the effort of reconstructing the tree of life.

Description: A simple diagnosis of the presence or absence of an infection is an uninformative metric when individuals differ considerably in their tolerance to different infection loads or resistance to rates of disease progression. Models that incorporate the relationship between the progression of the infection with the potential alternate outcomes provide a far more powerful predictive tool than diagnosis alone. The global decline of amphibians has been amplified by Batrachochytrium dendrobatidis , a pathogen that can cause the fatal disease chytridiomycosis. We measured the infection load and observed signs of disease in Litoria aurea . Receiver operating characteristic curves were used to quantify the dissimilarity between the infection loads of L. aurea that showed signs associated with chytridiomycosis and those that did not. Litoria aurea had a 78% probability of developing chytridiomycosis past a threshold of 68 zoospore equivalents (ZE) per swab and chytridiomycosis occurred within a variable range of 0.5–490 ZE. Studies should incorporate a species-specific threshold as a predictor of chytridiomycosis, rather than a binary diagnosis. Measures of susceptibility to chytridiomycosis must account not only for the ability of B. dendrobatidis to increase its abundance on the skin of amphibians but also to determine how each species tolerates these infection loads.

Description: Visualization of chromosomal dynamics is important for understanding many fundamental intra-nuclear processes. Efficient and reliable live-cell multicolor labeling of chromosomal loci can realize this goal. However, the current methods are constrained mainly by insufficient labeling throughput, efficiency, flexibility as well as photostability. Here we have developed a new approach to realize dual-color chromosomal loci imaging based on a modified single-guide RNA (sgRNA) of the CRISPR/Cas9 system. The modification of sgRNA was optimized by structure-guided engineering of the original sgRNA, consisting of RNA aptamer insertions that bind fluorescent protein-tagged effectors. By labeling and tracking telomeres, centromeres and genomic loci, we demonstrate that the new approach is easy to implement and enables robust dual-color imaging of genomic elements. Importantly, our data also indicate that the fast exchange rate of RNA aptamer binding effectors makes our sgRNA-based labeling method much more tolerant to photobleaching than the Cas9-based labeling method. This is crucial for continuous, long-term tracking of chromosomal dynamics. Lastly, as our method is complementary to other live-cell genomic labeling systems, it is therefore possible to combine them into a plentiful palette for the study of native chromatin organization and genome ultrastructure dynamics in living cells.

Description: The control and prevention of biofilm-related infections is an important public healthcare issue. Given the increasing antibiotic resistance among bacteria and fungi that cause serious infections in humans, promotion of new strategies combating microorganisms has been essential. One attractive approach to inactivate microorganisms is the use of semiconductor photo-catalysis, which has become the subject of extensive research. In this study, the bactericidal properties of four photo-catalysts, TiO 2 , TiO 2 -S, TiO 2 -Eu and TiO 2 -Eu-S, were investigated against established 24, 48, 72 and 96 h biofilms of Enterococcus . The exposure of biofilms to the catalysts induced the production of superoxide radical anions. The best photo-catalytic inactivation was achieved with the TiO 2 -Eu-S and TiO 2 -S nanopowders and 24 h biofilms. Transmission electron microscopy images showed significant changes in the structure of the biofilm cells following photo-inactivation. The results suggest that doping with europium and modifying the surface with sulphate groups enhanced the bactericidal activity of the TiO 2 nanoparticles against enterococcal biofilms.

Description: The Pseudomonas aeruginosa Chp chemosensory system regulates twitching motility, intracellular adenosine 3 '' 5 ' -cyclic monophosphate (cAMP) levels and is postulated to be involved in directional twitching towards phosphatidylethanolamine (PE). Because PilJ is the only methyl-accepting chemotaxis protein (MCP) identified in the Chp system, we determined the role of PilJ in mediating signal transduction for the distinct outputs of this system. Mutants that lack the periplasmic domain of PilJ ( pilJ 74-273 ) showed lower levels of cAMP but retained directional twitching towards PE. While initial studies revealed reduced twitching motility by PilJ 74-273 , this was due to decreased cAMP levels. Our data illustrate the importance of the periplasmic domain of PilJ in regulating cAMP. This is the first time a defined domain within PilJ has been identified as having a distinct role in signal transduction.

Description: Histone-like nucleoid-structuring (H-NS) proteins, which are conserved in Gram-negative bacteria, bind DNA and act as the global transcriptional repressors. In this study, we identified and characterized the xrvC gene encoding a H-NS protein in Xathomonas oryzae pv. oryzae ( Xoo ) Philippines strain PXO99 A . Compared with the wild type, the xrvC -deficient mutant of PXO99 A (named PXO99 xrvC ) showed a reduced growth rate in both nutrient-rich and nutrient-limited media. Interestingly, PXO99 xrvC exhibited significantly reduced virulence on rice cultivar IRBB214, but its virulence on 31 other rice cultivars was not affected. Transcriptional analysis revealed that the expression of hrpG , hrpX and hpa1 and of 15 out of 18 tested non-TAL (transcription activator-like) effector genes was decreased significantly in the xrvC mutant compared with that in the wild type. In addition, loss of xrvC also impaired the induction of the rice susceptibility gene Os8N3 in IRBB214 by PXO99 A . Our results suggest that the xrvC gene is involved in bacterial growth, and it plays a vital role in virulence by positively regulating the expression of hrp genes and non-TAL effector genes in PXO99 A and the susceptibility gene Os8N3 in rice.

Description: In Xanthomonas oryzae pv. oryzae , the pathogen of bacterial leaf blight of rice, hrp gene expression is regulated by the key hrp regulators HrpG and HrpX. HrpG regulates hrpX and hrpA , and HrpX regulates the other hrp genes on hrpB–hrpF operons. We previously examined the expression of the HrpX-regulated hrp gene hrcU and demonstrated that hrp gene expression is highly induced in a certain nutrient-poor medium containing xylose. In the present study, we found that the induction level of HrpX-regulated hrp genes was higher in medium with xylose than in media with any other sugar sources (glucose, sucrose and fructose), but that expression of hrpG , hrpX and hrpA was independent of the sugar sources. In western blot analysis, the accumulation of HrpX was reduced in media with a sugar other than xylose, probably as a result of proteolysis, but the addition of xylose canceled this reduced accumulation of the protein. The results suggest that proteolysis of HrpX is an important hrp regulatory mechanism and that xylose specifically suppresses this proteolysis, resulting in active hrp gene expression in X. oryzae pv. oryzae .

Description: Antibiotic therapy has long term consequences in the intestinal microbiome. Clostridium difficile has a well-known role in antibiotic-associated diarrhea, but in addition, persistent infection with this organism may increase the risk for developing inflammatory bowel disease. Here, recent literature on how the intestinal microbiome is altered by antibiotic therapy is presented.

Description: The aim of the present study was to verify whether penicillin-resistant, ampicillin-susceptible Enterococcus faecalis (PRASEF) occurred in Brazil prior to the beginning of the 21st century, and to verify whether ampicillin susceptibility can predict susceptibility to other β-lactams in E. faecalis with this inconsistent phenotype. The presence of polymorphisms in the pbp4 gene and genetic diversity among the isolates were investigated. Of 21 PRASEF analyzed, 5 (23.8%) and 4 (19.0%) were imipenem and piperacillin resistant simultaneously by disk diffusion and broth dilution respectively, contradicting the current internationally accepted standards of susceptibility testing. Sequencing of pbp4 gene revealed an amino acid substitution (Asp-573-〉Glu) in all PRASEF isolates but not in the penicillin-susceptible, ampicillin-susceptible E. faecalis . Most PRASEF (90.5%) had related pulsed-field gel electrophoresis profiles, but were different from other PRASEF described to date. Results demonstrate that penicillin-resistant, ampicillin-susceptible phenotype was already a reality in the 1990s in E. faecalis isolates in different Brazilian states, and some of these isolates were also imipenem- and piperacillin-resistant; therefore, internationally accepted susceptibility criteria cannot be applied to these isolates. According to pbp4 gene sequencing, this study suggests that a specific amino acid substitution in pbp4 gene found in all PRASEF analyzed is associated with penicillin resistance.

Description: The Mycobacterium tuberculosis Rv0679c protein is a surface protein that contributes to host cell invasion. We previously showed that a single nucleotide transition of the Rv0679c gene leads to a single amino acid substitution from asparagine to lysine at codon 142 in the Beijing genotype family. In this study, we examined the immunological effect of this substitution. Several recombinant proteins were expressed in Escherichia coli and Mycobacterium smegmatis and characterized with antisera and two monoclonal antibodies named 5D4-C2 and 8G10-H2. A significant reduction of antibody binding was detected by enzyme-linked immunosorbent assay (ELISA) and western blot analysis in the Lys142-type protein. This reduction of 8G10-H2 binding was more significant, with the disappearance of a signal in the proteins expressed by recombinant mycobacteria in western blot analysis. In addition, epitope mapping analysis of the recombinant proteins showed a linear epitope by 5D4-C2 and a discontinuous epitope by 8G10-H2. The antibody recognizing the conformational epitope detected only mycobacterial Asn142-type recombinant protein. Our results suggest that a single amino acid substitution of Rv0679c has potency for antigenic change in Beijing genotype strains.

Description: Listeria monocytogenes possesses the highest number of leucine-rich repeat (LRR)-containing proteins among all Gram-positive bacteria; these LRR-containing molecules are known as the ‘internalin’ family. To understand the functions of largely uncharacterized LRR-containing molecules, we constructed seven deletion mutants in the L. monocytogenes H7858 strain targeting genes in this family and tested their virulence. Among the seven mutants, the LMOh7858_0369 strain and the LMOh7858_2546 strain showed significantly impaired invasiveness of HepG2 cells. We further tested the virulence of these two strains in the intravascular sepsis model using BALB/c mice. Interestingly, the LMOh7858_0369 strain showed significant reduction in organ colonization, bacteremia and invasion of the brain compared with the parental wild-type strain. Host immune responses to listerial intravascular infection were measured at 24 and 72 h post-infection. Transcript levels of several proinflammatory cytokines and chemokines were significantly lower when induced by the lmOh7858_0369 strain than when induced by the wild type. These results suggest that the putative LRR-containing protein encoded by LMOh7858_0369 might be a novel virulence factor of the L. monocytogenes H7858 strain.

Description: Several Gram-positive pathogens scavenge host-derived heme to satisfy their nutritional iron requirement. However, heme is a toxic molecule capable of damaging the bacterial cell. Gram-positive pathogens within the phylum Firmicutes overcome heme toxicity by sensing heme through HssRS, a two-component system that regulates the heme detoxification transporter HrtAB. Here we show that heme sensing by HssRS and heme detoxification by HrtAB occur in the insect pathogen Bacillus thuringiensis . We find that in B. thuringiensis , HssRS directly regulates an operon, hrmXY , encoding hypothetical membrane proteins that are not found in other Firmicutes with characterized HssRS and HrtAB systems. This novel HssRS-regulated operon or its orthologs BMB171_c3178 and BMB171_c3330 are required for maximal heme resistance. Furthermore, the activity of HrmXY is not dependent on expression of HrtAB. These results suggest that B. thuringiensis senses heme through HssRS and induces expression of separate membrane-localized systems capable of overcoming different aspects of heme toxicity.

Description: The aim of this study was to understand the pathogenesis of motile aeromonas septicemia caused by an emergent, high virulent Aeromonas hydrophila (vAh) in channel catfish, Ictalurus punctatus . Adipose fin clipped catfish were challenged with vAh using a waterborne challenge method, and the distribution of vAh over a time course was detected and quantified using real-time polymerase chain reaction. The results showed that 77.8% of fish died within 48 h post challenge with mean day to death of 1.5 days. At 2 h post challenge, vAh (inferred from genomic DNA copies or genome equivalents) was detected in all external and internal tissues sampled. Gill had the highest vAh cells at 1 h post challenge. Spleen harbored the most vAh cells among internal organs at 4 h post challenge. The tissues/organs with most vAh cells detected at 8 h post challenge were adipose fin, blood, intestine, kidney and skin, while liver showed the highest vAh cells at 24 h post challenge. These results suggest that vAh was able to rapidly proliferate and spread, following wound infection, through the fish blood circulation system and cause mortality within 8–24 h.

Description: The major Staphylococcus aureus autolysin, Atl, has been implicated in attachment to surfaces and release of extracellular DNA during biofilm formation under laboratory conditions. Consistent with this, polyclonal antibodies to the amidase and glucosaminidase domains of Atl inhibited in vitro biofilm formation. However, in a murine model of device-related infection the community-associated S. aureus strain USA300 LAC JE2 established a successful infection in the absence of atl . These data indicate that Atl activity is not required for biofilm production in this infection model and reveal the importance of characterizing the contribution of biofilm phenotypes to virulence under in vivo conditions.

Description: Environmental biosurveillance and microbial ecology studies use PCR-based assays to detect and quantify microbial taxa and gene sequences within a complex background of microorganisms. However, the fragmentary nature and growing quantity of DNA-sequence data make group-specific assay design challenging. We solved this problem by developing a software platform that enables PCR-assay design at an unprecedented scale. As a demonstration, we developed quantitative PCR assays for a globally widespread, ecologically important bacterial group in soil, Acidobacteria Group 1. A total of 33 684 Acidobacteria 16S rRNA gene sequences were used for assay design. Following 1 week of computation on a 376-core cluster, 83 assays were obtained. We validated the specificity of the top three assays, collectively predicted to detect 42% of the Acidobacteria Group 1 sequences, by PCR amplification and sequencing of DNA from soil. Based on previous analyses of 16S rRNA gene sequencing, Acidobacteria Group 1 species were expected to decrease in response to elevated atmospheric CO 2 . Quantitative PCR results, using the Acidobacteria Group 1-specific PCR assays, confirmed the expected decrease and provided higher statistical confidence than the 16S rRNA gene-sequencing data. These results demonstrate a powerful capacity to address previously intractable assay design challenges.

Description: Due to dramatic advances in DNA technology, quantitative measures of annotation data can now be obtained in continuous coordinates across the entire genome, allowing various heterogeneous ‘genomic landscapes’ to emerge. Although much effort has been devoted to comparing DNA sequences, not much attention has been given to comparing these large quantities of data comprehensively. In this article, we introduce a method for rapidly detecting local regions that show high correlations between genomic landscapes. We overcame the size problem for genome-wide data by converting the data into series of symbols and then carrying out sequence alignment. We also decomposed the oscillation of the landscape data into different frequency bands before analysis, since the real genomic landscape is a mixture of embedded and confounded biological processes working at different scales in the cell nucleus. To verify the usefulness and generality of our method, we applied our approach to well investigated landscapes from the human genome, including several histone modifications. Furthermore, by applying our method to over 20 genomic landscapes in human and 12 in mouse, we found that DNA replication timing and the density of Alu insertions are highly correlated genome-wide in both species, even though the Alu elements have amplified independently in the two genomes. To our knowledge, this is the first method to align genomic landscapes at multiple scales according to their shape.

Description: The numerous discovered cases of domesticated transposable element (TE) proteins led to the recognition that TEs are a significant source of evolutionary innovation. However, much less is known about the reverse process, whether and to what degree the evolution of TEs is influenced by the genome of their hosts. We addressed this issue by searching for cases of incorporation of host genes into the sequence of TEs and examined the systems-level properties of these genes using the Saccharomyces cerevisiae and Drosophila melanogaster genomes. We identified 51 cases where the evolutionary scenario was the incorporation of a host gene fragment into a TE consensus sequence, and we show that both the yeast and fly homologues of the incorporated protein sequences have central positions in the cellular networks. An analysis of selective pressure (Ka/Ks ratio) detected significant selection in 37% of the cases. Recent research on retrovirus-host interactions shows that virus proteins preferentially target hubs of the host interaction networks enabling them to take over the host cell using only a few proteins. We propose that TEs face a similar evolutionary pressure to evolve proteins with high interacting capacities and take some of the necessary protein domains directly from their hosts.

Description: Presence of excess unaltered, wild-type (WT) DNA providing no information of biological or clinical value often masks rare alterations containing diagnostic or therapeutic clues in cancer, prenatal diagnosis, infectious diseases or organ transplantation. With the surge of high-throughput technologies there is a growing demand for removing unaltered DNA over large pools-of-sequences. Here we present nuclease-assisted minor-allele enrichment with probe-overlap (NaME-PrO), a single-step approach with broad genome coverage that can remove WT-DNA from numerous sequences simultaneously, prior to genomic analysis. NaME-PrO employs a double-strand-DNA-specific nuclease and overlapping oligonucleotide-probes interrogating WT-DNA targets and guiding nuclease digestion to these sites. Mutation-containing DNA creates probe-DNA mismatches that inhibit digestion, thus subsequent DNA-amplification magnifies DNA-alterations at all selected targets. We demonstrate several-hundred-fold mutation enrichment in diverse human samples on multiple clinically relevant targets including tumor samples and circulating DNA in 50-plex reactions. Enrichment enables routine mutation detection at 0.01% abundance while by adjusting conditions it is possible to sequence mutations down to 0.00003% abundance, or to scan tumor-suppressor genes for rare mutations. NaME-PrO introduces a simple and highly parallel process to remove un-informative DNA sequences and unmask clinically and biologically useful alterations.

Description: Edwardsiella tarda is distributed widely in a variety of hosts. Eha has recently been found to be its virulence regulator. In order to explore the mechanism of its regulation, we investigated the survival rates of wild type strain ET13, and its eha mutant and complemented strains in RAW264.7 macrophages under light microscopic observation as well as by counting bacterial CFUs on the plates. All of the different strains could live within the macrophages; however, the intracellular numbers of the wild type were significantly higher than the mutant when the incubation time extended 4 h or 6 h ( P 〈 0.05). Furthermore, more ROS were produced by the mutant-infected cells, indicating that Eha may enhance ET13's capacity to detoxify ROS. In agreement with this, we found that the mutant exhibited more sensitivity by H 2 O 2 disk inhibitory assay and less survival ability with H 2 O 2 treatment. We further demonstrated that the bacterial antioxidant enzymes SodC and KatG were regulated by Eha with qRT-PCR and β -galactosidase assay. Collectively, our data show Eha is required for E. tarda to resist the oxidative stress from the macrophages.

Description: Direct interaction between pathogens and host cells often is a prerequisite for colonization, infection and dissemination. Regulated production of capsular polysaccharide (CPS), which is made of hyaluronic acid, is essential for the pathogenicity of Streptococcus equi subsp. Zooepidemicus (SEZ). Here, we constructed a CPS-deleted mutant and analyzed it along with the parental wild-type strain in attachment and invasion of mammalian epithelial and endothelial cell lines. The CPS-deleted mutant exhibited significant increase in adherence and invasion by several orders of magnitude compared with the wild-type strain through quantitative analysis and electron microscopy observation. After the wild-type strain was recovered from invaded cells, its morphology was analyzed by visual methods and scanning electron microscopy, which revealed that its capsule was almost completely absent. Capsule measurements showed a similar result in which CPS production was nearly attenuated to the same extent as in the CPS-deleted mutant. qPCR assays revealed a marked reduction in the transcriptional levels of the CPS biosynthesis genes, has operon. Moreover, the repression in capsular production was stable inheritance. Our findings indicate that SEZ is a facultative intracellular bacterium, capsule attenuation in SEZ contributes to attachment and invasion in interactions with host cells, and the active regulation of capsule breakdown is controlled by SEZ during internalization.

Description: Burkholderia pseudomallei causes melioidosis, a potentially fatal infectious disease in tropical and subtropical countries worldwide. The intracellular behaviour of this pathogen in host cells has been reported to impact the severity of melioidosis, including the development of septicaemia, a consequence of pneumonia melioidosis. We previously identified a predicted cation transporter protein, BPSS1228, that participates in the transitional stage of this intracellular pathogen. For further analysis, in this study B. pseudomallei bpss1228 mutant and complemented strains were constructed and bacterial infectivity on human lung epithelial cells, A549, investigated in vitro . Burkholderia pseudomallei bpss1228 mutant showed impaired bacterial adhesion and invasion into A549 cells compared with wild-type strain, while the deficient phenotypes were restored to wild-type levels by the complemented strain. Additionally, the inactivation of bpss1228 in the mutant strain affected flagella-based swimming on a semi-solid surface and resistance to acid stresses simulating intracellular environments. These observations of BPSS1228 relating to B. pseudomallei infection strategies shed a new light on its association with intracellular B. pseudomallei during the interaction with host cells.

Description: Formation of a transient sub-population of bacteria, referred to as persisters, is one of the most important and least understood mechanisms that bacteria employ to evade elimination. Persister cells appear to be slow-growing bacteria that are broadly protected from a wide range of antibiotics. Using both theoretical and experimental methods, we show that alternating the application and withdrawal of antibiotics can be an effective treatment—as long as the timing of the protocol is estimated with precision. More specifically, we demonstrate that timing the alternating treatment based on theoretical predictions is confirmed using experimental observations. These results support a large class of theoretical studies that show that, even without complete understanding of the biological mechanisms, these models can provide insight into properties of the system.

Description: The phytopathogen Pseudomonas syringae pv. tabaci 11528 ( P. syringae 11528), causing wild-fire disease in soybean and tobacco plants, processes PsyI-PsyR quorum-sensing (QS) system, in which PsyI is the N -(3-oxo-hexanoyl)-homoserine lactone (3OC6-HSL) synthase. In comparison to P. syringae 11528 AHL-deficient mutant, 845 3OC6-HSL-dependent genes were identified using RNA sequencing (RNA-seq) in the AHL-deficient mutant grown with exogenous 3OC6-HSL in the transition from the exponential to the stationary phase, and many of them were associated with virulence, which were negatively regulated. The gene ontology and KEGG pathway enrichment analysis of those genes presented that the most pronounced regulation was involved in bacterial motility. Moreover, similar expression profiles of genes during growth phases were observed in both the wild type and the AHL-deficient mutant with exogenous 3OC6-HSL compared with the AHL-deficient mutant. These findings imply that 3OC6-HSL has a critical contribution to the QS-dependent regulation on gene expression, and 3OC6-HSL-dependent regulation may play a significant role in plant infection.

Description: Helicobacter pylori commonly infects the epithelial layer of the human stomach and in some individuals causes peptic ulcers, gastric adenocarcinoma or gastric lymphoma. Helicobacter pylori is a genetically diverse species, and the most important bacterial virulence factor that increases the risk of developing disease, versus asymptomatic colonization, is the cytotoxin associated gene pathogenicity island ( cag PAI). Socially housed rhesus macaques are often naturally infected with H. pylori similar to that which colonizes humans, but little is known about the cag PAI. Here we show that H. pylori strains isolated from naturally infected rhesus macaques have a cag PAI very similar to that found in human clinical isolates, and like human isolates, it encodes a functional type IV secretion system. These results provide further support for the relevance of rhesus macaques as a valid experimental model for H. pylori infection in humans.

Description: Sporisorium scitamineum is the fungus that causes sugarcane smut disease. Despite of the importance of sugarcane for Brazilian agribusiness and the persistence of the pathogen in most cropping areas, genetic variation studies are still missing for Brazilian isolates. In this study, sets of isolates were analyzed using two molecular markers (AFLP and telRFLP) and ITS sequencing. Twenty-two whips were collected from symptomatic plants in cultivated sugarcane fields of Brazil. A total of 41 haploid strains of compatible mating types were selected from individual teliospores and used for molecular genetic analyses. telRFLP and ITS analyses were expanded to six Argentine isolates, where the sugarcane smut was first recorded in America. Genetic relationship among strains suggests the human-mediated dispersal of S. scitamineum within the Brazilian territory and between the two neighboring countries. Two genetically distinct groups were defined by the combined analysis of AFLP and telRFLP. The opposite mating-type strains derived from single teliospores were clustered together into these main groups, but had not always identical haplotypes. telRFLP markers analyzed over two generations of selfing and controlled outcrossing confirmed the potential for emergence of new variants and occurrence of recombination, which are relevant events for evolution of virulence and environmental adaptation.

Description: A number of studies have shown that transcriptome analysis in terms of chromosomal location can reveal regions of non-random transcriptional activity within the genome. Genomic clusters of differentially expressed genes can identify genomic patterns of structural organization, underlying copy number variations or long-range epigenetic regulation such as X-chromosome inactivation. Here we apply an integrative bioinformatics analysis to a collection of 315 freely available mouse pluripotent stem cell samples to discover transcriptional clusters in the genome. We show that over half of the analysed samples (56.83%) carry whole or partial-chromosome spanning clusters which recur in genomic regions previously implicated in chromosomal imbalances. Strikingly, we found that the presence of such large-clusters is linked to the differential expression of a limited number of genes, common to all samples carrying clusters irrespectively of the chromosome where the cluster is found. We have used these genes to train and test classification models that can predict samples that carry large-scale clusters on any chromosome with over 90% accuracy. Our findings suggest that there is a common downstream activation in these cells that affects a limited number of nodes. We propose that this effect is linked to selective advantage and identify potential driver genes.

Description: The rapidly growing amount of genomic sequence data being generated and made publicly available necessitate the development of new data storage and archiving methods. The vast amount of data being shared and manipulated also create new challenges for network resources. Thus, developing advanced data compression techniques is becoming an integral part of data production and analysis. The HapMap project is one of the largest public resources of human single-nucleotide polymorphisms (SNPs), characterizing over 3 million SNPs genotyped in over 1000 individuals. The standard format and biological properties of HapMap data suggest that a dedicated genetic compression method can outperform generic compression tools. We propose a compression methodology for genetic data by introducing H ap Z ipper , a lossless compression tool tailored to compress HapMap data beyond benchmarks defined by generic tools such as gzip , bzip2 and lzma . We demonstrate the usefulness of H ap Z ipper by compressing HapMap 3 populations to 〈5% of their original sizes. H ap Z ipper is freely downloadable from https://bitbucket.org/pchanda/hapzipper/downloads/HapZipper.tar.bz2 .

Description: Spliced alignment plays a central role in the precise identification of eukaryotic gene structures. Even though many spliced alignment programs have been developed, recent rapid progress in DNA sequencing technologies demands further improvements in software tools. Benchmarking algorithms under various conditions is an indispensable task for the development of better software; however, there is a dire lack of appropriate datasets usable for benchmarking spliced alignment programs. In this study, we have constructed two types of datasets: simulated sequence datasets and actual cross-species datasets. The datasets are designed to correspond to various real situations, i.e. divergent eukaryotic species, different types of reference sequences, and the wide divergence between query and target sequences. In addition, we have developed an extended version of our program Spaln , which incorporates two additional features to the scoring scheme of the original version, and examined this extended version, Spaln2, together with the original Spaln and other representative aligners based on our benchmark datasets. Although the effects of the modifications are not individually striking, Spaln2 is consistently most accurate and reasonably fast in most practical cases, especially for plants and fungi and for increasingly divergent pairs of target and query sequences.

Description: Various ‘omics’ technologies, including microarrays and gas chromatography mass spectrometry, can be used to identify hundreds of interesting genes, proteins and metabolites, such as differential genes, proteins and metabolites associated with diseases. Identifying metabolic pathways has become an invaluable aid to understanding the genes and metabolites associated with studying conditions. However, the classical methods used to identify pathways fail to accurately consider joint power of interesting gene/metabolite and the key regions impacted by them within metabolic pathways. In this study, we propose a powerful analytical method referred to as Subpathway-GM for the identification of metabolic subpathways. This provides a more accurate level of pathway analysis by integrating information from genes and metabolites, and their positions and cascade regions within the given pathway. We analyzed two colorectal cancer and one metastatic prostate cancer data sets and demonstrated that Subpathway-GM was able to identify disease-relevant subpathways whose corresponding entire pathways might be ignored using classical entire pathway identification methods. Further analysis indicated that the power of a joint genes/metabolites and subpathway strategy based on their topologies may play a key role in reliably recalling disease-relevant subpathways and finding novel subpathways.

Description: Multiplex analytical systems that allow detection of multiple nucleic acid targets in one assay can provide rapid characterization of a sample while still saving cost and resources. However, few systems have proven to offer a solution for mid-plex (e.g. 10- to 50-plex) analysis that is high throughput and cost effective. Here we describe the combined use of fluorescence color and melting temperature (T m ) as a virtual 2D label that enables homogenous detection of one order of magnitude more targets than current strategies on real-time polymerase chain reaction platform. The target was first hybridized with a pair of ligation oligonucleotides, one of which harbored an artificial sequence that had a unique T m when hybridized with a reporter fluorogenic probe. The ligated products were then amplified by a universal primer pair and denatured by a melting curve analysis procedure. The targets were identified by their respective T m values in the corresponding fluorescence detection channels. The proof-of-principle of this approach was validated by genotyping 15 high-risk human papillomaviruses and 48 human single-nucleotide polymorphisms. The robustness of this method was demonstrated by analyzing a large number of clinical samples in both cases. The combined merits of multiplexity, flexibility and simplicity should make this approach suitable for a variety of applications.

Description: The chromatin structure of eukaryotic telomeres plays an essential role in telomere functions. However, their study might be impaired by the presence of interstitial telomeric sequences (ITSs), which have a widespread distribution in different model systems. We have developed a simple approach to study the chromatin structure of Arabidopsis telomeres independently of ITSs by analyzing ChIP-seq data. This approach could be used to study the chromatin structure of telomeres in some other eukaryotes. The analysis of ChIP-seq experiments revealed that Arabidopsis telomeres have higher density of histone H3 than centromeres, which might reflects their short nucleosomal organization. These experiments also revealed that Arabidopsis telomeres have lower levels of heterochromatic marks than centromeres (H3K9 Me2 and H3K27 Me ), higher levels of some euchromatic marks (H3K4 Me2 and H3K9Ac) and similar or lower levels of other euchromatic marks (H3K4 Me3 , H3K36 Me2 , H3K36 Me3 and H3K18Ac). Interestingly, the ChIP-seq experiments also revealed that Arabidopsis telomeres exhibit high levels of H3K27 Me3 , a repressive mark that associates with many euchromatic genes. The epigenetic profile of Arabidopsis telomeres is closely related to the previously defined chromatin state 2. This chromatin state is found in 23% of Arabidopsis genes, many of which are repressed or lowly expressed. At least, in part, this scenario is similar in rice.

Description: DExD/H box RNA helicases play essential roles in various biological processes in prokaryotes and eukaryotes. By screening Pseudomonas aeruginosa strains with mutations in various DExD/H box helicase genes, we identified that deaD was required for bacterial cytotoxicity and virulence in a mouse acute pneumonia model. Compared to a wild-type strain and its complementation strain, the deaD mutant induced less production of proinflammatory cytokines, neutrophil infiltration and lung damage during infection. We further found that the RNA helicase activity of DeaD was required for the expression of type III secretion system (T3SS) genes. Overexpression of ExsA, a master activator of the T3SS, restored the expression of T3SS genes as well as the virulence of the deaD mutant, suggesting that the attenuated virulence of the deaD mutant was mainly due to the defective T3SS. Overall, our results reveal a role of DeaD in the virulence of P. aeruginosa .

Description: Klebsiella pneumoniae is an opportunistic pathogen that commonly causes nosocomial infections in the urinary tract, respiratory tract, lung, wound sites and blood in individuals with debilitating diseases. Klebsiella pneumoniae is still a cause of severe pneumonia in alcoholics in Africa and Asia, and the predominant primary pathogen of primary liver abscess in Taiwan and Southeast Asia, particularly in Asian and Hispanic patients, and individuals with diabetes mellitus. In the United States and Europe, K. pneumoniae infections are most frequently associated with nosocomial infections. The emergence of antibiotic-resistant strains of K. pneumoniae worldwide has become a cause of concern where extended-spectrum β-lactamases (ESBLs) and carbapenemase-producing strains have been isolated with increasing frequency. The pathogen's ability to form biofilms on inserted devices such as urinary catheter has been proposed as one of the important mechanisms in nosocomially acquired and persistent infections, adding to the increased resistance to currently used antibiotics. In this review, infections caused by K. pneumoniae , antibiotic resistance and formation of biofilm will be discussed.

Description: Copy-number variants (CNVs) are a major form of genetic variation and a risk factor for various human diseases, so it is crucial to accurately detect and characterize them. It is conceivable that allele-specific reads from high-throughput sequencing data could be leveraged to both enhance CNV detection and produce allele-specific copy number (ASCN) calls. Although statistical methods have been developed to detect CNVs using whole-genome sequence (WGS) and/or whole-exome sequence (WES) data, information from allele-specific read counts has not yet been adequately exploited. In this paper, we develop an integrated method, called AS-GENSENG, which incorporates allele-specific read counts in CNV detection and estimates ASCN using either WGS or WES data. To evaluate the performance of AS-GENSENG, we conducted extensive simulations, generated empirical data using existing WGS and WES data sets and validated predicted CNVs using an independent methodology. We conclude that AS-GENSENG not only predicts accurate ASCN calls but also improves the accuracy of total copy number calls, owing to its unique ability to exploit information from both total and allele-specific read counts while accounting for various experimental biases in sequence data. Our novel, user-friendly and computationally efficient method and a complete analytic protocol is freely available at https://sourceforge.net/projects/asgenseng/ .

Description: Increased sequencing of microbial genomes has revealed that prevailing prokaryotic species assignments can be inconsistent with whole genome information for a significant number of species. The long-standing need for a systematic and scalable species assignment technique can be met by the genome-wide Average Nucleotide Identity (gANI) metric, which is widely acknowledged as a robust measure of genomic relatedness. In this work, we demonstrate that the combination of gANI and the alignment fraction (AF) between two genomes accurately reflects their genomic relatedness. We introduce an efficient implementation of AF,gANI and discuss its successful application to 86.5M genome pairs between 13,151 prokaryotic genomes assigned to 3032 species. Subsequently, by comparing the genome clusters obtained from complete linkage clustering of these pairs to existing taxonomy, we observed that nearly 18% of all prokaryotic species suffer from anomalies in species definition. Our results can be used to explore central questions such as whether microorganisms form a continuum of genetic diversity or distinct species represented by distinct genetic signatures. We propose that this precise and objective AF,gANI-based species definition: the MiSI (Microbial Species Identifier) method, be used to address previous inconsistencies in species classification and as the primary guide for new taxonomic species assignment, supplemented by the traditional polyphasic approach, as required.

Description: Classically or alternatively activated macrophages (M1 and M2, respectively) play distinct and important roles for microbiocidal activity, regulation of inflammation and tissue homeostasis. Despite this, their transcriptional regulatory dynamics are poorly understood. Using promoter-level expression profiling by non-biased deepCAGE we have studied the transcriptional dynamics of classically and alternatively activated macrophages. Transcription factor (TF) binding motif activity analysis revealed four motifs, NFKB1_REL_RELA, IRF1,2, IRF7 and TBP that are commonly activated but have distinct activity dynamics in M1 and M2 activation. We observe matching changes in the expression profiles of the corresponding TFs and show that only a restricted set of TFs change expression. There is an overall drastic and transient up-regulation in M1 and a weaker and more sustainable up-regulation in M2. Novel TFs, such as Thap6, Maff , (M1) and Hivep1, Nfil3, Prdm1 , (M2) among others, were suggested to be involved in the activation processes. Additionally, 52 (M1) and 67 (M2) novel differentially expressed genes and, for the first time, several differentially expressed long non-coding RNA (lncRNA) transcriptome markers were identified. In conclusion, the finding of novel motifs, TFs and protein-coding and lncRNA genes is an important step forward to fully understand the transcriptional machinery of macrophage activation.

Description: In recent years, extraintestinal pathogenic Escherichia coli (ExPEC) has been found to pose a great threat to human and animal health, but its pathogenic mechanism is not fully understood yet. Capsular polysaccharide, an essential virulence factor in these bacteria, can damage the host immune system, and kpsM is a member of the gene cluster responsible for capsular polysaccharide synthesis. In this study, whole sequence alignment of the virulent strain PCN033 and the attenuated strain PCN061 revealed that kpsM exists in PCN033 but not in PCN061. To determine its function and biological characteristics, we deleted kpsM from PCN033 by homologous recombination. The results of adhesion assays, phagocytosis assays and serum bactericidal assays together with the results of colonization assays in mice indicate that the deletion of kpsM decreases the virulence of porcine ExPEC. Our findings about the biological characteristics of kpsM help to elucidate the complex pathogenic mechanism of ExPEC.

Description: The RAG1/RAG2 endonuclease initiates V(D)J recombination at antigen receptor loci but also binds to thousands of places outside of these loci. RAG2 localizes directly to lysine 4 trimethylated histone 3 (H3K4me3) through a plant homeodomain (PHD) finger. The relative contribution of RAG2-dependent and RAG1-intrinsic mechanisms in determining RAG1 binding patterns is not known. Through analysis of deep RAG1 ChIP-seq data, we provide a quantitative description of the forces underlying genome-wide targeting of RAG1. Surprisingly, sequence-specific DNA binding contributes minimally to RAG1 targeting outside of antigen receptor loci. Instead, RAG1 binding is driven by two distinct modes of interaction with chromatin: the first is driven by H3K4me3, promoter-focused and dependent on the RAG2 PHD, and the second is defined by H3K27Ac, enhancer-focused and dependent on ‘non-core’ portions of RAG1. Based on this and additional chromatin and genomic features, we formulated a predictive model of RAG1 targeting to the genome. RAG1 binding sites predicted by our model correlate well with observed patterns of RAG1-mediated breaks in human pro-B acute lymphoblastic leukemia. Overall, this study provides an integrative model for RAG1 genome-wide binding and off-target activity and reveals a novel role for the RAG1 non-core region in RAG1 targeting.

Description: Asexual development of phytopathogenic fungi such as Magnaporthe oryzae involves morphological changes that require spatiotemporal regulation of polarized growth. ADP-ribosylation factor 6 (Arf6) is a small GTPase known to regulate membrane trafficking and organization of the actin cytoskeleton at the cell surface, and consequently has an impact on cell morphology and polarity. In this study, we have functionally characterized the Arf6 homolog in M. oryzae , showing that arf6 exhibits hyperbranching at hyphal tips and morphologically abnormal conidia as a result of defective polarized growth. arf6 hyphae are also defective in endocytosis as evidenced by a significant delay of FM4-64 uptake. Most arf6 conidia display reduced conidial length, and have defects in conidial septum formation and nuclear distribution. Furthermore, arf6 conidia show a disorganized actin cytoskeleton with random distribution of actin patches at the cell cortex and reduced accumulation of tropomyosin. Arf6-GFP is found to concentrate at the septum area and possibly in endocytic vesicles. Taken together, our data indicate that Arf6 plays an essential role in endocytosis and polarity establishment during asexual development of M. oryzae.

Description: As a leading pathogen, Edwardsiella piscicida can cause hemorrhagic septicemia in fish and gastro-intestinal infections in humans. The two-component regulatory system EsrA-EsrB plays essential roles in pathogenesis through the type III and type VI secretion systems, and hemolysin production in E. piscicida . It is unclear whether other virulence- or stress response-associated genes are regulated by EsrA-EsrB. In this study, the proteomes of wild-type E. piscicida EIB202 and esrB mutant strains were compared to reveal EsrB regulon components after growth in Luria–Bertani broth (LB). Overall, the expression levels of nine genes exhibited significant changes, and five of them required the presence of EsrB, while others exhibited higher levels in the esrB mutant. The diverse functions of these proteins were identified, including amino acid metabolism, oxidative stress defense and energy production. Interestingly, superoxidase dismutase and thiol peroxidase were the most significantly down-regulated by EsrB. Furthermore, other reported reactive oxygen species (ROS) resistance-related genes were also down-regulated by EsrB as revealed by quantitative real-time. Compared with the wild-type and the complement strain esrB + , esrB displayed a significantly enhanced ROS resistance. These results demonstrated that EsrB plays important roles in the ROS resistance pathway in E. piscicida grown in LB conditions.

Description: The sepsis caused by Vibrio vulnificus is characterized by an average incubation period of 26 h and a high mortality rate exceeding 50%. The fast growth and dissemination of V. vulnificus in vivo lead to poor clinical outcomes in patients. Therefore, elucidation of the proliferation mechanisms of this organism in vivo may lead to the development of an effective therapeutic strategy. In this study, we focused on the low oxygen concentration in the intestinal milieu because of its drastic difference from that in air. Fumarate and nitrate reduction regulatory protein (FNR) is known to be a global transcriptional regulator for adaptation to anaerobic conditions in various bacteria. We generated a strain of V. vulnificus in which the fnr gene was replaced with an erythromycin resistance gene ( fnr :: erm mutant). When the fnr :: erm mutant was tested in a growth competition assay against the wild-type (WT) in vivo , the competitive index of fnr :: erm mutant to WT in the intestinal loop and liver was 0.378 ± 0.192 (mean ± SD) and 0.243 ± 0.123, respectively. These data suggested that FNR is important for the proliferation of V. vulnificus in the intestine to achieve a critical mass to be able to invade the systemic circulation.

Description: The transduction mediated by bacteriophages is considered to be one of the primary driving forces in horizontal gene transfer in staphylococci, which is crucial to their adaptation and successful evolution. For a transduction to be effective, it is generally accepted that the recipient strain should be susceptible to the transducing phage. In this study, we demonstrate that the plasmid DNAs are effectively transduced into the recipient Staphylococcus aureus strains in spite of their insensitivity to the lytic action of the transducing phage, provided that these phages adsorb effectively to the bacterial cells. The tetracycline and penicillinase plasmids were transduced to insensitive laboratory and clinical strains by bacteriophages 29, 52A and 80α as well as by prophage 53 and naturally occurring prophages induced from donor lysogenic strains. Comparable frequencies of transduction were achieved in both phage-sensitive and phage-insensitive recipient strains. We have demonstrated that such mechanisms as the restriction of DNA and lysogenic immunity which are responsible for insensitivity of cells to phages may not be a barrier to the transfer, maintenance and effective spread of plasmids to a wider range of potential recipients in the staphylococcal population.

Description: Helicobacter pylori is an important cause of gastric pathologies and persistent infection can lead to stomach cancer. Virulent H. pylori strains encode a type IV secretion system responsible for translocation of the oncogenic CagA protein into cells of the gastric mucosa. Gene HP0522 encodes the essential component Cag (Cag3), and we show by gel filtration and cross-linking that purified Cag forms high molecular mass complexes. In contrast, its interaction partner CagT is mostly monomeric, but co-fractionates after gel filtration. Analysis by transmission electron microscopy revealed that purified Cag complexes can self-assemble ring-like structures. Cag-overexpressing Escherichia coli exhibits membrane-associated circular profiles in regions of the cell envelope with intense immunogold labelling with a Cag-specific antiserum. Our results suggest that Cag has the capacity to form macromolecular structures contributing to the assembly of the type IV secretion system.