Genome-wide Distribution Research Articles

High-throughput capture of transcription factor-driven epigenome dynamics using PHILO ChIP-seq.

Assessing the dynamics of chromatin features and transcription factor (TF) binding at scale remains a significant challenge in plants. Here, we present PHILO (Plant HIgh-throughput LOw input) ChIP-seq, a high-throughput ChIP-seq platform that enables the cost-effective and extensive capture of TF binding and genome-wide distributions of histone modifications. The PHILO ChIP-seq pipeline is adaptable to many plant species, requires very little starting material (1mg), and provides the option to use MNase (micrococcal nuclease) for chromatin fragmentation. By employing H3K9ac PHILO ChIP-seq on eight Arabidopsis thaliana jasmonic acid (JA) pathway mutants, with the simultaneous processing of over 100 samples, we not only recapitulated but also expanded the current understanding of the intricate interplay between the master TFs MYC2/3/4 and various chromatin regulators. Additionally, our analyses brought to light previously unknown histone acetylation patterns within the regulatory regions of MYC2 target genes in Arabidopsis, which is also conserved in tomato (Solanum lycopersicum). In summary, our PHILO ChIP-seq platform demonstrates its high effectiveness in investigating TF binding and chromatin dynamics on a large scale in plants, paving the way for the cost-efficient realization of complex experimental setups.

MalKinID: A Classification Model for Identifying Malaria Parasite Genealogical Relationships Using Identity-by-Descent.

Pathogen genomics is a powerful tool for tracking infectious disease transmission. In malaria, identity-by-descent (IBD) is used to assess the genetic relatedness between parasites and has been used to study transmission and importation. In theory, IBD can be used to distinguish genealogical relationships to reconstruct transmission history or identify parasites for quantitative-trait-locus experiments. MalKinID (Malaria Kinship Identifier) is a new classification model designed to identify genealogical relationships among malaria parasites based on genome-wide IBD proportions and IBD segment distributions. MalKinID was calibrated to the genomic data from three laboratory-based genetic crosses (yielding 440 parent-child [PC] and 9060 full-sibling [FS] comparisons). MalKinID identified lab generated F1 progeny with >80% sensitivity and showed that 0.39 (95% CI 0.28, 0.49) of the second-generation progeny of a NF54 and NHP4026 cross were F1s and 0.56 (0.45, 0.67) were backcrosses of an F1 with the parental NF54 strain. In simulated outcrossed importations, MalKinID reconstructs genealogy history with high precision and sensitivity, with F1-scores exceeding 0.84. However, when importation involves inbreeding, such as during serial co-transmission, the precision and sensitivity of MalKinID declined, with F1-scores (the harmonic mean of precision and sensitivity) of 0.76 (0.56, 0.92) and 0.23 (0.0, 0.4) for PC and FS and <0.05 for second-degree and third-degree relatives. Disentangling inbred relationships required adapting MalKinID to perform multi-sample comparisons. Genealogical inference is most powered when 1) outcrossing is the norm or 2) multi-sample comparisons based on a predefined pedigree are used. MalKinID lays the foundations for using IBD to track parasite transmission history and for separating progeny for quantitative-trait-locus experiments.

ChIPmentation for Epigenomic Analysis in Fission Yeast.

Histone modifications and transcription factor-DNA interactions regulate vital processes such as transcription, recombination, repair, and accurate chromosome segregation. Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) has been instrumental in studying genome-wide distribution of DNA-bound or chromatin-associated factors and histone posttranslational modifications (PTMs). Here, we describe a ChIPmentation protocol adapted for fission yeast, Schizosaccharomyces pombe. This method merges Tn5 mediated tagmentation with existing ChIP protocols, resulting in lower sample input requirements with significant reduction in hands-on time and sample preparation costs.

Dynamic atlas of histone modifications and gene regulatory networks in endosperm of bread wheat.

