CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer.

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

Exome Sequencing Reveals a Phenotype Modifying Variant in ZNF528 in Primary Osteoporosis With a COL1A2 Deletion.

We studied a family with severe primary osteoporosis carrying a heterozygous p.Arg8Phefs*14 deletion in COL1A2, leading to haploinsufficiency. Three affected individuals carried the mutation and presented nearly identical spinal fractures but lacked other typical features of either osteogenesis imperfecta or Ehlers-Danlos syndrome. Although mutations leading to haploinsufficiency in COL1A2 are rare, mutations in COL1A1 that lead to less protein typically result in a milder phenotype. We hypothesized that other genetic factors may contribute to the severe phenotype in this family. We performed whole-exome sequencing in five family members and identified in all three affected individuals a rare nonsense variant (c.1282C > T/p.Arg428*, rs150257846) in ZNF528. We studied the effect of the variant using qPCR and Western blot and its subcellular localization with immunofluorescence. Our results indicate production of a truncated ZNF528 protein that locates in the cell nucleus as per the wild-type protein. ChIP and RNA sequencing analyses on ZNF528 and ZNF528-c.1282C > T indicated that ZNF528 binding sites are linked to pathways and genes regulating bone morphology. Compared with the wild type, ZNF528-c.1282C > T showed a global shift in genomic binding profile and pathway enrichment, possibly contributing to the pathophysiology of primary osteoporosis. We identified five putative target genes for ZNF528 and showed that the expression of these genes is altered in patient cells. In conclusion, the variant leads to expression of truncated ZNF528 and a global change of its genomic occupancy, which in turn may lead to altered expression of target genes. ZNF528 is a novel candidate gene for bone disorders and may function as a transcriptional regulator in pathways affecting bone morphology and contribute to the phenotype of primary osteoporosis in this family together with the COL1A2 deletion. © 2020 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Enhancer Dysfunction in 3D Genome and Disease.

Spatiotemporal patterns of gene expression depend on enhancer elements and other factors during individual development and disease progression. The rapid progress of high-throughput techniques has led to well-defined enhancer chromatin properties. Various genome-wide methods have revealed a large number of enhancers and the discovery of three-dimensional (3D) genome architecture showing the distant interacting mechanisms of enhancers that loop to target gene promoters. Whole genome sequencing projects directed at cancer have led to the discovery of substantial enhancer dysfunction in misregulating gene expression and in tumor initiation and progression. Results from genome-wide association studies (GWAS) combined with functional genomics analyses have elucidated the functional impacts of many cancer risk-associated variants that are enriched within the enhancer regions of chromatin. Risk variants dysregulate the expression of enhancer variant-associated genes via 3D genomic interactions. Moreover, these enhancer variants often alter the chromatin binding affinity for cancer-relevant transcription factors, which in turn leads to aberrant expression of the genes associated with cancer susceptibility. In this review, we investigate the extent to which these genetic regulatory circuits affect cancer predisposition and how the recent development of genome-editing methods have enabled the determination of the impacts of genomic variation and alteration on cancer phenotype, which will eventually lead to better management plans and treatment responses to human cancer in the clinic.

Oncogenic regulatory circuits driven by 19q13 rs11672691 underlies prostate cancer aggressiveness.

The 19q13 allele rs11672691 has been reproducibly found in association with aggressive form of prostate cancer, yet the underlying mechanism remains totally unknown. We have recently uncovered a mechanism by which rs11672691 influenced a novel oncogenic regulatory circuit, including HOXA2, PCAT19 and CEACAM21, thereby contributing to prostate cancer aggressiveness.

Biology and Clinical Implications of the 19q13 Aggressive Prostate Cancer Susceptibility Locus.

