Expression Of Nearby Genes Research Articles

Integrating Multi-Omics Data to Uncover Prostate Tissue DNA Methylation Biomarkers and Target Genes for Prostate Cancer Risk.

Previous studies have indicated that specific CpG sites may be linked to the risk of prostate cancer (PCa) by regulating the expression of PCa target genes. However, most existing studies aim to identify DNA methylation (DNAm) biomarkers through blood tissue genetic instruments, which impedes the identification of relevant biomarkers in prostate tissue. To identify PCa risk-associated CpG sites in prostate tissue, we established genetic prediction models of DNAm levels using data from normal prostate samples in the GTEx (N = 108) and assessed associations between genetically predicted DNAm in prostate and PCa risk by studying 122,188 cases and 604,640 controls. We observed significant associations for 3879 CpG sites, including 926 at novel genomic loci. Among them, DNAm levels of 80 CpG sites located at novel loci are significantly associated with expression levels of 45 neighboring genes in normal prostate tissue. Of these genes, 11 further exhibit significant associations with PCa risk for their predicted expression levels in prostate tissue. Intriguingly, a total of 31 CpG sites demonstrate consistent association patterns across the methylation-gene expression-PCa risk pathway. Our findings suggest that specific CpG sites may be related to PCa risk by modulating the expression of nearby target genes.

Abstract C092: Transcriptome analysis of human endogenous retroviral elements in tumor microenvironment of prostate cancer patients of different ancestry

Abstract Background: Prostate cancer (PCa) poses unique challenges for Non-Hispanic Black Americans (BA), who often face poorer outcomes than Non-Hispanic White Americans (WA). Even after adjusting for socioeconomic factors, BAs are diagnosed younger and experience more severe clinical effects, potentially due to the tumor microenvironment. Within PCa, Endogenous Retroviruses (ERVs) have been shown to regulate immune response genes and have prognostic value in PCa. Despite constituting 8% of the human, ERVs' role in the tumor microenvironment remains incompletely understood. These findings highlight the need to target microenvironment dysregulation in solid tumor therapies, especially for ERV expression in prostate cancer. Methods: Using RNA-sequencing, we investigated the expression of ERVs in BA (n=18) and WA (n=14) tumor adjacent stroma (TAS) samples. We analyzed reads by mapping them to the human genome (hg38). We then used bedtools and data from the Human Endogenous Retrovirus database (HERVd) to calculate the number of reads at each ERV locus (n=519,060). Differentially expressed ERV elements were then identified in an Audic Claverie Test (fold- change cutoff of 2 and padj < 0.05). The expression of protein-coding genes in the same BA (n=18) and WA (n=14) TAS samples was investigated using RNA-sequencing. Regulation of significant ERV elements were compared to the regulation of nearby immune response genes between the two races. We developed primary stroma (carcinoma associated fibroblast, CAFs) culture of different races with PCa using fresh prostatectomy tissues. To investigate the biological functions of ERVs in CAFs, we developed FANA oligonucleotides (synthetic single- stranded nucleic acid analogs that can modulate gene expression by enzymatic degradation of a target RNA) and corresponding primers to knock down specific ERV sequences near key immune response genes in prostate cancer. Results: We examined 32 prostatectomy specimens of BA (n=18) and WA (n=14) PCa. Of the 5,786 statistically significant differentially expressed ERV elements between BA and WA TAS that were identified, 3,274 elements had an increase in ERV expression for BA TAS, and 2,512 elements had a decrease in ERV expression for BA TAS. 61 MER41 elements were found to be differentially expressed between the two races. The MER family sequences contain binding sites for transcription factors that are crucial in cancer, such as STAT1. This transcription factor can regulate genes involved in cell proliferation, apoptosis, and DNA repair. MER41 is in close proximity (within 5000 base pairs) to IFNA6, an immune response gene. We found both MER41 and IFNA6 to be downregulated in BA TAS, as compared to their WA counterparts. We developed FANA oligonucleotides to knockdown MER41E, MER61-Int, HERVL-Int, MLT1A, and LTR61E1 to investigate their effect on nearby gene expression. Conclusion: The CAF models from PCa patients of different races can be used to explore the biological function of ERV inactivation on anti-tumor immune response genes within the tumor microenvironment. Citation Format: Tara S. K. Jennings, Vinay Kumar, Anton N. Nguyen, Michael M. Ittmann, Patricia Castro, Farah Rahmatpanah. Transcriptome analysis of human endogenous retroviral elements in tumor microenvironment of prostate cancer patients of different ancestry [abstract]. In: Proceedings of the 17th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2024 Sep 21-24; Los Angeles, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2024;33(9 Suppl):Abstract nr C092.

