Developmental Gene Regulation Research Articles

A Novel Solid Media-Free In-Planta Soybean (Glycine max. (L) Merr.) Transformation Approach

Soybean’s lengthy protocols for transgenic plant production are a bottleneck in the transgenic breeding of this crop. Explants cultured on a medium for an extended duration exhibit unanticipated modifications. Stress-induced somaclonal variations and in vitro contaminations also cause substantial losses of transgenic plants. This effect could potentially be mitigated by direct shoot regeneration without solid media or in-planta transformation. The current study focused primarily on developing a rapid and effective media-free in-planta transformation technique for three soybean genotypes (Wm82) and our newly developed two hybrids, designated as ZX-16 and ZX-3. The whole procedure for a transgenic plant takes the same time as a stable grown seedling. Multiple axillary shoots were regenerated on stable-grown soybean seedlings without the ectopic expression of developmental regulatory genes. An approximate amount of 200 µL medium with a growth regulator was employed for shoot organogenesis and growth. The maximal shoot regeneration percentages in the Wm82 and ZX-3 genotypes were 87.1% and 84.5%, respectively. The stable transformation ranged from 3% to 8.0%, with an average of 5.5%. This approach seems to be the opposite of the hairy root transformation method, which allowed transgenic shoots to be regenerated on normal roots. Further improvement regarding an increase in the transformation efficiency and of this technique for a broad range of soybean genotypes and other dicot species would be extremely beneficial in achieving increased stable transformation.

Number of human protein interactions correlates with structural, but not regulatory conservation of the respective genes.

The differential ratio of nonsynonymous to synonymous nucleotide substitutions (dN/dS) is a common measure of the rate of structural evolution in proteincoding genes. In addition, we recently suggested that the proportion of transposable elements in gene promoters that host functional genomic sites serves as a marker of the rate of regulatory evolution of genes. Such functional genomic regions may include transcription factor binding sites and modified histone binding loci. Here, we constructed a model of the human interactome based on 600,136 documented molecular interactions and investigated the overall relationship between the number of interactions of each protein and the rate of structural and regulatory evolution of the corresponding genes. By evaluating a total of 4,505 human genes and 1,936 molecular pathways we found a general correlation between structural and regulatory evolution rate metrics (Spearman 0.08-0.16 and 0.25-0.37 for gene and pathway levels, respectively, p < 0.01). Further exploration revealed in the established human interactome model lack of correlation between the rate of gene regulatory evolution and the number of protein interactions on gene level, and weak negative correlation (∼0.15) on pathway level. We also found a statistically significant negative correlation between the rate of gene structural evolution and the number of protein interactions (Spearman -0.11 and -0.3 for gene and pathway levels, respectively, p < 0.01). Our result suggests stronger structural rather than regulatory conservation of genes whose protein products have multiple interaction partners.

Understanding developmental system drift.

Developmental system drift (DSD) occurs when the genetic basis for homologous traits diverges over time despite conservation of the phenotype. In this Review, we examine the key ideas, evidence and open problems arising from studies of DSD. Recent work suggests that DSD may be pervasive, having been detected across a range of different organisms and developmental processes. Although developmental research remains heavily reliant on model organisms, extrapolation of findings to non-model organisms can be error-prone if the lineages have undergone DSD. We suggest how existing data and modelling approaches may be used to detect DSD and estimate its frequency. More direct study of DSD, we propose, can inform null hypotheses for how much genetic divergence to expect on the basis of phylogenetic distance, while also contributing to principles of gene regulatory evolution.

Exploring an economic and highly efficient genetic transformation and genome-editing system for radish through developmental regulators and visible reporter.

