Metagenomic binning and association of plasmids with bacterial host genomes using DNA methylation.

Shotgun metagenomics methods enable characterization of microbial communities in human microbiome and environmental samples. Assembly of metagenome sequences does not output whole genomes, so computational binning methods have been developed to cluster sequences into genome 'bins'. These methods exploit sequence composition, species abundance, or chromosome organization but cannot fully distinguish closely related species and strains. We present a binning method that incorporates bacterial DNA methylation signatures, which are detected using single-molecule real-time sequencing. Our method takes advantage of these endogenous epigenetic barcodes to resolve individual reads and assembled contigs into species- and strain-level bins. We validate our method using synthetic and real microbiome sequences. In addition to genome binning, we show that our method links plasmids and other mobile genetic elements to their host species in a real microbiome sample. Incorporation of DNA methylation information into shotgun metagenomics analyses will complement existing methods to enable more accurate sequence binning.

Colorectal cancer stages transcriptome analysis.

Colorectal cancer (CRC) is the third most common cancer and the second leading cause of cancer-related deaths in the United States. The purpose of this study was to evaluate the gene expression differences in different stages of CRC. Gene expression data on 433 CRC patient samples were obtained from The Cancer Genome Atlas (TCGA). Gene expression differences were evaluated across CRC stages using linear regression. Genes with p≤0.001 in expression differences were evaluated further in principal component analysis and genes with p≤0.0001 were evaluated further in gene set enrichment analysis. A total of 377 patients with gene expression data in 20,532 genes were included in the final analysis. The numbers of patients in stage I through IV were 59, 147, 116 and 55, respectively. NEK4 gene, which encodes for NIMA related kinase 4, was differentially expressed across the four stages of CRC. The stage I patients had the highest expression of NEK4 genes, while the stage IV patients had the lowest expressions (p = 9*10-6). Ten other genes (RNF34, HIST3H2BB, NUDT6, LRCh4, GLB1L, HIST2H4A, TMEM79, AMIGO2, C20orf135 and SPSB3) had p value of 0.0001 in the differential expression analysis. Principal component analysis indicated that the patients from the 4 clinical stages do not appear to have distinct gene expression pattern. Network-based and pathway-based gene set enrichment analyses showed that these 11 genes map to multiple pathways such as meiotic synapsis and packaging of telomere ends, etc. Ten of these 11 genes were linked to Gene Ontology terms such as nucleosome, DNA packaging complex and protein-DNA interactions. The protein complex-based gene set analysis showed that four genes were involved in H2AX complex II. This study identified a small number of genes that might be associated with clinical stages of CRC. Our analysis was not able to find a molecular basis for the current clinical staging for CRC based on the gene expression patterns.

Association Between Early-Life Antibiotic Use and the Risk of Islet or Celiac Disease Autoimmunity.

Importance: Evidence is lacking regarding the consequences of antibiotic use in early life and the risk of certain autoimmune diseases. Objective: To test the association between early-life antibiotic use and islet or celiac disease (CD) autoimmunity in genetically at-risk children prospectively followed up for type 1 diabetes (T1D) or CD. Design, Setting, and Participants: HLA-genotyped newborns from Finland, Germany, Sweden, and the United States were enrolled in the prospective birth cohort of The Environmental Determinants of Diabetes in the Young (TEDDY) study between November 20, 2004, and July 8, 2010. The dates of analysis were November 20, 2004, to August 31, 2014. Individuals from the general population and those having a first-degree relative with T1D were enrolled if they had 1 of 9 HLA genotypes associated with a risk for T1D. Exposures: Parental reports of the most common antibiotics (cephalosporins, penicillins, and macrolides) used between age 3 months and age 4 years were recorded prospectively. Main Outcomes and Measures: Islet autoimmunity and CD autoimmunity were defined as being positive for islet or tissue transglutaminase autoantibodies at 2 consecutive clinic visits at least 3 months apart. Hazard ratios and 95% CIs calculated from Cox proportional hazards regression models were used to assess the relationship between antibiotic use in early life before seroconversion and the development of autoimmunity. Results: Participants were 8495 children (49.0% female) and 6558 children (48.7% female) enrolled in the TEDDY study who were tested for islet and tissue transglutaminase autoantibodies, respectively. Exposure to and frequency of use of any antibiotic assessed in this study in early life or before seroconversion did not influence the risk of developing islet autoimmunity or CD autoimmunity. Cumulative use of any antibiotic during the first 4 years of life was not associated with the appearance of any autoantibody (hazard ratio [HR], 0.98; 95% CI, 0.95-1.01), multiple islet autoantibodies (HR, 0.99; 95% CI, 0.95-1.03), or the transglutaminase autoantibody (HR, 1.00; 95% CI, 0.98-1.02). Conclusions and Relevance: The use of the most prescribed antibiotics during the first 4 years of life, regardless of geographic region, was not associated with the development of autoimmunity for T1D or CD. These results suggest that a risk of islet or tissue transglutaminase autoimmunity need not influence the recommendations for clinical use of antibiotics in young children at risk for T1D or CD.

Corrigendum: Integrating DNA Methylation and Gene Expression Data in the Development of the Soybean-Bradyrhizobium N2-Fixing Symbiosis.

[This corrects the article on p. 518 in vol. 7, PMID: 27148207.].

Integrating DNA Methylation and Gene Expression Data in the Development of the Soybean-Bradyrhizobium N2-Fixing Symbiosis.

Very little is known about the role of epigenetics in the differentiation of a bacterium from the free-living to the symbiotic state. Here genome-wide analysis of DNA methylation changes between these states is described using the model of symbiosis between soybean and its root nodule-forming, nitrogen-fixing symbiont, Bradyrhizobium diazoefficiens. PacBio resequencing of the B. diazoefficiens genome from both states revealed 43,061 sites recognized by five motifs with the potential to be methylated genome-wide. Of those sites, 3276 changed methylation states in 2921 genes or 35.5% of all genes in the genome. Over 10% of the methylation changes occurred within the symbiosis island that comprises 7.4% of the genome. The CCTTGAG motif was methylated only during symbiosis with 1361 adenosines methylated among the 1700 possible sites. Another 89 genes within the symbiotic island and 768 genes throughout the genome were found to have methylation and significant expression changes during symbiotic development. Of those, nine known symbiosis genes involved in all phases of symbiotic development including early infection events, nodule development, and nitrogenase production. These associations between methylation and expression changes in many B. diazoefficiens genes suggest an important role of the epigenome in bacterial differentiation to the symbiotic state.