Role of the X Chromosome in Alzheimer Disease Genetics

The X chromosome has remained enigmatic in Alzheimer disease (AD), yet it makes up 5% of the genome and carries a high proportion of genes expressed in the brain, making it particularly appealing as a potential source of unexplored genetic variation in AD. To perform the first large-scale X chromosome-wide association study (XWAS) of AD. This was a meta-analysis of genetic association studies in case-control, family-based, population-based, and longitudinal AD-related cohorts from the US Alzheimer's Disease Genetics Consortium, the Alzheimer's Disease Sequencing Project, the UK Biobank, the Finnish health registry, and the US Million Veterans Program. Risk of AD was evaluated through case-control logistic regression analyses. Data were analyzed between January 2023 and March 2024. Genetic data available from high-density single-nucleotide variant microarrays and whole-genome sequencing and summary statistics for multitissue expression and protein quantitative trait loci available from published studies were included, enabling follow-up genetic colocalization analyses. A total of 1 629 863 eligible participants were selected from referred and volunteer samples, 477 596 of whom were excluded for analysis exclusion criteria. The number of participants who declined to participate in original studies was not available. Risk of AD, reported as odds ratios (ORs) with 95% CIs. Associations were considered at X chromosome-wide (P < 1 × 10-5) and genome-wide (P < 5 × 10-8) significance. Primary analyses are nonstratified, while secondary analyses evaluate sex-stratified effects. Analyses included 1 152 284 participants of non-Hispanic White, European ancestry (664 403 [57.7%] female and 487 881 [42.3%] male), including 138 558 individuals with AD. Six independent genetic loci passed X chromosome-wide significance, with 4 showing support for links between the genetic signal for AD and expression of nearby genes in brain and nonbrain tissues. One of these 4 loci passed conservative genome-wide significance, with its lead variant centered on an intron of SLC9A7 (OR, 1.03; 95% CI, 1.02-1.04) and colocalization analyses prioritizing both the SLC9A7 and nearby CHST7 genes. Of these 6 loci, 4 displayed evidence for escape from X chromosome inactivation with regard to AD risk. This large-scale XWAS of AD identified the novel SLC9A7 locus. SLC9A7 regulates pH homeostasis in Golgi secretory compartments and is anticipated to have downstream effects on amyloid β accumulation. Overall, this study advances our knowledge of AD genetics and may provide novel biological drug targets. The results further provide initial insights into elucidating the role of the X chromosome in sex-based differences in AD.

The pan-tandem repeat map highlights multiallelic variants underlying gene expression and agronomic traits in rice

Tandem repeats (TRs) are genomic regions that tandemly change in repeat number, which are often multiallelic. Their characteristics and contributions to gene expression and quantitative traits in rice are largely unknown. Here, we survey rice TR variations based on 231 genome assemblies and the rice pan-genome graph. We identify 227,391 multiallelic TR loci, including 54,416 TR variations that are absent from the Nipponbare reference genome. Only 1/3 TR variations show strong linkage with nearby bi-allelic variants (SNPs, Indels and PAVs). Using 193 panicle and 202 leaf transcriptomic data, we reveal 485 and 511 TRs act as QTLs independently of other bi-allelic variations to nearby gene expression, respectively. Using plant height and grain width as examples, we identify and validate TRs contributions to rice agronomic trait variations. These findings would enhance our understanding of the functions of multiallelic variants and facilitate rice molecular breeding.

Biased Retention of Environment-Responsive Genes Following Genome Fractionation.

The molecular underpinnings and consequences of cycles of whole-genome duplication (WGD) and subsequent gene loss through subgenome fractionation remain largely elusive. Endogenous drivers, such as transposable elements (TEs), have been postulated to shape genome-wide dominance and biased fractionation, leading to a conserved least-fractionated (LF) subgenome and a degenerated most-fractionated (MF) subgenome. In contrast, the role of exogenous factors, such as those induced by environmental stresses, has been overlooked. In this study, a chromosome-scale assembly of the alpine buckler mustard (Biscutella laevigata; Brassicaceae) that underwent a WGD event about 11 million years ago is coupled with transcriptional responses to heat, cold, drought, and herbivory to assess how gene expression is associated with differential gene retention across the MF and LF subgenomes. Counteracting the impact of TEs in reducing the expression and retention of nearby genes across the MF subgenome, dosage balance is highlighted as a main endogenous promoter of the retention of duplicated gene products under purifying selection. Consistent with the "turn a hobby into a job" model, about one-third of environment-responsive duplicates exhibit novel expression patterns, with one copy typically remaining conditionally expressed, whereas the other copy has evolved constitutive expression, highlighting exogenous factors as a major driver of gene retention. Showing uneven patterns of fractionation, with regions remaining unbiased, but with others showing high bias and significant enrichment in environment-responsive genes, this mesopolyploid genome presents evolutionary signatures consistent with an interplay of endogenous and exogenous factors having driven gene content following WGD-fractionation cycles.