Dissecting the genetic basis of seed traits in wheat is impeded by limited genetic polymorphisms and significant variations caused by environmental conditions and seed position in a spikelet. Seed performance is largely determined by endosperm development controlled by spatiotemporal variation in gene activities, which is greatly affected by chromatin status. Here, we map genome-wide dynamic distributions of H3K27me3, H3K4me3 and H3K9ac modifications and profile gene transcription across wheat endosperm development. The combinatorial effects of active and repressive marks ensure spatiotemporal dynamic gene expression, especially for starch biosynthesis. By scanning the transcription factor binding motifs in the ATAC-seq peaks, hub regulators are identified from the regulatory network. In addition, we observe significant correlations between sequence polymorphisms of hub regulators and variations in seed traits in a germplasm population. Thus, the analysis of genomic regulatory activities together with genetic variation provides a robust approach to dissect seed traits in bread wheat.

Heterochromatin dynamics during the initial stages of sexual development in Plasmodium falciparum

Asexual replication of Plasmodium falciparum in the human blood results in exponential parasite growth and causes all clinical symptoms of malaria. However, at each round of the replicative cycle, some parasites convert into sexual precursors called gametocytes, which develop through different stages until they become infective to mosquito vectors. The genome-wide distribution of heterochromatin, a type of chromatin generally refractory to gene expression, is identical at all asexual blood stages, but is altered in stage II/III and more mature gametocytes. However, it is not known if these changes occur concomitantly with sexual conversion or at a later time during gametocyte development. Using a transgenic line in which massive sexual conversion can be conditionally induced, we show that the genome-wide distribution of heterochromatin at the initial stages of sexual development (i.e., sexual rings and stage I gametocytes) is almost identical to asexual blood stages, and major changes do not occur until stage II/III. However, we found that at loci with heterochromatin alterations, transcriptional changes associated with sexual development typically precede, rather than follow, changes in heterochromatin occupancy.

Genome-wide profiling the epigenetic landscape of histone variant TH2B in murine oocytes and pre-implantation embryos.

The histone variant TH2B, enriched in oocytes, sperm, and early embryos, decreases as embryos differentiate into pre-gastrula stages. Despite its presence, the role of TH2B in epigenetic reprogramming during early embryonic development remains largely under-researched. Our study employed ultra-low-input ChIP-seq (ULI-ChIP) to analyze the genome-wide distribution of TH2B in MII oocytes and early embryos. We found that TH2B is enriched in the chromatin of oocytes and 2-cell stage embryos but becomes less prevalent after the 2-cell stage. Correlation analysis revealed that the TH2B chromatin patterns in sperm and preimplantation embryos are more similar to each other than to those in MII oocytes. Gene ontology (GO) analysis of TH2B-occupied loci linked them to various developmental processes, including oogenesis, fertilization, chromatin modification, and transcription regulation. The study also identified a strong association of TH2B with specific transposable elements (TEs), particularly long terminal repeats (LTRs), which are known to regulate preimplantation development. Additionally, early embryos showed H3K9me3 marks at TH2B-bound loci. TH2B exhibited strong correlations with H2A.Z and H3.3 in the 2-cell and 8-cell stages, a positive association with H3K27Ac and H3K4me3, and a negative correlation with H3K27me3. Allelic reprogramming analysis of TH2B in embryos from C57BL/6J and DBA/2J crosses revealed differential dynamics between maternal and paternal alleles, with a notable paternal bias at the promoter in 2-cell embryos. Thus, TH2B's enrichment in early embryonic stages and its association with key regulatory regions and histone modifications underscore its importance in ZGA and subsequent developmental processes.

Untargeted CUT&Tag and BG4 CUT&Tag are both enriched at G-quadruplexes and accessible chromatin.

