|Home | About | Journals | Submit | Contact Us | Français|
Translocations are chromosomal rearrangements that are frequently associated with a variety of disease states and developmental disorders. We identified 2 families with brachydactyly type E (BDE) resulting from different translocations affecting chromosome 12p. Both translocations caused downregulation of the parathyroid hormone-like hormone (PTHLH) gene by disrupting the cis-regulatory landscape. Using chromosome conformation capturing, we identified a regulator on chromosome 12q that interacts in cis with PTHLH over a 24.4-megabase distance and in trans with the sex-determining region Y-box 9 (SOX9) gene on chromosome 17q. The element also harbored a long noncoding RNA (lncRNA). Silencing of the lncRNA, PTHLH, or SOX9 revealed a feedback mechanism involving an expression-dependent network in humans. In the BDE patients, the human lncRNA was upregulated by the disrupted chromosomal association. Moreover, the lncRNA occupancy at the PTHLH locus was reduced. Our results document what we believe to be a novel in cis– and in trans–acting DNA and lncRNA regulatory feedback element that is reciprocally regulated by coding genes. Furthermore, our findings provide a systematic and combinatorial view of how enhancers encoding lncRNAs may affect gene expression in normal development.
Translocations are chromosomal rearrangements that can be inherited and are often associated with cancers and congenital disorders. One of the great difficulties in determining mechanistic causes is the fact that translocations can affect gene expression and regulatory elements and alter chromosomal architecture by changing the relative genomic position of genes and their regulators. Thus, structural genomic changes may effect the regulatory landscape. Such perturbations may be involved in the well-known observations of position effects (1, 2). We combined several methods to resolve some of the complex translocation-related position effects on gene regulation, based on 2 families with brachydactyly type E (BDE). BDE is characterized by a shortening of metacarpals and metatarsals (3). The morphogenesis gene PTHLH, encoding parathyroid hormone–like hormone on chromosome 12p11.2, has been linked to BDE with short stature (4, 5). PTHLH is involved in the regulation of cartilage differentiation and digit condensation (6). During digit condensation, PTHLH and other factors act in concert to balance an intricate signaling network, ultimately impinging on BMP-mediated and sex-determining region Y-box 9–directed (SOX9-directed) chondrogenesis (7). Dysregulated expressions of the haploinsufficient genes PTHLH or SOX9 can cause skeletal phenotypes (2, 4, 5, 8–10). In previous work with a BDE family with balanced translocation t(8;12)(q13;p11.2), we discovered a breakpoint-dependent PTHLH downregulation that elucidated new mechanisms resulting in a Mendelian phenotype (4). The cis regulation of gene expression is accompanied by distinct tissue- and cell type–specific histone modifications that enable or disable the binding of transcription factors to evolutionarily highly conserved DNA elements (11, 12). cis-regulatory elements (CREs) can be intra- or intergenic (13, 14). They are typically located up to 1 megabase (Mb) from their coregulated genes (8, 15–17), but may be farther away. Enhancers may be transcribed into lncRNAs that increase gene expression in cis (18). Furthermore, lncRNAs can themselves act as enhancers; their expression is then correlated with the expression of associated genes (19). The lncRNA transcripts have more than 200 nt, display high interspecies conservation, and are multiexonic (20). Finally, lncRNAs are essential for the development of complex organisms (21). The most accepted model for CRE communication involves direct interaction of CRE with a gene promoter through spatial proximity in chromatin loops (22).
We here investigated another family with autosomal-dominant BDE and a balanced t(4;12)(q13.2-13.3;p11.2) translocation. The translocation-mediated genomic disruption of a CRE caused PTHLH repression and led to BDE in both translocation families. We studied the native PTHLH CREs and determined how they were altered. We discovered what we believe to be a novel CRE that acts both in cis and in trans. Within this CRE exists a lncRNA that was pivotal to its function. lncRNAs are known to be important to developmental programs and tumor suppressor gene networks (23). We observed reciprocal gene and lncRNA regulation that provide a systematic and combinatorial view of how enhancers encoding lncRNAs may affect gene expression in normal development. Dysfunction of the CRE and its lncRNA caused a Mendelian disease.
We analyzed a pedigree demonstrating autosomal-dominant inheritance. The roentgenograms were diagnostic of BDE (ref. 3, Figure Figure1A,1A, and Supplemental Figure 1, A and B; supplemental material available online with this article; doi: 10.1172/JCI65508DS1). Affected individuals exhibited short stature (less than third percentile). A relative shortening of the extremities in comparison to the trunk was evident. The arm span/height ratio was reduced (Supplemental Table 1). High-resolution karyotyping identified the balanced translocation t(4;12)(q13.2-13.3;p11.2) (Figure (Figure1A).1A). BDE-causative HOXD13 and PTHLH mutations were excluded by sequencing (5, 24). Breakpoint allele-specific sequencing identified the derivative chromosome 4 [der(4)] breakpoint 145,460 bp downstream of PTHLH at position 28,147,744 bp (UCSC assembly hg18). The breakpoint was 45,433 bp downstream of the previously described der(8) breakpoint in t(8;12)(q13;p11.2) (4). Microdeletions and 7SL pseudogene insertions at the breakpoints did not eliminate or create transcription factor binding sites (TFBSs; Supplemental Figure 1, E and F). To focus on the chondrogenic function of PTHLH, we induced chondrogenic differentiation in fibroblasts of an affected patient and 2 nonaffected controls. Using quantitative RT-PCR, we observed downregulation of PTHLH (P = 0.04) and its downstream targets, ADAMTS-7 and ADAMTS-12 (Supplemental Figure 1, G and H). PTHLH was downregulated, and its insufficiency correlated with BDE in affected family members.
Chromosome conformation capture (3C) techniques can identify gene-associated cis or trans regulators (25, 26). The 3C method was extended to include ChIP, circularization, and cloning (6C technique). We purposely used a non–genome-wide approach to investigate the cis-regulatory landscape of the candidate gene PTHLH. The ChIP was done on chromatin of the chondrocyte cell line C28/I2 with anti-RNA polymerase II (recognized phosphorylated S5) and anti-H3K4me1 antibodies (27). We searched for H3K4me1-enriched CREs that were linked by chromatin loops to the active PTHLH or SOX9 promoters, our bait sequences. SOX9, located on chromosome 17 and expressed in C28/I2, is the major transcription factor commanding chondrogenesis in the limb, known for long-range cis regulation (9, 28), and served as positive control. The intramolecular ultra-high dilution ligation resulted in circularized chromatin (29). Using inverse-orientated oligos within the baits, we PCR amplified any CREs in the circular chromatin, subcloned the nested PCR amplicons, and sequenced the plasmids (Figure (Figure11B).
We restrictively classified the reads: the bait sequence had to be adjacent to a CRE. These specific sequences of approximately 50 nt of either the promoter or the promoter plus CRE were interspersed by several nt due to the blunt-end repair of chromatin ends (Figure (Figure1C).1C). Randomly distributed positions at all chromosomes without the baits were classified as nonspecific. The 2% interaction frequency of SOX9 reads identified 4 CREs on chromosome 17, whereas 16.2% showed 9 CREs for PTHLH on chromosome 12 (Figure (Figure1D1D and Figure Figure2,2, A and B). The nonspecific sequences of 92.6% in SOX9 6C compared with 47.9% in PTHLH 6C could be explained by the lower SOX9 expression observed (Figure (Figure1E).1E). Surprisingly, we detected an interaction frequency of 1.4% in SOX9 6C for an in trans element on chromosome 12 (Figure (Figure1D).1D). The same element was also detected in PTHLH 6C (Figure (Figure11C).
We then used ChIP to quantify the H3K4me1 and H3K4me2 enrichment of the 6C-identified CREs by quantitative PCR (qPCR; Supplemental Figure 2 and ref. 30). To check for regulatory features of the putative 6C-identified CREs, we analyzed the interspecies PhastCons conservation, the ENCODE H3K4me1 enrichment, and the ESPERR regulatory potential in UCSC. We performed comparative analysis of TFBSs on homologous DNA sequences from different species (rVISTA 2.0; zPicture) and identified highly conserved TFBSs (Supplemental Table 2 and Supplemental Figures 3 and 4). SOX5, SOX9, PBX1, AP1, and HOXA motifs were among those detected. All are known to regulate limb bud development and chondrogenesis (reviewed in ref. 31). Based on significant H3K4me1 enrichment, the most abundantly found 6C interactions, and the in silico analysis, we selected 3 CREs for further experiments: re70373 for SOX9 and re18527 and re52431 for PTHLH (Figure (Figure2,2, A and B). re52431 on chromosome 12q interacted with SOX9 in trans. Chromosomal aberrations (dup and del) in 3 patients with short stature and shortened phalanges in the DECIPHER database of chromosomal imbalances and phenotypes (patient nos. 248901, 248785, and 250426) included the re52431 region, which indicates that the region is important during chondrogenesis.
New 6C libraries validated the reproducibility of the 6C technique. To determine bias, a human lymphoblastoid cell line (LCL), marginally expressing SOX9 and PTHLH, was the control (Figure (Figure1E).1E). In a direct PCR approach, oligos were aligned in each of the 3 selected CREs and within the PTHLH or SOX9 promoter. We reproduced the chondrogenic in cis and in trans interactions (Figure (Figure2C).2C). In the LCL, re18527 did not interact, and re52431 showed a significantly lesser interaction by χ2 test, with PTHLH or SOX9. re70373 interacted with SOX9 in both cell lines, 1.64 Mb more downstream than a previous reported CRE for SOX9 (9). The higher interaction frequencies in C28/I2 compared with those in LCL argued for the detection of chondrogenic-specific CREs. Due to the cis- and trans-chromosomal interactions of re52431, we termed this regulator CISTR-ACT.
The distances between the 6C-identified CREs and the associated gene promoter were far greater than previously reported (17). re70373 was 2.74 Mb downstream of SOX9. re18527 was 9.48 Mb upstream, and CISTR-ACT 24.43 Mb downstream, of PTHLH (Figure (Figure2D).2D). We constructed luciferase reporters with the specific SOX9 and PTHLH promoter sequences downstream of re70373, re18527, and CISTR-ACT to investigate their enhancing potential. Since the latter were found combined with PTHLH in 6C, we constructed composite plasmids: CISTR-ACTF represented all identified TFBSs, whereas CISTR-ACTS was from the most conserved part (Supplemental Figure 4). Significant enhancement of luciferase reporter activity for all CRE-PTHLH and CRE-SOX9 constructs was observed. re70373 combined with CISTR-ACTS and CISTR-ACTF increased SOX9 promoter activity 2.5- and 3.8-fold, respectively (Figure (Figure2E).2E). CISTR-ACTS and CISTR-ACTF activated the PTHLH promoter less effectively than they did SOX9. The combination of re18527 with either CISTR-ACTS or CISTR-ACTF showed an additive activation effect on the PTHLH promoter (Figure (Figure22F).
To address spatial chromatin loops that might form between the selected CREs and the promoters as a prerequisite of gene regulation, we performed colocalization-FISH on C28/I2 nuclei (8). We captured 80 nuclei in each experiment and analyzed signal distances less than 1 μm. Bona fide colocalization signals merged below 0.27 μm (Mann-Whitney rank-sum test; Figure Figure3A).3A). For negative controls, we chose intervening probes between the gene and the CREs and equidistant loci from the CREs in the other direction on chromosome 12 or 17 (Supplemental Figure 5).
The measured distances between re18527 and PTHLH were significantly shorter than the negative controls (P = 9.31 × 10–3; Figure Figure3B3B and Supplemental Figure 5A). For investigation of the cis distance of 24.43 Mb between PTHLH and CISTR-ACT, probe pairs with accordant genomic distances on chromosome 15 around the protein tyrosine phosphatase–like A domain containing 1 (PTPLAD1) locus were used as further controls. PTPLAD1 was not expressed in C28/I2. Compared with all controls, CISTR-ACT was in spatial proximity to PTHLH (Figure (Figure3C3C and Supplemental Figure 5A). The colocalization of the re70373 element and SOX9 was barely significant; the intervening probe detected closer proximity (P = 0.028; Figure Figure3D3D and Supplemental Figure 5B). The SOX9–CISTR-ACT in trans contact was compared with SOX9–chromosome 12q and SOX9–chromosome 15q probe pairs (Figure (Figure3E).3E). Moreover, 2 probe pairs of SOX9 with positions on the exterior nucleus–located chromosome 1 showed no spatial distance less than 1 μm. In all instances, CISTR-ACT had significant contact with SOX9. Putative chromosome 12 with 17 translocations were excluded in metaphase FISH (Supplemental Figure 5C). The statistical analysis showed undetermined signal distances, but specific cis or trans interactions of CISTR-ACT with PTHLH or SOX9 happened in a fraction of cells. The frequency distributions of colocalized or very close signals were significantly higher in the experimental measurements than in the negative controls (Table (Table1).1). Overall, the in cis genomic distances between promoter and CRE were inversely proportional to the interaction frequencies in cells. Consistent with the 6C results, the FISH data excluded the possibility that incidental intrachromosomal folding or interchromosomal territorial contacts occurred.
Interestingly, we found an expressed sequence tag, DA125942, with 2 exons overlapping with CISTR-ACT on 12q in the UCSC Genome Browser. DA125942 represented a nonpolyadenylated lncRNA of 581 bp (Supplemental Figure 4). Rapid amplification of cDNA ends (RACE) experiments confirmed the full-length transcript in C28/I2 cells. To determine whether DA125942 transcription was independent of an external promoter and regulator, we constructed a vector lacking regulatory sequences (pGL-luci) and added re18527, CISTR-ACTS, or CISTR-ACTF alone or in combination. Using quantitative RT-PCR, we observed highly upregulated DA125942 in C28/I2 cells transfected with CISTR-ACTF (Figure (Figure4A).4A). This result was suggestive of 2 concepts: first, CISTR-ACTF harbors an internal promoter driving the transcription; and second, mature DA125942 results from splicing. Kim and colleagues reported similar observations (18). In disrupted CISTR-ACTS plasmids, DA125942 was not expressed, and re18527 negatively influenced DA125942 expression (Figure (Figure44A).
lncRNAs can enhance and specifically influence gene regulation in cis over several Mb (19). We proposed that DA125942 could conduct in cis regulation of PTHLH expression and simultaneously could also be involved in the in trans regulation of SOX9. To test this issue, we used siRNA knockdown in C28/I2 cells to deplete DA125942, PTHLH, and SOX9. Depletion of DA125942 resulted in downregulation of PTHLH and SOX9. Surprisingly, knocking down PTHLH and SOX9 significantly increased DA125942 expression by 4- and 6.8-fold, respectively, suggestive of reverse action of these gene products on lncRNA transcription. Silencing of SOX9 led to significantly reduced PTHLH expression, and vice versa (Figure (Figure4B).4B). These results implied an interactive influence of DA125942, PTHLH, and SOX9 on transcript levels. Overexpression of PTHLH or SOX9 led to upregulation of DA125942, SOX9, or PTHLH. DA125942 overexpression repressed the expression of the coding genes (Figure (Figure4C).4C). These data confirmed the expression dependency of the network.
The HOXC cluster is localized 187 kb downstream of DA125942, and several lncRNAs are encoded by, and interact with, HOX loci (32, 33). We tested whether DA125942 regulated HOXC or other genes either in cis or in trans by 3′ expression profiling. We overexpressed and depleted DA125942 in C28/I2 cells. The combined ANOVA resulted in 1,842 differentially expressed probes (5% false discovery rate [FDR]; Supplemental Table 4). To distinguish up- and downregulation patterns, k-means clustering was performed (Supplemental Figure 6A). Chromatin-remodeling genes, i.e., protein arginine methyltransferase 1 (PRMT1), high mobility group nucleosomal binding domain 2 (HMGN2), polycomb group ring-finger 6 (PCGF6), and enhancer of zeste homolog 2 (EZH2), and several transcription factors were differentially expressed in clusters 1, 3, 4, and 5 (Supplemental Figure 6B). In cluster 2, 33 genes had annotations linking them to organ morphogenesis (GO:0009887; P = 7.51 × 10–7); 11 genes as singletons were excluded in STRING network (Figure (Figure4D).4D). In addition to PTHLH and SOX9, we found essential morphogenesis genes for early limb patterning and chondrogenesis. Hypoxia inducible factor 1 (HIF1A) regulates SOX9 (34). Noggin (NOG) and follistatin (FST) regulate outgrowth and chondrogenesis in the chicken limb bud (35, 36). ETS variant 4 (ETV4) is involved in limb outgrowth and in the anterior-posterior patterning of the autopode (37). The HOXB genes HOXB3, HOXB5–HOXB7, and HOXB8 (singleton; data not shown) and the HOXC genes code for transcription factors involved in diverse morphogenetic processes, including limb formation (38, 39). None of the HOXC genes was differentially expressed. Osteoprotegerin (TNFRSF11B) and the vitamin D receptor (VDR) maintain mineral bone density in osteoblasts (40). The HOXB cluster, ETV4, NOG, and SOX9 were all located on chromosome 17, indicative of several DA125942-mediated gene regulations.
We investigated whether the 2 independent translocations on 12p influenced, next to the PTHLH downregulation, the distant DA125942 and SOX9 expressions in BDE patients. Significant DA125942 upregulation in chondrogenically differentiated fibroblasts of BDE patients with t(8;12) or t(4;12) was found (Figure (Figure4,4, E and F). SOX9 was marginally, but not significantly, downregulated (P = 0.059; Supplemental Figure 7A). Next, we addressed whether DA125942 occupied the PTHLH and SOX9 loci and whether the interaction was altered in the translocations t(8;12) and t(4;12). We used the recently described chromatin isolation by RNA purification (ChIRP) method with split probe sets 1 and 2, tiling DA125942 (~80% tiled) on chromatin of C28/I2 and chondrogenically differentiated BDE fibroblasts (41). For qPCR, we eluted RNA to measure lncRNA enrichment and DNA to detect DA125942 occupancy at the PTHLH or SOX9 loci; approximately 12 kb of each gene was analyzed. A probe set against lacZ, the GAPDH and CISTR-ACT loci, and a healthy proband served as controls; nonprecipitated RNA or chromatin was the input control. DA125942 was 1- to 5.6-fold enriched compared with the input. As expected, in chromatin from BDE patients, more lncRNA was retrieved; GAPDH and lacZ were not detected (Figure (Figure4G).4G). Since the lncRNA occupancy at genomic loci is not well understood, we analyzed low and highly conserved, intronic, exonic, and promoter sequences especially in regions of ENCODE-annotated TFBSs. We detected DA125942 binding in the eluted C28/I2 DNA at the PTHLH and SOX9 loci (Figure (Figure4,4, H and I). In qPCR of the restricted chromatin from BDE patients, we analyzed the highest and lowest enriched positions that we had found in C28/I2 cells. Compared with the healthy proband, binding of DA125942 at the PTHLH locus was significantly reduced in the BDE patients’ chromatin (Figure (Figure4J).4J). Binding at the SOX9 locus was not significantly altered (Figure (Figure4K);4K); the GAPDH or DA125942 loci were not detected (Supplemental Figure 7, B and C).
The homologous regions to the human CISTR-ACT in mouse and in rat were detected on the basis of the DNA sequence, including 1-exon lncRNA transcripts in both species (mouse, AW491522; rat, BM385392). By phylogenetic and comparative analysis of TFBSs by Mulan (42, 43), we found that 121 multiconserved TFBSs overlapped in humans, mice, and rats. The lncRNA sequences were evolutionarily more divergent (Figure (Figure5,5, A and B). In contrast to the human genome, the lncRNAs AW491522 and BM385392 were not located on the chromosome encoding Pthlh (Figure (Figure5C),5C), raising the possibility of solely in trans–acting regulation. The approved mouse ATDC5 cells and rat chondrocytes (RCSs) were used to study chondrogenesis (44, 45). ATDC5 cells were chondrogenically differentiated for 10 days to induce major chondrogenic markers (Supplemental Figure 8, A and B). We again used siRNA transfections to knock down Pthlh, Sox9, and each lncRNA to investigate whether there is an expression-associated mechanism corresponding with that found in humans. AW491522 or BM385392 knockdown led to Pthlh and Sox9 downregulation. In contrast to the observations in human C28/I2 cells, Pthlh or Sox9 depletion caused lncRNA downregulation in ATDC5 cells and RCSs (Figure (Figure5,5, D and E), indicative of an evolutionarily more complex gene regulation in humans. By ChIRP experiments with 2 probe sets targeting the homologous regions of the rodent lncRNAs (~80% tiled), we retrieved approximately 60% (Figure (Figure6A).6A). In qPCRs on eluted DNA, homologous regions at the rodent Pthlh or Sox9 loci were occupied by either AW491522 or BM385392 (Figure (Figure6,6, B–E). The regions corresponded to the human regions occupied by DA125942; enrichment of the lncRNA loci and Gapdh was not detected (Supplemental Figure 8, C–F). Further investigation is required to determine whether the exposed CISTR-ACT locations on the outer chromosomal q arms support the trans regulations. Overall, these data indicate that the lncRNAs act in a contact-based manner and that the network is conserved among the mammals studied herein.
To investigate whether the human and mouse CISTR-ACTs act as classical enhancers in the developing limb buds, we generated transgenic embryos using pronuclei microinjection. The CISTR-ACTs were cloned in front of a lacZ reporter. All stages analyzed showed consistent lacZ expression patterns during limb development. Human transgene-driven lacZ expression was detected in the distal limb bud mesenchyme at E10.5 and E11.5, with a recognizable shift toward the anterior part (Figure (Figure7A).7A). At E12.5, lacZ was expressed overlapping with condensations of digits 1, 2 and 4 and in the first and fourth interdigital mesenchyme (Figure (Figure7A).7A). At E13.5, the strongest expression was maintained overlapping with the condensations of digits 1 and 2 and the first interdigit. Weaker staining was visible overlapping the condensations of digit 3 and the distal tip of digit 4 as well as in interdigits 2, 3, and 4. Thus, lacZ expression overlapped with domains known to be involved in skeletal patterning, as well as those undergoing and/or regulating chondrogenesis (46). Mouse transgene-driven lacZ expression analyzed at E10.5 and E12.5 showed a strong overlap with the staining of the human counterparts in both cases. At E12.5, staining from the mouse transgene showed restriction to the first interdigit and regions overlapping condensations of digits 1 and 2. The domains, in part, correlated with previous documentation of the early limb patterning and chondrogenesis genes identified in our expression profiling, including Hoxb8, Etv4, Pthlh, and Sox9 (37, 47–49).
These results showed that the translocation-mediated disruption of the PTHLH–CISTR-ACT chromosomal association affected the chromosomes’ feedback control and caused PTHLH and DA125942 dysregulation. The findings in the BDE patients were analogous to those from the in vitro and in vivo studies described herein and expand our knowledge about remote-controlled gene and lncRNA regulation in altered chromatin architecture (Figure (Figure7,7, B and C).
Our results elucidate how the translocation-mediated disruption of the PTHLH–CISTR-ACT association caused disease and shed light on the complexity of a molecular pathogenic mechanism that we believe to be novel. These basic chromatin-related epigenetic findings could be of great translational value, especially in other diseases with chromosomal rearrangements. We demonstrated consistent results in patient material and in both in vitro and in vivo models. We documented that chromatin loops involving conserved DNA elements and lncRNA existed for intra- and inter-chromosomal communication, resulting in the regulation of gene and lncRNA expression. We elucidated CISTR-ACT to be a chondrogenic regulator, encoding lncRNA.
The human CISTR-ACT pinpointed SOX9 in trans on chromosome 17. We observed spatial nuclear proximity of the cis elements with distances of 9.48 and 24.43 Mb to PTHLH with colocalization-FISH. Proximal CRE promoter signals occurred. These temporary interactions happened in a limited number of cells, which indicates that adequate target gene expression and maintenance of tissue specificity were effectively fulfilled. In studies of a human ectopic β-globin locus control region (LCR), similar observations were made. Rare ultra-long range regulations and endogenous trans β-globin activation happened only in “jackpot” cells with enriched interchromosomal interaction (50). Amano and colleagues studied murine limbs and observed that sonic hedgehog (Shh) cis regulation was only active in some fractions within the area of Shh expression, but not in others. These findings confirm that chromatin looping is also tissue and developmental stage specific (51). Based on our observations using FISH and 6C, we suggest that tissue-specific fine regulation mediated by chromatin loops occurs infrequently and unpredictably.
We characterized the function of the CISTR-ACT–encoded lncRNAs in siRNA and overexpression experiments. We recapitulated the expression-associated network in 3 species. lncRNA silencing led to PTHLH or SOX9 repressions in humans and rodents, whereas their altered expression coincided with lncRNA dysregulation. Some lncRNAs have enhancer-like functions and are associated with gene expression (19). In contrast to this prior report, we here present evidence for a regulatory feedback expression network. Coding genes could regulate lncRNAs, and vice versa (Figure (Figure7B).7B). lncRNAs may be necessary to establish transcription by recruiting secondary mediators, as previously shown for the epigenetic state regulation of rRNA genes (52). In ChIRP experiments, we documented that the lncRNAs bound at the PTHLH and SOX9 loci. This finding indicates that the lncRNAs act directly on the coding genes’ loci. The translocation-altered chromatin architecture seemed to suppress the looping of CISTR-ACT necessary to find and adequately regulate PTHLH by its lncRNA or vice versa. PTHLH has a nonclassical nuclear localization signal and can bind RNA in the nucleus (53, 54).
Wang and colleagues observed that the long intergenic ncRNA (lincRNA) HOTTIP activates several HOXA genes by proximity and bound protein mediators, thereby maintaining active chromatin. Ectopic HOTTIP did not regulate HOXA genes in their study and generated a dominant-negative effect (33). Supported by their findings, we believe that regular chromosomal contacts bring the CISTR-ACT DNA locus in proximity of its target genes to bring down its lncRNA as an intermediate for targeting the regulated coding genes. Hence, the translocation-mediated altered chromosomal architecture disrupted the chromatin looping and the ability of CISTR-ACT’s lncRNA to find PTHLH. Moreover, DA125942 overexpression led to PTHLH and SOX9 repression. The lncRNA’s overexpression may lead to the binding of proteins that can generate a dominant-negative effect, as the repressed PTHLH leading to autosomal-dominant BDE.
lincRNAs interact with chromatin-remodeling complexes in a well-documented manner for HOTAIR at the HOXC gene cluster and XIST for the X chromosome inactivation (32, 55, 56). Chromatin-remodeling genes and ncRNAs determine tissue-specific epigenetic states and gene expression (57). In studies of HOTAIR or XIST ncRNA, others have proposed that lncRNAs mediate genetic traffic to establish specific epigenetic marks, resulting in guided gene expression (58). The differentially expressed genes suggest that DA125942 has several functions. The preparation of transcriptionally active euchromatin through the lncRNA might go hand in hand with specific morphogene regulation, mediated by the lncRNA and by DNA- or RNA-bound factors. The underlying highly conserved CISTR-ACT DNA sequences, with the predicted chondrogenic TFBSs, showed tissue-specific regulation in the transgenic embryos. The CISTR-ACT–driven lacZ expression domains partially overlapped with the expression patterns of Hoxb8, Etv4, Pthlh, and Sox9. These morphogenesis genes were also differentially expressed in expression profiling. Moreover, the domains anatomically covered the most apparent human BDE regions.
Translocation-mediated chromosome 12 disruption led to the dysregulation of PTHLH and DA125942. The zinc-finger protein CTCF is the best-described regulator mediating either intra- or interchromosomal genome organization (59, 60). Aberrant genomic architecture results in disorganization and failed recruitment of secondary mediators to genomic locations. Thus, the tissue-specific feedback expression of PTHLH and DA125942 does not occur. We suggest that SOX9 and PTHLH share common sites of transcription and belong to a chondrogenic transcription factory (61, 62). As far as we know from studying chromosomal aberration syndromes and intact chromosomes, certain chromosomal territories must be maintained to ensure proper gene regulation (63, 64). In a 3D map of human fibroblasts, the chromosome 12 and 17 territories are adjacent, although unambiguous deterministic chromosomal neighborhoods are not observed (65). We hypothesize that the trans-mediated regulation of SOX9 can occur because the major der(12)s carrying CISTR-ACT remain in their chromosomal territory. CISTR-ACT is still able to communicate with SOX9 on 17q, although the PTHLH and lncRNA dysregulations marginally affect SOX9 expression. PTHLH on 12p appears to be translocated to the chromosomal territories of chromosomes 4 and 8, respectively. Access to the CISTR-ACT is prevented by the translocation (Figure (Figure77C).
Taken together, these observations broaden our awareness of regulatory distances across the centromere, further elucidate the incidence of trans-regulatory DNA-DNA and RNA-lncRNA communication, and suggest that molecular dysfunctions that result in clinical phenotypes are highly complex. We provided genetic evidence that chromosomal communicators work synergistically at both the DNA and RNA levels in a combinatorial and responsive fashion. The reciprocal fine regulation of gene and lncRNA expression is a systems biological process, accounting for effective development or tissue perpetuation. We believe that nucleus-based networks probably serve to maintain temporary, consecutive, and tissue-specific expression of certain gene products. The recruitment mechanism supporting the locus-specific contact must be investigated further. We argue that biophysical and biochemical processes work together stochastically and dynamically. Chromosomes own their territories; although mitotic and meiotic conformational changes occur, they are not randomly organized. Nuclear organization is evolutionarily conserved and directed (64, 66). Brownian motion, directed biochemical mechanisms, and effects that can be described using chaos theory could regulate the diverse interactions (67). We propose that determination meets variability.
After obtaining written informed consent, we studied 2 affected (III:1 and III:4) and 2 nonaffected (III:3 and IV:1) individuals. Skin fibroblasts from 1 affected individual (III:4) were obtained.
Fibroblasts were cultured in MEM 199, C28/I2 and ATDC5 cells were cultured in DMEM/F12 at a 1:1 ratio (68), and RCSs were cultured in DMEM; all were supplemented with 10% FCS, 100 U/ml penicillin, and 100 μg/ml streptomycin. The chondrogenic differentiation of fibroblasts was done as previously reported (4); ITS-X (Life Technologies) induced ATDC5 chondrogenic differentiation.
Metaphase FISH analysis was performed using standard protocols. BACs and fosmids were ordered at RZPD ( www.imagenes-bio.de). 2D colocalization-FISH was performed as previously reported (69). The Olympus-BX61 microscope with U-TV1X2 and PlanApo 60×/1.4 oil lens in combination with the ASI CCD-camera 1300QDS and FISHView software (version 220.127.116.11) were used for signal detection and measurements. Measurements were made independently by 2 different investigators.
See Supplemental Methods for the detailed 6C protocol. ChIRP was done according to a previously published protocol (41); see Supplemental Table 5 for probes.
5 × 104 cells were seeded 12–24 hours prior to transfection. Equal molarities in transiently transfected C28/I2 cells compensated for the diverse plasmid sizes. For reporter assays, 75 fmol each of the CRE promoter luciferase constructs were transfected with Fugene HD according to the manufacturer’s recommendations (Roche). The SOX9 promoter sequence was on chromosome 17: 67,626,313–67,629,147 bp (UCSC assembly hg18). The vector backbone was pGL2-basic (Promega); 12.5 ng pRL-TK (Promega) was added to control transfection efficiency. Lysates were analyzed using the Dual-Glo Assay (Promega). 100 nM siRNA were transfected with Dharmafect (Dharmacon) for 48 hours in C28/I2 and ATDC5 cells, 120 nM siRNA were used for RCS transfections. Scrambled, GAPDH, SOX9, and PTHLH siRNAs were validated (Ambion). Different custom-designed siRNAs for the lncRNAs were obtained. For visual simplicity, results using only 1 siRNA are shown.
DA125942 siRNAs were as follows: siRNA-1 sense, GCACAUGACCACAUGGAAATT; siRNA-1 antisense, UUUCCAUGUGGUCAUGUGCGA; siRNA-2 sense, GUUUCUCUCAGCAAAUUCATT; siRNA-2 antisense, UGAAUUUGCUGAGAGAAACTG; siRNA-3 sense, CAAGGACACUGAAAAAGCUTT; siRNA-3 antisense, AGCUUUUUCAGUGUCCUUGAG. AW491522 siRNAs were as follows: siRNA-1 sense, AGGCCCUGAUGCUAGUGAATT; siRNA-1 antisense, UUCACUAGCAUCAGGGCCUTG; siRNA-2 sense, GCACGCAUCUCUUCCACCATT; siRNA-2 antisense, UGGUAAGAGAUGCGUGCGG. BM385392 siRNAs were as follows: siRNA-1 sense, GAGAUCUGGUCAGAGCAAATT; siRNA-1 antisense, UUUGCUCUGACCAGAUCUCTA; siRNA-2 sense, CAAUACUUCUGGACUCUCUTT; siRNA-2 antisense, AGAGAGUCCAGAAGUAUUGCA. Negative control was scrambled siRNA. SOX9, DA125942, and PTHLH were significantly downregulated (>75%). The siRNA for GAPDH was a reference for transfection efficiency and downregulation. Overexpression transfections were done as previously described by Fugene X-tremeGENE HP (Roche). PTHLH (SC307833) and SOX9 (SC321884) expression plasmids were obtained from OriGene.
Genomic sequences were analyzed by NCBI BLAST, UCSC BLAT, repeat masking at GIRI Institute, and in the UCSC Genome Browser (March 2006; NCBI36/hg18) (70, 71). rVISTA version 2.0 and Transfac Professional version 10.2 were used to determine mammalian conserved TFBSs (72, 73). Functional experiments were done at least 3 times. Number of experiments (n) are shown in figure legends. Significance was determined by Student’s t test. Signal detections in colocalization-FISH were statistically analyzed by Mann-Whitney rank-sum test; data represent mean ± SD. Expression profiling was done on ht12 arrays (Illumina) in biological triplicates, followed by ANOVA, multiple testing (FDR 5%), cluster, GProfiler enrichment, and STRING analysis. The microarray dataset was deposited in GEO (accession no. GSE39896; http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE39896).
On the BeadStudio 18.104.22.16806 Platform array, data were quantile normalized at the probe level without background correction, and consistently low-expressed transcripts (minimum detection P > 0.05) were rejected. Data were log2 transformed after an offset addition (16) and analyzed on significant expression differences among the 4 groups (siRNA-depleted and overexpressed DA125942 C28/I2 RNA, controls with empty expression vector, and scrambled siRNA) in triplicate using ANOVA followed by a FDR (74) multiple testing correction. Probes that underwent 5% FDR were k means clustered with average Euclidean distance function and k = 5, according to a local minimum in the Davies-Bouldin k estimation (75). Clusters 1, 2, and 4 were selected according to those outstanding expression statuses in DA125942 silencing and further investigated by functional enrichment and interaction. Enrichment analysis was performed using G:Profiler (76), with a simulation-based analytical threshold for significance estimation. Physical and functional interactions between proteins were determined using the String platform (77) and a medium confidence score of 0.4. Resulting singletons were rejected before network visualization.
Human (chromosome 12; 52,430,503–52,431,992 bp; assembly hg18) and mouse (chromosome 15; 102,576,412–102,577,508 bp; assembly mm9) CISTR-ACT were cloned in the Sfi-Hsp68LacZ reporter construct (78). Gel-extracted linearized DNA was injected in fertilized mouse oocytes according to standard procedures. At each time point, 25–42 embryos were analyzed by whole-mount lacZ staining.
Data are presented as mean ± SD unless otherwise indicated. Differences between experimental groups were analyzed by 2-tailed Student’s t test unless otherwise specified. A P value less than 0.05 was considered statistically significant.
The local ethics committee (Charité-Universitätsmedizin Berlin, Ethikausschuss 1; EA3/009/07) approved the study. All subjects provided written informed consent. The local regulatory animal authority (Landesamt für Gesundheit und Soziales Berlin — LAGeSo; 0143/07) approved the pronuclear injections and oviduct transfers in mice.
We thank all family members for their cooperation. We thank I. Hollfinger, E. Bartels-Klein, Y. Wefeld-Neuenfeld, A. Mühl, A. Käbner-Frensel, and E.M. Gerlach for technical assistance. The Deutsche Forschungsgemeinschaft (DFG; grants BA1773/5-1, MA5028/1-2, and RU1508/1-1) supported the study. P.G. Maass was supported by a Max-Delbrück-Center and a DFG scholarship. This study made use of data generated by the DECIPHER Consortium. A full list of centers that contributed to the generation of the data is available from http://decipher.sanger.ac.uk and via email from ku.ca.regnas@rehpiced. The Wellcome Trust provided funding.
Conflict of interest: The authors have declared that no conflict of interest exists.
Citation for this article: J Clin Invest. 2012;122(11):3990–4002. doi:10.1172/JCI65508.
See the related Commentary beginning on page 3837.