Targeted next-generation sequencing panel to investigate antiplatelet adverse reactions in acute coronary syndrome patients undergoing percutaneous coronary intervention with stenting

Antiplatelet therapy, the gold standard of care for patients with acute coronary syndrome (ACS) undergoing percutaneous coronary intervention (PCI), is one of the therapeutic approaches most associated with the development of adverse drug reactions (ADRs). Although numerous studies have shown that pharmacological intervention based on a limited number of high-evidence variants (primarily CYP2C19*2 and *3) can reduce the incidence of major adverse cardiovascular events (MACEs), ADRs still occur at variable rates (10.1 % in our case) despite personalized therapy.This study aimed to identify novel genetic variants associated with the endpoint of MACEs 12 months after PCI by designing and analyzing a targeted gene panel. We sequenced 244 ACS-PCI-stent patients (109 with event and 135 without event) and 99 controls without structural cardiovascular disease and performed an association analysis to search for unexpected genetic variants.No single nucleotide polymorphisms reached genomic significance after correction, but three novel variants, including ABCA1 (rs2472434), KLB (rs17618244), and ZNF335 (rs3827066), may play a role in MACEs in ACS patients. These genetic variants are involved in regulating high-density lipoprotein levels and cholesterol deposition, and as they are regulatory variants, they may affect the expression of nearby lipid metabolism-related genes. Our findings suggest new targets (both at the gene and pathway levels) that may increase susceptibility to MACEs, but further research is needed to clarify the role and impact of the identified variants before these findings can be incorporated into the therapeutic decision-making process.

Mapping parental DMRs predictive of local and distal methylome remodeling in epigenetic F1 hybrids.

F1 hybrids derived from a cross between two inbred parental lines often display widespread changes in DNA methylation and gene expression patterns relative to their parents. An emerging challenge is to understand how parental epigenomic differences contribute to these events. Here, we generated a large mapping panel of F1 epigenetic hybrids, whose parents are isogenic but variable in their DNA methylation patterns. Using a combination of multi-omic profiling and epigenetic mapping strategies we show that differentially methylated regions in parental pericentromeres act as major reorganizers of hybrid methylomes and transcriptomes, even in the absence of genetic variation. These parental differentially methylated regions are associated with hybrid methylation remodeling events at thousands of target regions throughout the genome, both locally (in cis) and distally (in trans). Many of these distally-induced methylation changes lead to nonadditive expression of nearby genes and associate with phenotypic heterosis. Our study highlights the pleiotropic potential of parental pericentromeres in the functional remodeling of hybrid genomes and phenotypes.

Identification and characterization of structural variants related to meat quality in pigs using chromosome-level genome assemblies

BackgroundMany studies have been performed to identify various genomic loci and genes associated with the meat quality in pigs. However, the full genetic architecture of the trait still remains unclear in part because of the lack of accurate identification of related structural variations (SVs) which resulted from the shortage of target breeds, the limitations of sequencing data, and the incompleteness of genome assemblies. The recent generation of a new pig breed with superior meat quality, called Nanchukmacdon, and its chromosome-level genome assembly (the NCMD assembly) has provided new opportunities.ResultsBy applying assembly-based SV calling approaches to various genome assemblies of pigs including Nanchukmacdon, the impact of SVs on meat quality was investigated. Especially, by checking the commonality of SVs with other pig breeds, a total of 13,819 Nanchukmacdon-specific SVs (NSVs) were identified, which have a potential effect on the unique meat quality of Nanchukmacdon. The regulatory potentials of NSVs for the expression of nearby genes were further examined using transcriptome- and epigenome-based analyses in different tissues.ConclusionsWhole-genome comparisons based on chromosome-level genome assemblies have led to the discovery of SVs affecting meat quality in pigs, and their regulatory potentials were analyzed. The identified NSVs will provide new insights regarding genetic architectures underlying the meat quality in pigs. Finally, this study confirms the utility of chromosome-level genome assemblies and multi-omics analysis to enhance the understanding of unique phenotypes.

