Behaviors of nucleosomes with mutant histone H4s in euchromatic domains of living human cells.

Since Robert Feulgen first stained DNA in the cell, visualizing genome chromatin has been a central issue in cell biology to uncover how chromatin is organized and behaves in the cell. To approach this issue, we have developed single-molecule imaging of nucleosomes, a basic unit of chromatin, to unveil local nucleosome behavior in living cells. In this study, we investigated behaviors of nucleosomes with various histone H4 mutants in living HeLa cells to address the role of H4 tail acetylation, including H4K16Ac and others, which are generally associated with more transcriptionally active chromatin regions. We ectopically expressed wild-type (wt) or mutated H4s (H4K16 point; H4K5,8,12,16 quadruple; and H4 tail deletion) fused with HaloTag in HeLa cells. Cells that expressed wtH4-Halo, H4K16-Halo mutants, and multiple H4-Halo mutants had euchromatin-concentrated distribution. Consistently, the genomic regions of the wtH4-Halo nucleosomes corresponded to Hi-C contact domains (or topologically associating domains, TADs) with active chromatin marks (A-compartment). Utilizing single-nucleosome imaging, we found that none of the H4 deacetylation or acetylation mimicked H4 mutants altered the overall local nucleosome motion. This finding suggests that H4 mutant nucleosomes embedded in the condensed euchromatic domains with excess endogenous H4 nucleosomes cannot cause an observable change in the local motion. Interestingly, H4 with four lysine-to-arginine mutations displayed a substantial freely diffusing fraction in the nucleoplasm, whereas H4 with a truncated N-terminal tail was incorporated in heterochromatic regions as well as euchromatin. Our study indicates the power of single-nucleosome imaging to understand individual histone/nucleosome behavior reflecting chromatin environments in living cells.

Chromatin organization and behavior in HRAS-transformed mouse fibroblasts.

In higher eukaryotic cells, a string of nucleosomes, where long genomic DNA is wrapped around core histones, are rather irregularly folded into a number of condensed chromatin domains, which have been revealed by super-resolution imaging and Hi-C technologies. Inside these domains, nucleosomes fluctuate and locally behave like a liquid. The behavior of chromatin may be highly related to DNA transaction activities such as transcription and repair, which are often upregulated in cancer cells. To investigate chromatin behavior in cancer cells and compare those of cancer and non-cancer cells, we focused on oncogenic-HRAS (Gly12Val)-transformed mouse fibroblasts CIRAS-3 cells and their parental 10T1/2 cells. CIRAS-3 cells are tumorigenic and highly metastatic. First, we found that HRAS-induced transformation altered not only chromosome structure, but also nuclear morphology in the cell. Using single-nucleosome imaging/tracking in live cells, we demonstrated that nucleosomes are locally more constrained in CIRAS-3 cells than in 10T1/2 cells. Consistently, heterochromatin marked with H3K27me3 was upregulated in CIRAS-3 cells. Finally, Hi-C analysis showed enriched interactions of the B-B compartment in CIRAS-3 cells, which likely represents transcriptionally inactive chromatin. Increased heterochromatin may play an important role in cell migration, as they have been reported to increase during metastasis. Our study also suggests that single-nucleosome imaging provides new insights into how local chromatin is structured in living cells.

Transcriptome profile of subsynaptic myonuclei at the neuromuscular junction in embryogenesis.

Skeletal muscle fiber is a large syncytium with multiple and evenly distributed nuclei. Adult subsynaptic myonuclei beneath the neuromuscular junction (NMJ) express specific genes, the products of which coordinately function in the maintenance of the pre- and post-synaptic regions. However, the gene expression profiles that promote the NMJ formation during embryogenesis remain largely unexplored. We performed single-nucleus RNA sequencing (snRNA-seq) analysis of embryonic and neonatal mouse diaphragms, and found that each myonucleus had a distinct transcriptome pattern during the NMJ formation. Among the previously reported NMJ-constituting genes, Dok7, Chrna1, and Chrnd are specifically expressed in subsynaptic myonuclei at E18.5. In the E18.5 diaphragm, ca. 10.7% of the myonuclei express genes for the NMJ formation (Dok7, Chrna1, and Chrnd) together with four representative ß-catenin regulators (Amotl2, Ptprk, Fam53b, and Tcf7l2). Additionally, the temporal gene expression patterns of these seven genes are synchronized in differentiating C2C12 myoblasts. Amotl2 and Ptprk are expressed in the sarcoplasm, where ß-catenin serves as a structural protein to organize the membrane-anchored NMJ structure. In contrast, Fam53b and Tcf7l2 are expressed in the myonucleus, where ß-catenin serves as a transcriptional coactivator in Wnt/ß-catenin signaling at the NMJ. In C2C12 myotubes, knockdown of Amotl2 or Ptprk markedly, and that of Fam53b and Tcf7l2 less efficiently, impair the clustering of acetylcholine receptors. In contrast, knockdown of Fam53b and Tcf7l2, but not of Amotl2 or Ptprk, impairs the gene expression of Slit2 encoding an axonal attractant for motor neurons, which is required for the maturation of motor nerve terminal. Thus, Amotl2 and Ptprk exert different roles at the NM compared to Fam53b and Tcf7l2. Additionally, Wnt ligands originating from the spinal motor neurons and the perichondrium/chondrocyte are likely to work remotely on the subsynaptic nuclei and the myotendinous junctional nuclei, respectively. We conclude that snRNA-seq analysis of embryonic/neonatal diaphragms reveal a novel coordinated expression profile especially in the Wnt/ß-catenin signaling that regulate the formation of the embryonic NMJ.

