Chromatin insulation orchestrates matrix metalloproteinase gene cluster expression reprogramming in aggressive breast cancer tumors.

BACKGROUND: Triple-negative breast cancer (TNBC) is an aggressive subtype that exhibits a high incidence of distant metastases and lacks targeted therapeutic options. Here we explored how the epigenome contributes to matrix metalloprotease (MMP) dysregulation impacting tumor invasion, which is the first step of the metastatic process. METHODS: We combined RNA expression and chromatin interaction data to identify insulator elements potentially associated with MMP gene expression and invasion. We employed CRISPR/Cas9 to disrupt the CCCTC-Binding Factor (CTCF) binding site on an insulator element downstream of the MMP8 gene (IE8) in two TNBC cellular models. We characterized these models by combining Hi-C, ATAC-seq, and RNA-seq with functional experiments to determine invasive ability. The potential of our findings to predict the progression of ductal carcinoma in situ (DCIS), was tested in data from clinical specimens. RESULTS: We explored the clinical relevance of an insulator element located within the Chr11q22.2 locus, downstream of the MMP8 gene (IE8). This regulatory element resulted in a topologically associating domain (TAD) boundary that isolated nine MMP genes into two anti-correlated expression clusters. This expression pattern was associated with worse relapse-free (HR = 1.57 [1.06 - 2.33]; p = 0.023) and overall (HR = 2.65 [1.31 - 5.37], p = 0.005) survival of TNBC patients. After CRISPR/Cas9-mediated disruption of IE8, cancer cells showed a switch in the MMP expression signature, specifically downregulating the pro-invasive MMP1 gene and upregulating the antitumorigenic MMP8 gene, resulting in reduced invasive ability and collagen degradation. We observed that the MMP expression pattern predicts DCIS that eventually progresses into invasive ductal carcinomas (AUC = 0.77, p < 0.01). CONCLUSION: Our study demonstrates how the activation of an IE near the MMP8 gene determines the regional transcriptional regulation of MMP genes with opposing functional activity, ultimately influencing the invasive properties of aggressive forms of breast cancer.

3-D chromatin conformation, accessibility, and gene expression profiling of triple-negative breast cancer.

OBJECTIVES: Triple-negative breast cancer (TNBC) is a highly aggressive breast cancer subtype with limited treatment options. Unlike other breast cancer subtypes, the scarcity of specific therapies and greater frequencies of distant metastases contribute to its aggressiveness. We aimed to find epigenetic changes that aid in the understanding of the dissemination process of these cancers. DATA DESCRIPTION: Using CRISPR/Cas9, our experimental approach led us to identify and disrupt an insulator element, IE8, whose activity seemed relevant for cell invasion. The experiments were performed in two well-established TNBC cellular models, the MDA-MB-231 and the MDA-MB-436. To gain insights into the underlying molecular mechanisms of TNBC invasion ability, we generated and characterized high-resolution chromatin interaction (Hi-C) and chromatin accessibility (ATAC-seq) maps in both cell models and complemented these datasets with gene expression profiling (RNA-seq) in MDA-MB-231, the cell line that showed more significant changes in chromatin accessibility. Altogether, our data provide a comprehensive resource for understanding the spatial organization of the genome in TNBC cells, which may contribute to accelerating the discovery of TNBC-specific alterations triggering advances for this devastating disease.

Index hopping on the Illumina HiseqX platform and its consequences for ancient DNA studies.

The high-throughput capacities of the Illumina sequencing platforms and the possibility to label samples individually have encouraged wide use of sample multiplexing. However, this practice results in read misassignment (usually <1%) across samples sequenced on the same lane. Alarmingly high rates of read misassignment of up to 10% were reported for lllumina sequencing machines with exclusion amplification chemistry. This may make use of these platforms prohibitive, particularly in studies that rely on low-quantity and low-quality samples, such as historical and archaeological specimens. Here, we use barcodes, short sequences that are ligated to both ends of the DNA insert, to directly quantify the rate of index hopping in 100-year old museum-preserved gorilla (Gorilla beringei) samples. Correcting for multiple sources of noise, we identify on average 0.470% of reads containing a hopped index. We show that sample-specific quantity of misassigned reads depends on the number of reads that any given sample contributes to the total sequencing pool, so that samples with few sequenced reads receive the greatest proportion of misassigned reads. This particularly affects ancient DNA samples, as these frequently differ in their DNA quantity and endogenous content. Through simulations we show that even low rates of index hopping, as reported here, can lead to biases in ancient DNA studies when multiplexing samples with vastly different quantities of endogenous material.