Genome-wide association studies (GWAS) have identified rs11672691 at 19q13 associated with aggressive prostate cancer (PCa). Here, we independently confirmed the finding in a cohort of 2,738 PCa patients and discovered the biological mechanism underlying this association. We found an association of the aggressive PCa-associated allele G of rs11672691 with elevated transcript levels of two biologically plausible candidate genes, PCAT19 and CEACAM21, implicated in PCa cell growth and tumor progression. Mechanistically, rs11672691 resides in an enhancer element and alters the binding site of HOXA2, a novel oncogenic transcription factor with prognostic potential in PCa. Remarkably, CRISPR/Cas9-mediated single-nucleotide editing showed the direct effect of rs11672691 on PCAT19 and CEACAM21 expression and PCa cellular aggressive phenotype. Clinical data demonstrated synergistic effects of rs11672691 genotype and PCAT19/CEACAM21 gene expression on PCa prognosis. These results provide a plausible mechanism for rs11672691 associated with aggressive PCa and thus lay the ground work for translating this finding to the clinic.

High-throughput screening of prostate cancer risk loci by single nucleotide polymorphisms sequencing.

Functional characterization of disease-causing variants at risk loci has been a significant challenge. Here we report a high-throughput single-nucleotide polymorphisms sequencing (SNPs-seq) technology to simultaneously screen hundreds to thousands of SNPs for their allele-dependent protein-binding differences. This technology takes advantage of higher retention rate of protein-bound DNA oligos in protein purification column to quantitatively sequence these SNP-containing oligos. We apply this technology to test prostate cancer-risk loci and observe differential allelic protein binding in a significant number of selected SNPs. We also test a unique application of self-transcribing active regulatory region sequencing (STARR-seq) in characterizing allele-dependent transcriptional regulation and provide detailed functional analysis at two risk loci (RGS17 and ASCL2). Together, we introduce a powerful high-throughput pipeline for large-scale screening of functional SNPs at disease risk loci.

A Large-Scale, Exome-Wide Association Study of Han Chinese Women Identifies Three Novel Loci Predisposing to Breast Cancer.

Genome-wide association studies have identified more than 90 susceptibility loci for breast cancer. However, the missing heritability is evident, and the contributions of coding variants to breast cancer susceptibility have not yet been systematically evaluated. Here, we present a large-scale whole-exome association study for breast cancer consisting of 24,162 individuals (10,055 cases and 14,107 controls). In addition to replicating known susceptibility loci (e.g., ESR1, FGFR2, and TOX3), we identify two novel missense variants in C21orf58 (rs13047478, Pmeta = 4.52 × 10-8) and ZNF526 (rs3810151, Pmeta = 7.60 × 10-9) and one new noncoding variant at 7q21.11 (P < 5 × 10-8). C21orf58 and ZNF526 possessed functional roles in the control of breast cancer cell growth, and the two coding variants were found to be the eQTL for several nearby genes. rs13047478 was significantly (P < 5.00 × 10-8) associated with the expression of genes MCM3AP and YBEY in breast mammary tissues. rs3810151 was found to be significantly associated with the expression of genes PAFAH1B3 (P = 8.39 × 10-8) and CNFN (P = 3.77 × 10-4) in human blood samples. C21orf58 and ZNF526, together with these eQTL genes, were differentially expressed in breast tumors versus normal breast. Our study reveals additional loci and novel genes for genetic predisposition to breast cancer and highlights a polygenic basis of disease development.Significance: Large-scale genetic screening identifies novel missense variants and a noncoding variant as predisposing factors for breast cancer. Cancer Res; 78(11); 3087-97. ©2018 AACR.

Hyperinsulinemia shifted energy supply from glucose to ketone bodies in early nonalcoholic steatohepatitis from high-fat high-sucrose diet induced Bama minipigs.

