Identification of an alternative short ARID5B isoform associated with B-ALL survival.

Utilizing RNA sequence (RNA-Seq) splice junction data from a cohort of 1841 B-cell acute lymphoblastic leukemia (B-ALL) patients we define transcriptionally distinct isoforms of ARID5B, a risk-associated gene identified in genome wide association studies (GWAS), which associate with disease survival. Short (S) and long (L) ARID5B transcripts, which differ in an encoded BAH-like chromatin interaction domain, show remarkable correlation to the isoform splicing pattern. Testing of the ARID5B proximal promoter of the S & L isoforms indicated that both are functionally independent in luciferase reporter assays. Increased short isoform expression is associated with decreased event-free and overall survival. The abundance of short and long transcripts strongly correlates to B-ALL prognostic stratification, where B-ALL subtypes with poor outcomes express a higher proportion of the S-isoform. These data demonstrate that the analysis of independent promoters and alternative splicing events are essential for improved risk stratification and a more complete understanding of disease pathology.

ARID5B regulates fatty acid metabolism and proliferation at the Pre-B cell stage during B cell development.

During B cell development in bone marrow, large precursor B cells (large Pre-B cells) proliferate rapidly, exit the cell cycle, and differentiate into non-proliferative (quiescent) small Pre-B cells. Dysregulation of this process may result in the failure to produce functional B cells and pose a risk of leukemic transformation. Here, we report that AT rich interacting domain 5B (ARID5B), a B cell acute lymphoblastic leukemia (B-ALL) risk gene, regulates B cell development at the Pre-B stage. In both mice and humans, we observed a significant upregulation of ARID5B expression that initiates at the Pre-B stage and is maintained throughout later stages of B cell development. In mice, deletion of Arid5b in vivo and ex vivo exhibited a significant reduction in the proportion of immature B cells but an increase in large and small Pre-B cells. Arid5b inhibition ex vivo also led to an increase in proliferation of both Pre-B cell populations. Metabolic studies in mouse and human bone marrow revealed that fatty acid uptake peaked in proliferative B cells then decreased during non-proliferative stages. We showed that Arid5b ablation enhanced fatty acid uptake and oxidation in Pre-B cells. Furthermore, decreased ARID5B expression was observed in tumor cells from B-ALL patients when compared to B cells from non-leukemic individuals. In B-ALL patients, ARID5B expression below the median was associated with decreased survival particularly in subtypes originating from Pre-B cells. Collectively, our data indicated that Arid5b regulates fatty acid metabolism and proliferation of Pre-B cells in mice, and reduced expression of ARID5B in humans is a risk factor for B cell leukemia.

Coordinated allele-specific histone acetylation at the differentially methylated regions of imprinted genes.

Genomic imprinting is an epigenetic inheritance system characterized by parental allele-specific gene expression. Allele-specific DNA methylation and chromatin composition are two epigenetic modification systems that control imprinted gene expression. To get a general assessment of histone lysine acetylation at imprinted genes we measured allele-specific acetylation of a wide range of lysine residues, H3K4, H3K18, H3K27, H3K36, H3K79, H3K64, H4K5, H4K8, H4K12, H2AK5, H2BK12, H2BK16 and H2BK46 at 11 differentially methylated regions (DMRs) in reciprocal mouse crosses using multiplex chromatin immunoprecipitation SNuPE assays. Histone acetylation marks generally distinguished the methylation-free alleles from methylated alleles at DMRs in mouse embryo fibroblasts and embryos. Acetylated lysines that are typically found at transcription start sites exhibited stronger allelic bias than acetylated histone residues in general. Maternally methylated DMRs, that usually overlap with promoters exhibited higher levels of acetylation and a 10% stronger allele-specific bias than paternally methylated DMRs that reside in intergenic regions. Along the H19/Igf2 imprinted domain, allele-specific acetylation at each lysine residue depended on functional CTCF binding sites in the imprinting control region. Our results suggest that many different histone acetyltransferase and histone deacetylase enzymes must act in concert in setting up and maintaining reciprocal parental allelic histone acetylation at DMRs.

Mice with Fabp4-Cre ablation of Arid5b are resistant to diet-induced obesity and hepatic steatosis.

