Genome-wide screening and functional validation of methylation barriers near promoters.

CpG islands near promoters are normally unmethylated despite being surrounded by densely methylated regions. Aberrant hypermethylation of these CpG islands has been associated with the development of various human diseases. Although local genetic elements have been speculated to play a role in protecting promoters from methylation, only a limited number of methylation barriers have been identified. In this study, we conducted an integrated computational and experimental investigation of colorectal cancer methylomes. Our study revealed 610 genes with disrupted methylation barriers. Genomic sequences of these barriers shared a common 41-bp sequence motif (MB-41) that displayed homology to the chicken HS4 methylation barrier. Using the CDKN2A (P16) tumor suppressor gene promoter, we validated the protective function of MB-41 and showed that loss of such protection led to aberrant hypermethylation. Our findings highlight a novel sequence signature of cis-acting methylation barriers in the human genome that safeguard promoters from silencing.

A Mouse-Specific Model to Detect Genes under Selection in Tumors.

The mouse is a widely used model organism in cancer research. However, no computational methods exist to identify cancer driver genes in mice due to a lack of labeled training data. To address this knowledge gap, we adapted the GUST (Genes Under Selection in Tumors) model, originally trained on human exomes, to mouse exomes via transfer learning. The resulting tool, called GUST-mouse, can estimate long-term and short-term evolutionary selection in mouse tumors, and distinguish between oncogenes, tumor suppressor genes, and passenger genes using high-throughput sequencing data. We applied GUST-mouse to analyze 65 exomes of mouse primary breast cancer models and 17 exomes of mouse leukemia models. Comparing the predictions between cancer types and between human and mouse tumors revealed common and unique driver genes. The GUST-mouse method is available as an open-source R package on github.

Evolution of COVID-19 Health Disparities in Arizona.

COVID-19 burdens are disproportionally high in underserved and vulnerable communities in Arizona. As the pandemic progressed, it is unclear if the initial associated health disparities have changed. This study aims to elicit the dynamic landscape of COVID-19 disparities at the community level and identify newly emerging vulnerable subpopulations. Findings from this study can inform interventions to increase health equity among minoritized communities in the Southwest, other regions of the US, and globally. We compiled biweekly COVID-19 case counts of 274 zip code tabulation areas (ZCTAs) in Arizona from October 21, 2020, to November 25, 2021, a time spanning multiple waves of COVID-19 case growth. Within each biweekly period, we tested the associations between the growth rate of COVID-19 cases and the population composition in a ZCTA including race/ethnicity, income, employment, and age using multiple regression analysis. We then compared the associations across time periods to discover temporal patterns of health disparities. The association between the percentage of Latinx population and the COVID-19 growth rate was positive before April 2021 but gradually converted to negative afterwards. The percentage of Black population was not associated with the COVID-19 growth rate at the beginning of the study but became positive after January 2021 which persisted till the end of the study period. Young median age and high unemployment rate emerged as new risk factors around mid-August 2021. Based on these findings, we identified 37 ZCTAs that were highly vulnerable to future fast escalation of COVID-19 cases. As the pandemic progresses, vulnerabilities associated with Latinx ethnicity improved gradually, possibly bolstered by culturally responsive programs in Arizona to support Latinx. Still communities with disadvantaged social determinants of health continued to struggle. Our findings inform the need to adjust current resource allocations to support the design and implementation of new interventions addressing the emerging vulnerabilities at the community level.

Sleeping Beauty Screen Identifies RREB1 and Other Genetic Drivers in Human B-cell Lymphoma.

Follicular lymphoma and diffuse large B-cell lymphoma (DLBCL) are the most common non-Hodgkin lymphomas distinguishable by unique mutations, chromosomal rearrangements, and gene expression patterns. Here, it is demonstrated that early B-cell progenitors express 2',3'-cyclic-nucleotide 3' phosphodiesterase (CNP) and that when targeted with Sleeping Beauty (SB) mutagenesis, Trp53R270H mutation or Pten loss gave rise to highly penetrant lymphoid diseases, predominantly follicular lymphoma and DLBCL. In efforts to identify the genetic drivers and signaling pathways that are functionally important in lymphomagenesis, SB transposon insertions were analyzed from splenomegaly specimens of SB-mutagenized mice (n = 23) and SB-mutagenized mice on a Trp53R270H background (n = 7) and identified 48 and 12 sites with statistically recurrent transposon insertion events, respectively. Comparison with human data sets revealed novel and known driver genes for B-cell development, disease, and signaling pathways: PI3K-AKT-mTOR, MAPK, NFκB, and B-cell receptor (BCR). Finally, functional data indicate that modulating Ras-responsive element-binding protein 1 (RREB1) expression in human DLBCL cell lines in vitro alters KRAS expression, signaling, and proliferation; thus, suggesting that this proto-oncogene is a common mechanism of RAS/MAPK hyperactivation in human DLBCL. IMPLICATIONS: A forward genetic screen identified new genetic drivers of human B-cell lymphoma and uncovered a RAS/MAPK-activating mechanism not previously appreciated in human lymphoid disease. Overall, these data support targeting the RAS/MAPK pathway as a viable therapeutic target in a subset of human patients with DLBCL.