G-quadruplex DNA structures (G4s) form within single-stranded DNA in nucleosome-free chromatin. As G4s modulate gene expression and genomic stability, genome-wide mapping of G4s has generated strong research interest. Recently, the Cleavage Under Targets and Tagmentation (CUT&Tag) method was performed with the G4-specific BG4 antibody to target Tn5 transposase to G4s. While this method generated a novel high-resolution map of G4s, we unexpectedly observed a strong correlation between the genome-wide signal distribution of BG4 CUT&Tag and accessible chromatin. To examine whether untargeted Tn5 cutting at accessible chromatin contributes to BG4 CUT&Tag signal, we examined the genome-wide distribution of signal from untargeted (i.e. negative control) CUT&Tag datasets. We observed that untargeted CUT&Tag signal distribution was highly similar to both that of accessible chromatin and of BG4 CUT&Tag. We also observed that BG4 CUT&Tag signal increased at mapped G4s, but this increase was accompanied by a concomitant increase in untargeted CUT&Tag at the same loci. Consequently, enrichment of BG4 CUT&Tag over untargeted CUT&Tag was not increased at mapped G4s. These results imply that either the vast majority of accessible chromatin regions contain mappable G4s or that the presence of G4s within accessible chromatin cannot reliably be determined using BG4 CUT&Tag alone.

H3K4me3 changes occur in cell wall genes during the development of Fagopyrum tataricum morphogenic and non-morphogenic calli.

Epigenetic changes accompany the dynamic changes in the cell wall composition during the development of callus cells. H3K4me3 is responsible for active gene expression and reaction to environmental cues. Chromatin immunoprecipitation (ChIP) is a powerful technique for studying the interplay between epigenetic modifications and the DNA regions of interest. In combination with sequencing, it can provide the genome-wide enrichment of the specific epigenetic mark, providing vital information on its involvement in the plethora of cellular processes. Here, we describe the genome-wide distribution of H3K4me3 in morphogenic and non-morphogenic callus of Fagopyrum tataricum. Levels of H3K4me3 were higher around the transcription start site, in agreement with the role of this mark in transcriptional activation. The global levels of methylation were higher in the non-morphogenic callus, which indicated increased gene activation compared to the morphogenic callus. We also employed ChIP to analyse the changes in the enrichment of this epigenetic mark on the cell wall-related genes in both calli types during the course of the passage. Enrichment of H3K4me3 on cell wall genes was specific for callus type, suggesting that the role of this mark in cell-wall remodelling is complex and involved in many processes related to dedifferentiation and redifferentiation. This intricacy of the cell wall composition was supported by the immunohistochemical analysis of the cell wall epitopes' distribution of pectins and extensins. Together, these data give a novel insight into the involvement of H3K4me3 in the regeneration processes in F. tataricum in vitro callus tissue culture.

Chromosome-level genome assembly of the morabine grasshopper Vandiemenella viatica19

Morabine grasshoppers in the Vandiemenella viatica species group, which show karyotype diversity, have been studied for their ecological distribution and speciation in relation to their genetic and chromosomal diversity. They are good models for studying sex chromosome evolution as “old” and newly emerged sex chromosomes co-exist within the group. Here we present a reference genome for the viatica19 chromosomal race, that possesses the ancestral karyotype within the group. Using PacBio HiFi and Hi-C sequencing, we generated a chromosome-level assembly of 4.09 Gb in span, scaffold N50 of 429 Mb, and complete BUSCO score of 98.1%, containing 10 pseudo-chromosomes. We provide Illumina datasets of males and females, used to identify the X chromosome. The assembly contains 19,034 predicted protein-coding genes, and a total of 75.21% of repetitive DNA sequences. By leveraging HiFi reads, we mapped the genome-wide distribution of methylated bases (5mC and 6 mA). This comprehensive assembly offers a robust reference for morabine grasshoppers and supports further research into speciation and sex chromosome diversification within the group and its related species.

5-Hydroxymethylcytosines in circulating cell-free DNA as a diagnostic biomarker for nasopharyngeal carcinoma