Metagenomic Thermometer.

Various microorganisms exist in environments, and each of them has its optimal growth temperature (OGT). The relationship between genomic information and OGT of each species has long been studied, and one such study revealed that OGT of prokaryotes can be accurately predicted based on the fraction of seven amino acids (IVYWREL) among all encoded amino-acid sequences in its genome. Extending this discovery, we developed a 'Metagenomic Thermometer' as a means of predicting environmental temperature based on metagenomic sequences. Temperature prediction of diverse environments using publicly available metagenomic data revealed that the Metagenomic Thermometer can predict environmental temperatures with small temperature changes and little influx of microorganisms from other environments. The accuracy of the Metagenomic Thermometer was also confirmed by a demonstration experiment using an artificial hot water canal. The Metagenomic Thermometer was also applied to human gut metagenomic samples, yielding a reasonably accurate value for human body temperature. The result further suggests that deep body temperature determines the dominant lineage of the gut community. Metagenomic Thermometer provides a new insight into temperature-driven community assembly based on amino-acid composition rather than microbial taxa.

Diversity and compositional differences of the airborne microbiome in a biophilic indoor environment.

Biophilic design based on indoor planting plays an important role in human physical and mental well-being. To investigate and assess the effects of indoor planting on air quality, we sequenced 16S rRNA gene amplicons to compare the airborne bacterial microbiomes of three planting rooms before and after installing natural materials (plants, soil, water, etc.) with distinct biophilic attributes. Incorporation of indoor plantings significantly increased the taxonomic diversity of the airborne microbiome in each room, and we observed different microbiome compositions in each room. The proportional contribution of each bacterial source to the airborne microbiome in the indoor planting rooms was estimated by SourceTracker2. This analysis revealed that the proportion of airborne microbial sources (e.g., plants and soil) varied depending on the natural materials installed. Our results have important implications for indoor planting with biophilic design to control the indoor airborne microbiome.

Condensed but liquid-like domain organization of active chromatin regions in living human cells.

In eukaryotes, higher-order chromatin organization is spatiotemporally regulated as domains, for various cellular functions. However, their physical nature in living cells remains unclear (e.g., condensed domains or extended fiber loops; liquid-like or solid-like). Using novel approaches combining genomics, single-nucleosome imaging, and computational modeling, we investigated the physical organization and behavior of early DNA replicated regions in human cells, which correspond to Hi-C contact domains with active chromatin marks. Motion correlation analysis of two neighbor nucleosomes shows that nucleosomes form physically condensed domains with ~150-nm diameters, even in active chromatin regions. The mean-square displacement analysis between two neighbor nucleosomes demonstrates that nucleosomes behave like a liquid in the condensed domain on the ~150 nm/~0.5 s spatiotemporal scale, which facilitates chromatin accessibility. Beyond the micrometers/minutes scale, chromatin seems solid-like, which may contribute to maintaining genome integrity. Our study reveals the viscoelastic principle of the chromatin polymer; chromatin is locally dynamic and reactive but globally stable.

Integrated gut microbiome and metabolome analyses identified fecal biomarkers for bowel movement regulation by Bifidobacterium longum BB536 supplementation: A RCT.