Imprinting defects at human 14q32 locus alters gene expression and is associated with the pathobiology of osteosarcoma.

Osteosarcoma is the most common primary bone malignancy affecting children and adolescents. Although several genetic predisposing conditions have been associated with osteosarcoma, our understanding of its pathobiology is rather limited. Here we show that, first, an imprinting defect at human 14q32-locus is highly prevalent (87%) and specifically associated with osteosarcoma patients < 30 years of age. Second, the average demethylation at differentially methylated regions (DMRs) in the 14q32-locus varied significantly compared to genome-wide demethylation. Third, the 14q32-locus was enriched in both H3K4-me3 and H3K27-me3 histone modifications that affected expression of all imprinted genes and miRNAs in this region. Fourth, imprinting defects at 14q32 - DMRs are present in triad DNA samples from affected children and their biological parents. Finally, imprinting defects at 14q32-DMRs were also observed at higher frequencies in an Rb1/Trp53 mutation-induced osteosarcoma mouse model. Further analysis of normal and tumor tissues from a Sleeping Beauty mouse model of spontaneous osteosarcoma supported the notion that these imprinting defects may be a key factor in osteosarcoma pathobiology. In conclusion, we demonstrate that imprinting defects at the 14q32 locus significantly alter gene expression, may contribute to the pathogenesis of osteosarcoma, and could be predictive of survival outcomes.

Epigenetic synergy between decitabine and platinum derivatives.

BACKGROUND: Aberrant epigenetic silencing of tumor suppressor genes has been recognized as a driving force in cancer. Epigenetic drugs such as the DNA methylation inhibitor decitabine reactivate genes and are effective in myeloid leukemia, but resistance often develops and efficacy in solid tumors is limited. To improve their clinical efficacy, we searched among approved anti-cancer drugs for an epigenetic synergistic combination with decitabine. RESULTS: We used the YB5 cell line, a clonal derivative of the SW48 colon cancer cell line that contains a single copy of a hypermethylated cytomegalovirus (CMV) promoter driving green fluorescent protein (GFP) to screen for drug-induced gene reactivation and synergy with decitabine. None of the 16 anti-cancer drugs tested had effects on their own. However, in combination with decitabine, platinum compounds showed striking synergy in activating GFP. This was dose dependent, observed both in concurrent and sequential combinations, and also seen with other alkylating agents. Clinically achievable concentrations of carboplatin at (25 µM) and decitabine reactivated GFP in 28 % of the YB5 cells as compared to 15 % with decitabine alone. Epigenetic synergy was also seen at endogenously hypermethylated tumor suppressor genes such as MLH1 and PDLIM4. Genome-wide studies showed that reactivation of hypermethylated genes by the combination was significantly better than that induced by decitabine alone or carboplatin alone. Platinum compounds did not enhance decitabine-induced hypomethylation. Rather, we found significantly inhibited HP1α expression by carboplatin and the combination. This was accompanied by increased histone H3 lysine 4 (H3K4) trimethylation and histone H3 lysine 9 (H3K9) acetylation at reactivated genes (P < 0.0001) and reduced occupancy by methyl-binding proteins including MeCP2 and methyl-CpG-binding domain protein 2 (MBD2) (P < 0.0001). CONCLUSIONS: Our results suggest that the combination of decitabine with platinum analogs shows epigenetic synergy that might be exploited in the treatment of different cancers.

HomeRun Vector Assembly System: a flexible and standardized cloning system for assembly of multi-modular DNA constructs.