Mice with global deletion of Arid5b expression are lean and resistant to diet-induced obesity, and Arid5b is required for adipogenesis in a variety of in vitro models. To determine whether the lean phenotype of Arid5b-/- mice can be explained by its absence in adipose tissues, we generated mice with Fabp4-mediated ablation of Arid5b. Arid5b expression was ablated in adipocytes, from Fabp4-CREpos; Arid5bFLOX/FLOX (FSKO) mice. FSKO mice were not lean when maintained on standard chow, but males were resistant to weight gains when placed on high-fat diets (HFD). This was mainly due to decreased lipid accumulation in subcutaneous (inguinal) white adipose tissue (IWAT), and the liver. Lipid accumulation proceeded normally in gonadal WAT (GWAT) and glucose intolerance developed to the same degree in FSKO and WT controls when subjected to HFD. CD68-positive macrophages were also significantly reduced in both inguinal and gonadal fat depots. RNA-Seq analysis of IWAT adipocytes from FSKO mice on HFD revealed significant decreases in the expression of genes associated with inflammation. Although Arid5b expression was normal in livers of FSKO mice, tissue weight gains and triglyceride accumulation, and expression of genes involved in lipid metabolism were markedly reduced in livers of FSKO mice on HFD. These results suggest that Arid5b plays a critical role in lipid accumulation in specific WAT depots, and in the inflammatory signaling from WAT depots to liver that lead to lipid accumulation and hepatic steatosis.

Meta-analysis of breast cancer microarray studies in conjunction with conserved cis-elements suggest patterns for coordinate regulation.

BACKGROUND: Gene expression measurements from breast cancer (BrCa) tumors are established clinical predictive tools to identify tumor subtypes, identify patients showing poor/good prognosis, and identify patients likely to have disease recurrence. However, diverse breast cancer datasets in conjunction with diagnostic clinical arrays show little overlap in the sets of genes identified. One approach to identify a set of consistently dysregulated candidate genes in these tumors is to employ meta-analysis of multiple independent microarray datasets. This allows one to compare expression data from a diverse collection of breast tumor array datasets generated on either cDNA or oligonucleotide arrays. RESULTS: We gathered expression data from 9 published microarray studies examining estrogen receptor positive (ER+) and estrogen receptor negative (ER-) BrCa tumor cases from the Oncomine database. We performed a meta-analysis and identified genes that were universally up or down regulated with respect to ER+ versus ER- tumor status. We surveyed both the proximal promoter and 3' untranslated regions (3'UTR) of our top-ranking genes in each expression group to test whether common sequence elements may contribute to the observed expression patterns. Utilizing a combination of known transcription factor binding sites (TFBS), evolutionarily conserved mammalian promoter and 3'UTR motifs, and microRNA (miRNA) seed sequences, we identified numerous motifs that were disproportionately represented between the two gene classes suggesting a common regulatory network for the observed gene expression patterns. CONCLUSION: Some of the genes we identified distinguish key transcripts previously seen in array studies, while others are newly defined. Many of the genes identified as overexpressed in ER- tumors were previously identified as expression markers for neoplastic transformation in multiple human cancers. Moreover, our motif analysis identified a collection of specific cis-acting target sites which may collectively play a role in the differential gene expression patterns observed in ER+ versus ER- breast cancer tumors. Importantly, the gene sets and associated DNA motifs provide a starting point with which to explore the mechanistic basis for the observed expression patterns in breast tumors.

Prevalence and type of BRCA mutations in Hispanics undergoing genetic cancer risk assessment in the southwestern United States: a report from the Clinical Cancer Genetics Community Research Network.

