Dynamic changes in mRNA nucleocytoplasmic localization in the nitrate response of Arabidopsis roots.

Nitrate is a nutrient and signal that regulates gene expression. The nitrate response has been extensively characterized at the organism, organ, and cell-type-specific levels, but intracellular mRNA dynamics remain unexplored. To characterize nuclear and cytoplasmic transcriptome dynamics in response to nitrate, we performed a time-course expression analysis after nitrate treatment in isolated nuclei, cytoplasm, and whole roots. We identified 402 differentially localized transcripts (DLTs) in response to nitrate treatment. Induced DLT genes showed rapid and transient recruitment of the RNA polymerase II, together with an increase in the mRNA turnover rates. DLTs code for genes involved in metabolic processes, localization, and response to stimulus indicating DLTs include genes with relevant functions for the nitrate response that have not been previously identified. Using single-molecule RNA FISH, we observed early nuclear accumulation of the NITRATE REDUCTASE 1 (NIA1) transcripts in their transcription sites. We found that transcription of NIA1, a gene showing delayed cytoplasmic accumulation, is rapidly and transiently activated; however, its transcripts become unstable when they reach the cytoplasm. Our study reveals the dynamic localization of mRNAs between the nucleus and cytoplasm as an emerging feature in the temporal control of gene expression in response to nitrate treatment in Arabidopsis roots.

Widespread changes to the translational landscape in a maize microRNA biogenesis mutant.

MicroRNAs are short, non-coding RNAs that repress gene expression in both plants and animals and have diverse functions related to growth, development, and stress responses. The ribonuclease, DICER-LIKE1 (DCL1) is required for two steps in plant miRNA biogenesis: cleavage of the primary miRNAs (pri-miRNAs) to release a hairpin structure, called the precursor miRNA (pre-miRNA) and cleavage of the pre-miRNA to generate the miRNA/miRNA* duplex. The mature miRNA guides the RNA-induced silencing complex to target RNAs with complementary sequences, resulting in translational repression and/or RNA cleavage of target mRNAs. However, the relative contribution of translational repression versus mRNA degradation by miRNAs remains unknown at the genome-level in crops, especially in maize. The maize fuzzy tassel (fzt) mutant contains a hypomorphic mutation in DCL1 resulting in broad developmental defects. While most miRNAs are reduced in fzt, the levels of miRNA-targeted mRNAs are not dramatically increased, suggesting that translational regulation by miRNAs may be common. To gain insight into the repression mechanism of plant miRNAs, we combined ribosome profiling and RNA-sequencing to globally survey miRNA activities in maize. Our data indicate that translational repression contributes significantly to regulation of most miRNA targets and that approximately one-third of miRNA targets are regulated primarily at the translational level. Surprisingly, ribosomes appear altered in fzt mutants suggesting that DCL1 may also have a role in ribosome biogenesis. Thus, DICER-LIKE1 shapes the translational landscape in plants through both miRNA-dependent and miRNA-independent mechanisms.

Small-molecule inhibition of the METTL3/METTL14 complex suppresses neuroblastoma tumor growth and promotes differentiation.

The N6-methyladenosine (m6A) RNA modification is an important regulator of gene expression. m6A is deposited by a methyltransferase complex that includes methyltransferase-like 3 (METTL3) and methyltransferase-like 14 (METTL14). High levels of METTL3/METTL14 drive the growth of many types of adult cancer, and METTL3/METTL14 inhibitors are emerging as new anticancer agents. However, little is known about the m6A epitranscriptome or the role of the METTL3/METTL14 complex in neuroblastoma, a common pediatric cancer. Here, we show that METTL3 knockdown or pharmacologic inhibition with the small molecule STM2457 leads to reduced neuroblastoma cell proliferation and increased differentiation. These changes in neuroblastoma phenotype are associated with decreased m6A deposition on transcripts involved in nervous system development and neuronal differentiation, with increased stability of target mRNAs. In preclinical studies, STM2457 treatment suppresses the growth of neuroblastoma tumors in vivo. Together, these results support the potential of METTL3/METTL14 complex inhibition as a therapeutic strategy against neuroblastoma.