A screen for regeneration-associated silencer regulatory elements in zebrafish

Regeneration involves gene expression changes explained in part by context-dependent recruitment of transcriptional activators to distal enhancers. Silencers that engage repressive transcriptional complexes are less studied than enhancers and more technically challenging to validate, but they potentially have profound biological importance for regeneration. Here, we identified candidate silencers through a screening process that examined the ability of DNA sequences to limit injury-induced gene expression in larval zebrafish after fin amputation. A short sequence (s1) on chromosome 5 near several genes that reduce expression during adult fin regeneration could suppress promoter activity in stable transgenic lines and diminish nearby gene expression in knockin lines. High-resolution analysis of chromatin organization identified physical associations of s1 with gene promoters occurring preferentially during fin regeneration, and genomic deletion of s1 elevated the expression of these genes after fin amputation. Our study provides methods to identify "tissue regeneration silencer elements" (TRSEs) with the potential to reduce unnecessary or deleterious gene expression during regeneration.

CRISPR/Cas9-mediated deletion of a GA-repeat in human GPM6B leads to disruption of neural cell differentiation from NT2 cells

The human neuron-specific gene, GPM6B (Glycoprotein membrane 6B), is considered a key gene in neural cell functionality. This gene contains an exceptionally long and strictly monomorphic short tandem repeat (STR) of 9-repeats, (GA)9. STRs in regulatory regions, may impact on the expression of nearby genes. We used CRISPR-based tool to delete this GA-repeat in NT2 cells, and analyzed the consequence of this deletion on GPM6B expression. Subsequently, the edited cells were induced to differentiate into neural cells, using retinoic acid (RA) treatment. Deletion of the GA-repeat significantly decreased the expression of GPM6B at the RNA (p < 0.05) and protein (40%) levels. Compared to the control cells, the edited cells showed dramatic decrease of the astrocyte and neural cell markers, including GFAP (0.77-fold), TUBB3 (0.57-fold), and MAP2 (0.2-fold). Subsequent sorting of the edited cells showed an increased number of NES (p < 0.01), but a decreased number of GFAP (p < 0.001), TUBB3 (p < 0.05), and MAP2 (p < 0.01), compared to the control cells. In conclusion, CRISPR/Cas9-mediated deletion of a GA-repeat in human GPM6B, led to decreased expression of this gene, which in turn, disrupted differentiation of NT2 cells into neural cells.

Transposable elements mediate genetic effects altering the expression of nearby genes in colorectal cancer

Transposable elements (TEs) are prevalent repeats in the human genome, play a significant role in the regulome, and their disruption can contribute to tumorigenesis. However, TE influence on gene expression in cancer remains unclear. Here, we analyze 275 normal colon and 276 colorectal cancer samples from the SYSCOL cohort, discovering 10,231 and 5,199 TE-expression quantitative trait loci (eQTLs) in normal and tumor tissues, respectively, of which 376 are colorectal cancer specific eQTLs, likely due to methylation changes. Tumor-specific TE-eQTLs show greater enrichment of transcription factors, compared to shared TE-eQTLs suggesting specific regulation of their expression in tumor. Bayesian networks reveal 1,766 TEs as mediators of genetic effects, altering the expression of 1,558 genes, including 55 known cancer driver genes and show that tumor-specific TE-eQTLs trigger the driver capability of TEs. These insights expand our knowledge of cancer drivers, deepening our understanding of tumorigenesis and presenting potential avenues for therapeutic interventions.

OP08 Multi-ancestry genome-wide association study of inflammatory bowel disease identifies 125 novel loci and directly implicates new genes in disease susceptibility