Radish (Raphanus sativus L.) is one of the most important root vegetable crops worldwide. However, gene function exploration and germplasm innovation still face tremendous challenges due to its extremely low transformation efficiency. Here, an economic and highly efficient genetic transformation method for radish was explored by Agrobacterium rhizogenes-mediated transformation with the help of combining special developmental regulator (DR) genes and the visual identification reporter. Firstly, the RUBY gene, a betalain biosynthesis system, could result in a visual red-violet color used as a convenient and effective reporter for monitoring transgenic hairy roots screening of radish. However, the hairy roots-to-shoots conversion system of radish still stands as a barrier to the obtainment of whole transgenic plants, although different hormone combinations and various culture conditions were tried. Following, two DR genes including Wuschel2 (Wus2) and isopentenyl transferase (ipt), as well as their combination Wus2-ipt were introduced for the shoot regeneration capacity improvement. The results showed that the transgenic shoots could be directly generated without externally supplying any hormones in the presence of a Wus2-ipt combination. Then, Wus2-ipt along with the RUBY reporter was employed to establish an efficient genetic transformation system of radish. Moreover, this system was applied in generating gene-edited radish plants and the phytoene desaturase (RsPDS) gene was effectively knockout through albino phenotype observation and sequencing analysis. These findings have the potential to be widely applied in genetic transformation and genome-editing genetic improvement of other vegetable species.

Comparative single-cell multiome identifies evolutionary changes in neural progenitor cells during primate brain development.

Understanding the cellular and genetic mechanisms driving human-specific features of cortical development remains a challenge. We generated a cell-type resolved atlas of transcriptome and chromatin accessibility in the developing macaque and mouse prefrontal cortex (PFC). Comparing with published human data, our findings demonstrate that although the cortex cellular composition is overall conserved across species, progenitor cells show significant evolutionary divergence in cellular properties. Specifically, human neural progenitors exhibit extensive transcriptional rewiring in growth factor and extracellular matrix (ECM) pathways. Expression of the human-specific progenitor marker ITGA2 in the fetal mouse cortex increases the progenitor proliferation and the proportion of upper-layer neurons. These transcriptional divergences are primarily driven by altered activity in the distal regulatory elements. The chromatin regions with human-gained accessibility are enriched with human-specific sequence changes and polymorphisms linked to intelligence and neuropsychiatric disorders. Our results identify evolutionary changes in neural progenitors and putative gene regulatory mechanisms shaping primate brain evolution.

Mono-allelic epigenetic regulation of bi-directional silencing of RNA Polymerase II polycistronic transcription initiation in Trypanosoma brucei.

The single-cell parasite Trypanosoma brucei causes lethal diseases in both humans and livestock. T. brucei undergoes multiple developmental changes to adapt in different environments during its digenetic life cycle. With protein-coding genes organized as polycistronic transcription and apparent absence of promoter-mediated regulation of transcription initiation, it is believed that developmental gene regulation in trypanosomes is essentially post-transcriptional. In this study, we found reversible Pol II transcriptional silencing of two adjacent polycistronic gene arrays that correlates with the novel DNA base J modification of the shared promoter region. Our findings support epigenetic regulation of Pol II transcription initiation as a viable mechanism of gene expression control in T. brucei . This has implications for our understanding how trypanosomes utilize polycistronic genome organization to regulate gene expression during its life cycle.

A PRE loop at the dac locus acts as a topological chromatin structure that restricts and specifies enhancer-promoter communication.

Three-dimensional (3D) genome folding has a fundamental role in the regulation of developmental genes by facilitating or constraining chromatin interactions between cis-regulatory elements (CREs). Polycomb response elements (PREs) are a specific kind of CRE involved in the memory of transcriptional states in Drosophila melanogaster. PREs act as nucleation sites for Polycomb group (PcG) proteins, which deposit the repressive histone mark H3K27me3, leading to the formation of a class of topologically associating domain (TAD) called a Polycomb domain. PREs can establish looping contacts that stabilize the gene repression of key developmental genes during development. However, the mechanism by which PRE loops fine-tune gene expression is unknown. Using clustered regularly interspaced short palindromic repeats and Cas9 genome engineering, we specifically perturbed PRE contacts or enhancer function and used complementary approaches including 4C-seq, Hi-C and Hi-M to analyze how chromatin architecture perturbation affects gene expression. Our results suggest that the PRE loop at the dac gene locus acts as a constitutive 3D chromatin scaffold during Drosophila development that forms independently of gene expression states and has a versatile function; it restricts enhancer-promoter communication and contributes to enhancer specificity.

