Gene Regulatory Regions Research Articles

Transcriptome and DNA methylome divergence of inflorescence development between 2 ecotypes in Panicum hallii.

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the 2 major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMR-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

Abstract 4744: Epigenetic reprogramming induced by lipids fosters mammary cell plasticity in non-transformed breast epithelial cells

Abstract Introduction. Lipid exposure of non-transformed breast epithelial cells results in increased flux modifications, and profound changes in gene expression. The upregulated genes are involved in stemness, neural development, and neural crest pathways. Neural genes are highly expressed in an estrogen receptor negative breast cancer (ERnegBC): the triple negative subtype. We aim to identify mechanisms that link lipids and epigenetic reprogramming to the genesis of ERnegBC. Methods. Non-transformed MCF-10A cells exposed to octanoic acid (OA) were utilized for U13C-glucose tracing. S-adenosylmethionine (SAM), 2-HG and other metabolites were analyzed following treatment with OA ± PHGDH inhibitor. CUT&RUN for H3K4me3 and H3K27me3 was performed, and genes affected by OA (PMID: 28263391) were compared with OA-responsive peaks. OA-induced genes with CUT&RUN peaks were measured in breast organoids derived from reduction mammoplasty samples exposed to OA. The Aldefluour assay was used to identify stem-like (ALDH+) cells. Results. U13C-glucose tracing in presence of OA showed decreased glucose uptake, glycolysis and pentose phosphate pathway flux, with little effect on the TCA cycle. One-carbon (1C)-THF was redirected to the methionine cycle increasing flux to methylation. The methyl donor SAM and the demethylase inhibitor 2-HG increased after 15- and 30-min OA exposure, respectively; PHGDH inhibitor blocked these increases likely via effects on 1C metabolism. H3K4me3 CUT&RUN revealed 661 differential peaks (FDR < 0.05) comparing OA to control. The genes associated with these peaks are involved in neuron differentiation, neural crest, EMT and neural genes expressed in ERnegBC. 73% of H3K4me3 OA-associated peaks were in regulatory regions of OA-induced genes (FDR < 0.01), including NGFR (log2FC = 11.7), NGF (log2FC = 8), NTRK1 (log2FC = 2.1). OA-induced genes with enriched H3K4me3 peaks were also upregulated in breast organoids exposed to OA, such as NGFR (FC = 2.2), NGF (FC = 4.7), NTRK1 (FC = 6.6). Motif analysis revealed an overrepresentation of transcription factors (p < 0.05) associated with EMT (Zeb1, Slug), neural functions (E2A, AP1, JunB), neuronal-injury (Atf3) and stress response (Chop, Atf4). Twelve H3K27me3 peaks were enriched in control (FDR < 0.05) and associated with increased gene expression in OA, among them were LGR6 (log2FC = 1.9), a gene associated with ERnegBC and PLAG1 (log2FC = 2.8) a stem cell marker. ALDH+ cells’ percentage increased by at least 10% after OA exposure. Conclusions. Lipid exposure affects the production of SAM and 2-HG, which results in epigenetic fostered plasticity, leading to reprogramming/selecting cells with a multi-potential embryonic or stem-like state and/or differentiation to a neural/neural crest-like state which may facilitate ERnegBC genesis. Citation Format: Mariana Bustamante Eduardo, Gannon Cottone, Shiyu Liu, Seema Khan, Susan Clare. Epigenetic reprogramming induced by lipids fosters mammary cell plasticity in non-transformed breast epithelial cells. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4744.

NAGS, CPS1, and SLC25A13 (Citrin) at the Crossroads of Arginine and Pyrimidines Metabolism in Tumor Cells.