Background: Bifidobacterium longum BB536 supplementation can be used to regulate bowel movements in various people, including healthy subjects and patients with irritable bowel syndrome (IBS); however, individuals vary in their responses to B. longum BB536 treatment. One putative factor is the gut microbiota; recent studies have reported that the gut microbiota mediates the effects of diet or drugs on the host. Here, we investigated intestinal features, such as the microbiome and metabolome, related to B. longum BB536 effectiveness in increasing bowel movement frequency. Results: A randomized, double-blind controlled crossover trial was conducted with 24 adults who mainly tended to be constipated. The subjects received a two-week dietary intervention consisting of B. longum BB536 in acid-resistant seamless capsules or similarly encapsulated starch powder as the placebo control. Bowel movement frequency was recorded daily, and fecal samples were collected at several time points, and analyzed by metabologenomic approach that consists of an integrated analysis of metabolome data obtained using mass spectrometry and microbiome data obtained using high-throughput sequencing. There were differences among subjects in B. longum intake-induced bowel movement frequency. The responders were predicted by machine learning based on the microbiome and metabolome features of the fecal samples collected before B. longum intake. The abundances of eight bacterial genera were significantly different between responders and nonresponders. Conclusions: Intestinal microbiome and metabolome profiles might be utilized as potential markers of improved bowel movement after B. longum BB536 supplementation. These findings have implications for the development of personalized probiotic treatments.

Leaf Bleaching in Rice: A New Disease in Vietnam Caused by Methylobacterium indicum, Its Genomic Characterization and the Development of a Suitable Detection Technique.

A new disease in rice that is characterized by leaf bleaching was recently identified in some fields in the Mekong Delta region of Vietnam. The present study was the first to isolate and identify the pathogen of this disease. We confirmed that leaf bleaching symptoms were due to infection with Methylobacterium indicum bacteria using molecular biology approaches. A full-length genome analysis of pathogenic Methylobacterium strain VL1 revealed that it comprises a single chromosome and six plasmids, with a total size of 7.05| |Mbp and GC content of 70.5%. The genomic features of VL1 were similar to those of the non-pathogenic M. indicum strain SE2.11T; however, VL1 possessed additional unique genes, including those related to homoserine lactone biosynthesis. We established a loop-mediated isothermal amplification (LAMP) assay using the unique sequences of VL1 as target sequences for the rapid and simple detection of pathogenic M. indicum strains. Our initial evaluation demonstrated that the LAMP assay successfully distinguished between pathogenic and non-pathogenic strains infecting rice plants in a rapid and sensitive manner. The present results provide insights into the pathogenesis and development of control measures for novel rice diseases.

PZLAST: an ultra-fast amino acid sequence similarity search server against public metagenomes.

SUMMARY: : Similarity searches of amino acid sequences against the public metagenomic data can provide users insights about the function of sequences based on the environmental distribution of similar sequences. However, a considerable reduction in the amount of data or the accuracy of the result was necessary to conduct sequence similarity searches against public metagenomic data, because of the vast data size more than Terabytes. Here, we present an ultra-fast service for the highly accurate amino acid sequence similarity search, called PZLAST, which can search the user's amino acid sequences to several Terabytes of public metagenomic sequences in ∼10-20 min. PZLAST accomplishes its search speed by using PEZY-SC2, which is a Multiple Instruction Multiple Data many-core processor. Results of PZLAST are summarized by the ontology-based environmental distribution of similar sequences. PZLAST can be used to predict the function of sequences and mine for homologs of functionally important gene sequences. AVAILABILITY AND IMPLEMENTATION: PZLAST is freely accessible at https://pzlast.riken.jp/meta without requiring registration. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Key bacterial taxa and metabolic pathways affecting gut short-chain fatty acid profiles in early life.

Infant gut microbiota development affects the host physiology throughout life, and short-chain fatty acids (SCFAs) are promising key metabolites mediating microbiota-host relationships. Here, we investigated dense longitudinally collected faecal samples from 12 subjects during the first 2 years (n = 1048) to identify early life gut SCFA patterns and their relationships with the microbiota. Our results revealed three distinct phases of progression in the SCFA profiles: early phase characterised by low acetate and high succinate, middle-phase characterised by high lactate and formate and late-phase characterised by high propionate and butyrate. Assessment of the SCFA-microbiota relationships revealed that faecal butyrate is associated with increased Clostridiales and breastfeeding cessation, and that diverse and personalised assemblage of Clostridiales species possessing the acetyl-CoA pathway play major roles in gut butyrate production. We also found an association between gut formate and some infant-type bifidobacterial species, and that human milk oligosaccharides (HMO)-derived fucose is the substrate for formate production during breastfeeding. We identified genes upregulated in fucose and fucosylated HMO utilisation in infant-type bifidobacteria. Notably, bifidobacteria showed interspecific and intraspecific variation in the gene repertoires, and cross-feeding of fucose contributed to gut formate production. This study provides an insight into early life SCFA-microbiota relationships, which is an important step for developing strategies for modulating lifelong health.