ObjectiveTo evaluate the diagnostic value of 5-hydroxymethylcytosines (5hmC) in circulating cell-free DNA (cfDNA) for nasopharyngeal carcinoma (NPC) and to develop a diagnostic model. MethodsGenome-wide 5hmC profiles in cfDNA from 174 NPC patients and 146 non-cancer individuals were analyzed using the 5hmC-Seal technique. A cfDNA 5hmC-based diagnostic model to identify NPC patients was developed using least absolute shrinkage and selection operator (LASSO) logistic regression, and performance was evaluated with receiver operating characteristic (ROC) curves and confusion matrices. ResultsThe 5hmC-Seal data from patients with NPC showed a different genome-wide distribution than non-tumor samples. Our initial analysis revealed a 12-gene-based 5hmC marker panel to be an accurate diagnostic model effectively distinguishing between NPC samples and non-cancerous samples (training set: area under curve (AUC)= 0.97 [95 % CI: 0.94–0.99]; and test set: AUC= 0.93 [95 % CI: 0.88–0.98]) superior to EBV DNA testing. The diagnostic score performed well in differentiating the non-cancer subjects from early-stage NPC (training set: AUC=0.99 [95 % CI: 0.98–1]; test set: AUC=0.98 [95 % CI: 0.95–1]), and advanced-stage NPC (training set: AUC=0.96 [95 % CI: 0.93–0.99]; test set: AUC=0.93 [95 % CI: 0.88–0.98]). Notably, in EBV-negative patients, the diagnostic scores showed excellent capacity for distinguishing EBV-negative patients with NPC from non-cancer subjects in both the training set (AUC= 0.94 [95 % CI: 0.88–1]) and test set (AUC=0.91 [95 % CI: 0.81–1]). Conclusion5hmC modifications in cfDNA are promising noninvasive biomarkers for NPC, offering high sensitivity and specificity, particularly for early-stage and EBV-negative NPC.

Calibration of variant effect predictors on genome-wide data masks heterogeneous performance across genes

In silico variant effect predictions are available for nearly all missense variants but played a minimal role in clinical variant classification because they were deemed to provide only supporting evidence. Recently, the ClinGen Sequence Variant Interpretation (SVI) Working Group updated recommendations for variant effect prediction use. By analyzing control pathogenic and benign variants across all genes, they were able to compute evidence strength for predictor score intervals with some intervals generating moderate, strong, or even very strong evidence. However, this genome-wide approach could obscure heterogeneous predictor performance in different genes. We quantified the gene-by-gene performance of two top predictors, REVEL and BayesDel, by analyzing control variants in each predictor score interval in 3,668 disease-relevant genes. Approximately 10% of intervals had sufficient control variants for analysis, and ∼70% of these intervals exceeded the maximum number of incorrect predictions implied by the SVI recommendations. These trending discordant intervals arose owing to the divergence of the gene-specific distribution of predictions from the genome-wide distribution, suggesting that gene-specific calibration is needed in many cases. Approximately 22% of ClinVar missense variants of uncertain significance in genes we analyzed (REVEL= 100,629, BayesDel= 71,928) had predictions in trending discordant intervals. Thus, genome-wide calibrations could result in many variants receiving inappropriate evidence strength. To facilitate a review of the SVI's calibrations, we developed a web application enabling visualization of gene-specific predictions and trending concordant and discordant intervals.

ERCC2 mutations alter the genomic distribution pattern of somatic mutations and are independently prognostic in bladder cancer

Excision repair cross-complementation group 2 (ERCC2) encodes the DNA helicase xeroderma pigmentosum group D, which functions in transcription and nucleotide excision repair. Point mutations in ERCC2 are putative drivers in around 10% of bladder cancers (BLCAs) and a potential positive biomarker for cisplatin therapy response. Nevertheless, the prognostic significance directly attributed to ERCC2 mutations and its pathogenic role in genome instability remain poorly understood. We first demonstrated that mutant ERCC2 is an independent predictor of prognosis in BLCA. We then examined its impact on the somatic mutational landscape using a cohort of ERCC2 wild-type (n= 343) and mutant (n= 39) BLCA whole genomes. The genome-wide distribution of somatic mutations is significantly altered in ERCC2 mutants, including T[C>T]N enrichment, altered replication time correlations, and CTCF-cohesin binding site mutation hotspots. We leverage these alterations to develop a machine learning model for predicting pathogenic ERCC2 mutations, which may be useful to inform treatment of patients with BLCA.

