Looping mediated interaction between the promoter and 3' UTR regulates type II collagen expression in chondrocytes

Cell layer-specific expression of the homeotic MADS-box transcription factor PhDEF contributes to modular petal morphogenesis in petunia

Floral homeotic MADS-box transcription factors ensure the correct morphogenesis of floral organs, which are organized in different cell layers deriving from distinct meristematic layers. How cells from these distinct layers acquire their respective identities and coordinate their growth to ensure normal floral organ morphogenesis is unresolved. Here, we studied petunia (Petunia × hybrida) petals that form a limb and tube through congenital fusion. We identified petunia mutants (periclinal chimeras) expressing the B-class MADS-box gene DEFICIENS in the petal epidermis or in the petal mesophyll, called wico and star, respectively. Strikingly, wico flowers form a strongly reduced tube while their limbs are almost normal, while star flowers form a normal tube but greatly reduced and unpigmented limbs, showing that petunia petal morphogenesis is highly modular. These mutants highlight the layer-specific roles of PhDEF during petal development. We explored the link between PhDEF and petal pigmentation, a well-characterized limb epidermal trait. The anthocyanin biosynthesis pathway was strongly downregulated in star petals, including its major regulator ANTHOCYANIN2 (AN2). We established that PhDEF directly binds to the AN2 terminator in vitro and in vivo, suggesting that PhDEF might regulate AN2 expression and therefore petal epidermis pigmentation. Altogether, we show that cell layer-specific homeotic activity in petunia petals differently impacts tube and limb development, revealing the relative importance of the different cell layers in the modular architecture of petunia petals.

Towards a comprehensive regulatory map of Mammalian Genomes

Genome mapping studies have generated a nearly complete collection of genes for the human genome, but we still lack an equivalently vetted inventory of human regulatory sequences. Cis-regulatory modules (CRMs) play important roles in controlling when, where, and how much a gene is expressed. We developed a training data-free CRM-prediction algorithm, the Mammalian Regulatory MOdule Detector (MrMOD) for accurate CRM prediction in mammalian genomes. MrMOD provides genome position-fixed CRM models similar to the fixed gene models for the mouse and human genomes using only genomic sequences as the inputs with one adjustable parameter - the significance p-value. Importantly, MrMOD predicts a comprehensive set of high-resolution CRMs in the mouse and human genomes including all types of regulatory modules not limited to any tissue, cell type, developmental stage, or condition. We computationally validated MrMOD predictions used a compendium of 21 orthogonal experimental data sets including thousands of experimentally defined CRMs and millions of putative regulatory elements derived from hundreds of different tissues, cell types, and stimulus conditions obtained from multiple databases. In ovo transgenic reporter assay demonstrates the power of our prediction in guiding experimental design. We analyzed CRMs located in the chromosome 17 using unsupervised machine learning and identified groups of CRMs with multiple lines of evidence supporting their functionality, linking CRMs with upstream binding transcription factors and downstream target genes. Our work provides a comprehensive base pair resolution annotation of the functional regulatory elements and non-functional regions in the mammalian genomes.

PLK2 Single Nucleotide Variant in Gastric Cancer Patients Affects miR-23b-5p Binding

PURPOSE

Chromosomal instability is a hallmark of gastric cancer (GC). It can be driven by single nucleotide variants (SNVs) in cell cycle genes. We investigated the associations between SNVs in candidate genes, PLK2, PLK3, and ATM, and GC risk and clinicopathological features.

MATERIALS AND METHODS

The genotyping study included 542 patients with GC and healthy controls. Generalized linear models were used for the risk and clinicopathological association analyses. Survival analysis was performed using the Kaplan-Meier method. The binding of candidate miRs was analyzed using a luciferase reporter assay.

RESULTS

The PLK2 Crs15009-Crs963615 haplotype was under-represented in the GC group compared to that in the control group (Pcorr=0.050). Male patients with the PLK2 rs963615 CT genotype had a lower risk of GC, whereas female patients had a higher risk (P=0.023; P=0.026). The PLK2 rs963615 CT genotype was associated with the absence of vascular invasion (P=0.012). The PLK3 rs12404160 AA genotype was associated with a higher risk of GC in the male population (P=0.015). The ATM Trs228589-Ars189037-Grs4585 haplotype was associated with a higher risk of GC (P<0.001). The ATM rs228589, rs189037, and rs4585 genotypes TA+AA, AG+GG, and TG+GG were associated with the absence of perineural invasion (P=0.034). In vitro analysis showed that the cancer-associated miR-23b-5p mimic specifically bound to the PLK2 rs15009 G allele (P=0.0097). Moreover, low miR-23b expression predicted longer 10-year survival (P=0.0066) in patients with GC.