Aberrant mucosal immunoreaction to tonsillar microbiota in immunoglobulin A nephropathy.

BACKGROUND: Immunoglobulin A nephropathy (IgAN) is the most common glomerulonephritis worldwide, characterized by mesangial polymeric IgA1 deposition. IgAN is believed to develop owing to aberrant mucosal immunoreaction against commensals in the tonsils. However, the exact interrelation between pathogenic IgA and mucosal microbiota in IgAN patients is unclear. METHODS: Biopsy-proven IgAN or recurrent tonsillitis (RT) patients who had undergone tonsillectomy were enrolled. We used 16S ribosomal RNA gene amplicon sequencing with a flow cytometry-based bacterial cell sorting technique) and immunoglobulin repertoire sequencing of the IgA heavy chain to characterize IgA-coated bacteria of the tonsillar microbiota (IgA-SEQ) and their corresponding IgA repertoire. Furthermore, we fractionated patient serum using gel-filtration chromatography and performed flow cytometry-based analysis of IgA binding to bacteria cultured from incised tonsils. RESULTS: Tonsillar proliferation-inducing ligand and B-cell activating factor levels were significantly higher in IgAN than in RT patients. IgA-SEQ for tonsillar microbiota revealed the preferential binding ability of IgA to Bacteroidetes in IgAN tonsils compared with those from RT patients. Expression of immunoglobulin heavy (IGH) constant alpha 1 with IGH variable 3-30 was significantly higher in IgAN than that in RT, and positively correlated with the IgA-coated enrichment score of Bacteroidetes. Serum polymeric IgA, comprising high levels of GdIgA1, exhibited considerable binding to Bacteroidetes strains cultured from the tonsils of IgAN patients. CONCLUSIONS: These findings provide evidence that aberrant mucosal immune responses to tonsillar anaerobic microbiota, primarily consisting of members of the phylum Bacteroidetes, are involved in IgAN pathophysiology.

RNA-binding motifs of hnRNP K are critical for induction of antibody diversification by activation-induced cytidine deaminase.

Activation-induced cytidine deaminase (AID) is the key enzyme for class switch recombination (CSR) and somatic hypermutation (SHM) to generate antibody memory. Previously, heterogeneous nuclear ribonucleoprotein K (hnRNP K) was shown to be required for AID-dependent DNA breaks. Here, we defined the function of major RNA-binding motifs of hnRNP K, GXXGs and RGGs in the K-homology (KH) and the K-protein-interaction (KI) domains, respectively. Mutation of GXXG, RGG, or both impaired CSR, SHM, and cMyc/IgH translocation equally, showing that these motifs were necessary for AID-dependent DNA breaks. AID-hnRNP K interaction is dependent on RNA; hence, mutation of these RNA-binding motifs abolished the interaction with AID, as expected. Some of the polypyrimidine sequence-carrying prototypical hnRNP K-binding RNAs, which participate in DNA breaks or repair bound to hnRNP K in a GXXG and RGG motif-dependent manner. Mutation of the GXXG and RGG motifs decreased nuclear retention of hnRNP K. Together with the previous finding that nuclear localization of AID is necessary for its function, lower nuclear retention of these mutants may worsen their functional deficiency, which is also caused by their decreased RNA-binding capacity. In summary, hnRNP K contributed to AID-dependent DNA breaks with all of its major RNA-binding motifs.

Latent environment allocation of microbial community data.

As data for microbial community structures found in various environments has increased, studies have examined the relationship between environmental labels given to retrieved microbial samples and their community structures. However, because environments continuously change over time and space, mixed states of some environments and its effects on community formation should be considered, instead of evaluating effects of discrete environmental categories. Here we applied a hierarchical Bayesian model to paired datasets containing more than 30,000 samples of microbial community structures and sample description documents. From the training results, we extracted latent environmental topics that associate co-occurring microbes with co-occurring word sets among samples. Topics are the core elements of environmental mixtures and the visualization of topic-based samples clarifies the connections of various environments. Based on the model training results, we developed a web application, LEA (Latent Environment Allocation), which provides the way to evaluate typicality and heterogeneity of microbial communities in newly obtained samples without confining environmental categories to be compared. Because topics link words and microbes, LEA also enables to search samples semantically related to the query out of 30,000 microbiome samples.

The Relationship Between Microbial Community Structures and Environmental Parameters Revealed by Metagenomic Analysis of Hot Spring Water in the Kirishima Area, Japan.