Urea cycle enzymes and transporters collectively convert ammonia into urea in the liver. Aberrant overexpression of carbamylphosphate synthetase 1 (CPS1) and SLC25A13 (citrin) genes has been associated with faster proliferation of tumor cells due to metabolic reprogramming that increases the activity of the CAD complex and pyrimidine biosynthesis. N-acetylglutamate (NAG), produced by NAG synthase (NAGS), is an essential activator of CPS1. Although NAGS is expressed in lung cancer derived cell lines, expression of the NAGS gene and its product was not evaluated in tumors with aberrant expression of CPS1 and citrin. We used data mining approaches to identify tumor types that exhibit aberrant overexpression of NAGS, CPS1, and citrin genes, and evaluated factors that may contribute to increased expression of the three genes and their products in tumors. Median expression of NAGS, CPS1, and citrin mRNA was higher in glioblastoma multiforme (GBM), glioma, and stomach adenocarcinoma (STAD) samples compared to the matched normal tissue. Median expression of CPS1 and citrin mRNA was higher in the lung adenocarcinoma (LUAD) sample while expression of NAGS mRNA did not differ. High NAGS expression was associated with an unfavorable outcome in patients with glioblastoma and GBM. Low NAGS expression was associated with an unfavorable outcome in patients with LUAD. Patterns of DNase hypersensitive sites and histone modifications in the upstream regulatory regions of NAGS, CPS1, and citrin genes were similar in liver tissue, lung tissue, and A549 lung adenocarcinoma cells despite different expression levels of the three genes in the liver and lung. Citrin gene copy numbers correlated with its mRNA expression in glioblastoma, GBM, LUAD, and STAD samples. There was little overlap between NAGS, CPS1, and citrin sequence variants found in patients with respective deficiencies, tumor samples, and individuals without known rare genetic diseases. The correlation between NAGS, CPS1, and citrin mRNA expression in the individual glioblastoma, GBM, LUAD, and STAD samples was very weak. These results suggest that the increased cytoplasmic supply of either carbamylphosphate, produced by CPS1, or aspartate may be sufficient to promote tumorigenesis, as well as the need for an alternative explanation of CPS1 activity in the absence of NAGS expression and NAG.

PASTE, Don't Cut: Genome Editing Tool Looks Beyond CRISPR and Prime

PASTE, Don't Cut: Genome Editing Tool Looks Beyond CRISPR and Prime

An in Silico Approach to Identifying TF Binding Sites: Analysis of the Regulatory Regions of BUSCO Genes from Fungal Species in the Ceratocystidaceae Family.

Transcriptional regulation controls gene expression through regulatory promoter regions that contain conserved sequence motifs. These motifs, also known as regulatory elements, are critically important to expression, which is driving research efforts to identify and characterize them. Yeasts have been the focus of such studies in fungi, including in several in silico approaches. This study aimed to determine whether in silico approaches could be used to identify motifs in the Ceratocystidaceae family, and if present, to evaluate whether these correspond to known transcription factors. This study targeted the 1000 base-pair region upstream of the start codon of 20 single-copy genes from the BUSCO dataset for motif discovery. Using the MEME and Tomtom analysis tools, conserved motifs at the family level were identified. The results show that such in silico approaches could identify known regulatory motifs in the Ceratocystidaceae and other unrelated species. This study provides support to ongoing efforts to use in silico analyses for motif discovery.

Abstract P6-11-03: Lipid-rich environment induces epigenetic reprogramming in non-transformed breast epithelial cells enhancing a mammary cell plasticity