CONCLUSIONS

PLK2, PLK3, and ATM SNVs could potentially be helpful for the prediction of GC risk and clinicopathological features. PLK2 rs15009 affects the binding of miR-23b-5p. MiR-23b-5p expression status could serve as a prognostic marker for survival in patients with GC.

An Unbiased Molecular Approach Using 3'-UTRs Resolves the Avian Family-Level Tree of Life

Presumably, due to a rapid early diversification, major parts of the higher-level phylogeny of birds are still resolved controversially in different analyses or are considered unresolvable. To address this problem, we produced an avian tree of life, which includes molecular sequences of one or several species of ∼90% of the currently recognized family-level taxa (429 species, 379 genera) including all 106 family-level taxa of the nonpasserines and 115 of the passerines (Passeriformes). The unconstrained analyses of noncoding 3-prime untranslated region (3′-UTR) sequences and those of coding sequences yielded different trees. In contrast to the coding sequences, the 3′-UTR sequences resulted in a well-resolved and stable tree topology. The 3′-UTR contained, unexpectedly, transcription factor binding motifs that were specific for different higher-level taxa. In this tree, grebes and flamingos are the sister clade of all other Neoaves, which are subdivided into five major clades. All nonpasserine taxa were placed with robust statistical support including the long-time enigmatic hoatzin (Opisthocomiformes), which was found being the sister taxon of the Caprimulgiformes. The comparatively late radiation of family-level clades of the songbirds (oscine Passeriformes) contrasts with the attenuated diversification of nonpasseriform taxa since the early Miocene. This correlates with the evolution of vocal production learning, an important speciation factor, which is ancestral for songbirds and evolved convergent only in hummingbirds and parrots. As 3′-UTR-based phylotranscriptomics resolved the avian family-level tree of life, we suggest that this procedure will also resolve the all-species avian tree of life

Crosstalk of promoter and terminator during RNA polymerase II transcription cycle

The transcription cycle of RNAPII is comprised of three consecutive steps; initiation, elongation and termination. It has been assumed that the initiation and termination steps occur in spatial isolation, essentially as independent events. A growing body of evidence, however, has challenged this dogma. First, factors involved in initiation and termination exhibit both a genetic and a physical interaction during transcription. Second, the initiation and termination factors have been found to occupy both ends of a transcribing gene. Third, physical interaction of initiation and termination factors occupying distal ends of a gene sometime results in the entire terminator region of a genes looping back and contact its cognate promoter, thereby forming a looped gene architecture during transcription. A logical interpretation of these findings is that the initiation and termination steps of transcription do not occur in isolation. There is extensive communication of factors occupying promoter and terminator ends of a gene during transcription cycle. This review entails a discussion of the promoter-terminator crosstalk and its implication in the context of transcription.

Enhancer Function of MicroRNA-3681 Derived from Long Terminal Repeats Represses the Activity of Variable Number Tandem Repeats in the 3' UTR of SHISA7

microRNAs (miRNAs) are non-coding RNA molecules involved in the regulation of gene expression. miRNAs inhibit gene expression by binding to the 3' untranslated region (UTR) of their target gene. miRNAs can originate from transposable elements (TEs), which comprise approximately half of the eukaryotic genome and one type of TE, called the long terminal repeat (LTR) is found in class of retrotransposons. Amongst the miRNAs derived from LTR, hsa-miR-3681 was chosen and analyzed using bioinformatics tools and experimental analysis. Studies on hsa-miR-3681 have been scarce and this study provides the relative expression analysis of hsa-miR-3681-5p from humans, chimpanzees, crab-eating monkeys, and mice. Luciferase assay for hsa-miR-3681-5p and its target gene SHISA7 supports our hypothesis that the number of miRNA binding sites affects target gene expression. Especially, the variable number tandem repeat (VNTR) and hsa-miR-3681-5p share the binding sites in the 3' UTR of SHISA7, which leads the enhancer function of hsa-miR-3681-5p to inhibit the activity of VNTR. In conclusion, hsa-miR-3681-5p acts as a super-enhancer and the enhancer function of hsa-miR-3681-5p acts as a repressor of VNTR activity in the 3' UTR of SHISA7.