Genome-Wide CRISPR Screen Identifies an NF2-Adherens Junction Mechanistic Dependency for Cardiac Lineage.

BACKGROUND: Cardiomyocyte differentiation involves a stepwise clearance of repressors and fate-restricting regulators through the modulation of BMP (bone morphogenic protein)/Wnt-signaling pathways. However, the mechanisms and how regulatory roadblocks are removed with specific developmental signaling pathways remain unclear. METHODS: We conducted a genome-wide CRISPR screen to uncover essential regulators of cardiomyocyte specification in human embryonic stem cells using a myosin heavy chain 6 (MYH6)-GFP (green fluorescence protein) reporter system. After an independent secondary single guide ribonucleic acid validation of 25 candidates, we identified NF2 (neurofibromin 2), a moesin-ezrin-radixin like (MERLIN) tumor suppressor, as an upstream driver of early cardiomyocyte lineage specification. Independent monoclonal NF2 knockouts were generated using CRISPR-Cas9, and cell states were inferred through bulk RNA sequencing and protein expression analysis across differentiation time points. Terminal lineage differentiation was assessed by using an in vitro 2-dimensional-micropatterned gastruloid model, trilineage differentiation, and cardiomyocyte differentiation. Protein interaction and post-translation modification of NF2 with its interacting partners were assessed using site-directed mutagenesis, coimmunoprecipitation, and proximity ligation assays. RESULTS: Transcriptional regulation and trajectory inference from NF2-null cells reveal the loss of cardiomyocyte identity and the acquisition of nonmesodermal identity. Sustained elevation of early mesoderm lineage repressor SOX2 and upregulation of late anticardiac regulators CDX2 and MSX1 in NF2 knockout cells reflect a necessary role for NF2 in removing regulatory roadblocks. Furthermore, we found that NF2 and AMOT (angiomotin) cooperatively bind to YAP (yes-associated protein) during mesendoderm formation, thereby preventing YAP activation, independent of canonical MST (mammalian sterile 20-like serine-threonine protein kinase)-LATS (large tumor suppressor serine-threonine protein kinase) signaling. Mechanistically, cardiomyocyte lineage identity was rescued by wild-type and NF2 serine-518 phosphomutants, but not NF2 FERM (ezrin-radixin-meosin homology protein) domain blue-box mutants, demonstrating that the critical FERM domain-dependent formation of the AMOT-NF2-YAP scaffold complex at the adherens junction is required for early cardiomyocyte lineage differentiation. CONCLUSIONS: These results provide mechanistic insight into the essential role of NF2 during early epithelial-mesenchymal transition by sequestering the repressive effect of YAP and relieving regulatory roadblocks en route to cardiomyocytes.

Quantitative RUBY reporter assay for gene regulation analysis.

Various reporter genes have been developed to study gene expression pattern and gene regulation. The RUBY reporter gene was recently developed and widely used, because of its visible and noninvasive advantages. However, quantitative analysis of RUBY gene expression levels was lacking. In this study, we introduce a novel betalain quantification method in combination with the tobacco transient expression system. The betalain produced in tobacco leaves was extracted and purified, and its concentration was quantitatively measured. We successfully applied this approach in studying the transcriptional regulation of ARC5 gene by transcription factors CPD25 and CPD45. Furthermore, with this method, we showed that the gene expression of RCA and Rbcs1A gene were regulated by light, transcription factors HY5 and PIFs through G-box and I-box elements. The development of this betalain quantification approach with the tobacco transient expression system offers a cost-effective and intuitive strategy for studying the regulatory mechanism of gene expression.

Tandem gene duplication selected by activation of horizontally transferred gene in bacteria.