Advances in molecular and synthetic biology call for efficient assembly of multi-modular DNA constructs. We hereby present a novel modular cloning method that obviates the need for restriction endonucleases and significantly improves the efficiency in the design and construction of complex DNA molecules by standardizing all DNA elements and cloning reactions. Our system, named HomeRun Vector Assembly System (HVAS), employs a three-tiered vector series that utilizes both multisite gateway cloning and homing endonucleases, with the former building individual functional modules and the latter linking modules into the final construct. As a proof-of-principle, we first built a two-module construct that supported doxycycline-induced expression of green fluorescent protein (GFP). Further, with a three-module construct we showed quantitatively that there was minimal promoter leakage between neighbouring modules. Finally, we developed a method, in vitro Cre recombinase-mediated cassette exchange (RMCE) cloning, to regenerate a gateway destination vector from a previous multisite gateway cloning reaction, allowing access to existing DNA element libraries in conventional gateway entry clones, and simple creation of constructs ready for in vivo RMCE. We believe these methods constitute a useful addition to the standard molecular cloning techniques that could potentially support industrial scale synthesis of DNA constructs.

Dose-dependent differential mRNA target selection and regulation by let-7a-7f and miR-17-92 cluster microRNAs.

MicroRNAs (miRNAs) are important players of post-transcriptional gene regulation. Individual miRNAs can target multiple mRNAs and a single mRNA can be targeted by many miRNAs. We hypothesized that miRNAs select and regulate their targets based on their own expression levels, those of their target mRNAs and triggered feedback loops. We studied the effects of varying concentrations of let-7a-7f and the miR-17-92 cluster plasmids on the reporter genes carrying either DICER- or cMYC -3'UTR in Huh-7 cells. We showed that let-7 significantly downregulated expression of DICER 3'UTR reporter at lower concentrations, but selectively downregulated expression of a cMYC 3'UTR reporter at higher dose. This miRNA dose-dependent target selection was also confirmed in other target genes, including CCND1, CDKN1 and E2F1. After overexpressing let-7a-7f or the miR-17-92 clusters at wide-ranging doses, the target genes displayed a nonlinear correlation to the transfected miRNA. Further, by comparing the expression levels of let-7a and miR-17-5p, along with their selected target genes in 3 different cell lines, we found that the knockdown dose of each miRNA was directly related to their baseline expression level, that of the target gene and feedback loops. These findings were supported by gene modulation studies using endogenous levels of miR-29, -1 and -206 and a luciferase reporter system in multiple cell lines. Finally, we determined that the miR-17-92 cluster affected cell viability in a dose-dependent manner. In conclusion, we have shown that miRNAs potentially select their targets in a dose-dependent and nonlinear fashion that affects biological function; and this represents a novel mechanism by which miRNAs orchestrate the finely tuned balance of cell function.

Repetitive elements and enforced transcriptional repression co-operate to enhance DNA methylation spreading into a promoter CpG-island.

Repression of many tumor suppressor genes in cancer is concurrent with aberrantly increased DNA methylation levels at promoter CpG islands (CGIs). About one-fourth of empirically defined human promoters are surrounded by or contain clustered repetitive elements. It was previously observed that a sharp transition of methylation exists between highly methylated repetitive elements and unmethylated promoter-CGIs in normal tissues. The factors that lead to aberrant CGI hypermethylation in cancer remain poorly understood. Here, we established a site-specific integration system with enforced local transcriptional repression in colorectal cancer cells and monitored the occurrence of initial de novo methylation at specific CG sites adjacent to the CGI of the INSL6 promoter, which could be accelerated by binding a KRAB-containing transcriptional factor. Additional repetitive elements from P16 and RIL (PDLIM4), if situated adjacent to the promoter of INSL6, could confer DNA methylation spreading into the CGI particularly in the setting of KRAB-factor binding. However, a repressive chromatin alone was not sufficient to initiate DNA methylation, which required specific DNA sequences and was integration-site (and/or cell-line) specific. Overall, these results demonstrate a requirement for specific DNA sequences to trigger de novo DNA methylation, and repetitive elements as cis-regulatory factors to cooperate with advanced transcriptional repression in promoting methylation spreading.

Mechanisms of resistance to decitabine in the myelodysplastic syndrome.