Life stage–specific poly(A) site selection regulated by Trypanosoma brucei DRBD18

The kinetoplastid parasite, Trypanosoma brucei, undergoes a complex life cycle entailing slender and stumpy bloodstream forms in mammals and procyclic and metacyclic forms (MFs) in tsetse fly hosts. The numerous gene regulatory events that underlie T. brucei differentiation between hosts, as well as between active and quiescent stages within each host, take place in the near absence of transcriptional control. Rather, differentiation is controlled by RNA-binding proteins (RBPs) that associate with mRNA 3' untranslated regions (3'UTRs) to impact RNA stability and translational efficiency. DRBD18 is a multifunctional T. brucei RBP, shown to impact mRNA stability, translation, export, and processing. Here, we use single-cell RNAseq to characterize transcriptomic changes in cell populations that arise upon DRBD18 depletion, as well as to visualize transcriptome-wide alterations to 3'UTR length. We show that in procyclic insect stages, DRBD18 represses expression of stumpy bloodstream form and MF transcripts. Additionally, DRBD18 regulates the 3'UTR lengths of over 1,500 transcripts, typically promoting the use of distal polyadenylation sites, and thus the inclusion of 3'UTR regulatory elements. Remarkably, comparison of polyadenylation patterns in DRBD18 knockdowns with polyadenylation patterns in stumpy bloodstream forms shows numerous similarities, revealing a role for poly(A) site selection in developmental gene regulation, and indicating that DRBD18 controls this process for a set of transcripts. RNA immunoprecipitation supports a direct role for DRBD18 in poly(A) site selection. This report highlights the importance of alternative polyadenylation in T. brucei developmental control and identifies a critical RBP in this process.

Emergence of the ZNF675 Gene During Primate Evolution-Influenced Human Neurodevelopment Through Changing HES1 Autoregulation.

In this study, we investigated recurrent copy number variations (CNVs) in the 19p12 locus, which are associated with neurodevelopmental disorders. The two genes in this locus, ZNF675 and ZNF681, arose via gene duplication in primates, and their presence in several pathological CNVs in the human population suggests that either or both of these genes are required for normal human brain development. ZNF675 and ZNF681 are members of the Krüppel-associated box zinc finger (KZNF) protein family, a class of transcriptional repressors important for epigenetic silencing of specific genomic regions. About 170 primate-specific KZNFs are present in the human genome. Although KZNFs are primarily associated with repressing retrotransposon-derived DNA, evidence is emerging that they can be co-opted for other gene regulatory processes. We show that genetic deletion of ZNF675 causes developmental defects in cortical organoids, and our data suggest that part of the observed neurodevelopmental phenotype is mediated by a gene regulatory role of ZNF675 on the promoter of the neurodevelopmental gene Hes family BHLH transcription factor 1 (HES1). We also find evidence for the recently evolved regulation of genes involved in neurological disorders, microcephalin 1 and sestrin 3. We show that ZNF675 interferes with HES1 auto-inhibition, a process essential for the maintenance of neural progenitors. As a striking example of how some KZNFs have integrated into preexisting gene expression networks, these findings suggest the emergence of ZNF675 has caused a change in the balance of HES1 autoregulation. The association of ZNF675 CNV with human developmental disorders and ZNF675-mediated regulation of neurodevelopmental genes suggests that it evolved into an important factor for human brain development.

Comprehensive Insights into the Remarkable Function and Regulatory Mechanism of FluG during Asexual Development in Beauveria bassiana.