Analysis of genomic copy number variations through whole-genome scan in Yunling cattle.

Yunling cattle is a new breed of beef cattle bred in Yunnan Province, China, which has the advantages of fast growth, excellent meat quality, improved tolerance ability, and important landscape value. Copy number variation (CNV) is a significant source of gene structural variation and plays a crucial role in evolution and phenotypic diversity. Based on the latest reference genome ARS-UCD2.0, this study analyzed the genome-wide distribution of CNVs in Yunling cattle using short-read whole-genome sequencing data (n = 129) and single-molecule long-read sequencing data (n = 1), and a total of 16,507 CNVs were detected. After merging CNVs with overlapping genomic positions, 3,728 CNV regions (CNVRs) were obtained, accounting for 0.61% of the reference genome. The functional analysis indicated significant enrichment of CNVRs in 96 GO terms and 57 KEGG pathways, primarily related to cell adhesion, signal transduction, neuromodulation, and nutritional metabolism. Additionally, 111 CNVRs overlapped with 76 quantitative trait loci (QTLs), including Subcutaneous fat thickness QTL, Longissimus muscle area QTL, and Marbling score QTL. Several CNVR-overlapping genes, including BZW1, AOX1, and LOC100138449, overlap with regions associated with meat color and quality QTLs. Furthermore, Vst analysis showed that PSMB4, ERICH1, SMC2, and PPP4R3A were highly divergent between Yunling and Brahman cattle. In summary, we have constructed the genomic CNV map of Yunling cattle for the first time using whole-genome resequencing. This provides valuable genetic variation resources for the study of the Yunling cattle genome and contributes to the study of economic traits in Yunling cattle.

MalKinID: A Likelihood-Based Model for Identifying Malaria Parasite Genealogical Relationships Using Identity-by-Descent.

Pathogen genomics is a powerful tool for tracking infectious disease transmission. In malaria, identity-by-descent (IBD) is used to assess the genetic relatedness between parasites and has been used to study transmission and importation. In theory, IBD can be used to distinguish genealogical relationships to reconstruct transmission history or identify parasites for genotype-to-phenotype quantitative-trait-locus experiments. MalKinID (Malaria Kinship Identifier) is a new likelihood-based classification model designed to identify genealogical relationships among malaria parasites based on genome-wide IBD proportions and IBD segment distributions. MalKinID was calibrated to the genomic data from three laboratory-based genetic crosses (yielding 440 parent-child and 9060 full-sibling comparisons). MalKinID identified lab generated F1 progeny with >80% sensitivity and showed that 0.39 (95% CI 0.28, 0.49) of the second-generation progeny of a NF54 and NHP4026 cross were F1s and 0.56 (0.45, 0.67) were backcrosses of an F1 with the parental NF54 strain. In simulated outcrossed importations, MalKinID accurately reconstructs genealogy history with high precision and sensitivity, with F1-scores exceeding 0.84. However, when importation involves inbreeding, such as during serial co-transmission, the precision and sensitivity of MalKinID declined, with F1-scores of 0.76 (0.56, 0.92) and 0.23 (0.0, 0.4) for PC and FS and <0.05 for second-degree and third-degree relatives. Genealogical inference is most powered 1) when outcrossing is the norm or 2) when multi-sample comparisons based on a predefined pedigree are used. MalKinID lays the foundations for using IBD to track parasite transmission history and for separating progeny for quantitative-trait-locus experiments.

Metapresence: a tool for accurate species detection in metagenomics based on the genome-wide distribution of mapping reads.