Regulation of the muscle-specific AMP-activated protein kinase alpha2beta2gamma3 complexes by AMP and implications of the mutations in the gamma3-subunit for the AMP dependence of the enzyme.

The AMP-activated protein kinase is an evolutionarily conserved heterotrimer that is important for metabolic sensing in all eukaryotes. The muscle-specific isoform of the regulatory gamma-subunit of the kinase, AMP-activated protein kinase gamma3, has a key role in glucose and fat metabolism in skeletal muscle, as suggested by metabolic characterization of humans, pigs and mice harboring substitutions in the AMP-binding Bateman domains of gamma3. We demonstrate that AMP-activated protein kinase alpha2beta2gamma3 trimers are allosterically activated approximately three-fold by AMP with a half-maximal stimulation (A(0.5)) at 1.9 +/- 0.5 or 2.6 +/- 0.3 microm, as measured for complexes expressed in Escherichia coli or mammalian cells, respectively. We show that mutations in the N-terminal Bateman domain of gamma3 (R225Q, H306R and R307G) increased the A(0.5) values for AMP, whereas the fold activation of the enzyme by 200 microm AMP remained unchanged in comparison to the wild-type complex. The mutations in the C-terminal Bateman domain of gamma3 (H453R and R454G), on the other hand, substantially reduced the fold stimulation of the complex by 200 microm AMP, and resulted in AMP dependence curves similar to those of the double mutant, R225Q/R454G. A V224I mutation in gamma3, known to result in a reduced glycogen content in pigs, did not affect the fold activation or the A(0.5) values for AMP. Importantly, we did not detect any increase in phosphorylation of Thr172 of alpha2 by the upstream kinases in the presence of increasing concentrations of AMP. Taken together, the data show that different mutations in gamma3 exert different effects on the allosteric regulation of the alpha2beta2gamma3 complex by AMP, whereas we find no evidence for their role in regulating the level of phosphorylation of alpha2 by upstream kinases.

Foxe3 haploinsufficiency in mice: a model for Peters' anomaly.

PURPOSE: To evaluate the importance in anterior segment dysgenesis of genetic variation in Foxe3, a gene encoding a forkhead transcription factor specifically expressed in the lens. METHODS: The phenotype of mice heterozygous for a mutation in the DNA-binding domain of Foxe3 was examined from histologic sections, and DNA binding by the encoded protein was investigated by gel-shift assay. FOXE3 from human patients with Peters' anomaly was PCR amplified and sequenced. RESULTS: The dysgenetic lens (dyl) allele of Foxe3 was found to encode a protein unable to bind DNA. Approximately 40% of mice heterozygous for Foxe3(dyl) have corneal and lenticular defects. The phenotype is variable but typically consists of the equivalent of Peters' anomaly in humans, with central corneal opacity, keratolenticular adhesion, and, in some cases, anterior polar cataract. In a small cohort (n = 13) of patients with Peters' anomaly, shown to be normal in the PAX6 locus, one individual was found to be heterozygous for a nonconservative missense mutation in FOXE3. The mutation, which does not occur in 116 chromosomes from a control population, substitutes leucine for arginine 90 at a highly conserved position in the forkhead domain. CONCLUSIONS: Haploinsufficiency of Foxe3 in a mouse model causes anterior segment dysgenesis similar to Peters' anomaly. Although causality could not be shown in the human case, the presence of a rare, nonconservative substitution in FOXE3 of a patient with Peters' anomaly is interesting, in light of the phenotypic similarities with the mutant mice.

A comparative gene expression database for invertebrates.

BACKGROUND: As whole genome and transcriptome sequencing gets cheaper and faster, a great number of 'exotic' animal models are emerging, rapidly adding valuable data to the ever-expanding Evo-Devo field. All these new organisms serve as a fantastic resource for the research community, but the sheer amount of data, some published, some not, makes detailed comparison of gene expression patterns very difficult to summarize - a problem sometimes even noticeable within a single lab. The need to merge existing data with new information in an organized manner that is publicly available to the research community is now more necessary than ever. DESCRIPTION: In order to offer a homogenous way of storing and handling gene expression patterns from a variety of organisms, we have developed the first web-based comparative gene expression database for invertebrates that allows species-specific as well as cross-species gene expression comparisons. The database can be queried by gene name, developmental stage and/or expression domains. CONCLUSIONS: This database provides a unique tool for the Evo-Devo research community that allows the retrieval, analysis and comparison of gene expression patterns within or among species. In addition, this database enables a quick identification of putative syn-expression groups that can be used to initiate, among other things, gene regulatory network (GRN) projects.

Relationship between AMPK and the transcriptional balance of clock-related genes in skeletal muscle.