Abstract Background Genome-wide association studies (GWASs) have identified 320 loci associated with inflammatory bowel disease (IBD). About 15% of these loci include genes that lie in pathways targeted by drugs currently approved to treat IBD. Thus, identifying new loci can not only provide further insights into causal biology but also reveal novel therapeutic targets. Methods We performed a multi-ancestry GWAS meta-analysis of 67 cohorts, from European, East Asian and South Asian populations. The final dataset included 125,992 IBD patients (59,734 with CD; 57,565 with UC) and more than 1.2 million population controls (Fig 1). Genotype imputation was undertaken using the TOPMed diverse population panel and association tests were performed using REGENIE. The results were meta-analysed using a fixed effects meta-analysis via Metal. To prioritise effector genes and generate therapeutic hypotheses, we integrated our results with bulk and single cell transcriptome data from relevant tissues and cell types. Through colocalization analyses we identified those signals driving the expression of nearby genes (cis-eQTL), suggesting that aberrant expression of these genes underpins the association. Results We identified 290 new independent disease susceptibility association signals. Out of these, 125 new signals are located &amp;gt;1Mb from any previously reported loci. Six new loci contain genes implicated in monogenic syndromes that include colitis GUCY2C, G6PC3, NFKBIA, PIK3CD, PIK3R1, ZBTB24. Up to 40% of all disease susceptibility association signals colocalize with a cis-eQTL. For instance, a new variant associated with increased risk of UC (P=8x10-14) decreases the expression of the transmembrane mitochondrial protein TSPO in T cells, B cells, macrophages, and transverse colon. Pro-inflammatory macrophages show consistent TSPO downregulation. Loss of TSPO in mice causes exacerbated DSS-induced colitis. Ongoing clinical trials are evaluating whether targeting this gene could ameliorate the immune response in multiple sclerosis patients. We also identified a new CD susceptibility variant (P=3x10-16) associated with increased expression of WNK1 in monocytes and macrophages. Whereas loss of WNK1 in mice show increased expression of IL-1B in response to NLRP3 inflammasome activation; it also stimulates autophagy through the upregulation of the PI3KC3 complex. Impaired autophagy has shown to be key in promoting CD susceptibility. Conclusion Integrating molecular data from multiple cell types and conditions not only expands our ability to determine novel candidate causal genes, it also provides additional evidence of the utility of genetics as an instrument to guide the development of new therapeutic drugs or the repurposing of existing molecules.

Transposable elements acquire time- and sex-specific transcriptional and epigenetic signatures along mouse fetal gonad development.

Gonadal sex determination in mice is a complex and dynamic process, which is crucial for the development of functional reproductive organs. The expression of genes involved in this process is regulated by a variety of genetic and epigenetic mechanisms. Recently, there has been increasing evidence that transposable elements (TEs), which are a class of mobile genetic elements, play a significant role in regulating gene expression during embryogenesis and organ development. In this study, we aimed to investigate the involvement of TEs in the regulation of gene expression during mouse embryonic gonadal development. Through bioinformatics analysis, we aimed to identify and characterize specific TEs that operate as regulatory elements for sex-specific genes, as well as their potential mechanisms of regulation. We identified TE loci expressed in a time- and sex-specific manner along fetal gonad development that correlate positively and negatively with nearby gene expression, suggesting that their expression is integrated to the gonadal regulatory network. Moreover, chromatin accessibility and histone post-transcriptional modification analyses in differentiating supporting cells revealed that TEs are acquiring a sex-specific signature for promoter-, enhancer-, and silencer-like elements, with some of them being proximal to critical sex-determining genes. Altogether, our study introduces TEs as the new potential players in the gene regulatory network that controls gonadal development in mammals.

Genetic variation is a key determinant of chromatin accessibility and drives differences in the regulatory landscape of C57BL/6J and 129S1/SvImJ mice.

Most common genetic variants associated with disease are located in non-coding regions of the genome. One mechanism by which they function is through altering transcription factor (TF) binding. In this study, we explore how genetic variation is connected to differences in the regulatory landscape of livers from C57BL/6J and 129S1/SvImJ mice fed either chow or a high-fat diet. To identify sites where regulatory variation affects TF binding and nearby gene expression, we employed an integrative analysis of H3K27ac ChIP-seq (active enhancers), ATAC-seq (chromatin accessibility) and RNA-seq (gene expression). We show that, across all these assays, the genetically driven (i.e. strain-specific) differences in the regulatory landscape are more pronounced than those modified by diet. Most notably, our analysis revealed that differentially accessible regions (DARs, N=29635, FDR<0.01 and fold change>50%) are almost always strain-specific and enriched with genetic variation. Moreover, proximal DARs are highly correlated with differentially expressed genes. We also show that TF binding is affected by genetic variation, which we validate experimentally using ChIP-seq for TCF7L2 and CTCF. This study provides detailed insights into how non-coding genetic variation alters the gene regulatory landscape, and demonstrates how this can be used to study the regulatory variation influencing TF binding.

Genome‐wide association study for Alzheimer disease in APOE ɛ4/ɛ4 carriers

