A spatially resolved multiomic single-cell atlas of soybean development.

Cis-regulatory elements (CREs) precisely control spatiotemporal gene expression in cells. Using a spatially resolved single-cell atlas of gene expression with chromatin accessibility across ten soybean tissues, we identified 103 distinct cell types and 303,199 accessible chromatin regions (ACRs). Nearly 40% of the ACRs showed cell-type-specific patterns and were enriched for transcription factor (TF) motifs defining diverse cell identities. We identified de novo enriched TF motifs and explored conservation of gene regulatory networks underpinning legume symbiotic nitrogen fixation. With comprehensive developmental trajectories for endosperm and embryo, we uncovered the functional transition of the three sub-cell types of endosperm, identified 13 sucrose transporters sharing the DOF11 motif that were co-up-regulated in late peripheral endosperm and identified key embryo cell-type specification regulators during embryogenesis, including a homeobox TF that promotes cotyledon parenchyma identity. This resource provides a valuable foundation for analyzing gene regulatory programs in soybean cell types across tissues and life stages.

Dynamic evolution of the heterochromatin sensing histone demethylase IBM1.

Heterochromatin is critical for maintaining genome stability, especially in flowering plants, where it relies on a feedback loop involving the H3K9 methyltransferase, KRYPTONITE (KYP), and the DNA methyltransferase CHROMOMETHYLASE3 (CMT3). The H3K9 demethylase INCREASED IN BONSAI METHYLATION 1 (IBM1) counteracts the detrimental consequences of KYP-CMT3 activity in transcribed genes. IBM1 expression in Arabidopsis is uniquely regulated by methylation of the 7th intron, allowing it to monitor global H3K9me2 levels. We show the methylated intron is prevalent across flowering plants and its underlying sequence exhibits dynamic evolution. We also find extensive genetic and expression variations in KYP, CMT3, and IBM1 across flowering plants. We identify Arabidopsis accessions resembling weak ibm1 mutants and Brassicaceae species with reduced IBM1 expression or deletions. Evolution towards reduced IBM1 activity in some flowering plants could explain the frequent natural occurrence of diminished or lost CMT3 activity and loss of gene body DNA methylation, as cmt3 mutants in A. thaliana mitigate the deleterious effects of IBM1.

Haplotype-resolved genome assembly of Populus tremula × P. alba reveals aspen-specific megabase satellite DNA.

Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717-1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence-sensitive analyses. Here we report a telomere-to-telomere genome for this hybrid aspen with two chromosome-scale, haplotype-resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array-based and array-less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus) but PtaM147-like sequences occur in LTR-retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis-driven research to probe into the function of the abundant and aspen-specific PtaM147 satellite DNA.

DNMT1 mutant ants develop normally but have disrupted oogenesis.

Although DNA methylation is an important gene regulatory mechanism in mammals, its function in arthropods remains poorly understood. Studies in eusocial insects have argued for its role in caste development by regulating gene expression and splicing. However, such findings are not always consistent across studies, and have therefore remained controversial. Here we use CRISPR/Cas9 to mutate the maintenance DNA methyltransferase DNMT1 in the clonal raider ant, Ooceraea biroi. Mutants have greatly reduced DNA methylation, but no obvious developmental phenotypes, demonstrating that, unlike mammals, ants can undergo normal development without DNMT1 or DNA methylation. Additionally, we find no evidence of DNA methylation regulating caste development. However, mutants are sterile, whereas in wild-type ants, DNMT1 is localized to the ovaries and maternally provisioned into nascent oocytes. This supports the idea that DNMT1 plays a crucial but unknown role in the insect germline.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae.

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus. We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration. Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis-elements of duplicated genes are associated with convergent expression. Genome-wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa-originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation. Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis-trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

Heterochromatin is a quantitative trait associated with spontaneous epiallele formation.

Epialleles are meiotically heritable variations in expression states that are independent from changes in DNA sequence. Although they are common in plant genomes, their molecular origins are unknown. Here we show, using mutant and experimental populations, that epialleles in Arabidopsis thaliana that result from ectopic hypermethylation are due to feedback regulation of pathways that primarily function to maintain DNA methylation at heterochromatin. Perturbations to maintenance of heterochromatin methylation leads to feedback regulation of DNA methylation in genes. Using single base resolution methylomes from epigenetic recombinant inbred lines (epiRIL), we show that epiallelic variation is abundant in euchromatin, yet, associates with QTL primarily in heterochromatin regions. Mapping three-dimensional chromatin contacts shows that genes that are hotspots for ectopic hypermethylation have increases in contact frequencies with regions possessing H3K9me2. Altogether, these data show that feedback regulation of pathways that have evolved to maintain heterochromatin silencing leads to the origins of spontaneous hypermethylated epialleles.

The catalytic core of DEMETER guides active DNA demethylation in Arabidopsis.