Horizontal gene transfer occurs frequently in bacteria, but the mechanism driving activation and optimization of the expression of horizontally transferred genes (HTGs) in new recipient strains is not clear. Our previous study found that spontaneous tandem DNA duplication resulted in rapid activation of HTGs. Here, we took advantage of this finding to develop a novel technique for tandem gene duplication, named tandem gene duplication selected by activation of horizontally transferred gene in bacteria (TDAH), in which tandem duplication was selected by the activation of horizontally transferred selectable marker gene. TDAH construction does not contain any reported functional elements based on homologous or site-specific recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection and 9-bp-microhomology border sequences for precise illegitimate recombination. One transformation and 3 days were enough to produce a high-copy strain, so its procedure is simple and fast. Without subsequent knockout of the endogenous recombination system, TDAH could also generate the relatively stable high-copy tandem duplication for plasmid-carried and genome-integrated DNA. TDAH also showed an excellent capacity for increase gene expression and worked well in different industrial bacteria. We also applied TDAH to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5 times and remained stable even after 12 subcultures. TDAH is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. KEY POINTS: • We develop a novel and efficient technique (TDAH) for tandem gene duplication in bacterium. TDAH is based on the mechanism of HTG rapid activation. TDAH does not contain any reported functional elements based on homologous recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection, so its construction and procedure are very simple and fast. • TDAH is the first reported selected and stable tandem-gene-duplication technique in which the selected high-copy plasmid-carried and genome-integrated DNA could remain stable without the subsequent knockout of recombination system. • TDAH showed an excellent capacity for regulating gene expression and worked well in different industrial bacteria, indicating it is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. • TDAH was applied to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5-fold and remained stable even after 12 subcultures.

Size-dependent enhancement of gene expression by Plasmodium 5'UTR introns.

BACKGROUND: Eukaryotic genes contain introns that are removed by the spliceosomal machinery during mRNA maturation. Introns impose a huge energetic burden on a cell; therefore, they must play an essential role in maintaining genome stability and/or regulating gene expression. Many genes (> 50%) in Plasmodium parasites contain predicted introns, including introns in 5' and 3' untranslated regions (UTR). However, the roles of UTR introns in the gene expression of malaria parasites remain unknown. METHODS: In this study, an episomal dual-luciferase assay was developed to evaluate gene expression driven by promoters with or without a 5'UTR intron from four Plasmodium yoelii genes. To investigate the effect of the 5'UTR intron on endogenous gene expression, the pytctp gene was tagged with 3xHA at the N-terminal of the coding region, and parasites with or without the 5'UTR intron were generated using the CRISPR/Cas9 system. RESULTS: We showed that promoters with 5'UTR introns had higher activities in driving gene expression than those without 5'UTR introns. The results were confirmed in recombinant parasites expressing an HA-tagged gene (pytctp) driven by promoter with or without 5'UTR intron. The enhancement of gene expression was intron size dependent, but not the DNA sequence, e.g. the longer the intron, the higher levels of expression. Similar results were observed when a promoter from one strain of P. yoelii was introduced into different parasite strains. Finally, the 5'UTR introns were alternatively spliced in different parasite development stages, suggesting an active mechanism employed by the parasites to regulate gene expression in various developmental stages. CONCLUSIONS: Plasmodium 5'UTR introns enhance gene expression in a size-dependent manner; the presence of alternatively spliced mRNAs in different parasite developmental stages suggests that alternative slicing of 5'UTR introns is one of the key mechanisms in regulating parasite gene expression and differentiation.

Regulatory networks of senescence-associated gene-transcription factors promote degradation in Moso bamboo shoots.

Bamboo cultivation, particularly Moso bamboo (Phyllostachys edulis), holds significant economic importance in various regions worldwide. Bamboo shoot degradation (BSD) severely affects productivity and economic viability. However, despite its agricultural consequences, the molecular mechanisms underlying BSD remain unclear. Consequently, we explored the dynamic changes of BSD through anatomy, physiology and the transcriptome. Our findings reveal ruptured protoxylem cells, reduced cell wall thickness and the accumulation of sucrose and reactive oxygen species (ROS) during BSD. Transcriptomic analysis underscored the importance of genes related to plant hormone signal transduction, sugar metabolism and ROS homoeostasis in this process. Furthermore, BSD appears to be driven by the coexpression regulatory network of senescence-associated gene transcription factors (SAG-TFs), specifically PeSAG39, PeWRKY22 and PeWRKY75, primarily located in the protoxylem of vascular bundles. Yeast one-hybrid and dual-luciferase assays demonstrated that PeWRKY22 and PeWRKY75 activate PeSAG39 expression by binding to its promoter. This study advanced our understanding of the molecular regulatory mechanisms governing BSD, offering a valuable reference for enhancing Moso bamboo forest productivity.

