Dataset Information

Comparative genomic analysis of the human genome and six bat genomes using unsupervised machine learning: Mb-level CpG and TFBS islands.

ABSTRACT:

Background

Emerging infectious disease-causing RNA viruses, such as the SARS-CoV-2 and Ebola viruses, are thought to rely on bats as natural reservoir hosts. Since these zoonotic viruses pose a great threat to humans, it is important to characterize the bat genome from multiple perspectives. Unsupervised machine learning methods for extracting novel information from big sequence data without prior knowledge or particular models are highly desirable for obtaining unexpected insights. We previously established a batch-learning self-organizing map (BLSOM) of the oligonucleotide composition that reveals novel genome characteristics from big sequence data.

Results

In this study, using the oligonucleotide BLSOM, we conducted a comparative genomic study of humans and six bat species. BLSOM is an explainable-type machine learning algorithm that reveals the diagnostic oligonucleotides contributing to sequence clustering (self-organization). When unsupervised machine learning reveals unexpected and/or characteristic features, these features can be studied in more detail via the much simpler and more direct standard distribution map method. Based on this combined strategy, we identified the Mb-level enrichment of CG dinucleotide (Mb-level CpG islands) around the termini of bat long-scaffold sequences. In addition, a class of CG-containing oligonucleotides were enriched in the centromeric and pericentromeric regions of human chromosomes. Oligonucleotides longer than tetranucleotides often represent binding motifs for a wide variety of proteins (e.g., transcription factor binding sequences (TFBSs)). By analyzing the penta- and hexanucleotide composition, we observed the evident enrichment of a wide range of hexanucleotide TFBSs in centromeric and pericentromeric heterochromatin regions on all human chromosomes.

Conclusion

Function of transcription factors (TFs) beyond their known regulation of gene expression (e.g., TF-mediated looping interactions between two different genomic regions) has received wide attention. The Mb-level TFBS and CpG islands are thought to be involved in the large-scale nuclear organization, such as centromere and telomere clustering. TFBSs, which are enriched in centromeric and pericentromeric heterochromatin regions, are thought to play an important role in the formation of nuclear 3D structures. Our machine learning-based analysis will help us to understand the differential features of nuclear 3D structures in the human and bat genomes.

SUBMITTER: Iwasaki Y

PROVIDER: S-EPMC9264310 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background?-Glutamyl hydrolase (GGH) regulates intracellular folates and antifolates such as methotrexate (MTX) for proper nucleotide biosynthesis and antifolate-induced cytotoxicity, respectively. In addition to genetic polymorphism and karyotypic abnormalities, methylation of CpG island 1 (CpG1) in the promoter region is found to modulate GGH activity by reducing GGH mRNA expression in acute lymphoblastic leukemia (ALL) cells. We aim to investigate methylation status of two CpG islands (CpG1 and CpG2) in the GGH promoter region in pediatric patients with ALL and acute myelogenous leukemia (AML).Methods70B-ALL, 29 AML, 10 ITP (idiopathic thrombocytopenic purpura) and 40 healthy children are recruited in the present study. MS-HRM (methylation-sensitive high-resolution melting) and bisulfite sequencing PCR (BSP) are used to detect methylation change and its level in CpG1 and CpG2 in the GGH promoter region. GGH mRNA expression is quantified by real-time PCR. Correlation between CpG island methylation and GGH mRNA expression is assessed by statistical software.ResultsMethylations of CpG1 are detected in leukemia cells samples obtained from 30.9% (21/68) of patients with ALL and 20.7% (6/29) of patients with AML. These methylations are not detected in the controls. Methylations of CpG2 are detected in leukemia cell samples obtained from 44.1% (30/68) of the ALL patients and 37.9% (11/29) of the AML patients. These percentages are significantly higher than that observed in the control cell samples: 6.0% (3/50) (Fisher's exact test, P = 0.000). The abundance of CpG1 methylation in all leukemia cell samples is classified as Grade I (methylation level is less than 10%) and the abundance of CpG2 methylation in leukemia cell samples is classified in separate grades. Our results indicate that methylation of CpG1 or hypermethylation (the methylation level is greater than 10%) of CpG2 could significantly reduce GGH mRNA expression in leukemia cells from the ALL and AML patients (ALL-CpG1: t = 4.632, P = 0.000; ALL-CpG2: t = 3.250, P = 0.006; AML-CpG1: t = -2.254, P = 0.037; AML-CpG2: t = 1.328, P = 0.202).ConclusionEither methylation of CpG1 or hypermethylation of CpG2 in GGH promoter region can significantly reduce GGH mRNA expression in pediatric patients with acute leukemia, which can improve the response to treatment.

Project description:BackgroundMammalian CpG islands (CGIs) normally escape DNA methylation in all adult tissues and developmental stages. However, in our previous study we unexpectedly identified many methylated CGIs in human peripheral blood leukocytes. Methylated CpG dinucleotides convert to TpG dinucleotides through deaminization of their cytosine bases more frequently than hypomethylated CpG dinucleotides. Therefore, we wondered how methylated CGIs in germline or non-germline cells maintain their CpG-rich sequences. It is known that events such as germline hypomethylation, CpG selection, biased gene conversion (BGC), and frequent CpG fixation can contribute to the maintenance of CpG-rich sequences in methylated CGIs in germline or non-germline cells. However, it has not been investigated which of the processes maintain CpG-rich sequences of methylated CGIs in each genomic position.ResultsIn this study, we comprehensively examined the contribution of the processes described above to the maintenance of CpG-rich sequences in methylated CGIs in germline and non-germline cells which were classified by genomic positions. Approximately 60-80% of CGIs with high methylation in H1 cell line (H1-HM) in all the genomic positions showed a low average CpG→TpG/CpA substitution rate. In contrast, fewer than half the numbers of CGIs with H1-HM in all the genomic positions showed a low average CpG→TpG/CpA substitution rate and low levels of methylation in sperm cells (SPM-LM). Furthermore, a small fraction of CGIs with a low average CpG→TpG/CpA substitution rate and high levels of methylation in sperm cells (SPM-HM) showed CpG selection. On the other hand, independent of the positions in genes, most CGIs with SPM-HM showed a slightly higher average TpG/CpA→CpG substitution rate compared with those with SPM-LM.ConclusionsRelatively high numbers (approximately 60-80%) of CGIs with H1-HM in all the genomic positions preserve their CpG-rich sequences by a low CpG→TpG/CpA substitution rate caused mainly by their SPM-LM, and for those with SPM-HM partly by CpG selection and TpG/CpA→CpG fixation. BGC has little contribution to the maintenance of CpG-rich sequences of CGIs with SPM-HM which were classified by genomic positions.