Abstract Introduction. The identification of women specifically at risk for estrogen receptor negative breast cancer (ER- BC) and the targeted treatment of this disease are significantly unmet clinical needs. We analyzed the gene expression profiles of epithelial cells from the contralateral unaffected breasts (CUBs) of BC patients and identified a lipid metabolism gene signature enriched in the CUBs of women with ER- BC (PMID: 28263391). Subsequently, we observed that exposure of non-transformed breast epithelial cells to lipids results in significant changes in gene expression and histone posttranslational modifications; and increased flux through multiple metabolic reactions, including those of serine and methionine (PMID: 35508495). Interestingly, the upregulated genes are involved in neural development, axon guidance, neural crest pathways and stemness. Neural genes are highly expressed in Triple Negative Breast Cancers (TNBCs) especially in the C2 cluster. Given that mammary stem/progenitor cells have distinct metabolic properties compared to other mammary cell subsets, we hypothesized that upon lipid exposure, stem-like cells have a survival advantage, and that lipid induces epigenetic reprogramming into neural-like cells which may foster a malignant transformation. Methods. To interrogate potential mechanisms linking lipids and epigenetic reprogramming, we performed CUT&RUN for H3K4me3 and H3K27me3 in non-transformed, estrogen and progesterone receptor negative MCF-10A cells cells exposed to vehicle or octanoic acid (OA) for 24 hrs. Peaks were called using MACS2 and differential peaks identified with DiffBind. Differentially expressed genes affected by OA (PMID: 28263391) were compared with target genes from the CUT&RUN. To determine the contribution in OA exposed cells of serine and methionine metabolism to S-adenosyl methionine (SAM), 13C-glucose tracing was performed. The Aldefluour assay was used to identify stem-like (ALDH+) cells in lipid-exposed MCF-10A cells. To determine if lipid-exposed cells adopt a neural-like phenotype, MCF-10A cells were grown on Poly-D-Lysine/Laminin (PDL/LM) coated plates. Results. A total of 661 differential peaks were identified for H3K4me3 (FDR < 0.05) encompassing TNBC C2 markers (NRTN, CHRM3) and genes involved in neuron differentiation, axonogenesis (NGFR, NGF, FOXA2), neural crest migration (NTRK2, MMP2) and EMT (DLL4, MMP15). Approximately 73% of H3K4me3 OA-associated peaks encompass regulatory regions of genes that experienced a significant increase in gene expression with OA exposure (FDR < 0.01), including NGFR (log2 FC = 11.7), FOXA2 (log2 FC = 11.6), NGF (log2 FC = 8). Twelve H3K27me3 peaks were significantly enriched in vehicle (FDR < 0.05) and associated with increased gene expression in OA, among them were the stem cell markers LGR6 (log2 FC = 1.9) and PLAG1 (log2 FC = 2.8). 13C stable isotope tracing revealed that in presence of OA, glucose contributes to increased fractional abundance of the SAM M+1 isotopologue (adj p = 0.003) indicating that carbons derived from the serine synthesis pathway are used for re-methylation of homocysteine to methionine. Vehicle treated cells growing on PDL/LM plates assumed an epithelial phenotype while OA-treated cells adopted a neurite-like phenotype. Upon OA treatment the percentage of ALDH+ cells increased by a minimum of 10%. Conclusions. The increase of fractional abundance of SAM and the upregulation of neural genes regulated by H3K4me3 as well as the enrichment of ALDH+ cells and the development of a neural-like phenotype in a rich lipid environment, suggests that lipid exposure affects the production of SAM, which results in epigenetic fostered plasticity, leading to reprogramming/selecting cells with a multi-potential embryonic or stem-like state and/or differentiation to a neural/neural crest-like state. Citation Format: Mariana Bustamante Eduardo, Gannon Cottone, Seema Khan, Susan Clare. Lipid-rich environment induces epigenetic reprogramming in non-transformed breast epithelial cells enhancing a mammary cell plasticity [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P6-11-03.

The Histone Variant H2A.Z C-Terminal Domain Has Locus-Specific Differential Effects on H2A.Z Occupancy and Nucleosome Localization.

The incorporation of histone variant H2A.Z into nucleosomes creates specialized chromatin domains that regulate DNA-templated processes, such as gene transcription. In Saccharomyces cerevisiae, the diverging H2A.Z C terminus is thought to provide the H2A.Z exclusive functions. To elucidate the roles of this H2A.Z C terminus genome-wide, we used derivatives in which the C terminus was replaced with the corresponding region of H2A (ZA protein), or the H2A region plus a transcriptional activating peptide (ZA-rII'), with the intent of regenerating the H2A.Z-dependent regulation globally. The distribution of these H2A.Z derivatives indicates that the H2A.Z C-terminal region is crucial for both maintaining the occupation level of H2A.Z and the proper positioning of targeted nucleosomes. Interestingly, the specific contribution on incorporation efficiency versus nucleosome positioning varies enormously depending on the locus analyzed. Specifically, the role of H2A.Z in global transcription regulation relies on its C-terminal region. Remarkably, however, this mostly involves genes without a H2A.Z nucleosome in the promoter. Lastly, we demonstrate that the main chaperone complex which deposits H2A.Z to gene regulatory region (SWR1-C) is necessary to localize all H2A.Z derivatives at their specific loci, indicating that the differential association of these derivatives is not due to impaired interaction with SWR1-C. IMPORTANCE We provide evidence that the Saccharomyces cerevisiae C-terminal region of histone variant H2A.Z can mediate its special function in performing gene regulation by interacting with effector proteins and chaperones. These functional interactions allow H2A.Z not only to incorporate to very specific gene regulatory regions, but also to facilitate the gene expression process. To achieve this, we used a chimeric protein which lacks the native H2A.Z C-terminal region but contains an acidic activating region, a module that is known to interact with components of chromatin-remodeling entities and/or transcription modulators. We reasoned that because this activating region can fulfill the role of the H2A.Z C-terminal region, at least in part, the role of the latter would be to interact with these activating region targets.

The association between genetic variants at 3'-UTR and 5'-URR of HLA-G gene and the clinical outcomes of patients with leukemia receiving hematopoietic stem cell transplantation.

In addition to the classical human leukocyte antigen (HLA) genes, the outcomes of post-hematopoietic stem cell transplantation (HSCT) are associated with human leukocyte antigen (HLA)-related genes and non-HLA genes involved in immune regulation. HLA-G gene plays an important role in immune tolerance, assisting immune escape of tumor cells, and decrease of transplant rejection. In this study, we explored the association of genetic variants at the 3'-untranslated region (3'-UTR) and 5'-upstream regulatory region (5'-URR) of HLA-G gene with the adverse outcomes of patients with leukemia receiving HSCT. The genomic DNAs of 164 patients who had acute leukemia and received HSCT were collected for analysis. Nine single nucleotide polymorphisms (SNPs) and six haplotypes in the 3'-UTR and 27 SNPs and 6 haplotypes in the 5'-URR were selected to investigate their relationship with the development of adverse outcomes for patients receiving HSCT, including mortality, relapse, and graft-versus-host disease. Our results revealed that two SNPs (rs371194629 and rs9380142) and one haplotype (UTR-3) located in the 3'-UTR and two SNPs (rs3823321 and rs1736934) and one haplotype (G0104a) located in the 5'-URR of HLA-G were associated with the occurrence of chronic GVHD or development of any forms of GVHD. No SNP was found to associate with the occurrence of mortality and relapse for patients receiving HSCT. These SNPs and haplotypes may play important roles in regulating immune tolerance of allografts post-HSCT that can be used to predict the risk of poor outcomes after receiving HSCT and giving preventive treatment to patients on time.

