Deletion of CTCF sites in the SHH locus alters enhancer-promoter interactions and leads to acheiropodia - PubMed (original) (raw)

. 2021 Apr 16;12(1):2282.

doi: 10.1038/s41467-021-22470-z.

Yichi Zhang 1 2 3, Chenling Xiong 1 2, Jingjing Zhao 1 2, Ilias Georgakopoulos-Soares 1 2, Lauren Kane 4, Kirsty Jamieson 2, Michael J Bamshad 5 6 7, Deborah A Nickerson 6 7; University of Washington Center for Mendelian Genomics; Yin Shen 2 8, Laura A Lettice 4, Elizabeth Lemos Silveira-Lucas 9, Florence Petit 10, Nadav Ahituv 11 12

Affiliations

Deletion of CTCF sites in the SHH locus alters enhancer-promoter interactions and leads to acheiropodia

Aki Ushiki et al. Nat Commun. 2021.

Abstract

Acheiropodia, congenital limb truncation, is associated with homozygous deletions in the LMBR1 gene around ZRS, an enhancer regulating SHH during limb development. How these deletions lead to this phenotype is unknown. Using whole-genome sequencing, we fine-mapped the acheiropodia-associated region to 12 kb and show that it does not function as an enhancer. CTCF and RAD21 ChIP-seq together with 4C-seq and DNA FISH identify three CTCF sites within the acheiropodia-deleted region that mediate the interaction between the ZRS and the SHH promoter. This interaction is substituted with other CTCF sites centromeric to the ZRS in the disease state. Mouse knockouts of the orthologous 12 kb sequence have no apparent abnormalities, showcasing the challenges in modelling CTCF alterations in animal models due to inherent motif differences between species. Our results show that alterations in CTCF motifs can lead to a Mendelian condition due to altered enhancer-promoter interactions.

PubMed Disclaimer

Conflict of interest statement

The authors declare no competing interests.

Figures

Fig. 1

Fig. 1. Fine-mapping of the acheiropodia-associated deletion.

a Pedigree of acheiropodia family with proband indicated via the arrow. Squares and circles represent males and females, respectively. b WGS alignments showing a homozygous 12 kb deletion in the acheiropodia proband. The Y-axis is the read depth (number of reads for each nucleotide). The deletion appears in a heterozygous manner in both parents. BP: breakpoint; P: proband; M: mother; F: father. c PCR amplification using three different primers pairs, whose location is indicated in b, further confirming the breakpoint in the proband (P) and mother (M) and father (F). PCR was performed several times using different primer sets to validate the deletion. d Sanger sequencing of the acheiropodia patient showing the breakpoint sequence which also has a CA insertion.

Fig. 2

Fig. 2. Mouse transgenic enhancer assay for the 12 kb acheiropodia-associated sequence.

Schematic representation of the mouse transgenic enhancer assay (upper panel) showing the 12 kb acheiropodia-associated sequence cloned upstream of the Hsp68 promoter-LacZ gene. Enhancer activity, as visualized by LacZ staining, was not observed in the ZPA for the six PCR positive E11.5 mouse embryos (lower panel).

Fig. 3

Fig. 3. CTCF and RAD21 distribution in the LMBR1-SHH locus.

a CTCF and RAD21 ChIP-seq enrichment in lymphoblastoid cells from wildtype (WT) and proband (Mut) at the LMBR1-SHH locus. GM12878 (lymphoblastoid cell line) TAD boundaries are shown in orange and gray horizontal bar. ZRS and the acheiropodia-associated deleted region are shown in orange and blue vertical lines respectively. CTCF orientations are shown as red triangles. The _Y_-axis is the signal _P_-value to reject the null hypothesis that the signal at that location is present in the control. b Zoom in of the region around the LMBR1 gene. c CTCF motif from JASPAR [

http://jaspar.genereg.net/

] and CTCF motif scores, as assigned by FIMO, overlapping CTCF peaks in the LMBR1 locus.

Fig. 4

Fig. 4. Chromatin interactions with the SHH promoter.

a 4C contact profiles in lymphoblastoid cells from wildtype (WT) and proband (Mut) at the _LMBR1_-SHH locus. The viewpoint is depicted by a black arrowhead. The median and 20th and 80th percentiles of sliding 2–50 kb windows determine the main trend line. The color scale represents enrichment relative to the maximum medium value attainable at 12 kb resolution. CTCF and RAD21 ChIP-seq peaks are shown as black and blue vertical line respectively. The ZRS and the acheiropodia-associated deleted region are shown as orange and blue vertical lines respectively. CTCF orientations are shown as red triangles. b Zoom in of the region around the LMBR1 gene.

