Abstract
DNA 6mA modification, an important newly discovered epigenetic mark, plays a crucial role in organisms and has been attracting more and more attention in recent years. The soybean is economically the most important bean in the world, providing vegetable protein for millions of people. However, the distribution pattern and function of 6mA in soybean are still unknown. In this study, we decoded 6mA modification to single-nucleotide resolution in wild and cultivated soybeans, and compared the 6mA differences between cytoplasmic and nuclear genomes and between wild and cultivated soybeans. The motif of 6mA in the nuclear genome was conserved across the two kinds of soybeans, and ANHGA was the most dominant motif in wild and cultivated soybeans. Genes with 6mA modification in the nucleus had higher expression than those without modification. Interestingly, 6mA distribution patterns in cytoplasm for each soybean were significantly different from those in nucleus, which was reported for the first time in soybean. Our research provides a new insight in the deep analysis of cytoplasmic genomic DNA modification in plants.