Methylation of nonessential genes in cutaneous melanoma - Rule Out hypothesis.

Differential methylation plays an important role in melanoma development and is associated with survival, progression and response to treatment. However, the mechanisms by which methylation promotes melanoma development are poorly understood. The traditional explanation of selective advantage provided by differential methylation postulates that hypermethylation of regulatory 5'-cytosine-phosphate-guanine-3' dinucleotides (CpGs) downregulates the expression of tumor suppressor genes and therefore promotes tumorigenesis. We believe that other (not necessarily alternative) explanations of the selective advantages of methylation are also possible. Here, we hypothesize that melanoma cells use methylation to shut down transcription of nonessential genes - those not required for cell survival and proliferation. Suppression of nonessential genes allows tumor cells to be more efficient in terms of energy and resource usage, providing them with a selective advantage over the tumor cells that transcribe and subsequently translate genes they do not need. We named the hypothesis the Rule Out (RO) hypothesis. The RO hypothesis predicts higher methylation of CpGs located in regulatory regions (CpG islands) of nonessential genes. It also predicts the higher methylation of regulatory CpGs linked to nonessential genes in melanomas compared to nevi and lower expression of nonessential genes in malignant (derived from melanoma) versus normal (derived from nonaffected skin) melanocytes. The analyses conducted using in-house and publicly available data found that all predictions derived from the RO hypothesis hold, providing observational support for the hypothesis.

The Effects of Peripubertal THC Exposure in Neurodevelopmental Rat Models of Psychopathology.

Adolescent exposure to cannabinoids as a postnatal environmental insult may increase the risk of psychosis in subjects exposed to perinatal insult, as suggested by the two-hit hypothesis of schizophrenia. Here, we hypothesized that peripubertal Δ9-tetrahydrocannabinol (aTHC) may affect the impact of prenatal methylazoxymethanol acetate (MAM) or perinatal THC (pTHC) exposure in adult rats. We found that MAM and pTHC-exposed rats, when compared to the control group (CNT), were characterized by adult phenotype relevant to schizophrenia, including social withdrawal and cognitive impairment, as revealed by social interaction test and novel object recognition test, respectively. At the molecular level, we observed an increase in cannabinoid CB1 receptor (Cnr1) and/or dopamine D2/D3 receptor (Drd2, Drd3) gene expression in the prefrontal cortex of adult MAM or pTHC-exposed rats, which we attributed to changes in DNA methylation at key regulatory gene regions. Interestingly, aTHC treatment significantly impaired social behavior, but not cognitive performance in CNT groups. In pTHC rats, aTHC did not exacerbate the altered phenotype nor dopaminergic signaling, while it reversed cognitive deficit in MAM rats by modulating Drd2 and Drd3 gene expression. In conclusion, our results suggest that the effects of peripubertal THC exposure may depend on individual differences related to dopaminergic neurotransmission.

Transcription factor CgPOU3F4-like regulates expression of pheomelanin synthesis related gene CgB-aat1 in the Pacific oyster (Crassostrea gigas)

