Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation ofcis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.
more »
« less
This content will become publicly available on December 1, 2026
From duplication to divergence: Single-cell insights into transcriptional and cis-regulatory landscapes in soybean
Abstract Gene duplication is a major source of evolutionary innovation, enabling the emergence of novel expression patterns and functions. Leveraging single-cell genomics, we investigated the transcriptional and cis-regulatory landscapes of duplicated genes in cultivated soybean (Glycine max), which has undergone 2 rounds of whole-genome duplication. Our analysis revealed extensive diversity of transcriptional profiles within and across tissues among duplicated gene pairs. Within-tissue divergence was largely attributable to genetic variation in their associated accessible chromatin regions (ACRs), where cis-regulatory elements reside, whereas cross-tissue divergence was more likely shaped by dynamics in ACR chromatin accessibility profiles across tissues. Distinct duplication mechanisms also likely give rise to different types of cis-regulatory variants, contributing variably to transcriptional divergence. By comparing ACRs associated with gene sets derived from 2 rounds of whole-genome duplication and sharing a common ancestral gene, we found that most ACRs retained one or multiple corresponding duplicated sequences in which mutations gradually accumulated over time, while a subset likely arose de novo. Finally, we traced the evolution of cell-type-specific expression and cell-type-specific ACRs within duplicated gene sets, illustrating a powerful framework for identifying candidate regulatory regions driving cell-type-specific expression. Collectively, our findings highlight the important role of cis-regulatory evolution in shaping transcriptional divergence in a spatiotemporal manner, uncovered with the resolution of single-cell genomics.
more »
« less
- Award ID(s):
- 1856627
- PAR ID:
- 10687878
- Publisher / Repository:
- American Society of Plant Biologists
- Date Published:
- Journal Name:
- The Plant Cell
- Volume:
- 37
- Issue:
- 12
- ISSN:
- 1040-4651
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Purugganan, Michael (Ed.)Subgenome dominance after whole-genome duplication (WGD) has been observed in many plant species. However, the degree to which the chromatin environment affects this bias has not been explored. Here, we compared the dominant subgenome (maize1) and the recessive subgenome (maize2) with respect to patterns of sequence substitutions, genes expression, transposable element accumulation, small interfering RNAs, DNA methylation, histone modifications, and accessible chromatin regions (ACRs). Our data show that the degree of bias between subgenomes for all the measured variables does not vary significantly when both of the WGD genes are located in pericentromeric regions. Our data further indicate that the location of maize1 genes in chromosomal arms is pivotal for maize1 to maintain its dominance, but location has a less effect on maize2 homoeologs. In addition to homoeologous genes, we compared ACRs, which often harbor cis-regulatory elements, between the two subgenomes and demonstrate that maize1 ACRs have a higher level of chromatin accessibility, a lower level of sequence substitution, and are enriched in chromosomal arms. Furthermore, we find that a loss of maize1 ACRs near their nearby genes is associated with a reduction in purifying selection and expression of maize1 genes relative to their maize2 homoeologs. Taken together, our data suggest that chromatin environment and cis-regulatory elements are important determinants shaping the divergence and evolution of duplicated genes.more » « less
-
Human-specific segmental duplications (HSDs) contain millions of base pairs of sequence unique to the human genome, including genes that shape neurodevelopment. Despite their young age (<6 million years), HSD genes exhibit widespread regulatory divergence, with paralog-specific expression patterns documented across a variety of tissues and cell types. Using long-read expression and epigenomic data, we show that human-specific paralogs tend to have lower activity than the shared, ancestral ones. To systematically characterize the cis-regulatory elements (CREs) within HSDs and understand patterns of regulatory change in recently evolved gene families, we conducted a massively parallel reporter assay of 7,160 human duplicated and chimpanzee orthologous sequences in lymphoblastoid (GM12878) and neuroblastoma (SH-SY5Y) cell lines. A large proportion (14–24%) of sequences exhibited differential activity relative to the chimpanzee ortholog (or between human paralogs), mostly with small fold-differences. Combining measured activity levels across all assayed sequences, predicted differences in cis-regulatory activity correlated with mRNA levels in SH-SY5Y. Differentially active CREs were validated for CHRFAM7A, HYDIN2, and SRGAP2C that may contribute to paralog-specific expression patterns and thereby to human-specific traits. While we find some changes in CRE activity shared between duplicate paralogs likely driving regulatory divergence in gene expression, consideration of non-shared adjacent sequences to duplications suggests a larger role for altered genome positional effects. In all, this work suggests that functional divergence of duplicated CREs contributes moderately to regulatory divergence of HSD genes and uncovers enhancers that are candidate drivers of human-specific regulatory patterns.more » « less
-
Abstract Motivation Genome-wide profiles of chromatin accessibility and gene expression in diverse cellular contexts are critical to decipher the dynamics of transcriptional regulation. Recently, convolutional neural networks have been used to learn predictive cis-regulatory DNA sequence models of context-specific chromatin accessibility landscapes. However, these context-specific regulatory sequence models cannot generalize predictions across cell types. Results We introduce multi-modal, residual neural network architectures that integrate cis-regulatory sequence and context-specific expression of trans-regulators to predict genome-wide chromatin accessibility profiles across cellular contexts. We show that the average accessibility of a genomic region across training contexts can be a surprisingly powerful predictor. We leverage this feature and employ novel strategies for training models to enhance genome-wide prediction of shared and context-specific chromatin accessible sites across cell types. We interpret the models to reveal insights into cis- and trans-regulation of chromatin dynamics across 123 diverse cellular contexts. Availability and implementation The code is available at https://github.com/kundajelab/ChromDragoNN. Supplementary information Supplementary data are available at Bioinformatics online.more » « less
-
Lineage-specific evolution of regulatory landscapes in a polyploid plant and its diploid progenitorsCis-regulatory elements are specific DNA sequences that control gene expression in a spatiotemporal manner, and variation within these elements represents a major source of phenotypic diversity and evolutionary innovation. Nevertheless, how regulatory elements evolve and shape gene expression remains poorly understood, particularly in plants. The well-resolved phylogeny of allopolyploid peanut (Arachis hypogaea) and its diploid progenitors,Arachis duranensisandArachis ipaensis, provides an ideal system to investigate the regulatory evolution at a lineage-specific level. By integrating comparative analyses of sequence similarity, chromatin accessibility, histone modifications, conserved noncoding sequences (CNSs), and gene expression, we reconstructed the evolutionary trajectories of Accessible Chromatin Regions (ACRs), where regulatory elements typically reside, and revealed their distinct contributions to homoeolog expression bias, unequal expressions between homoeologs. Most ACRs exhibited high sequence similarity, comparable chromatin accessibility, and conserved states for H3K4me3, H3K56ac, and H3K36me3, indicating regulatory stability after hybridization and polyploidization. However, a subset of novel ACRs emerged de novo from previously nonregulatory regions or through sequence mutations in preexisting ACRs, arising at different rates and evolutionary stages. Notably, even highly sequence-conserved ACRs exhibited substantial variation in chromatin accessibility, consistent with CNS composition differences and minor sequence variation, although causal relationships remain to be demonstrated. Our analyses further revealed a complex spectrum of CNS dynamics within the diploid-polyploid framework. Overall, our study provides empirical insights into the fine-scale evolution of plant regulatory landscapes and complements previous large-scale comparisons across distant lineages.more » « less
An official website of the United States government