Fig. 5

Fig. 5. DNA FISH showing the SHH promoter interaction with the acheiropodia-associated region.

a Schematic of the LMBR1-SHH locus showing the ZRS and the acheiropodia-associated deleted region via orange and blue vertical lines respectively. CTCF orientations are shown as red triangles and the locations to which the DNA FISH probes hybridize to are depicted by light blue bars. b Images of representative nuclei from DNA FISH analysis of parental and proband lymphoblastoid cells showing FISH signals for SHH, LSC2, LSC3-5, and 12 kb probes. Scale bars: 5 µm. c Violin plots showing the distribution of interprobe distances (µm) between SHH – LSC1, SHH – LSC2, and SHH – Deletion. The wild-type allele was distinguished from the mutant allele in the parental cell line using the 12 kb probe. SHH – deletion is measured from SHH to 12 kb probe on the wild-type allele and from SHH to LSC3-5 probe on the mutant allele. The statistical significance between datasets was examined by a two-sided Mann–Whitney U-test, ** = 0.004537 and **** = 2.083×10−11 (n = 75–150 alleles).

Fig. 6

Fig. 6. Human and mouse genomic comparisons and phenotype of mice where the orthologous region was deleted.

a Comparison of the LMBR1-SHH locus between human and mice. CTCF site deletions analyzed by Paliou et al. are marked by purple lines and those generated by #, Williamson et al. are denoted by yellow lines or gray rectangle, and marked by *. CTCF motif orientation is shown via red triangles. The acheiropodia-associated deletion and its mouse orthologous sequence are depicted by a blue rectangle. Mouse limb CTCF ChIP-seq data from ENCODE, Andrey et al., Paliou et al. and human CTCF ChIP-seq data from this study (WT = wild type; Mut=proband) are shown as black genomic tracks below the locus. The conservation track is adopted from the Ensembl Genome Browser with green lines indicating conserved sequences between humans and mice. b Shh gene expression levels dissected from E11.5 mouse autopods from wild type (WT) and knockout (KO) mice as determined by qRT-PCR. Each value represents the ratio of Shh gene expression to that of β-Actin, and values are mean ± standard deviation. The expression value of WT group was arbitrarily set at 1.0. Each dot represents one embryo and statistical differences were determined using a two-sided unpaired t test (P = 0.7796, N.S., not significant). Source data are provided as a Source Data file. c Whole-mount in situ hybridization for Shh of wild type (WT) and knockout (KO) E11.5 mouse embryos. Forelimbs (FL) and hindlimbs (HL) were dissected and shown in the lower panel. d Wild type (WT) and knockout (KO) E18.5 limb skeletal staining using alizarin red/alcian blue.

Fig. 7

Fig. 7. Proposed model for the aberrant chromatin structure of the _LMBR1_-SHH locus in the acheiropodia patient.

Model of chromatin structure in the _LMBR1_-SHH locus based on our ChIP-seq and 4C-seq data. CTCF sites are shown as red triangles and the cohesin complex is shown as a green ring. The ZRS is depicted as an orange oval.

Similar articles

Cited by

References

    1. Freire-Maia A, Freire-Maia N, Morton NE, Azevêdo ES, Quelce-Salgado A. Genetics of acheiropodia (the handless and footless families of Brazil). VI. Formal genetic analysis. Am. J. Hum. Genet. 1975;27:521–527. - PMC - PubMed
    1. Ianakiev P, et al. Acheiropodia is caused by a genomic deletion in C7orf2, the human orthologue of the Lmbr1 gene. Am. J. Hum. Genet. 2001;68:38–45. doi: 10.1086/316955. - DOI - PMC - PubMed
    1. Fagerberg L, et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol. Cell. Proteom. 2014;13:397–406. doi: 10.1074/mcp.M113.035600. - DOI - PMC - PubMed
    1. Williamson, I. et al. Developmentally regulated Shh expression is robust to TAD perturbations. Development146, dev179523 (2019). - PMC - PubMed
    1. Briscoe J, Thérond PP. The mechanisms of Hedgehog signalling and its roles in development and disease. Nat. Rev. Mol. Cell Biol. 2013;14:416–429. doi: 10.1038/nrm3598. - DOI - PubMed

Publication types

MeSH terms

Substances

Supplementary concepts

Grants and funding

LinkOut - more resources