Previous study has found that b (0, +) -type amino acid transporter 1 (CgB-aat1) plays an essential role on mantle pigmentation in the Pacific oyster Crassostrea gigas. However, the molecular regulation of CgB-aat1 gene expression remains unclear. Herein, three POU domain family members, CgPOU2F1, CgPOU3F4-like and CgPOU4F3-X1 were characterized and they all had POUs and HOX domains, respectively, which were important in transcriptional regulation. CgPOU3F4-like gene expression was the highest in mantle edge. Subsequently, the dual-luciferase reporter result showed that the core regulatory region of CgB-aat1 gene was from −632 to −350 bp of promoter. In transient co-transfection assays, the strongest activity was activated only by CgPOU3F4-like, suggesting CgPOU3F4-like was a valid transcriptional activator of CgB-aat1 gene promoter. And the structural integrity of CgPOU3F4-like was essential for its activation function. In addition, site directed mutagenesis assay was applied to detect three key binding sites between CgPOU3F4-like and core region of CgB-aat1 gene promoter, and this interaction was verified by ChIP test. Furthermore, CgPOU3F4-like knockdown by RNA interference led to obvious decreases in CgB-aat1 and cystathionine beta-synthase (CgCbs) expressions at both mRNA and protein levels. Collectively, these results indicate that CgPOU3F4-like positively regulate CgB-aat1 gene expression and it may be a critical upstream transcriptional regulation factor in pheomelanin synthesis in C. gigas.

A functional genomics pipeline to identify high-value asthma and allergy CpGs in the human methylome

DNA methylation of cytosines at cytosine-phosphate-guanine (CpG) dinucleotides (CpGs) is a widespread epigenetic mark, but genome-wide variation has been relatively unexplored due to the limited representation of variable CpGs on commercial high-throughput arrays. To explore this hidden portion of the epigenome, this study combined whole-genome bisulfite sequencing with in silico evidence of gene regulatory regions to design a custom array of high-value CpGs. This study focused on airway epithelial cells from children with and without allergic asthma because these cells mediate the effects of inhaled microbes, pollution, and allergens on asthma and allergic disease risk. This study identified differentially methylated regions from whole-genome bisulfite sequencing in nasal epithelial cell DNA from a total of 39 children with and without allergic asthma of both European and African ancestries. This study selected CpGs from differentially methylated regions, previous allergy or asthma epigenome-wide association studies (EWAS), or genome-wide association study loci, and overlapped them with functional annotations for inclusion on a custom Asthma&Allergy array. This study used both the custom and EPIC arrays to perform EWAS of allergic sensitization (AS) in nasal epithelial cell DNA from children in the URECA (Urban Environment and Childhood Asthma) birth cohort and using the custom array in the INSPIRE [Infant Susceptibility to Pulmonary Infections and Asthma Following RSV Exposure] birth cohort. Each CpG on the arrays was assigned to its nearest gene and its promotor capture Hi-C interacting gene and performed expression quantitative trait methylation (eQTM) studies for both sets of genes. Custom array CpGs were enriched for intermediate methylation levels compared to EPIC CpGs. Intermediate methylation CpGs were further enriched among those associated with AS and for eQTMs on both arrays. This study revealed signature features of high-value CpGs and evidence for epigenetic regulation of genes at AS EWAS loci that are robust to race/ethnicity, ascertainment, age, and geography.

Methylation of Regulatory Regions of DNA Repair Genes in Carotid Atherosclerosis

The status of DNA methylation in the human genome changes during the pathogenesis of common diseases and acts as a predictor of life expectancy. Therefore, it is of interest to investigate the methylation level of regulatory regions of genes responsible for general biological processes that are potentially significant for the development of age-associated diseases. Among them there are genes encoding proteins of DNA repair system, which are characterized by pleiotropic effects. Here, results of the targeted methylation analysis of two regions of the human genome (the promoter of the MLH1 gene and the enhancer near the ATM gene) in different tissues of patients with carotid atherosclerosis are present. Analysis of the methylation profiles of studied genes in various tissues of the same individuals demonstrated marked differences between leukocytes and tissues of the vascular wall. Differences in methylation levels between normal and atherosclerotic tissues of the carotid arteries were revealed only for two studied CpG sites (chr11:108089866 and chr11:108090020, GRCh37/hg19 assembly) in the ATM gene. Based on this, we can assume the involvement of ATM in the development of atherosclerosis. "Overload" of the studied regions with transcription factor binding sites (according to ReMapp2022 data) indicate that the tissue-specific nature of methylation of the regulatory regions of the MLH1 and ATM may be associated with expression levels of these genes in a particular tissue. It has been shown that inter-individual differences in the methylation levels of CpG sites are associated with sufficiently distant nucleotide substitutions.

Genetic Variants of SIRT1 Gene Promoter in Type 2 Diabetes.