PURPOSE: The DNA methylation inhibitor 5-aza-2'-deoxycytidine (DAC) is approved for the treatment of myelodysplastic syndromes (MDS), but resistance to DAC develops during treatment and mechanisms of resistance remain unknown. Therefore, we investigated mechanisms of primary and secondary resistance to DAC in MDS. PATIENTS AND METHODS: We performed Quantitative Real-Time PCR to examine expression of genes related to DAC metabolism prior to therapy in 32 responders and non-responders with MDS as well as 14 patients who achieved a complete remission and subsequently relapsed while on therapy (secondary resistance). We then performed quantitative methylation analyses by bisulfite pyrosequencing of 10 genes as well as Methylated CpG Island Amplification Microarray (MCAM) analysis of global methylation in secondary resistance. RESULTS: Most genes showed no differences by response, but the CDA/DCK ratio was 3 fold higher in non-responders than responders (P<.05), suggesting that this could be a mechanism of primary resistance. There were no significant differences at relapse in DAC metabolism genes, and no DCK mutations were detected. Global methylation measured by the LINE1 assay was lower at relapse than at diagnosis (P<.05). On average, the methylation of 10 genes was lower at relapse (16.1%) compared to diagnosis (18.1%) (P<.05). MCAM analysis showed decreased methylation of an average of 4.5% (range 0.6%-9.7%) of the genes at relapse. By contrast, new cytogenetic changes were found in 20% of patients. CONCLUSION: Pharmacological mechanisms are involved in primary resistance to DAC, whereas hypomethylation does not prevent a relapse for patients with DAC treatment.

Genomewide microRNA down-regulation as a negative feedback mechanism in the early phases of liver regeneration.

UNLABELLED: The liver is one of the few organs that have the capacity to regenerate in response to injury. We carried out genomewide microRNA (miRNA) microarray studies during liver regeneration in rats after 70% partial hepatectomy (PH) at early and mid time points to more thoroughly understand their role. At 3, 12, and 18 hours post-PH ∼40% of the miRNAs tested were up-regulated. Conversely, at 24 hours post-PH, ∼70% of miRNAs were down-regulated. Furthermore, we established that the genomewide down-regulation of miRNA expression at 24 hours was also correlated with decreased expression of genes, such as Rnasen, Dgcr8, Dicer, Tarbp2, and Prkra, associated with miRNA biogenesis. To determine whether a potential negative feedback loop between miRNAs and their regulatory genes exists, 11 candidate miRNAs predicted to target the above-mentioned genes were examined and found to be up-regulated at 3 hours post-PH. Using reporter and functional assays, we determined that expression of these miRNA-processing genes could be regulated by a subset of miRNAs and that some miRNAs could target multiple miRNA biogenesis genes simultaneously. We also demonstrated that overexpression of these miRNAs inhibited cell proliferation and modulated cell cycle in both Huh-7 human hepatoma cells and primary rat hepatocytes. From these observations, we postulated that selective up-regulation of miRNAs in the early phase after PH was involved in the priming and commitment to liver regeneration, whereas the subsequent genomewide down-regulation of miRNAs was required for efficient recovery of liver cell mass. CONCLUSION: Our data suggest that miRNA changes are regulated by negative feedback loops between miRNAs and their regulatory genes that may play an important role in the steady-state regulation of liver regeneration.

Chromatin remodeling is required for gene reactivation after decitabine-mediated DNA hypomethylation.

The DNA hypomethylating drug decitabine (DAC) reactivates silenced gene expression in cancer and is approved for the treatment of the myelodysplastic syndrome. Gene reactivation after DAC is variable and incompletely understood. Here, we established a cell line system (YB5) derived from the SW48 colon cancer cell line to study DAC-induced reactivation. YB5 contains a hypermethylated cytomegalovirus promoter driving green fluorescent protein (GFP), and the locus is transcriptionally silent. GFP reexpression can be achieved by DAC treatment, but the expression level of individual cells is heterogeneous. DAC-treated YB5 cells were separated into GFP-positive and GFP-negative subpopulations. By comparing DAC-treated sorted GFP-positive and GFP-negative cells, we found that their methylation levels were similarly decreased but that histone modifications and histone H3 densities were remarkably different. Despite a similar degree of (incomplete) DNA hypomethylation, GFP-positive cells reverted to an active chromatin structure marked by higher H3K9 acetylation, lower H3K27 trimethylation, and lower promoter nucleosome density. GFP-negative cells had histone modifications and promoter nucleosome density, similar to parental cells. On DAC withdrawal, gradual resilencing and remethylation occurred in both GFP-positive and GFP-negative cells, and the resilencing correlated with a gradual increase in nucleosome occupancy in GFP-positive cells. These data show that hypomethylation alone after DAC is insufficient for gene expression induction, and that chromatin resetting to an active state including nucleosome eviction is required for activation of protein expression. Our findings suggest that gene expression is the key in optimizing DAC treatment strategies in the clinic.