Fresh Insights Into SLC25A26: Potential New Therapeutic Target for Cancers: A Review.

SLC25A26 is the only known human mitochondrial S-adenosylmethionine carrier encoding gene. Recent studies have shown that SLC25A26 is abnormally expressed in some cancers, such as cervical cancer, low-grade glioma, non-small cell lung cancer, and liver cancer, which suggests SLC25A26 can affect the occurrence and development of some cancers. This article in brief briefly reviewed mitochondrial S-adenosylmethionine carrier in different species and its encoding gene, focused on the association of SLC25A26 aberrant expression and some cancers as well as potential mechanisms, summarized its potential for cancer prognosis, and characteristics of mitochondrial diseases caused by SLC25A26 mutation. Finally, we provide a brief expectation that needs to be further investigated. We speculate that SLC25A26 will be a potential new therapeutic target for some cancers.

Gene-environmental regulation of the postnatal post-mitotic neuronal maturation.

Embryonic neurodevelopment, particularly neural progenitor differentiation into post-mitotic neurons, has been extensively studied. While the number and composition of post-mitotic neurons remain relatively constant from birth to adulthood, the brain undergoes significant postnatal maturation marked by major property changes frequently disrupted in neural diseases. This review first summarizes recent characterizations of the functional and molecular maturation of the postnatal nervous system. We then review regulatory mechanisms controlling the precise gene expression changes crucial for the intricate sequence of maturation events, highlighting experience-dependent versus cell-intrinsic genetic timer mechanisms. Despite significant advances in understanding of the gene-environmental regulation of postnatal neuronal maturation, many aspects remain unknown. The review concludes with our perspective on exciting future research directions in the next decade.

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data.

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

Transcription Factor Binding Site in Promoter Determines the Pattern of Plasmid-Based Transgene Expression In Vivo.

Understanding the regulation of transgene expression is critical for the success of plasmid-based gene therapy and vaccine development. In this study, we used two sets of plasmid vectors containing secreted embryonic alkaline phosphatase or the mouse IL-10 gene as a reporter and investigated the role of promoter elements in regulating transgene expression in vivo. We demonstrated in mice that hydrodynamic transfer of plasmids with the CMV promoter resulted in a high level of reporter gene expression that declined rapidly over time. In contrast, when plasmids with albumin promoters were used, a lower but sustained gene expression pattern was observed. We also found that plasmids containing a shorter CMV promoter sequence with fewer transcription factor binding sites showed a decrease in the peak level of gene expression without changing the overall pattern of reporter gene expression. The replacement of regulatory elements in the CMV promoter with a single regulatory element of the albumin promoter changed the pattern of transient gene expression seen in the CMV promoter to a pattern of sustained gene expression identical to that of a full albumin promoter. ChIP analyses demonstrated an elevated binding of acetylated histones and TATA box-binding protein to the promoter carrying regulatory elements of the albumin promoter. These results suggest that the strength of a promoter is determined by the number of appropriate transcription factor binding sites, while gene expression persistence is determined by the presence of regulatory elements capable of recruiting epigenetic modifying complexes that make the promoter accessible for transcription. This study provides important insights into the mechanisms underlying gene expression regulation in vivo, which can be used to improve plasmid-based gene therapy and vaccine development.

Partial root-zone drying combined with nitrogen treatments mitigates drought responses in rice.