Asexual development is the main propagation and transmission mode of Beauveria bassiana and the basis of its pathogenicity. The regulation mechanism of conidiation and the key gene resources for utilization are key links to improving the conidia yield and quality of Beauveria bassiana. Their clarification may promote the industrialization of fungal pesticides. Here, we compared the regulation of morphology, resistance to external stress, virulence, and nutrient utilization capacity between the upstream developmental regulatory gene fluG and the key genes brlA, abaA, and wetA in the central growth and development pathway. The results showed that the ΔbrlA and ΔabaA mutants completely lost the capacity to conidiate and that the ΔwetA mutant had seriously reduced conidiation capacity. Although the deletion of fluG did not reduce the conidiation ability as much as deletions of brlA, abaA, and wetA, it significantly reduced the fungal response to external stress, virulence, and nutrient utilization, while the deletion of the three other genes had little effect. Via transcriptome analysis and screening the yeast nuclear system library, we found that the differentially expressed genes in the ΔfluG mutants were concentrated in the signaling pathways of ABC transporters, propionate metabolism, tryptophan metabolism, DNA replication, mismatch repair, and fatty acid metabolism. FluG directly acted on 40 proteins that were involved in various signaling pathways such as metabolism, oxidative stress, and cell homeostasis. The analysis indicated that the regulatory function of fluG was mainly involved in DNA replication, cell homeostasis, fungal growth and metabolism, and the response to external stress. Our results revealed the biological function of fluG in asexual development and the responses to several environmental stresses as well as its influence on the asexual development regulatory network in B. bassiana.

Organ transformation by environmental disruption of protein integrity and epigenetic memory in Drosophila.

Despite significant progress in understanding epigenetic reprogramming of cells, the mechanistic basis of "organ reprogramming" by (epi-)gene-environment interactions remained largely obscure. Here, we use the ether-induced haltere-to-wing transformations in Drosophila as a model for epigenetic "reprogramming" at the whole organism level. Our findings support a mechanistic chain of events explaining why and how brief embryonic exposure to ether leads to haltere-to-wing transformations manifested at the larval stage and on. We show that ether interferes with protein integrity in the egg, leading to altered deployment of Hsp90 and widespread repression of Trithorax-mediated establishment of active H3K4me3 chromatin marks throughout the genome. Despite this global reduction, Ubx targets and wing development genes preferentially retain higher levels of H3K4me3 that predisposes these genes for later up-regulation in the larval haltere disc, hence the wing-like outcome. Consistent with compromised protein integrity during the exposure, the penetrance of bithorax transformations increases by genetic or chemical reduction of Hsp90 function. Moreover, joint reduction in Hsp90 and trx gene dosage can cause bithorax transformations even without exposure to ether, supporting an underlying epistasis between Hsp90 and trx loss-of-functions. These findings implicate environmental disruption of protein integrity at the onset of histone methylation with altered epigenetic regulation of developmental patterning genes. The emerging picture provides a unique example wherein the alleviation of the Hsp90 "capacitor function" by the environment drives a morphogenetic shift towards an ancestral-like body plan. The morphogenetic impact of chaperone response during a major setup of epigenetic patterns may be a general scheme for organ reprogramming by environmental cues.

SVEN_5003 is a Major Developmental Regulator in Streptomyces venezuelae.

Many regulatory genes that affect cellular development in Streptomyces, such as the canonical bld genes, have already been identified. However, in this study, we identified sven_5003 in Streptomyces venezuelae as a major new developmental regulatory gene, the deletion of which leads to a bald phenotype, typical of bld mutants, under multiple growth conditions. Our data indicated that disruption of sven_5003 also has a differential impact on the production of the two antibiotics jadomycin and chloramphenicol. Enhanced production of jadomycin but reduced production of chloramphenicol were detected in our sven_5003 mutant strain (S. venezuelae D5003). RNA-Seq analysis indicated that SVEN_5003 impacts expression of hundreds of genes, including genes involved in development, primary and secondary metabolism, and genes of unknown function, a finding confirmed by real-time PCR analysis. Transcriptional analysis indicated that sven_5003 is an auto-regulatory gene, repressing its own expression. Despite the evidence indicating that SVEN_5003 is a regulatory factor, a putative DNA-binding domain was not predicted from its primary amino acid sequence, implying an unknown regulatory mechanism by SVEN_5003. Our findings revealed that SVEN_5003 is a pleiotropic regulator with a critical role in morphological development in S. venezuelae.