Step-Wise Chondrogenesis of Human Induced Pluripotent Stem Cells and Purification Via a Reporter Allele Generated by CRISPR-Cas9 Genome Editing

The differentiation of human induced pluripotent stem cells (hiPSCs) to prescribed cell fates enables the engineering of patient-specific tissue types, such as hyaline cartilage, for applications in regenerative medicine, disease modeling, and drug screening. In many cases, however, these differentiation approaches are poorly controlled and generate heterogeneous cell populations. Here, we demonstrate cartilaginous matrix production in three unique hiPSC lines using a robust and reproducible differentiation protocol. To purify chondroprogenitors (CPs) produced by this protocol, we engineered a COL2A1-GFP knock-in reporter hiPSC line by CRISPR-Cas9 genome editing. Purified CPs demonstrated an improved chondrogenic capacity compared with unselected populations. The ability to enrich for CPs and generate homogenous matrix without contaminating cell types will be essential for regenerative and disease modeling applications. Stem Cells 2019;37:65-76.

Plau and Tgfbr3 are YAP-regulated genes that promote keratinocyte proliferation

Yes-associated protein (YAP) is a mechanosensor protein and a downstream effector of the Hippo kinase pathway, which controls organ growth, cell proliferation, survival, maintenance and regeneration. Unphosphorylated YAP translocates to the nucleus where it acts as a cofactor of primarily the TEAD transcription factors to activate target gene transcription and cell proliferation. Perturbed YAP activation results in tumorigenesis. The pathways downstream of activated YAP that drive cell proliferation remain relatively unexplored. In this study, we employed YAP2-5SA-∆C transgenic mice, which overexpress a mildly activated YAP mutant protein in basal keratinocytes leading to increased proliferation of the epidermal stem/progenitor cell populations. We performed massively-parallel sequencing of skin biopsy mRNA (RNA-Seq) and found dysregulation of 1491 genes in YAP2-5SA-∆C skin, including many with roles in cell activation and proliferation. Furthermore, we found that 150 of these dysregulated genes harbored YAP/TEAD binding motifs in the 3′ UTR, suggesting that these may be direct YAP/TEAD target genes in the control of epidermal regeneration. Further validation and functional characterization assays identified Plau and Tgfbr3 as prime candidate genes that may be activated by epidermal YAP activity in the mouse skin in vivo to promote keratinocyte proliferation. This study provides novel insights into the mechanisms regulated by YAP that control tissue homeostasis, and in particular in conditions where YAP is aberrantly activated such as in neoplastic and regenerative skin disease.

Oncogenic p95HER2 regulates Na+-HCO3- cotransporter NBCn1 mRNA stability in breast cancer cells via 3'UTR-dependent processes

The Na⁺-HCO₃^- cotransporter NBCn1 (SLC4A7) is up-regulated in breast cancer, important for tumor growth, and a single nucleotide polymorphism (SNP), rs4973768, in its 3' untranslated region (3'UTR) correlates with increased breast cancer risk. We previously demonstrated that NBCn1 expression and promoter activity are strongly increased in breast cancer cells expressing a constitutively active oncogenic human epidermal growth factor receptor 2 (HER2) (p95HER2). Here, we address the roles of p95HER2 in regulating NBCn1 expression via post-transcriptional mechanisms. p95HER2 expression in MCF-7 cells reduced the rate of NBCn1 mRNA degradation. The NBCn1 3'UTR down-regulated luciferase reporter expression in control cells, and this was reversed by p95HER2, suggesting that p95HER2 counteracts 3'UTR-mediated suppression of NBCn1 expression. Truncation analyses identified three NBCn1 3'UTR regions of regulatory importance. Mutation of putative miRNA-binding sites (miR-374a/b, miR-200b/c, miR-29a/b/c, miR-488) in these regions did not have significant impact on 3'UTR activity. The NBCn1 3'UTR interacted directly with the RNA-binding protein human antigen R (HuR), and HuR knockdown reduced NBCn1 expression. Conversely, ablation of a distal AU-rich element increased 3'UTR-driven reporter activity, suggesting complex regulatory roles of these sites. The cancer-associated SNP variant decreased reporter expression in T-47D breast cancer cells, yet not in MCF-7, MDA-MB-231 and SK-BR-3 cells, arguing against a general role in regulating NBCn1 expression. Finally, p95HER2 expression increased total and plasma membrane NBCn1 protein levels and decreased the rate of NBCn1 protein degradation. Collectively, this is the first work to demonstrate 3'UTR-mediated NBCn1 regulation, shows that p95HER2 regulates NBCn1 expression at multiple levels, and substantiates the central position of p95HER2-NBCn1 signaling in breast cancer.