Shotgun metagenomics allows comprehensive sampling of the genomic information of microbes in a given environment and is a tool of choice for studying complex microbial systems. Mapping sequencing reads against a set of reference or metagenome-assembled genomes is in principle a simple and powerful approach to define the species-level composition of the microbial community under investigation. However, despite the widespread use of this approach, there is no established way to properly interpret the alignment results, with arbitrary relative abundance thresholds being routinely used to discriminate between present and absent species. Such an approach can be affected by significant biases, especially in the identification of rare species. Therefore, it is important to develop new metrics to overcome these biases. Here, we present Metapresence, a new tool to perform reliable identification of the species in metagenomic samples based on the distribution of mapped reads on the reference genomes. The analysis is based on two metrics describing the breadth of coverage and the genomic distance between consecutive reads. We demonstrate the high precision and wide applicability of the tool using data from various synthetic communities, a real mock community, and the gut microbiome of healthy individuals and antibiotic-associated-diarrhea patients. Overall, our results suggest that the proposed approach has a robust performance in hard-to-analyze microbial communities containing contaminated or closely related genomes in low abundance.IMPORTANCEDespite the prevalent use of genome-centric alignment-based methods to characterize microbial community composition, there lacks a standardized approach for accurately identifying the species within a sample. Currently, arbitrary relative abundance thresholds are commonly employed for this purpose. However, due to the inherent complexity of genome structure and biases associated with genome-centric approaches, this practice tends to be imprecise. Notably, it introduces significant biases, particularly in the identification of rare species. The method presented here addresses these limitations and contributes significantly to overcoming inaccuracies in precisely defining community composition, especially when dealing with rare members.

The genome-wide signature of short-term temporal selection

Despite evolutionary biology's obsession with natural selection, few studies have evaluated multigenerational series of patterns of selection on a genome-wide scale in natural populations. Here, we report on a 10-y population-genomic survey of the microcrustacean Daphnia pulex. The genome sequences of [Formula: see text]800 isolates provide insights into patterns of selection that cannot be obtained from long-term molecular-evolution studies, including the following: the pervasiveness of near quasi-neutrality across the genome (mean net selection coefficients near zero, but with significant temporal variance about the mean, and little evidence of positive covariance of selection across time intervals); the preponderance of weak positive selection operating on minor alleles; and a genome-wide distribution of numerous small linkage islands of observable selection influencing levels of nucleotide diversity. These results suggest that interannual fluctuating selection is a major determinant of standing levels of variation in natural populations, challenge the conventional paradigm for interpreting patterns of nucleotide diversity and divergence, and motivate the need for the further development of theoretical expressions for the interpretation of population-genomic data.

Genome-wide Methylation Dynamics and Context-dependent Gene Expression Variability in Differentiating Preadipocytes.

Obesity, characterized by the accumulation of excess fat, is a complex condition resulting from the combination of genetic and epigenetic factors. Recent studies have found correspondence between DNA methylation and cell differentiation, suggesting a role of the former in cell fate determination. There is a lack of comprehensive understanding concerning the underpinnings of preadipocyte differentiation, specifically when cells are undergoing terminal differentiation (TD). To gain insight into dynamic genome-wide methylation, 3T3 L1 preadipocyte cells were differentiated by a hormone cocktail. The genomic DNA was isolated from undifferentiated cells and 4 hours, 2 days postdifferentiated cells, and 15 days TD cells. We employed whole-genome bisulfite sequencing (WGBS) to ascertain global genomic DNA methylation alterations at single base resolution as preadipocyte cells differentiate. The genome-wide distribution of DNA methylation showed similar overall patterns in pre-, post-, and terminally differentiated adipocytes, according to WGBS analysis. DNA methylation decreases at 4 hours after differentiation initiation, followed by methylation gain as cells approach TD. Studies revealed novel differentially methylated regions (DMRs) associated with adipogenesis. DMR analysis suggested that though DNA methylation is global, noticeable changes are observed at specific sites known as "hotspots." Hotspots are genomic regions rich in transcription factor (TF) binding sites and exhibit methylation-dependent TF binding. Subsequent analysis indicated hotspots as part of DMRs. The gene expression profile of key adipogenic genes in differentiating adipocytes is context-dependent, as we found a direct and inverse relationship between promoter DNA methylation and gene expression.

Whole-genome landscape of histone H3K4me3 modification during sperm cell lineage development in tomato

