DNA Methylation Status of SHOX-Flanking CpG Islands in Healthy Individuals and Short Stature Patients with Pseudoautosomal Copy Number Variations.

SHOX resides in the short arm pseudoautosomal region (PAR1) of the sex chromosomes and escapes X inactivation. SHOX haploinsufficiency underlies idiopathic short stature (ISS) and Leri-Weill dyschondrosteosis (LWD). A substantial percentage of cases with SHOX haploinsufficiency arise from pseudoautosomal copy number variations (CNVs) involving putative enhancer regions of SHOX. Our previous study using peripheral blood samples showed that some CpG dinucleotides adjacent to SHOX exon 1 were hypomethylated in a healthy woman and methylated in a woman with gross X chromosomal rearrangements. However, it remains unknown whether submicroscopic pseudoautosomal CNVs cause aberrant DNA methylation of SHOX-flanking CpG islands. In this study, we examined the DNA methylation status of SHOX-flanking CpG islands in 50 healthy individuals and 10 ISS/LWD patients with pseudoautosomal CNVs. In silico analysis detected 3 CpG islands within the 20-kb region from the translation start site of SHOX. Pyrosequencing and bisulfite sequencing of genomic DNA samples revealed that these CpG islands were barely methylated in peripheral blood cells and cultured chondrocytes of healthy individuals, as well as in peripheral blood cells of ISS/LWD patients with pseudoautosomal CNVs. These results, in conjunction with our previous findings, indicate that the DNA methylation status of SHOX-flanking CpG islands can be affected by gross X-chromosomal abnormalities, but not by submicroscopic CNVs in PAR1. Such CNVs likely disturb SHOX expression through DNA methylation-independent mechanisms, which need to be determined in future studies.

SHOX far-downstream copy-number variations involving cis-regulatory nucleotide variants in two sisters with Leri-Weill dyschondrosteosis.

SHOX haploinsufficiency leading to Leri-Weill dyschondrosteosis (LWD) and idiopathic short stature typically results from intragenic mutations or copy-number variations (CNVs) involving SHOX and/or its putative enhancer regions that are distributed in the genomic interval between 400 kb and 840 kb from Xpter/Ypter. Here, we report two sisters with LWD, who carried a deletion in the far-downstream region of SHOX. The 0.62 Mb deletion contained 50 single nucleotide polymorphisms (SNPs) and short insertions and deletions (indels), whose genotypes were linked to SHOX expression levels in the Genotype-Tissue Expression portal. Notably, most of these SNPs/indels accumulated within a ~20 kb interval that was positioned ~900 kb away from Xpter/Ypter. These SNPs/indels showed similar minor allele frequencies, indicating that they reside within a haplotype block. The ~20 kb interval was not evolutionarily conserved; however, it was associated with the previously determined peak of chromosome conformation capture profiling (4C)-seq. Importantly, the deletion in the present cases partially overlapped with CNVs of three previous cases with skeletal deformity and/or short stature. The results indicate that far-downstream CNVs constitute rare genetic causes of SHOX haploinsufficiency. These CNVs possibly impair SHOX expression through copy-number changes of a human-specific cis-regulatory haplotype block. This notion awaits further validation.