AbstractBackgroundThe APOE ɛ4/ɛ4 genotype has the greatest risk for late‐onset Alzheimer disease (AD). However, modifying effects of other genes on AD risk in APOE ɛ4 homozygotes are not fully understood.MethodsWe conducted a genome‐wide association study (GWAS) for AD among 1,223 non‐Hispanic White APOE ɛ4/ɛ4 carriers assembled by the Alzheimer’s Disease Genetics Consortium (ADGC). Expression of genes within 50 kb of suggestively associated single nucleotide polymorphisms (SNPs, P&lt;10−6) were compared between AD and control brains from 576 Religious Orders Study and Memory and Aging Project (ROSMAP) subjects. We conducted survival analysis between carriers and non‐carriers of AD risk alleles for the suggestive SNPs in 4,287 Framingham Heart Study (FHS) participants. We compared expression of the top‐ranked genes between ɛ2/ɛ2 (APOE2) and ɛ4/ɛ4 (APOE4) human induced pluripotent stem cells (iPSC)‐derived neurons (iNeurons) and astrocytes (iAstrocytes) in a co‐culture system (Jun et al, 2022) and separately evaluated associations of their expression levels with neuropathological traits from 208 brains from FHS participants.ResultsWe identified association of AD with 11 SNPs at the suggestive level of which four near TRIB2, CTNNA2, ZNF443, and KLHL29 were genome‐wide significant (P&lt;5×10−8). Rs75983775 near TRIB2, the most significant finding (P = 9.4×10−11), was associated with decreased risk of AD in ADGC APOE ɛ4 homozygotes (odds ratio = 0.12) and increased TRIB2 expression in ROSMAP ɛ4‐carrier brains (P = 0.02). TRIB2 expression was increased in APOE2 compared to APOE4 iAstrocytes (log2 fold change (log2FC) = 0.41, adjusted P (adjP) = 0.03) and iNeurons (log2FC = 0.31, adjP = 0.002) and elevated the Iba1 level in FHS brains (P = 0.03). Association with TENM3 SNP rs56035906 was suggestively associated with reduced AD risk in ADGC APOE ɛ4 homozygotes (odds ratio = 0.30, P = 4.3×10−7) and with lower incidence of AD among FHS ɛ4 carriers (P = 0.03). TENM3 expression was decreased in APOE2 compared to APOE4 iNeurons (log2FC = 0.59, adjP = 2.4×10−4) and associated with increasing phosphorylated tau 181 (P = 0.007) and C4a protein expression (P = 9.6×10−4) in FHS brains.ConclusionWe identified variants that may protect against AD risk among APOE ɛ4/ɛ4 carriers, delay disease onset, modulate expression of nearby genes in iAstrocytes and iNeurons, and are associated with measures of tau and neuroinflammation in autopsied brains.

Functional TET2 gene polymorphisms increase the risk of neuroblastoma in Chinese children.

The 5-methylcytosine (m5C) is the key chemical modification in RNAs. As one of the demethylases in m5C, TET2 has been shown as a tumor suppressor. However, the impact of TET2 gene polymorphisms on neuroblastoma has not been elucidated. 402 neuroblastoma patients and 473 controls were genotyped for TET2 gene polymorphisms using the TaqMan method. The impact of TET2 gene polymorphisms on neuroblastoma susceptibility was determined using multivariate logistic regression analysis. We also adopted genotype-tissue expression database to explore the impact of TET2 gene polymorphisms on the expression of host and nearby genes. We used the R2 platform and Sangerbox tool to analyze the relationship between gene expression and neuroblastoma risk and prognosis through non-parametric testing and Kaplan-Meier analysis, respectively. We found the TET2 gene polymorphisms (rs10007915 G > C and rs7670522 A > C) and the combined 2-5 risk genotypes can significantly increase neuroblastoma risk. Stratification analysis showed that these significant associations were more prominent in certain subgroups. TET2 rs10007915 G > C and rs7670522 A > C are significantly associated with reduced expression of TET2 mRNA. Moreover, lower expression of TET2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastoma. The rs10007915 G > C and rs7670522 A > C are significantly related to the increased expression of inorganic pyrophosphatase 2 mRNA, and higher expression of PPA2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastomas. In summary, TET2 rs10007915 G > C and rs7670522 A > C significantly confer neuroblastoma susceptibility, and further research is needed to investigate the underlying mechanisms.

Identification of a Novel Gene Regulatory Element in Human Erythroid Progenitor Cells