BackgroundDuring male gametogenesis of flowering plants, sperm cell lineage (microspores, generative cells, and sperm cells) differentiated from somatic cells and acquired different cell fates. Trimethylation of histone H3 on lysine 4 (H3K4me3) epigenetically contributes to this process, however, it remained unclear how H3K4me3 influences the gene expression in each cell type. Here, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) to obtain a genome-wide landscape of H3K4me3 during sperm cell lineage development in tomato (Solanum lycopersicum).ResultsWe show that H3K4me3 peaks were mainly enriched in the promoter regions, and intergenic H3K4me3 peaks expanded as sperm cell lineage differentiated from somatic cells. H3K4me3 was generally positively associated with transcript abundance and served as a better indicator of gene expression in somatic and vegetative cells, compared to sperm cell lineage. H3K4me3 was mutually exclusive with DNA methylation at 3’ proximal of the transcription start sites. The microspore maintained the H3K4me3 features of somatic cells, while generative cells and sperm cells shared an almost identical H3K4me3 pattern which differed from that of the vegetative cell. After microspore division, significant loss of H3K4me3 in genes related to brassinosteroid and cytokinin signaling was observed in generative cells and vegetative cells, respectively.ConclusionsOur results suggest the asymmetric division of the microspore significantly reshapes the genome-wide distribution of H3K4me3. Selective loss of H3K4me3 in genes related to hormone signaling may contribute to functional differentiation of sperm cell lineage. This work provides new resource data for the epigenetic studies of gametogenesis in plants.

Microbiota regulates the TET1-mediated DNA hydroxymethylation program in innate lymphoid cell differentiation

Innate lymphoid cell precursors (ILCPs) develop into distinct subsets of innate lymphoid cells (ILCs) with specific functions. The epigenetic program underlying the differentiation of ILCPs into ILC subsets remains poorly understood. Here, we reveal the genome-wide distribution and dynamics of the DNA methylation and hydroxymethylation in ILC subsets and their respective precursors. Additionally, we find that the DNA hydroxymethyltransferase TET1 suppresses ILC1 but not ILC2 or ILC3 differentiation. TET1 deficiency promotes ILC1 differentiation by inhibiting TGF-β signaling. Throughout ILCP differentiation at postnatal stage, gut microbiota contributes to the downregulation of TET1 level. Microbiota decreases the level of cholic acid in the gut, impairs TET1 expression and suppresses DNA hydroxymethylation, ultimately resulting in an expansion of ILC1s. In adult mice, TET1 suppresses the hyperactivation of ILC1s to maintain intestinal homeostasis. Our findings provide insights into the microbiota-mediated epigenetic programming of ILCs, which links microbiota-DNA methylation crosstalk to ILC differentiation.

Genome-Wide Computational Prediction and Analysis of Noncoding RNAs in Oleidesulfovibrio alaskensis G20.

Noncoding RNAs (ncRNAs) play key roles in the regulation of important pathways, including cellular growth, stress management, signaling, and biofilm formation. Sulfate-reducing bacteria (SRB) contribute to huge economic losses causing microbial-induced corrosion through biofilms on metal surfaces. To effectively combat the challenges posed by SRB, it is essential to understand their molecular mechanisms of biofilm formation. This study aimed to identify ncRNAs in the genome of a model SRB, Oleidesulfovibrio alaskensis G20 (OA G20). Three in silico approaches revealed genome-wide distribution of 37 ncRNAs excluding tRNAs in the OA G20. These ncRNAs belonged to 18 different Rfam families. This study identified riboswitches, sRNAs, RNP, and SRP. The analysis revealed that these ncRNAs could play key roles in the regulation of several pathways of biosynthesis and transport involved in biofilm formation by OA G20. Three sRNAs, Pseudomonas P10, Hammerhead type II, and sX4, which were found in OA G20, are rare and their roles have not been determined in SRB. These results suggest that applying various computational methods could enrich the results and lead to the discovery of additional novel ncRNAs, which could lead to understanding the "rules of life of OA G20" during biofilm formation.