Type 2 diabetes (T2D) is a highly heterogeneous and polygenic disease. To date, genetic causes and underlying mechanisms for T2D remain unclear. SIRT1, one member of highly conserved NAD-dependent class III deacetylases, has been implicated in many human diseases. Accumulating evidence indicates that SIRT1 is involved in insulin resistance and impaired pancreatic β-cell function, the two hallmarks of T2D. Thus, we speculated that altered SIRT1 levels, resulting from the genetic variants within its regulatory region of SIRT1 gene, may contribute to the T2D development. In this study, the SIRT1 gene promoter was genetically analyzed in T2D patients (n = 218) and healthy controls (n = 358). A total of 20 genetic variants, including 7 single-nucleotide polymorphisms (SNPs), were identified. Five heterozygous genetic variants (g.4114-15InsA, g.4801G > A, g.4816G > C, g.4934G > T, and g.4963_64Ins17bp) and one SNP (g.4198A > C (rs35706870)) were identified in T2D patients, but in none of the controls. The frequencies of two SNPs (g.4540A > G (rs3740051) (OR: 1.75, 95% CI: 1.24-2.47, P < 0.001 in dominant genetic model) and g.4821G > T (rs35995735)) (OR: 3.58, 95% CI: 1.94-6.60, P < 0.001 in dominant genetic model) were significantly higher in T2D patients. Further association and haplotype analyses confirmed that these two SNPs were strongly linked, contributing to the T2D (OR: 1.442, 95% CI: 1.080-1.927, P < 0.05). Moreover, most of the genetic variants identified in T2D were disease-specific. Taken together, the genetic variants within SIRT1 gene promoter might contribute to the T2D development by altering SIRT1 levels. Underlying molecular mechanism needs to be further explored.

Dr.Nod: computational framework for discovery of regulatory non-coding drivers in tissue-matched distal regulatory elements.

The discovery of cancer driver mutations is a fundamental goal in cancer research. While many cancer driver mutations have been discovered in the protein-coding genome, research into potential cancer drivers in the non-coding regions showed limited success so far. Here, we present a novel comprehensive framework Dr.Nod for detection of non-coding cis-regulatory candidate driver mutations that are associated with dysregulated gene expression using tissue-matched enhancer-gene annotations. Applying the framework to data from over 1500 tumours across eight tissues revealed a 4.4-fold enrichment of candidate driver mutations in regulatory regions of known cancer driver genes. An overarching conclusion that emerges is that the non-coding driver mutations contribute to cancer by significantly altering transcription factor binding sites, leading to upregulation of tissue-matched oncogenes and down-regulation of tumour-suppressor genes. Interestingly, more than half of the detected cancer-promoting non-coding regulatory driver mutations are over 20 kb distant from the cancer-associated genes they regulate. Our results show the importance of tissue-matched enhancer-gene maps, functional impact of mutations, and complex background mutagenesis model for the prediction of non-coding regulatory drivers. In conclusion, our study demonstrates that non-coding mutations in enhancers play a previously underappreciated role in cancer and dysregulation of clinically relevant target genes.

In Silico Dissection of Regulatory Regions of PHT Genes from Saccharum spp. Hybrid and Sorghum bicolor and Expression Analysis of PHT Promoters under Osmotic Stress Conditions in Tobacco

Phosphorus (P) is the second-most essential macronutrient required for the growth and development of plants. It is involved in a number of cellular processes that contribute to the plant’s growth and development. This study investigated Saccharum spp. hybrid and Sorghum bicolor promoter regions of Phosphate transporters (PHT), viz., PHT1, PHT2, PHT3, PHT4, and PHO1, through in silico analysis. The transcription start sites (TSS), conserved motifs, and CpG islands were studied using various computational techniques. The distribution of TSSs indicated the highest promoter prediction scores (1.0). MSh2 and MSb4 were recognized as the common promoter motifs for PHT promoters, found in with 85 to 100% percentage of distribution. The CpG analysis revealed that the promoter regions of most PHT genes had low CpG density, indicating a possible tissue-specific expression. The PHT promoters were investigated for the presence of biotic- and abiotic-stress-associated transcription factor binding sites (TFbs) that revealed the presence of binding motifs for major transcription factors (TFs), namely, AP2/ERF, bHLH, bZIP, MYB, NAC, and WRKY. Therefore, the in-silico analysis of the promoter regions helps us to understand the regulation mechanism of phosphate transporter promoters and gene expression under stress management. The 5′ regulatory region of the EaPHT gene was isolated from Erianthus, a wild relative of the genus Saccharum. The promoter construct was prepared and transformed in tobacco wherein the promoter drove the expression of GUS. Analysis of GUS expression in transgenic tobacco revealed enhanced expression of GUS under salt-stress conditions. This is the first report of the isolation and characterization of a phosphate transporter gene promoter from Erianthus and is expected to be useful for the development of salt-stress transgenic crop plants.