Drought is a major stress affecting rice yields. Combining partial root-zone drying (PRD) and different nitrogen fertilizers reduces the damage caused by water stress in rice. However, the underlying molecular mechanisms remain unclear. In this study, we combined treatments with PRD and ammonia:nitrate nitrogen at 0:100 (PRD0:100) and 50:50 (PRD50:50) ratios or PEG and nitrate nitrogen at 0:100 (PEG0:100) ratios in rice. Physiological, transcriptomic, and metabolomic analyses were performed on rice leaves to identify key genes involved in water stress tolerance under different nitrogen forms and PRD pretreatments. Our results indicated that, in contrast to PRD0:100, PRD50:50 elevated the superoxide dismutase activity in leaves to accelerate the scavenging of ROS accumulated by osmotic stress, attenuated the degree of membrane lipid peroxidation, stabilized photosynthesis, and elevated the relative water content of leaves to alleviate the drought-induced osmotic stress. Moreover, the alleviation ability was better under PRD50:50 treatment than under PRD0:100. Integrated transcriptome and metabolome analyses of PRD0:100 vs PRD50:50 revealed that the differences in PRD involvement in water stress tolerance under different nitrogen pretreatments were mainly in photosynthesis, oxidative stress, nitrogen metabolism process, phytohormone signaling, and biosynthesis of other secondary metabolites. Some key genes may play an important role in these pathways, including OsGRX4, OsNDPK2, OsGS1;1, OsNR1.2, OsSUS7, and YGL8. Thus, the osmotic stress tolerance mediated by PRD and nitrogen cotreatment is influenced by different nitrogen forms. Our results provide new insights into osmotic stress tolerance mediated by PRD and nitrogen cotreatment, demonstrate the essential role of nitrogen morphology in PRD-induced molecular regulation, and identify genes that contribute to further improving stress tolerance in rice.

Ready for renascence in mosquito: The regulation of gene expression in Plasmodium sexual development.

Malaria remains one of the most perilous vector-borne infectious diseases for humans globally. Sexual gametocyte represents the exclusive stage at which malaria parasites are transmitted from the vertebrate to the Anopheles host. The feasible and effective approach to prevent malaria transmission is by addressing the sexual developmental processes, that is, gametocytogenesis and gametogenesis. Thus, this review will comprehensively cover advances in the regulation of gene expression surrounding the transmissible stages, including epigenetic, transcriptional, and post-transcriptional control.

Genome-wide identification, expression and function analysis of the MTP gene family in tulip (Tulipa gesneriana).

Currently, soil heavy metal contamination is a severe issue, particularly with Cd pollution. The metal tolerance protein (MTP) proteins, as plant divalent cation transporters, play a crucial role in the transport and tolerance of heavy metals in plants. This study conducted comprehensive identification and characterization of the MTP gene family in the tulip. A total of 11 TgMTP genes were identified and phylogenetically classified into three subfamilies. Conserved motif and gene structure analyses unveiled commonalities and variations among subfamily members. Expression profiling demonstrated several TgMTPs were markedly upregulated under Cd exposure, including the TgMTP7.1. Heterologous expression in yeast validated that TgMTP7.1 could ameliorate Cd sensitivity and enhance its tolerance. These results provide primary insights into the MTP gene family in tulip. Phylogenetic relationships and functional analyses establish a framework for elucidating the transporters and molecular mechanisms governing Cd accumulation and distribution in tulip. Key TgMTPs identified, exemplified by TgMTP7.1, may illuminate molecular breeding efforts aimed at developing Cd-tolerant cultivars for the remediation of soil Cd contamination.

Expression of Staphylococcus aureus translation elongation factor P is regulated by a stress-inducible promotor.

Translation elongation factor P, expressed by the efp gene, is a conserved protein closely related to bacterial virulence and environmental stress regulation responses, however, little is known about the efp gene expression regulations. Here, the strain of Staphylococcus aureus subsp. aureus NCTC 8325 was taken as the research object and cultured under different conditions, including different culture temperatures, pH, and antibiotics, to study the expression of the efp gene in S. aureus by qRT-PCR, the results showed that the expression of the efp gene is upregulated under high temperature (40 °C), acidic (pH 5.4) or alkaline (pH 9.4) culture conditions, but upregulated early and downregulated later under the conditions of 0.5 MIC antibiotics (chloramphenicol at the final concentration of 2 µg/mL and vancomycin at the final concentration of 0.25 µg/mL), indicating that the efp promoter in S. aureus is inducible. The efp promoter sequence and structure in S. aureus were predicted by bioinformatics methods, and the predicted promoter was validated by constructing a promoter-probe vector and a series of promoter mutants, the results showed that the efp promoter sequence in S. aureus, named Pro, located in 1,548,179-1,548,250 of the S. aureus genome (NC_007795.1), and the sequence of - 10 element is CCTTATAGT, - 35 element is TTTACT. The results above could lay a foundation for screening transcription factors involved in the expression of the efp gene and then exploring the transcriptional regulation mechanism of EF-P in S. aureus.