PURPOSE: To determine the prevalence and type of BRCA1 and BRCA2 (BRCA) mutations among Hispanics in the Southwestern United States and their potential impact on genetic cancer risk assessment (GCRA). PATIENTS AND METHODS: Hispanics (n = 746) with a personal or family history of breast and/or ovarian cancer were enrolled in an institutional review board-approved registry and received GCRA and BRCA testing within a consortium of 14 clinics. Population-based Hispanic breast cancer cases (n = 492) enrolled in the Northern California Breast Cancer Family Registry, negative by sequencing for BRCA mutations, were analyzed for the presence of the BRCA1 ex9-12del large rearrangement. RESULTS: Deleterious BRCA mutations were detected in 189 (25%) of 746 familial clinic patients (124 BRCA1, 65 BRCA2); 21 (11%) of 189 were large rearrangement mutations, of which 62% (13 of 21) were BRCA1 ex9-12del. Nine recurrent mutations accounted for 53% of the total. Among these, BRCA1 ex9-12del seems to be a Mexican founder mutation and represents 10% to 12% of all BRCA1 mutations in clinic- and population-based cohorts in the United States. CONCLUSION: BRCA mutations were prevalent in the largest study of Hispanic breast and/or ovarian cancer families in the United States to date, and a significant proportion were large rearrangement mutations. The high frequency of large rearrangement mutations warrants screening in every case. We document the first Mexican founder mutation (BRCA1 ex9-12del), which, along with other recurrent mutations, suggests the potential for a cost-effective panel approach to ancestry-informed GCRA.

Colitis locus on chromosome 2 impacting the severity of early-onset disease in mice deficient in GPX1 and GPX2.

BACKGROUND: Genetic background has a profound effect on inflammatory bowel disease. The Gpx1 and Gpx2 double knockout (GPX1/2-DKO) mice on a mixed C57BL/6 (B6) and 129S1/SvimJ (129) background exhibit spontaneous ileocolitis. The DKO mice on a B6 background have mild ileocolitis. We characterized the 129 DKO mice to identify a genetic locus affecting disease severity. METHODS: We backcrossed B6;129 DKO mice to 129 and analyzed for ileocolitis penetrance and severity at N5, N7, and N10. By correlating disease severity with single-nucleotide polymorphism (SNP) markers, we identified a colitis locus. RESULTS: As early as 9 days of age, 129 DKO N5 and N10 mice showed disease signs and morbidity. The N10 DKO mice had the severest colitis with nearly complete penetrance and high morbidity compared with other generations or backgrounds. 129 DKO mice had elevated colonic KC and SAA3 expression, shorter colon length, and cecal E. coli overgrowth compared to B6 DKO mice. Analysis of the B6 loci in 129 N5, N7, and N10 cohorts pointed to a region of chromosome 2: 119 Mbp contributing to mild symptoms. CONCLUSIONS: GPX1/2-DKO mice on 129 genetic background have the most aggressive colitis compared to B6;129 and B6 colonies. A B6 locus significantly contributing the resistance resides on chromosome 2: 119 Mbp. This region coincides with cytokine-deficiency-induced colitis susceptibility, Cdcs3, identified in the resistant B6 and sensitive C3H/HeJBir (C3Bir) with IL-10 deficiency. A three-way SNP analysis between 129, B6, and C3Bir locus points the major candidate genes to B2m, Dnajc17, Duox2, Pla2g4b, Pla2g4e, Pla2g4f and Slc30a4.

MIRA-SNuPE, a quantitative, multiplex method for measuring allele-specific DNA methylation.

5-methyl-C (5mC) and 5-hydroxymethyl-C (5hmC) are epigenetic marks with well known and putative roles in gene regulation, respectively. These two DNA covalent modifications cannot be distinguished by bisulfite sequencing or restriction digestion, the standard methods of 5mC detection. The methylated CpG island recovery assay (MIRA), however, specifically detects 5mC but not 5hmC. We further developed MIRA for the analysis of allele-specific CpG methylation at differentially methylated regions (DMRs) of imprinted genes. MIRA specifically distinguished between the parental alleles by capturing the paternally methylated H19/Igf2 DMR and maternally methylated KvDMR1 in mouse embryo fibroblasts (MEFs) carrying paternal and maternal duplication of mouse distal Chr7, respectively. MIRA in combination with multiplex single nucleotide primer extension (SNuPE) assays specifically captured the methylated parental allele from normal cells at a set of maternally and paternally methylated DMRs. The assay correctly recognized aberrant biallelic methylation in a case of loss-of imprinting. The MIRA-SNuPE assays revealed that placenta exhibited less DNA methylation bias at DMRs compared to yolk sac, amnion, brain, heart, kidney, liver and muscle. This method should be useful for the analysis of allele-specific methylation events related to genomic imprinting, X chromosome inactivation and for verifying and screening haplotype-associated methylation differences in the human population.