Early erythropoiesis involves commitment of multi-lineage progenitors to erythroid progenitor cells and differentiation into erythroid burst-forming unit (BFU-E) cells, then erythroid colony-forming unit (CFU-E) cells. In contrast to terminal erythroid differentiation, regulation of gene expression in early erythropoiesis is poorly defined. To address this knowledge gap, we identified active enhancers in erythroid progenitor cells linked to known disease genes and genes associated with erythrocyte traits in genome wide association studies. Human CD34+ hematopoietic stem and progenitor cells (HSPCs) isolated from mobilized peripheral blood were cultured towards the erythroid lineage. From these cultures, FACS-based methods were utilized to purify four stages of erythroid progenitor cells, termed EP1 to EP4, and five populations of terminally differentiated erythroblasts, from proerythroblasts to orthochromatic erythroblasts. ChIP-seq was performed using antibodies to H3K4me1, H3K27Ac, H3K27me3, and H3K9me2/3. Chromatin accessibility was determined by ATAC-seq. Patterns of gene expression were determined by RNA-seq. Active enhancers were identified using the ModHMM algorithm. Numerous erythroid progenitor cell-specific active enhancers were identified near genes linked to known disease genes and GWAS-associated erythrocyte traits. rs162375, a single nucleotide polymorphism (SNP) linked to mean corpuscular hemoglobin, was in a 475bp Long Terminal Repeat Retrotransposon (LTR) element upstream of the RRP1B gene. The LTR enhancer near the RRP1B gene is a member of the LTR15 family. Transposable elements are a source of novelty in primate gene regulation when they acquire binding sites for critical tissue-specific transcription factors, and then these tissue-specific regulatory regions direct tissue-specific expression of nearby genes. The SNP-containing LTR15 erythroid progenitor-specific enhancer upstream of RRP1B contains 2 conserved GATA motifs and a conserved RUNX motif. The enhancer region had a peak of open chromatin determined by ATAC-seq, and H3K27 acetylation and strong GATA1 binding determined by ChIP-seq. These findings were present only in erythroid progenitor cells, not in earlier hematopoietic cells, terminally differentiating erythroid cells, or other non-erythroid cells. CRISPR-based gene editing was utilized to delete the LTR15 enhancer in CD34+ HSPCs cultured towards the erythroid lineage into erythroid progenitor cells. mRNA expression of RRP1B was decreased to 70.2% of wild type. Luciferase-based reporter gene studies demonstrated the enhancer increased expression of a minimal promoter 3.2-fold over background. Mutation of the GATA and RUNX sites individually or together decreased reporter gene activity to background levels. Bioinformatic analyses identified 291 additional erythroid progenitor cell specific LTR15 element active enhancers, with most containing two GATA motifs and a RUNX motif. The majority of the novel enhancer regions were also marked by progenitor cell-specific ATAC peaks and strong GATA1 binding in ChIP-seq. GATA1 and RUNX, both present in early erythropoiesis, are required to activate the regulatory elements in erythroid progenitor cells. Decreasing levels of RUNX in terminal erythroid differentiation parallels loss of enhancer activity. Since LTR15 first appeared in Catarrhine, old-world monkeys, these data indicate this novel regulatory element has recently evolved to contribute to regulation of gene expression in early erythropoiesis in primates.

The Molecular Impacts of Retrotransposons in Development and Diseases.

Retrotransposons are invasive genetic elements that constitute substantial portions of mammalian genomes. They have the potential to influence nearby gene expression through their cis-regulatory sequences, reverse transcription machinery, and the ability to mold higher-order chromatin structures. Due to their multifaceted functions, it is crucial for host fitness to maintain strict regulation of these parasitic sequences to ensure proper growth and development. This review explores how subsets of retrotransposons have undergone evolutionary exaptation to enhance the complexity of mammalian genomes. It also highlights the significance of regulating these elements, drawing on recent studies conducted in human and murine systems.

Epigenetic Classification of Acute Myeloid Leukemia Revealed By Genome-Wide Chromatin Profiling