Abstract IA024: Single nuclei chromatin accessibility and transcriptomic map of breast tissues of women of diverse genetic ancestry

Abstract Single nuclei analysis is allowing robust classification of cell types in an organ that helps to establish relationships between cell-type specific gene expression and chromatin accessibility status of gene regulatory regions. Using the institutional resource of breast tissues of healthy donors of various genetic ancestry, we have developed a comprehensive chromatin accessibility (snATAC-seq) and gene expression (snRNA-seq) atlas of human breast tissues. Our analyses included 51,367 nuclei with average sequence coverage of 1195 genes per nuclei. These nuclei were derived from 22 donors of Ashkenazi descent, 20 of European non-Ashkenazi ancestry, 10 of Asian ancestry, 10 Hispanic/Latina, 10 Native American, 6 from BRCA1 mutation carriers, and 5 from BRCA2 mutation carriers. Although tissues from 20 women of African Ancestry were included in sequencing, poor quality sequencing reads prevented inclusion of data from those samples in the final analysis. Integrated analysis revealed 10 distinct cell types in the healthy breast, which included three major epithelial cell subtypes (mature luminal, luminal progenitor, basal), two endothelial subtypes, two adipocyte subtypes, fibroblasts, T-cells, and macrophages. Mature luminal cells could be further divided into two distinct hormone sensitive subtypes, HSa and HSb, with HSa subtype expressing higher levels of Estrogen Receptor alpha (ESR1) and NEK10, a tyrosine kinase that controls p53 activity and limits cell proliferation. The luminal progenitors could be broadly classified into alveolar progenitors (AP) and basal-luminal hybrid progenitors (BL). AP cells expressed higher levels of ELF5, a known marker of alveolar cells, compared to BL cells. Basal cells could be classified into two subtypes with Basal-BAa cells being the most dominant and expressing higher levels of TP63 and NFIB compared to Basal-BAb cells. ESR1 expression pattern was distinctly different in tissues from Native Americans compared to the rest, with a high level of ESR1 expression extending to AP cells. In fact, overall AP cell numbers were ~3-fold higher in Native Americans compared to others (18.9% versus 1.5-8.8%). Furthermore, Ingenuity pathway analysis of differentially expressed genes revealed elevated Estrogen Receptor signaling in ML and AP cell types of Native Americans compared to those of other genetic ancestry. Despite significant differences in expression and activity in Native Americans compared to others, the chromatin accessibility map of the ESR1 gene regulatory regions in ML and LP cells did not show any genetic ancestry dependent variability. In general, cell subtype-specific gene expression did not correlate with chromatin accessibility differences, suggesting that transcriptional regulation independent of chromatin accessibility governs cell type-specific gene expression in the breast. Collectively, these results reveal complexities in gene expression in different cell types of the breast Citation Format: Poornima Bhat-Nakshatri, Duojiao Chen, Aditi Khatpe, Hongyu Gao, Yunlong Liu, Anna Maria Storniolo, Harikrishna Nakshatri. Single nuclei chromatin accessibility and transcriptomic map of breast tissues of women of diverse genetic ancestry [abstract]. In: Proceedings of the 15th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2022 Sep 16-19; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2022;31(1 Suppl):Abstract nr IA024.

FHL3 promotes the formation of fast glycolytic muscle fibers by interacting with YY1 and muscle glycolytic metabolism.

The proportions of the various muscle fiber types are important in the regulation of skeletal muscle metabolism, as well as animal meat production. Four-and-a-half LIM domain protein 3 (FHL3) is highly expressed in fast glycolytic muscle fibers and differentially regulates the expression of myosin heavy chain (MyHC) isoforms at the cellular level. Whether FHL3 regulates the transformation of muscle fiber types in vivo and the regulatory mechanism is unclear. In this study, muscle-specific FHL3 transgenic mice were generated by random integration, and lentivirus-mediated gene knockdown or overexpression in muscles of mice or pigs was conducted. Functional analysis showed that overexpression of FHL3 in muscles significantly increased the proportion of fast-twitch myofibers and muscle mass but decreased muscle succinate dehydrogenase (SDH) activity and whole-body oxygen consumption. Lentivirus-mediated FHL3 knockdown in muscles significantly decreased muscle mass and the proportion of fast-twitch myofibers. Mechanistically, FHL3 directly interacted with the Yin yang 1 (YY1) DNA-binding domain, repressed the binding of YY1 to the fast glycolytic MyHC2b gene regulatory region, and thereby promoted MyHC2b expression. FHL3 also competed with EZH2 to bind the repression domain of YY1 and reduced H3K27me3 enrichment in the MyHC2b regulatory region. Moreover, FHL3 overexpression reduced glucose tolerance by affecting muscle glycolytic metabolism, and its mRNA expression in muscle was positively associated with hemoglobin A1c (HbA1c) in patients with type 2 diabetes. Therefore, FHL3 is a novel potential target gene for the treatment of muscle metabolism-related diseases and improvement of animal meat production.