Roles of Nuclear Orphan Receptors TR2 and TR4 during Hematopoiesis.

TR2 and TR4 (NR2C1 and NR2C2, respectively) are evolutionarily conserved nuclear orphan receptors capable of binding direct repeat sequences in a stage-specific manner. Like other nuclear receptors, TR2 and TR4 possess important roles in transcriptional activation or repression with developmental stage and tissue specificity. TR2 and TR4 bind DNA and possess the ability to complex with available cofactors mediating developmental stage-specific actions in primitive and definitive erythrocytes. In erythropoiesis, TR2 and TR4 are required for erythroid development, maturation, and key erythroid transcription factor regulation. TR2 and TR4 recruit and interact with transcriptional corepressors or coactivators to elicit developmental stage-specific gene regulation during hematopoiesis.

Transcriptional variations at perilesional site in vitiligo patients: probably a fight for cell-survival.

Exploration of gene expression variations is a potential source to unravel biological pathways involved in pathological changes in body and understand the mechanism underneath. Vitiligo patients were explored for gene expression changes transcriptionally at perilesional site in comparison to normal site of same patients for melanogenesis pathway (TYR, DCT & TYRP1) cell adhesion (MMPs & TIMP1), cell survival (BCL2 & BAX1) as well as proliferation, migration & development (SOX9, SOX10 & MITF) regulatory system, using skin biopsy samples. Results were also compared with changes in gene expression for melanocytes under stress after hydrogen peroxide treatment in-vitro. Gene amplification was carried out via real time PCR. We found increased expression of proliferation, migration & development regulatory genes as well as melanogenesis pathway genes at perilesional site of patients. In-vitro study also supports induced MITF expression and disturbed melanogenesis in melanocytes under stress. Expression level ratio of cell survival regulatory genes' (BCL2/BAX1) as well as cell adhesion regulatory genes (MMPs/TIMP1) was observed upregulated at patient's perilesional site however downregulated in hydrogen peroxide treated melanocytes in-vitro. Observed upregulated gene expression at perilesional site of patients may be via positive feedback loop in response to stress to increase cell tolerance power to survive against adverse conditions. Gene expression analysis suggests better cell survival and proliferation potential at perilesional site in vitiligo patients. It seems in-vivo conditions/growth factors supports cells to fight for survival to accommodate stressed conditions.

Interactions between U and V sex chromosomes during the life cycle of Ectocarpus.

In many animals and flowering plants, sex determination occurs in the diploid phase of the life cycle with XX/XY or ZW/ZZ sex chromosomes. However, in early diverging plants and most macroalgae, sex is determined by female (U) or male (V) sex chromosomes in a haploid phase called the gametophyte. Once the U and V chromosomes unite at fertilization to produce a diploid sporophyte, sex determination no longer occurs, raising key questions about the fate of the U and V sex chromosomes in the sporophyte phase. Here, we investigate genetic and molecular interactions of the UV sex chromosomes in both the haploid and diploid phases of the brown alga Ectocarpus. We reveal extensive developmental regulation of sex chromosome genes across its life cycle and implicate the TALE-HD transcription factor OUROBOROS in suppressing sex determination in the diploid phase. Small RNAs may also play a role in the repression of a female sex-linked gene, and transition to the diploid sporophyte coincides with major reconfiguration of histone H3K79me2, suggesting a more intricate role for this histone mark in Ectocarpus development than previously appreciated.

Developmentally dependent reprogramming of the Arabidopsis floral transcriptome under sufficient and limited water availability