Background Acute myeloid leukemia (AML) is a clinically and molecularly heterogeneous disease and is predominantly defined by genetic abnormalities and morphology. However, many lines of evidence suggest the role of epigenetic alterations in the pathogenesis and the classification of AML, whereas no large-scale comprehensive analyses of both genetic and epigenetic alterations in AML have been reported. Methods We enrolled a total of 1,564 primary AML cases (524 from Japanese and 1,040 from Swedish cohorts), whose diagnostic samples were subjected to ATAC sequencing together with targeted capture sequencing using an in-house gene panel including 331 known AML driver genes and an additional 1,158-1,317 probes for copy number detection. We also performed RNA (N=1,398) and H3K27ac-targeted ChIP sequencing (N=80). Normal bone marrow samples from 25 donors as well as 15 remission samples were analyzed by ATAC-seq as a control. Results ATAC-seq revealed approximately 300,000 reproducible peaks in primary AML. Most peaks were found in gene-distal regions, and explained a large variance across samples. Among all peaks, 72% were not detected in control samples and thus, specific to AML. Analysis of the correlations between ATAC-seq accessibility and gene expression revealed a subset of regulatory elements that regulate nearby gene expression; each peak was linked to a median of 1 gene (0-31), while each gene was associated with a median of 11 peaks (0-110). We next estimated the cellular composition of the bulk AML samples by deconvolution analysis using ATAC-seq data. The predominance of hematopoietic stem cell (HSC), monocyte, and erythroid signatures was associated with French-American-British (FAB) subtypes of minimal differentiation/without maturation subtypes (M0/M1), myelomonocytic/monoblastic/monocytic subtypes (M4/M5), and erythroid subtype (M6), respectively, suggesting bulk ATAC-seq of AML samples reflected the differentiation potential of leukemic cells. The ATAC-based clustering identified 11 distinct epigenetic subgroups, including three well-known genetic classes defined by gene fusions, such as PML::RARA, RUNX1::RUNX1T1, and CBFB::MYH11. By contrast, some genetic classes defined by a single mutation, such as NPM1 or CEBPA mutation, were further classified into multiple distinct epigenetic subgroups. For example, the NPM1-mutated subtype was separated into four subgroups based on the varying degrees of ATAC signatures associated with HSC and monocytes. Each subgroup of NPM1-mutated AML showed different clinical characteristics such as white blood cell counts, age, and prognosis, and co-mutation pattern in FLT3, DNMT3A, and IDH1/ IDH2 mutations. Gene ontology analysis revealed that the subgroup of NPM1-mutated AML with high HSC signature showed upregulated TGF-β signaling, while one with monocytic signature showed upregulated TNF-α signaling and interferon-γ response. Next, we evaluated transcription factor (TF) binding sites in a genome-wide manner, using TF motifs and foot printing in upregulated ATAC-peaks in each subgroup and revealed distinct profiles of activated TFs in each epigenetic subgroup of NPM1-mutated AML. In the subgroup characterized by high HSC signature, HSC-related TFs that belong to MECOM, HOX, and RUNX families, were more active. By contrast, other TFs, such as IRF, CEBP, and JUN/FOS families, were more enriched in the subgroup showing high monocyte signature. By developing a prediction model for these subgroups of NPM1-mutated AML based on the gene expression profile, we reproduced these subgroups and validated clinical and molecular features in an external Beat AML cohort. Utilizing ex vivo drug sensitivity dataset in the Beat AML, we identified candidate drugs showing subgroup-specific effectiveness such as panobinostat for NPM1-mutated AML with monocytic differentiation. Finally, ChIP-seq analysis identified a number of super-enhancers in the vicinity of leukemia-related TFs and oncogenes. Conclusion Through comprehensive profiling of chromatin accessibility and gene mutations/expression in a large cohort of AML patients, we demonstrate the epigenetic heterogeneity and identify novel AML subgroups that had unique epigenetic and clinical features. Our findings highlight the role of epigenetic alterations in AML pathogenesis and impact on molecular classification of AML.

Co-localization of clusters of TCR-regulated genes with TAD rearrangements

BackgroundGene expression has long been known to be influenced by the relative proximity of DNA regulatory elements. Topologically associating domains (TADs) are self-interacting genomic regions involved in regulating gene expression by controlling the proximity of these elements. Prior studies of TADs and their biological roles have revealed correlations between TAD changes and cellular differentiation. Here, we used Hi-C and RNA-seq data to correlate TCR-induced changes in TAD structure and gene expression in human CD4+ T cells.ResultsWe developed a pipeline, Differentially Expressed Gene Enrichment Finder (DEGEF), that identifies regions of differentially expressed gene enrichment. Using DEGEF, we found that TCR-regulated genes cluster non-uniformly across the genome and that these clusters preferentially localized in regions of TAD rearrangement. Interestingly, clusters of upregulated genes preferentially formed new Hi-C contacts compared to downregulated clusters, suggesting that TCR-activated CD4+ T cells may regulate genes by changing stimulatory contacts rather than inhibitory contacts.ConclusionsOur observations support a significant relationship between TAD rearrangements and changes in local gene expression. These findings indicate potentially important roles for TAD rearrangements in shaping their local regulatory environments and thus driving differential expression of nearby genes during CD4+ T cell activation. Moreover, they provide new insights into global mechanisms that regulate gene expression.