Purifying selection acts on coding and non-coding sequences of paralogous genes in Arabidopsis thaliana

Background

Whole-genome duplications in the ancestors of many diverse species provided the genetic material for evolutionary novelty. Several models explain the retention of paralogous genes. However, how these models are reflected in the evolution of coding and non-coding sequences of paralogous genes is unknown.

Results

Here, we analyzed the coding and non-coding sequences of paralogous genes in Arabidopsis thaliana and compared these sequences with those of orthologous genes in Arabidopsis lyrata. Paralogs with lower expression than their duplicate had more nonsynonymous substitutions, were more likely to fractionate, and exhibited less similar expression patterns with their orthologs in the other species. Also, lower-expressed genes had greater tissue specificity. Orthologous conserved non-coding sequences in the promoters, introns, and 3′ untranslated regions were less abundant at lower-expressed genes compared to their higher-expressed paralogs. A gene ontology (GO) term enrichment analysis showed that paralogs with similar expression levels were enriched in GO terms related to ribosomes, whereas paralogs with different expression levels were enriched in terms associated with stress responses.

Conclusions

Loss of conserved non-coding sequences in one gene of a paralogous gene pair correlates with reduced expression levels that are more tissue specific. Together with increased mutation rates in the coding sequences, this suggests that similar forces of purifying selection act on coding and non-coding sequences. We propose that coding and non-coding sequences evolve concurrently following gene duplication.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2803-2) contains supplementary material, which is available to authorized users.

Common Polymorphisms in the Solute Carrier SLC30A10 are Associated With Blood Manganese and Neurological Function

Manganese (Mn) is an essential nutrient in humans, but excessive exposure to Mn may cause neurotoxicity. Despite homeostatic regulation, Mn concentrations in blood vary considerably among individuals. We evaluated if common single-nucleotide polymorphisms (SNPs) in SLC30A10, which likely encodes an Mn transporter, influence blood Mn concentrations and neurological function. We measured blood Mn concentrations by ICP-MS or atomic absorption spectroscopy and genotyped 2 SLC30A10 non-coding SNPs (rs2275707 and rs12064812) by TaqMan PCR in cohorts from Bangladesh (N = 406), the Argentinean Andes (N = 198), and Italy (N = 238). We also measured SLC30A10 expression in whole blood by TaqMan PCR in a sub-group (N = 101) from the Andean cohort, and neurological parameters (sway velocity and finger-tapping speed) in the Italian cohort. The rs2275707 variant allele was associated with increased Mn concentrations in the Andes (8%, P = .027) and Italy (10.6%, P = .012), but not as clear in Bangladesh (3.4%, P = .21; linear regression analysis adjusted for age, gender, and plasma ferritin). This allele was also associated with increased sway velocity (15%, P = .033; adjusted for age and sex) and reduced SLC30A10 expression (−24.6%, P = .029). In contrast, the rs12064812 variant homozygous genotype was associated with reduced Mn concentrations, particularly in the Italian cohort (−18.4%, P = .04), and increased finger-tapping speed (8.7%, P = .025). We show that common SNPs in SLC30A10 are associated with blood Mn concentrations in 3 unrelated cohorts and that their influence may be mediated by altered SLC30A10 expression. Moreover, the SNPs appeared to influence neurological functions independent of blood Mn concentrations, suggesting that SLC30A10 could regulate brain Mn levels.

Chondrogenic potential of Slug-depleted human mesenchymal stem cells