Aortic Stress Activates an Adaptive Program in Thoracic Aortic Smooth Muscle Cells That Maintains Aortic Strength and Protects Against Aneurysm and Dissection in Mice.

When aortic cells are under stress, such as increased hemodynamic pressure, they adapt to the environment by modifying their functions, allowing the aorta to maintain its strength. To understand the regulation of this adaptive response, we examined transcriptomic and epigenomic programs in aortic smooth muscle cells (SMCs) during the adaptive response to AngII (angiotensin II) infusion and determined its importance in protecting against aortic aneurysm and dissection (AAD). We performed single-cell RNA sequencing and single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) analyses in a mouse model of sporadic AAD induced by AngII infusion. We also examined the direct effects of YAP (yes-associated protein) on the SMC adaptive response in vitro. The role of YAP in AAD development was further evaluated in AngII-infused mice with SMC-specific Yap deletion. In wild-type mice, AngII infusion increased medial thickness in the thoracic aorta. Single-cell RNA sequencing analysis revealed an adaptive response in thoracic SMCs characterized by upregulated genes with roles in wound healing, elastin and collagen production, proliferation, migration, cytoskeleton organization, cell-matrix focal adhesion, and PI3K-PKB/Akt (phosphoinositide-3-kinase-protein kinase B/Akt) and TGF-β (transforming growth factor beta) signaling. ScATAC-seq analysis showed increased chromatin accessibility at regulatory regions of adaptive genes and revealed the mechanical sensor YAP/transcriptional enhanced associate domains as a top candidate transcription complex driving the expression of these genes (eg, Lox, Col5a2, Tgfb2). In cultured human aortic SMCs, cyclic stretch activated YAP, which directly bound to adaptive gene regulatory regions (eg, Lox) and increased their transcript abundance. SMC-specific Yap deletion in mice compromised this adaptive response in SMCs, leading to an increased AAD incidence. Aortic stress triggers the systemic epigenetic induction of an adaptive response (eg, wound healing, proliferation, matrix organization) in thoracic aortic SMCs that depends on functional biomechanical signal transduction (eg, YAP signaling). Our study highlights the importance of the adaptive response in maintaining aortic homeostasis and preventing AAD in mice.

A human-relevant mixture of endocrine disrupting chemicals induces changes in hippocampal DNA methylation correlating with hyperactive behavior in male mice

Humans are ubiquitously exposed to endocrine disrupting chemicals (EDCs), substances that interfere with endogenous hormonal signaling. Exposure during early development is of particular concern due to the programming role of hormones during this period. A previous epidemiological study has shown association between prenatal co-exposure to 8 EDCs (Mixture N1) and language delay in children, suggesting an effect of this mixture on neurodevelopment. Furthermore, in utero exposure to Mixture N1 altered gene expression and behavior in adult mice. In this study, we investigated whether epigenetic mechanisms could underlie the long term effects of Mixture N1 on gene expression and behavior. To this end, we analyzed DNA methylation at regulatory regions of genes whose expression was affected by Mixture N1 in the hippocampus of in utero exposed mice using bisulfite-pyrosequencing. We show that Mixture N1 decreases DNA methylation in males at three genes that are part of the hypothalamus-pituitary-adrenal (HPA) axis: Nr3c1, Nr3c2, and Crhr1, coding for the glucocorticoid receptor, the mineralocorticoid receptor, and the corticotropin releasing hormone receptor 1, respectively. Furthermore, we show that the decrease in Nr3c1 methylation correlates with increased gene expression, and that Nr3c1, Nr3c2, and Crhr1 methylation correlates with hyperactivity and reduction in social behavior. These findings indicate that an EDC mixture corresponding to a human exposure scenario induces epigenetic changes, and thus programming effects, on the HPA axis that are reflected in the behavioral phenotypes of the adult male offspring.