BackgroundEnvironmental stresses negatively impact reproductive development and yield. Drought stress, in particular, has been examined during Arabidopsis reproductive development at morphological and transcriptomic levels. However, drought-responsive transcriptomic changes at different points in reproductive development remain unclear. Additionally, an investigation of the entire transcriptome at various stages during flower development is of great interest.ResultsHere, we treat Arabidopsis plants with well-watered and moderately and severely limiting water amounts when the first flowers reach maturity and generate RNA-seq datasets for early, middle, and late phases during flower development at 5, 6, and 7 days following treatment. Under different drought conditions, flowers in different developmental phases display differential sets of drought-responsive genes (DTGs), including those that are enriched in different GO functional categories, such as transcriptional regulation and response to stresses (early phase), lipid storage (middle phase), and pollen and seed development and metabolic processes (late phase). Some gene families have different members induced at different floral phases, suggesting that similar biochemical functions are carried out by distinct members. Developmentally-regulated genes (DVGs) with differential expression among the three floral phases belong to GO terms that are similar between water conditions, such as development and reproduction, metabolism and transport, and signaling and stress response. However, for different water conditions, such similar GO terms correspond to either distinct gene families or different members of a gene family, suggesting that drought affects the expression of distinct families or family members during reproductive development. A further comparison among transcriptomes of tissues collected on different days after treatment identifies differential gene expression, suggesting age-related genes (ARGs) might reflect the changes in the overall plant physiology in addition to drought response and development.ConclusionTogether, our study provides new insights into global transcriptome reprogramming and candidate genes for drought response, flower development, aging and coordination among these complex biological processes.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation

BackgroundRegulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging.ResultsHere we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes.We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters.ConclusionBased on these associations of (linearly) distal genes’ transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Morphological and dietary changes encoded in the genome of Beroe ovata, a ctenophore-eating ctenophore.

As the sister group to all other animals, ctenophores (comb jellies) are important for understanding the emergence and diversification of numerous animal traits. Efforts to explore the evolutionary processes that promoted diversification within Ctenophora are hindered by undersampling genomic diversity within this clade. To address this gap, we present the sequence, assembly and initial annotation of the genome of Beroe ovata. Beroe possess unique morphology, behavior, ecology and development. Unlike their generalist carnivorous kin, beroid ctenophores feed exclusively on other ctenophores. Accordingly, our analyses revealed a loss of chitinase, an enzyme critical for the digestion of most non-ctenophore prey, but superfluous for ctenophorivores. Broadly, our genomic analysis revealed that extensive gene loss and changes in gene regulation have shaped the unique biology of B. ovata. Despite the gene losses in B. ovata, our phylogenetic analyses on photosensitive opsins and several early developmental regulatory genes show that these genes are conserved in B. ovata. This additional sampling contributes to a more complete reconstruction of the ctenophore ancestor and points to the need for extensive comparisons within this ancient and diverse clade of animals. To promote further exploration of these data, we present BovaDB (http://ryanlab.whitney.ufl.edu/bovadb/), a portal for the B. ovata genome.

Autonomous differentiation of transgenic cells requiring no external hormone application: the endogenous gene expression and phytohormone behaviors.

The ectopic overexpression of developmental regulator (DR) genes has been reported to improve the transformation in recalcitrant plant species because of the promotion of cellular differentiation during cell culture processes. In other words, the external plant growth regulator (PGR) application during the tissue and cell culture process is still required in cases utilizing DR genes for plant regeneration. Here, the effect of Arabidopsis BABY BOOM (BBM) and WUSCHEL (WUS) on the differentiation of tobacco transgenic cells was examined. We found that the SRDX fusion to WUS, when co-expressed with the BBM-VP16 fusion gene, significantly influenced the induction of autonomous differentiation under PGR-free culture conditions, with similar effects in some other plant species. Furthermore, to understand the endogenous background underlying cell differentiation toward regeneration, phytohormone and RNA-seq analyses were performed using tobacco leaf explants in which transgenic cells were autonomously differentiating. The levels of active auxins, cytokinins, abscisic acid, and inactive gibberellins increased as cell differentiation proceeded toward organogenesis. Gene Ontology terms related to phytohormones and organogenesis were identified as differentially expressed genes, in addition to those related to polysaccharide and nitrate metabolism. The qRT-PCR four selected genes as DEGs supported the RNA-seq data. This differentiation induction system and the reported phytohormone and transcript profiles provide a foundation for the development of PGR-free tissue cultures of various plant species, facilitating future biotechnological breeding.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.