The Arabidopsis DEMETER (DME) DNA glycosylase demethylates the maternal genome in the central cell prior to fertilization and is essential for seed viability. DME preferentially targets small transposons that flank coding genes, influencing their expression and initiating plant gene imprinting. DME also targets intergenic and heterochromatic regions, but how it is recruited to these differing chromatin landscapes is unknown. The C-terminal half of DME consists of 3 conserved regions required for catalysis in vitro. We show that this catalytic core guides active demethylation at endogenous targets, rescuing dme developmental and genomic hypermethylation phenotypes. However, without the N terminus, heterochromatin demethylation is significantly impeded, and abundant CG-methylated genic sequences are ectopically demethylated. Comparative analysis revealed that the conserved DME N-terminal domains are present only in flowering plants, whereas the domain architecture of DME-like proteins in nonvascular plants mainly resembles the catalytic core, suggesting that it might represent the ancestral form of the 5mC DNA glycosylase found in plant lineages. We propose a bipartite model for DME protein action and suggest that the DME N terminus was acquired late during land plant evolution to improve specificity and facilitate demethylation at heterochromatin targets.

De Novo Transcriptome Analysis to Identify Anthocyanin Biosynthesis Genes Responsible for Tissue-Specific Pigmentation in Zoysiagrass (Zoysia japonica Steud.).

Zoysiagrass (Zoysia japonica Steud.) is commonly found in temperate climate regions and widely used for lawns, in part, owing to its uniform green color. However, some zoysiagrass cultivars accumulate red to purple pigments in their spike and stolon tissues, thereby decreasing the aesthetic value. Here we analyzed the anthocyanin contents of two zoysiagrass cultivars 'Anyang-jungji' (AJ) and 'Greenzoa' (GZ) that produce spikes and stolons with purple and green colors, respectively, and revealed that cyanidin and petunidin were primarily accumulated in the pigmented tissues. In parallel, we performed a de novo transcriptome assembly and identified differentially expressed genes between the two cultivars. We found that two anthocyanin biosynthesis genes encoding anthocyanidin synthase (ANS) and dihydroflavonol 4-reductase (DFR) were preferentially upregulated in the purple AJ spike upon pigmentation. Both ANS and DFR genes were also highly expressed in other zoysiagrass cultivars with purple spikes and stolons, but their expression levels were significantly low in the cultivars with green tissues. We observed that recombinant ZjDFR1 and ZjANS1 proteins successfully catalyze the conversions of dihydroflavonols into leucoanthocyanidins and leucoanthocyanidins into anthocyanidins, respectively. These findings strongly suggest that upregulation of ANS and DFR is responsible for tissue-specific anthocyanin biosynthesis and differential pigmentation in zoysiagrass. The present study also demonstrates the feasibility of a de novo transcriptome analysis to identify the key genes associated with specific traits, even in the absence of reference genome information.

AP endonucleases process 5-methylcytosine excision intermediates during active DNA demethylation in Arabidopsis.

DNA methylation is a primary epigenetic modification regulating gene expression and chromatin structure in many eukaryotes. Plants have a unique DNA demethylation system in that 5-methylcytosine (5mC) is directly removed by DNA demethylases, such as DME/ROS1 family proteins, but little is known about the downstream events. During 5mC excision, DME produces 3'-phosphor-α, ß-unsaturated aldehyde and 3'-phosphate by successive ß- and δ-eliminations, respectively. The kinetic studies revealed that these 3'-blocking lesions persist for a significant amount of time and at least two different enzyme activities are required to immediately process them. We demonstrate that Arabidopsis AP endonucleases APE1L, APE2 and ARP have distinct functions to process such harmful lesions to allow nucleotide extension. DME expression is toxic to E. coli due to excessive 5mC excision, but expression of APE1L or ARP significantly reduces DME-induced cytotoxicity. Finally, we propose a model of base excision repair and DNA demethylation pathway unique to plants.

Excision of 5-hydroxymethylcytosine by DEMETER family DNA glycosylases.

In plants and animals, 5-methylcytosine (5mC) serves as an epigenetic mark to repress gene expression, playing critical roles for cellular differentiation and transposon silencing. Mammals also have 5-hydroxymethylcytosine (5hmC), resulting from hydroxylation of 5mC by TET family-enzymes. 5hmC is abundant in mouse Purkinje neurons and embryonic stem cells, and regarded as an important intermediate for active DNA demethylation in mammals. However, the presence of 5hmC in plants has not been clearly demonstrated. In Arabidopsis, the DEMETER (DME) family DNA glycosylases efficiently remove 5mC, which results in DNA demethylation and transcriptional activation of target genes. Here we show that DME and ROS1 have a significant 5hmC excision activity in vitro, although we detected no 5hmC in Arabidopsis, suggesting that it is very unlikely for plants to utilize 5hmC as a DNA demethylation intermediate. Our results indicate that both plants and animals have 5mC in common but DNA demethylation systems have independently evolved with distinct mechanisms.