Circadian clocks coordinate physiological, behavioral, and biochemical events with predictable daily environmental changes by a self-sustained transcriptional feedback loop. CLOCK and ARNTL are transcriptional activators that regulate Per and Cry gene expression. PER and CRY inhibit their own transcription, and their turnover allows this cycle to restart. The transcription factors BHLHB2 and BHLHB3 repress Per activation, whereas orphan nuclear receptors of the NR1D and ROR families control Arntl expression. Here we show the AMP-activated protein kinase (AMPK)gamma(3) subunit is involved in the regulation of peripheral circadian clock function. AMPKgamma3 knockout (Prkag3(-/-)) mice or wild-type littermates were injected with saline or an AMPK activator, 5-amino-4-imidazole-carboxamide riboside (AICAR), and white glycolytic gastrocnemius muscle was removed for gene expression analysis. Genes involved in the regulation of circadian rhythms (Cry2, Nr1d1, and Bhlhb2) were differentially regulated in response to AICAR in wild-type mice but remained unaltered in Prkag3(-/-) mice. Basal expression of Per1 was higher in Prkag3(-/-) mice compared with wild-type mice. Distinct diurnal changes in the respiratory exchange ratio (RER) between the light and dark phase of the day were observed in wild-type mice but not Prkag3(-/-) mice. In summary, the expression profile of clock-related genes in skeletal muscle in response to AICAR, as well as the diurnal shift in energy utilization, is impaired in AMPKgamma(3) subunit knockout mice. Our results indicate AMPK heterotrimeric complexes containing the AMPKgamma(3) subunit may play a specific role in linking circadian oscillators and energy metabolism in skeletal muscle.

Foxf1 and Foxf2 control murine gut development by limiting mesenchymal Wnt signaling and promoting extracellular matrix production.

Development of the vertebrate gut is controlled by paracrine crosstalk between the endodermal epithelium and the associated splanchnic mesoderm. In the adult, the same types of signals control epithelial proliferation and survival, which account for the importance of the stroma in colon carcinoma progression. Here, we show that targeting murine Foxf1 and Foxf2, encoding forkhead transcription factors, has pleiotropic effects on intestinal paracrine signaling. Inactivation of both Foxf2 alleles, or one allele each of Foxf1 and Foxf2, cause a range of defects, including megacolon, colorectal muscle hypoplasia and agangliosis. Foxf expression in the splanchnic mesoderm is activated by Indian and sonic hedgehog secreted by the epithelium. In Foxf mutants, mesenchymal expression of Bmp4 is reduced, whereas Wnt5a expression is increased. Activation of the canonical Wnt pathway -- with nuclear localization of beta-catenin in epithelial cells -- is associated with over-proliferation and resistance to apoptosis. Extracellular matrix, particularly collagens, is severely reduced in Foxf mutant intestine, which causes epithelial depolarization and tissue disintegration. Thus, Foxf proteins are mesenchymal factors that control epithelial proliferation and survival, and link hedgehog to Bmp and Wnt signaling.

Differences in the embryonic expression patterns of mouse Foxf1 and -2 match their distinct mutant phenotypes.

Murine genes encoding the forkhead transcription factors Foxf1 and -2 are both expressed in derivatives of the splanchnic mesoderm, i.e., the mesenchyme of organs derived from the primitive gut. In addition, Foxf2 is also expressed in limbs and the central nervous system. Targeted mutagenesis of Foxf1 and -2 suggests that Foxf1 is the more important of the two mammalian FoxF genes with early embryonic lethality of null embryos and a haploinsufficiency phenotype affecting foregut-derived organs. In contrast, the only reported defect in Foxf2 null embryos is cleft palate. To investigate if the differences in mutant phenotype can be attributed to nonoverlapping expression patterns or if distinct functions of the encoded proteins have to be inferred, we analyzed the early embryonic expression of Foxf2 and compared it with that of the better investigated Foxf1. We find that in the early embryo, Foxf1 is completely dominating-in terms of expression-in extraembryonic and lateral plate mesoderm, consistent with the malformations and early lethality of Foxf1 null mutants. Along the developing gut, Foxf1 is highly expressed throughout, whereas Foxf2 expression is concentrated to the posterior part-fitting the foregut haploinsufficiency phenotypes of Foxf1 mutants. Foxf2, on the other hand, is more prominent than Foxf1 in mesenchyme around the oral cavity, as would be predicted from the cleft palate phenotype. The differences in expression pattern also highlight areas where defects should be sought for in the Foxf2 mutant, for example limbs, the posterior gut, genitalia, and derivatives of the neural crest mesenchyme.