Three ways of combining genotyping and resequencing in case-control association studies.

We describe three statistical results that we have found to be useful in case-control genetic association testing. All three involve combining the discovery of novel genetic variants, usually by sequencing, with genotyping methods that recognize previously discovered variants. We first consider expanding the list of known variants by concentrating variant-discovery in cases. Although the naive inclusion of cases-only sequencing data would create a bias, we show that some sequencing data may be retained, even if controls are not sequenced. Furthermore, for alleles of intermediate frequency, cases-only sequencing with bias-correction entails little if any loss of power, compared to dividing the same sequencing effort among cases and controls. Secondly, we investigate more strongly focused variant discovery to obtain a greater enrichment for disease-related variants. We show how case status, family history, and marker sharing enrich the discovery set by increments that are multiplicative with penetrance, enabling the preferential discovery of high-penetrance variants. A third result applies when sequencing is the primary means of counting alleles in both cases and controls, but a supplementary pooled genotyping sample is used to identify the variants that are very rare. We show that this raises no validity issues, and we evaluate a less expensive and more adaptive approach to judging rarity, based on group-specific variants. We demonstrate the important and unusual caveat that this method requires equal sample sizes for validity. These three results can be used to more efficiently detect the association of rare genetic variants with disease.

Allele-specific H3K79 Di- versus trimethylation distinguishes opposite parental alleles at imprinted regions.

Imprinted gene expression corresponds to parental allele-specific DNA CpG methylation and chromatin composition. Histone tail covalent modifications have been extensively studied, but it is not known whether modifications in the histone globular domains can also discriminate between the parental alleles. Using multiplex chromatin immunoprecipitation-single nucleotide primer extension (ChIP-SNuPE) assays, we measured the allele-specific enrichment of H3K79 methylation and H4K91 acetylation along the H19/Igf2 imprinted domain. Whereas H3K79me1, H3K79me2, and H4K91ac displayed a paternal-specific enrichment at the paternally expressed Igf2 locus, H3K79me3 was paternally biased at the maternally expressed H19 locus, including the paternally methylated imprinting control region (ICR). We found that these allele-specific differences depended on CTCF binding in the maternal ICR allele. We analyzed an additional 11 differentially methylated regions (DMRs) and found that, in general, H3K79me3 was associated with the CpG-methylated alleles, whereas H3K79me1, H3K79me2, and H4K91ac enrichment was specific to the unmethylated alleles. Our data suggest that allele-specific differences in the globular histone domains may constitute a layer of the "histone code" at imprinted genes.

A risk variant in an miR-125b binding site in BMPR1B is associated with breast cancer pathogenesis.

MicroRNAs regulate diverse cellular processes and play an integral role in cancer pathogenesis. Genomic variation within miRNA target sites may therefore be important sources for genetic differences in cancer risk. To investigate this possibility, we mapped HapMap single nucleotide polymorphisms (SNP) to putative miRNA recognition sites within genes dysregulated in estrogen receptor-stratified breast tumors and used local linkage disequilibrium patterns to identify high-ranking SNPs in the Cancer Genetic Markers of Susceptibility (CGEMS) breast cancer genome-wide association study for further testing. Two SNPs, rs1970801 and rs11097457, scoring in the top 100 from the CGEMS study, were in strong linkage disequilibrium with rs1434536, an SNP that resides within a miR-125b target site in the 3' untranslated region of the bone morphogenic receptor type 1B (BMPR1B) gene encoding a transmembrane serine/threonine kinase. We validated the CGEMS association findings for rs1970801 in an independent cohort of admixture-corrected cases identified from families with multiple case histories. Subsequent association testing of rs1434536 for these cases and CGEMS controls with imputed genotypes supported the association. Furthermore, luciferase reporter assays and overexpression of miR-125b-mimics combined with quantitative reverse transcription-PCR showed that BMPR1B transcript is a direct target of miR-125b and that miR-125b differentially regulates the C and T alleles of rs1434536. These results suggest that allele-specific regulation of BMPR1B by miR-125b explains the observed disease risk. Our approach is general and can help identify and explain the mechanisms behind disease association for alleles that affect miRNA regulation.