The use of short interfering RNA (siRNA) in combination with stem cells and biocompatible scaffolds is a promising strategy in regenerative medicine. Our experimental strategy was to explore the possibility of forcing or guiding the chondrogenic differentiation of human mesenchymal stem cells (hMSCs) by knocking down a negative regulator of chondrogenesis, Slug transcription factor (TF), thus altering cell behavior. We found that TGFβ-driven chondrogenic differentiation of hMSCs cultured onto a hyaluronan-based scaffold, HYAFF(®)-11, was strengthened after cell exposure to siRNA against Slug. Slug silencing was effective in promoting the expression of chondrogenic markers, including Col2A1, aggrecan, Sox9, LEF1, and TRPS1. In addition, we confirmed that HYAFF-11 is a good scaffold candidate for hMSC use in tissue engineering applications, and showed that it is effective in sustaining TGFβ3 treatment associated with a specific gene silencing. Interestingly, preliminary results from the experimental model described here suggested that, even in the absence of differentiation supplements, Slug silencing showed a pro-chondrogenic effect, highlighting both its potential use as an alternative to TGFβ treatment, and the critical role of the Slug TF in determining the fate of hMSCs.

A molecular basis for classic blond hair color in Europeans

Hair color differences are among the most obvious examples of phenotypic variation in humans. Although genome-wide association studies (GWAS) have implicated multiple loci in human pigment variation, the causative base-pair changes are still largely unknown. Here we dissect a regulatory region of the KITLG gene (encoding KIT ligand) that is significantly associated with common blond hair color in northern Europeans. Functional tests demonstrate that the region contains a regulatory enhancer that drives expression in developing hair follicles. This enhancer contains a common SNP (rs12821256) that alters a binding site for the lymphoid enhancer-binding factor 1 (LEF1) transcription factor, reducing LEF1 responsiveness and enhancer activity in cultured human keratinocytes. Mice carrying ancestral or derived variants of the human KITLG enhancer exhibit significant differences in hair pigmentation, confirming that altered regulation of an essential growth factor contributes to the classic blond hair phenotype found in northern Europeans.

The pattern and evolution of looped gene bendability

Gene looping, defined as the physical interaction between the promoter and terminator regions of a RNA polymerase II-transcribed gene, is widespread in yeast and mammalian cells. Gene looping has been shown to play important roles in transcription. Gene-loop formation is dependent on regulatory proteins localized at the 5' and 3' ends of genes, such as TFIIB. However, whether other factors contribute to gene looping remains to be elucidated. Here, we investigated the contribution of intrinsic DNA and chromatin structures to gene looping. We found that Saccharomyces cerevisiae looped genes show high DNA bendability around middle and 3/4 regions in open reading frames (ORFs). This bendability pattern is conserved between yeast species, whereas the position of bendability peak varies substantially among species. Looped genes in human cells also show high DNA bendability. Nucleosome positioning around looped ORF middle regions is unstable. We also present evidence indicating that this unstable nucleosome positioning is involved in gene looping. These results suggest a mechanism by which DNA bendability and unstable nucleosome positioning could assist in the formation of gene loops.

Chromatin loop organization of the junb locus in mouse dendritic cells

The junb gene behaves as an immediate early gene in bacterial lipopolysaccharide (LPS)-stimulated dendritic cells (DCs), where its transient transcriptional activation is necessary for the induction of inflammatory cytokines. junb is a short gene and its transcriptional activation by LPS depends on the binding of NF-κB to an enhancer located just downstream of its 3′ UTR. Here, we have addressed the mechanisms underlying the transcriptional hyper-reactivity of junb. Using transfection and pharmacological assays to complement chromatin immunoprecipitation analyses addressing the localization of histones, polymerase II, negative elongation factor (NELF)-, DRB sensitivity-inducing factor (DSIF)- and Positive Transcription Factor b complexes, we demonstrate that junb is a RNA Pol II-paused gene where Pol II is loaded in the transcription start site domain but poorly active. Moreover, High salt-Recovered Sequence, chromosome conformation capture (3C)- and gene transfer experiments show that (i) junb is organized in a nuclear chromatin loop bringing into close spatial proximity the upstream promoter region and the downstream enhancer and (ii) this configuration permits immediate Pol II release on the junb body on binding of LPS-activated NF-κB to the enhancer. Thus, our work unveils a novel topological framework underlying fast junb transcriptional response in DCs. Moreover, it also points to a novel layer of complexity in the modes of action of NF-κB.