Exploring SVA Insertion Polymorphisms in Shaping Differential Gene Expressions in the Central Nervous System.

Transposable elements (TEs) are repetitive elements which make up around 45% of the human genome. A class of TEs, known as SINE-VNTR-Alu (SVA), demonstrate the capacity to mobilise throughout the genome, resulting in SVA polymorphisms for their presence or absence within the population. Although studies have previously highlighted the involvement of TEs within neurodegenerative diseases, such as Parkinson's disease and amyotrophic lateral sclerosis (ALS), the exact mechanism has yet to be identified. In this study, we used whole-genome sequencing and RNA sequencing data of ALS patients and healthy controls from the New York Genome Centre ALS Consortium to elucidate the influence of reference SVA elements on gene expressions genome-wide within central nervous system (CNS) tissues. To investigate this, we applied a matrix expression quantitative trait loci analysis and demonstrate that reference SVA insertion polymorphisms can significantly modulate the expression of numerous genes, preferentially in the trans position and in a tissue-specific manner. We also highlight that SVAs significantly regulate mitochondrial genes as well as genes within the HLA and MAPT loci, previously associated within neurodegenerative diseases. In conclusion, this study continues to bring to light the effects of polymorphic SVAs on gene regulation and further highlights the importance of TEs within disease pathology.

The SATB1-MIR22-GBA axis mediates glucocerebroside accumulation inducing a cellular senescence-like phenotype in dopaminergic neurons.

Idiopathic Parkinson's disease (PD) is characterized by the loss of dopaminergic neurons in the substantia nigra pars compacta, which is associated with neuroinflammation and reactive gliosis. The underlying cause of PD and the concurrent neuroinflammation are not well understood. In this study, we utilize human and murine neuronal lines, stem cell-derived dopaminergic neurons, and mice to demonstrate that three previously identified genetic risk factors for PD, namely SATB1, MIR22HG, and GBA, are components of a single gene regulatory pathway. Our findings indicate that dysregulation of this pathway leads to the upregulation of glucocerebrosides (GluCer), which triggers a cellular senescence-like phenotype in dopaminergic neurons. Specifically, we discovered that downregulation of the transcriptional repressor SATB1 results in the derepression of the microRNA miR-22-3p, leading to decreased GBA expression and subsequent accumulation of GluCer. Furthermore, our results demonstrate that an increase in GluCer alone is sufficient to impair lysosomal and mitochondrial function, thereby inducing cellular senescence. Dysregulation of the SATB1-MIR22-GBA pathway, observed in both PD patients and normal aging, leads to lysosomal and mitochondrial dysfunction due to the GluCer accumulation, ultimately resulting in a cellular senescence-like phenotype in dopaminergic neurons. Therefore, our study highlights a novel pathway involving three genetic risk factors for PD and provides a potential mechanism for the senescence-induced neuroinflammation and reactive gliosis observed in both PD and normal aging.

Clonally heritable gene expression imparts a layer of diversity within cell types.

Cell types can be classified according to shared patterns of transcription. Non-genetic variability among individual cells of the same type has been ascribed to stochastic transcriptional bursting and transient cell states. Using high-coverage single-cell RNA profiling, we asked whether long-term, heritable differences in gene expression can impart diversity within cells of the same type. Studying clonal human lymphocytes and mouse brain cells, we uncovered a vast diversity of heritable gene expression patterns among different clones of cells of the same type in vivo. We combined chromatin accessibility and RNA profiling on different lymphocyte clones to reveal thousands of regulatory regions exhibiting interclonal variation, which could be directly linked to interclonal variation in gene expression. Our findings identify a source of cellular diversity, which may have important implications for how cellular populations are shaped by selective processes in development, aging, and disease. A record of this paper's transparent peer review process is included in the supplemental information.