Mechanisms of resistance to 5-aza-2'-deoxycytidine in human cancer cell lines.

5-aza-2'-deoxycytidine (DAC) is approved for the treatment of myelodysplastic syndromes, but resistance to this agent is common. In search for mechanisms of resistance, we measured the half maximal (50%) inhibitory concentration (IC(50)) of DAC and found it differed 1000-fold among a panel of cancer cell lines. The IC(50) was correlated with the doses of DAC that induced the most hypomethylation of long interspersed nuclear elements (LINE; R = 0.94, P < .001), but not with LINE methylation or DNA methyltransferase 1 (DNMT1), 3a, and 3b expression at baseline. Sensitivity to DAC showed a low correlation (R = 0.44, P = .11) to that of 5-azacytidine (AZA), but a good correlation to that of cytarabine (Ara-C; R = 0.89, P < .001). The 5 cell lines most resistant to DAC had a combination of low dCK, hENT1, and 2 transporters, and high cytosine deaminase. In an HL60 clone, resistance to DAC could be rapidly induced by drug exposure and was related to a switch from heterozygous to homozygous mutation of DCK. Transfection of wild-type DCK restored DAC sensitivity. DAC induced DNA breaks as evidenced by H2AX phosphorylation and increased homologous recombination rates by 7- to 10-fold. These results suggest that in vitro resistance to DAC can be explained by insufficient incorporation into DNA.

Genome-wide profiling of DNA methylation reveals a class of normally methylated CpG island promoters.

The role of CpG island methylation in normal development and cell differentiation is of keen interest, but remains poorly understood. We performed comprehensive DNA methylation profiling of promoter regions in normal peripheral blood by methylated CpG island amplification in combination with microarrays. This technique allowed us to simultaneously determine the methylation status of 6,177 genes, 92% of which include dense CpG islands. Among these 5,549 autosomal genes with dense CpG island promoters, we have identified 4.0% genes that are nearly completely methylated in normal blood, providing another exception to the general rule that CpG island methylation in normal tissue is limited to X inactivation and imprinted genes. We examined seven genes in detail, including ANKRD30A, FLJ40201, INSL6, SOHLH2, FTMT, C12orf12, and DPPA5. Dense promoter CpG island methylation and gene silencing were found in normal tissues studied except testis and sperm. In both tissues, bisulfite cloning and sequencing identified cells carrying unmethylated alleles. Interestingly, hypomethylation of several genes was associated with gene activation in cancer. Furthermore, reactivation of silenced genes could be induced after treatment with a DNA demethylating agent or in a cell line lacking DNMT1 and/or DNMT3b. Sequence analysis identified five motifs significantly enriched in this class of genes, suggesting that cis-regulatory elements may facilitate preferential methylation at these promoter CpG islands. We have identified a group of non-X-linked bona fide promoter CpG islands that are densely methylated in normal somatic tissues, escape methylation in germline cells, and for which DNA methylation is a primary mechanism of tissue-specific gene silencing.

Silencing of bidirectional promoters by DNA methylation in tumorigenesis.

CpG island methylation within promoters is known to silence individual genes in cancer. The involvement of this process in silencing gene pairs controlled by bidirectional promoters is unclear. In a screen for hypermethylated CpG islands in cancer, bidirectional promoters constituted 25.2% of all identified promoters, which matches with the genomic representation of bidirectional promoters. From the screen, we selected three bidirectional gene pairs for detailed analysis, WNT9A/CD558500, CTDSPL/BC040563, and KCNK15/BF195580. Levels of mRNA of all three pairs of genes were inversely correlated with the degree of promoter methylation in multiple cancer cell lines. Hypomethylation of these promoters induced by 5-aza-2'-deoxycytidine treatment reactivated or enhanced gene expression bidirectionally. The bidirectional nature of the WNT9A/CD558500 promoter was confirmed by luciferase assays, and hypermethylation down-regulated expression of both genes in the pair. Methylation of WNT9A/CD558500 and CTDSPL/BC040563 promoters occurs frequently in primary colon cancers and acute lymphoid leukemias (ALL), respectively, and methylation was correlated with decreased gene expression in ALL patient samples. Our study shows that hypermethylation of bidirectional promoter-associated CpG island silences two genes simultaneously, a property that should be taken into account when studying the functional consequences of hypermethylation in cancer.