Diverse microorganisms specifically inhabit extreme environments, such as hot springs and deep-sea hydrothermal vents. To test the hypothesis that the microbial community structure is predictable based on environmental factors characteristic of such extreme environments, we conducted correlation analyses of microbial taxa/functions and environmental factors using metagenomic and 61 types of physicochemical data of water samples from nine hot springs in the Kirishima area (Kyusyu, Japan), where hot springs with diverse chemical properties are distributed in a relatively narrow area. Our metagenomic analysis revealed that the samples can be classified into two major types dominated by either phylum Crenarchaeota or phylum Aquificae. The correlation analysis showed that Crenarchaeota dominated in nutrient-rich environments with high concentrations of ions and total carbons, whereas Aquificae dominated in nutrient-poor environments with low ion concentrations. These environmental factors were also important explanatory variables in the generalized linear models constructed to predict the abundances of Crenarchaeota or Aquificae. Functional enrichment analysis of genes also revealed that the separation of the two major types is primarily attributable to genes involved in autotrophic carbon fixation, sulfate metabolism and nitrate reduction. Our results suggested that Aquificae and Crenarchaeota play a vital role in the Kirishima hot spring water ecosystem through their metabolic pathways adapted to each environment. Our findings provide a basis to predict microbial community structures in hot springs from environmental parameters, and also provide clues for the exploration of biological resources in extreme environments.

Comprehensive microbiome analysis of tonsillar crypts in IgA nephropathy.

BACKGROUND: Immunoglobulin A nephropathy (IgAN) is the most prevalent primary chronic glomerular disease, in which the mucosal immune response elicited particularly in the tonsils or intestine has been estimated to be involved in the development of the disease. To explore the relationship between IgAN and bacterial flora in the tonsils, we conducted a comprehensive microbiome analysis. METHODS: We enrolled 48 IgAN patients, 21 recurrent tonsillitis (RT) patients without urine abnormalities and 30 children with tonsillar hyperplasia (TH) who had undergone tonsillectomy previously. Genomic DNA from tonsillar crypts of each patient was extracted, and V4 regions of the 16S ribosomal RNA gene were amplified and analysed using a high-throughput multiplexed sequencing approach. Differences in genus composition among the three study groups were statistically analysed by permutational multivariate analysis of variance and visualized by principal component analysis (PCA). RESULTS: Substantial diversity in bacterial composition was detected in each sample. Prevotella spp., Fusobacterium spp., Sphingomonas spp. and Treponema spp. were predominant in IgAN patients. The percentage of abundance of Prevotella spp., Haemophilus spp., Porphyromonas spp. and Treponema spp. in IgAN patients was significantly different from that in TH patients. However, there was no significant difference in the percentage of abundance of any bacterial genus between IgAN and RT patients. PCA did not distinguish IgAN from RT, although it discriminated TH. No significant differences in microbiome composition among the groups of IgAN patients according to clinicopathological parameters were observed. CONCLUSIONS: Similar patterns of bacteria are present in tonsillar crypts of both IgAN and RT patients, suggesting that the host response to these bacteria might be important in the development of IgAN.

H-NS Facilitates Sequence Diversification of Horizontally Transferred DNAs during Their Integration in Host Chromosomes.

Bacteria can acquire new traits through horizontal gene transfer. Inappropriate expression of transferred genes, however, can disrupt the physiology of the host bacteria. To reduce this risk, Escherichia coli expresses the nucleoid-associated protein, H-NS, which preferentially binds to horizontally transferred genes to control their expression. Once expression is optimized, the horizontally transferred genes may actually contribute to E. coli survival in new habitats. Therefore, we investigated whether and how H-NS contributes to this optimization process. A comparison of H-NS binding profiles on common chromosomal segments of three E. coli strains belonging to different phylogenetic groups indicated that the positions of H-NS-bound regions have been conserved in E. coli strains. The sequences of the H-NS-bound regions appear to have diverged more so than H-NS-unbound regions only when H-NS-bound regions are located upstream or in coding regions of genes. Because these regions generally contain regulatory elements for gene expression, sequence divergence in these regions may be associated with alteration of gene expression. Indeed, nucleotide substitutions in H-NS-bound regions of the ybdO promoter and coding regions have diversified the potential for H-NS-independent negative regulation among E. coli strains. The ybdO expression in these strains was still negatively regulated by H-NS, which reduced the effect of H-NS-independent regulation under normal growth conditions. Hence, we propose that, during E. coli evolution, the conservation of H-NS binding sites resulted in the diversification of the regulation of horizontally transferred genes, which may have facilitated E. coli adaptation to new ecological niches.