Polysome profiling shows extensive posttranscriptional regulation during human adipocyte stem cell differentiation into adipocytes

Adipocyte stem cells (hASCs) can proliferate and self-renew and, due to their multipotent nature, they can differentiate into several tissue-specific lineages, making them ideal candidates for use in cell therapy. Most attempts to determine the mRNA profile of self-renewing or differentiating stem cells have made use of total RNA for gene expression analysis. Several lines of evidence suggest that self-renewal and differentiation are also dependent on the control of protein synthesis by posttranscriptional mechanisms. We used adipogenic differentiation as a model, to investigate the extent to which posttranscriptional regulation controlled gene expression in hASCs. We focused on the initial steps of differentiation and isolated both the total mRNA fraction and the subpopulation of mRNAs associated with translating ribosomes. We observed that adipogenesis is committed in the first days of induction and three days appears as the minimum time of induction necessary for efficient differentiation. RNA-seq analysis showed that a significant percentage of regulated mRNAs were posttranscriptionally controlled. Part of this regulation involves massive changes in transcript untranslated regions (UTR) length, with differential extension/reduction of the 3'UTR after induction. A slight correlation can be observed between the expression levels of differentially expressed genes and the 3'UTR length. When we considered association to polysomes, this correlation values increased. Changes in the half lives were related to the extension of the 3'UTR, with longer UTRs mainly stabilizing the transcripts. Thus, changes in the length of these extensions may be associated with changes in the ability to associate with polysomes or in half-life.

Identification of thyroid hormone response elements in vivo using mice expressing a tagged thyroid hormone receptor α1

TRα1 (thyroid hormone receptor α1) is well recognized for its importance in brain development. However, due to the difficulties in predicting TREs (thyroid hormone response elements) in silico and the lack of suitable antibodies against TRα1 for ChIP (chromatin immunoprecipitation), only a few direct TRα1 target genes have been identified in the brain. Here we demonstrate that mice expressing a TRα1–GFP (green fluorescent protein) fusion protein from the endogenous TRα locus provide a valuable animal model to identify TRα1 target genes. To this end, we analysed DNA–TRα1 interactions in vivo using ChIP with an anti-GFP antibody. We validated our system using established TREs from neurogranin and hairless, and by verifying additional TREs from known TRα1 target genes in brain and heart. Moreover, our model system enabled the identification of novel TRα1 target genes such as RNF166 (ring finger protein 166). Our results demonstrate that transgenic mice expressing a tagged nuclear receptor constitute a feasible approach to study receptor–DNA interactions in vivo, circumventing the need for specific antibodies. Models like the TRα1–GFP mice may thus pave the way for genome-wide mapping of nuclear receptor-binding sites, and advance the identification of novel target genes in vivo.

The transcription factor activating enhancer-binding protein epsilon (AP-2ε) regulates the core promoter of type II collagen (COL2A1)

The transcription factor activating enhancer-binding protein epsilon (AP-2ε) was recently shown to be expressed during late chondrocyte differentiation, especially in hypertrophic chondrocytes. In this study, we were able to reveal that the promoter of the type II collagen (COL2A1) gene, encoding the extracellular matrix protein type II collagen, is specifically regulated by AP-2ε. Expression of COL2A1 is downregulated at the transition of chondroblasts into hypertrophic chondrocytes and our data provide evidence that AP-2ε is involved in this process. In reporter gene assays, overexpression of AP-2ε in cartilaginous cell lines resulted in a significant reduction in COL2A1 core promoter activity of ~ 45%. Furthermore, we found that this process is dose-dependent and highly specific for the epsilon isoform. Computational analysis offered only a single putative AP-2-binding motif within the chosen promoter fragment but site-directed mutagenesis revealed this motif to be regulatory inactive. After expanding our screening to motifs containing minor differences from the classical AP-2 consensus sequence (5'-GCCN3 GGC-3'), we determined the sequence 5'-GCCCAGGC-3' ranging from position -128 to -135 bp as an important regulatory site and responsible for COL2A1 downregulation through AP-2ε. Interaction of AP-2ε with the COL2A1 promoter at this site was confirmed by chromatin immunoprecipitation and electromobility shift assay. Further, our experiments suggest that at least one additional factor is involved in this process. This is the first study to prove regulation of COL2A1 by AP-2ε highlighting the role of the transcription factor as a modulator of cartilage development.