The minipig can serve as a good pharmacological model for human subjects. However, the long-term pathogenesis of high-calorie diet-induced metabolic syndromes, including NASH, has not been well described in minipigs. We examined the development of metabolic syndromes in Bama minipigs that were fed a high-fat, high-sucrose diet (HFHSD) for 23 months, by using histology and serum biochemistry and by profiling the gene expression patterns in the livers of HFHSD pigs compared to controls. The pathology findings revealed microvesicular steatosis, iron overload, arachidonic acid synthesis, lipid peroxidation, reduced antioxidant capacity, increased cellular damage, and inflammation in the liver. RNA-seq analysis revealed that 164 genes were differentially expressed between the livers of the HFHSD and control groups. The pathogenesis of early-stage NASH was characterized by hyperinsulinemia and by de novo synthesis of fatty acids and nascent triglycerides, which were deposited as lipid droplets in hepatocytes. Hyperinsulinemia shifted the energy supply from glucose to ketone bodies, and the high ketone body concentration induced the overexpression of cytochrome P450 2E1 (CYP2E1). The iron overload, CYP2E1 and alcohol dehydrogenase 4 overexpression promoted reactive oxygen species (ROS) production, which resulted in arachidonic and linoleic acid peroxidation and, in turn, led to malondialdehyde production and a cellular response to ROS-mediated DNA damage.

Long-term High-fat High-sucrose Diet Promotes Enlarged Islets and ß-Cell Damage by Oxidative Stress in Bama Minipigs.

OBJECTIVES: The effect of a long-term high-fat, high-caloric diet on the dysfunction of pancreas has not been clarified. We investigated the pancreatic histopathology and ß-cell apoptosis in Bama minipigs after 23 months on a high-fat high-sucrose diet (HFHSD). METHODS: Bama minipigs were randomly assigned to control (n = 6) and HFHSD groups (n = 6) for 23 months, and biochemical parameters were measured. Pancreata were subjected to histological analysis, followed by assessment with transmission electron microscopy. Lipid peroxidation was determined by the malondialdehyde concentration and antioxidant enzyme activity. Β-cell apoptosis was measured by an immunohistochemical method. RESULTS: In the HFHSD group, the islets were enlarged, and the pancreatic tissue had observed significant fatty infiltration. Moreover, the feeding program damaged the normal pancreatic tissue structure. The level of lipid peroxidation was increased, and the activities of pancreatic antioxidant enzymes were significantly decreased. The expression levels of caspase-3, Bax, and insulin were significantly increased (P < 0.05), and the expression levels of proliferating cell nuclear antigen and Bcl-2 were decreased (P < 0.05). CONCLUSIONS: The long-term HFHSD promotes pancreatic steatosis and oxidative stress, which increases ß-cell apoptosis as indicated by the activation of caspase-3 through the mitochondrial pathway (Bcl-2/Bax).

The study of regulatory effects of Pdx-1, MafA and NeuroD1 on the activity of porcine insulin promoter and the expression of human islet amyloid polypeptide.

The purpose of the present study was to determine the activation of porcine insulin promoter (PIP) by three transcription factors: pancreatic and duodenal homeobox 1 (Pdx-1), v-maf musculoaponeurotic fibrosarcoma oncogene (MafA) and neurogenic differentiation 1 (NeuroD1) in non-beta islet cells cultured in vitro. In addition, the expression of the exogenous human islet amyloid polypeptide (hIAPP) gene driving by PIP in porcine kidney 15 (PK15) cells co-transfected with these transcription factors was also examined. In the present study, a series of vectors for gene overexpression were constructed, including pGL3-Pdx-1, pGL3-MafA, pGL3-NeuroD1, pGL3-PIP-LUC and pcDNA3.1-PIP-hIAPP. The dual-luciferase reporter assay showed that the PIP activity was increased in PK15 cells when overexpressing the exogenous transcription factors Pdx-1, MafA and NeuroD1. Introducing the PIP-hIAPP expression vector into PK15 cells combined with exogenous Pdx-1, MafA and NeuroD1 resulted in the efficient expression of hIAPP at the gene level, but not the protein. The current systematic porcine insulin promoter analysis provided the basic information for future utilization of porcine insulin.