Novel biallelic PISD missense variants cause spondyloepimetaphyseal dysplasia with disproportionate short stature and fragmented mitochondrial morphology.

Biallelic variants in PISD cause a phenotypic spectrum ranging from short stature with spondyloepimetaphyseal dysplasia (SEMD) to a multisystem disorder affecting eyes, ears, bones, and brain. PISD encodes the mitochondrial-localized enzyme phosphatidylserine decarboxylase. The PISD precursor is self-cleaved to generate a heteromeric mature enzyme that converts phosphatidylserine to the phospholipid phosphatidylethanolamine. We describe a 17-year-old male patient, born to unrelated healthy parents, with disproportionate short stature and SEMD, featuring platyspondyly, prominent epiphyses, and metaphyseal dysplasia. Trio genome sequencing revealed compound heterozygous PISD variants c.569C>T; p.(Ser190Leu) and c.799C>T; p.(His267Tyr) in the patient. Investigation of fibroblasts showed similar levels of the PISD precursor protein in both patient and control cells. However, patient cells had a significantly higher proportion of fragmented mitochondria compared to control cells cultured under basal condition and after treatment with 2-deoxyglucose that represses glycolysis and stimulates respiration. Structural data from the PISD orthologue in Escherichia coli suggest that the amino acid substitutions Ser190Leu and His267Tyr likely impair PISD's autoprocessing activity and/or phosphatidylethanolamine biosynthesis. Based on the data, we propose that the novel PISD p.(Ser190Leu) and p.(His267Tyr) variants likely act as hypomorphs and underlie the pure skeletal phenotype in the patient.

Novel Alu insertion in the ZEB2 gene causing Mowat-Wilson syndrome.

Alu elements are short, interspersed elements located throughout the genome, playing a role in human diversity, and occasionally causing genetic diseases. Here, we report a novel Alu insertion causing Mowat-Wilson syndrome, a rare neurodevelopmental disorder, in an 8-year-old boy displaying the typical clinical features for Mowat-Wilson syndrome. The variant was not initially detected in genome sequencing data, but through deep phenotyping, which pointed to only one plausible candidate gene, manual inspection of genome sequencing alignment data enabled us to identify a de novo heterozygous Alu insertion in exon 8 of the ZEB2 gene. Nanopore long-read sequencing confirmed the Alu insertion, leading to the formation of a premature stop codon and likely haploinsufficiency of ZEB2. This underscores the importance of deep phenotyping and mobile element insertion analysis in uncovering genetic causes of monogenic disorders as these elements might be overlooked in standard next-generation sequencing protocols.

Clinical features and genotype-phenotype correlations in epilepsy patients with de novo DYNC1H1 variants.

OBJECTIVE: DYNC1H1 variants are involved on a disease spectrum from neuromuscular disorders to neurodevelopmental disorders. DYNC1H1-related epilepsy has been reported in small cohorts. We dissect the electroclinical features of 34 patients harboring de novo DYNC1H1 pathogenic variants, identify subphenotypes on the DYNC1H1-related epilepsy spectrum, and compare the genotype-phenotype correlations observed in our cohort with the literature. METHODS: Patients harboring de novo DYNC1H1 pathogenic variants were recruited through international collaborations. Clinical data were retrospectively collected. Latent class analysis was performed to identify subphenotypes. Multivariable binary logistic regression analysis was applied to investigate the association with DYNC1H1 protein domains. RESULTS: DYNC1H1-related epilepsy presented with infantile epileptic spasms syndrome (IESS) in 17 subjects (50%), and in 25% of these individuals the epileptic phenotype evolved into Lennox-Gastaut syndrome (LGS). In 12 patients (35%), focal onset epilepsy was defined. In two patients, the epileptic phenotype consisted of generalized myoclonic epilepsy, with a progressive phenotype in one individual harboring a frameshift variant. In approximately 60% of our cohort, seizures were drug-resistant. Malformations of cortical development were noticed in 79% of our patients, mostly on the lissencephaly-pachygyria spectrum, particularly with posterior predominance in a half of them. Midline and infratentorial abnormalities were additionally reported in 45% and 27% of subjects. We have identified three main classes of subphenotypes on the DYNC1H1-related epilepsy spectrum. SIGNIFICANCE: We propose a classification in which pathogenic de novo DYNC1H1 variants feature drug-resistant IESS in half of cases with potential evolution to LGS (Class 1), developmental and epileptic encephalopathy other than IESS and LGS (Class 2), or less severe focal or genetic generalized epilepsy including a progressive phenotype (Class 3). We observed an association between stalk domain variants and Class 1 phenotypes. The variants p.Arg309His and p.Arg1962His were common and associated with Class 1 subphenotype in our cohort. These findings may aid genetic counseling of patients with DYNC1H1-related epilepsy.

Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

DNA methylation signature classification of rare disorders using publicly available methylation data.

Disease-specific DNA methylation patterns (DNAm signatures) have been established for an increasing number of genetic disorders and represent a valuable tool for classification of genetic variants of uncertain significance (VUS). Sample size and batch effects are critical issues for establishing DNAm signatures, but their impact on the sensitivity and specificity of an already established DNAm signature has not previously been tested. Here, we assessed whether publicly available DNAm data can be employed to generate a binary machine learning classifier for VUS classification, and used variants in KMT2D, the gene associated with Kabuki syndrome, together with an existing DNAm signature as proof-of-concept. Using publicly available methylation data for training, a classifier for KMT2D variants was generated, and individuals with molecularly confirmed Kabuki syndrome and unaffected individuals could be correctly classified. The present study documents the clinical utility of a robust DNAm signature even for few affected individuals, and most importantly, underlines the importance of data sharing for improved diagnosis of rare genetic disorders.

DNA methylation episignature in Gabriele-de Vries syndrome.

PURPOSE: Gabriele-de Vries syndrome (GADEVS) is a rare genetic disorder characterized by developmental delay and/or intellectual disability, hypotonia, feeding difficulties, and distinct facial features. To refine the phenotype and to better understand the molecular basis of the syndrome, we analyzed clinical data and performed genome-wide DNA methylation analysis of a series of individuals carrying a YY1 variant. METHODS: Clinical data were collected for 13 individuals not yet reported through an international call for collaboration. DNA was collected for 11 of these individuals and 2 previously reported individuals in an attempt to delineate a specific DNA methylation signature in GADEVS. RESULTS: Phenotype in most individuals overlapped with the previously described features. We described 1 individual with atypical phenotype, heterozygous for a missense variant in a domain usually not involved in individuals with YY1 pathogenic missense variations. We also described a specific peripheral blood DNA methylation profile associated with YY1 variants. CONCLUSION: We reported a distinct DNA methylation episignature in GADEVS. We expanded the clinical profile of GADEVS to include thin/sparse hair and cryptorchidism. We also highlighted the utility of DNA methylation episignature analysis for classification of variants of unknown clinical significance.

Risks and Recommendations in Prenatally Detected De Novo Balanced Chromosomal Rearrangements from Assessment of Long-Term Outcomes.

The 6%-9% risk of an untoward outcome previously established by Warburton for prenatally detected de novo balanced chromosomal rearrangements (BCRs) does not account for long-term morbidity. We performed long-term follow-up (mean 17 years) of a registry-based nationwide cohort of 41 individuals carrying a prenatally detected de novo BCR with normal first trimester screening/ultrasound scan. We observed a significantly higher frequency of neurodevelopmental and/or neuropsychiatric disorders than in a matched control group (19.5% versus 8.3%, p = 0.04), which was increased to 26.8% upon clinical follow-up. Chromosomal microarray of 32 carriers revealed no pathogenic imbalances, illustrating a low prognostic value when fetal ultrasound scan is normal. In contrast, mate-pair sequencing revealed disrupted genes (ARID1B, NPAS3, CELF4), regulatory domains of known developmental genes (ZEB2, HOXC), and complex BCRs associated with adverse outcomes. Seven unmappable autosomal-autosomal BCRs with breakpoints involving pericentromeric/heterochromatic regions may represent a low-risk group. We performed independent phenotype-aware and blinded interpretation, which accurately predicted benign outcomes (specificity = 100%) but demonstrated relatively low sensitivity for prediction of the clinical outcome in affected carriers (sensitivity = 45%-55%). This sensitivity emphasizes the challenges associated with prenatal risk prediction for long-term morbidity in the absence of phenotypic data given the still immature annotation of the morbidity genome and poorly understood long-range regulatory mechanisms. In conclusion, we upwardly revise the previous estimates of Warburton to a morbidity risk of 27% and recommend sequencing of the chromosomal breakpoints as the first-tier diagnostic test in pregnancies with a de novo BCR.

Replicative and non-replicative mechanisms in the formation of clustered CNVs are indicated by whole genome characterization.

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both "simple" and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements.

Paroxysmal Cranial Dyskinesia and Nail-Patella Syndrome Caused by a Novel Variant in the LMX1B Gene.

BACKGROUND: In a Danish family, multiple individuals in five generations present with early-onset paroxysmal cranial dyskinesia, musculoskeletal abnormalities, and kidney dysfunction. OBJECTIVE: To demonstrate linkage and to identify the underlying genetic cause of disease. METHODS: Genome-wide single-nucleotide polymorphisms analysis, Sequence-Tagged-Site marker analyses, exome sequencing, and Sanger sequencing were performed. RESULTS: Linkage analyses identified a candidate locus on chromosome 9. Exome sequencing revealed a novel variant in LMX1B present in all affected individuals, logarithm of the odds (LOD) score of z = 6.54, predicted to be damaging. Nail-patella syndrome (NPS) is caused by pathogenic variants in LMX1B encoding a transcription factor essential to cytoskeletal and kidney growth and dopaminergic and serotonergic network development. NPS is characterized by abnormal musculoskeletal features and kidney dysfunction. Movement disorders have not previously been associated with NPS. CONCLUSIONS: Paroxysmal dyskinesia is a heretofore unrecognized feature of the NPS spectrum. The pathogenic mechanism might relate to aberrant dopaminergic circuits. © 2020 International Parkinson and Movement Disorder Society.

Multigenic truncation of the semaphorin-plexin pathway by a germline chromothriptic rearrangement associated with Moebius syndrome.

Moebius syndrome (MBS) is a congenital disorder caused by paralysis of the facial and abducens nerves. Although a number of candidate genes have been suspected, so far only mutations in PLXND1 and REV3L are confirmed to cause MBS. Here, we fine mapped the breakpoints of a complex chromosomal rearrangement (CCR) 46,XY,t(7;8;11;13) in a patient with MBS, which revealed 41 clustered breakpoints with typical hallmarks of chromothripsis. Among 12 truncated protein-coding genes, SEMA3A is known to bind to the MBS-associated PLXND1. Intriguingly, the CCR also truncated PIK3CG, which in silico interacts with REVL3 encoded by the other known MBS-gene REV3L, and with the SEMA3A/PLXND1 complex via FLT1. Additional studies of other complex rearrangements may reveal whether the multiple breakpoints in germline chromothripsis may predispose to complex multigenic disorders.

A splice-site variant in the lncRNA gene RP1-140A9.1 cosegregates in the large Volkmann cataract family.

Purpose: To identify the mutation for Volkmann cataract (CTRCT8) at 1p36.33. Methods: The genes in the candidate region 1p36.33 were Sanger and parallel deep sequenced, and informative single nucleotide polymorphisms (SNPs) were identified for linkage analysis. Expression analysis with reverse transcription polymerase chain reaction (RT-PCR) of the candidate gene was performed using RNA from different human tissues. Quantitative transcription polymerase chain reaction (qRT-PCR) analysis of the GNB1 gene was performed in affected and healthy individuals. Bioinformatic analysis of the linkage regions including the candidate gene was performed. Results: Linkage analysis of the 1p36.33 CCV locus applying new marker systems obtained with Sanger and deep sequencing reduced the candidate locus from 2.1 Mb to 0.389 Mb flanked by the markers STS-22AC and rs549772338 and resulted in an logarithm of the odds (LOD) score of Z = 21.67. The identified mutation, rs763295804, affects the donor splice site in the long non-coding RNA gene RP1-140A9.1 (ENSG00000231050). The gene including splice-site junctions is conserved in primates but not in other mammalian genomes, and two alternative transcripts were shown with RT-PCR. One of these transcripts represented a lens cell-specific transcript. Meta-analysis of the Cross-Linking-Immuno-Precipitation sequencing (CLIP-Seq) data suggested the RNA binding protein (RBP) eIF4AIII is an active counterpart for RP1-140A9.1, and several miRNA and transcription factors binding sites were predicted in the proximity of the mutation. ENCODE DNase I hypersensitivity and histone methylation and acetylation data suggest the genomic region may have regulatory functions. Conclusions: The mutation in RP1-140A9.1 suggests the long non-coding RNA as the candidate cataract gene associated with the autosomal dominant inherited congenital cataract from CCV. The mutation has the potential to destroy exon/intron splicing of both transcripts of RP1-140A9.1. Sanger and massive deep resequencing of the linkage region failed to identify alternative candidates suggesting the mutation in RP1-140A9.1 is causative for the CCV phenotype.

Very short DNA segments can be detected and handled by the repair machinery during germline chromothriptic chromosome reassembly.

Analyses at nucleotide resolution reveal unexpected complexity of seemingly simple and balanced chromosomal rearrangements. Chromothripsis is a rare complex aberration involving local shattering of one or more chromosomes and reassembly of the resulting DNA segments. This can influence gene expression and cause abnormal phenotypes. We studied the structure and mechanism of a seemingly balanced de novo complex rearrangement of four chromosomes in a boy with developmental and growth delay. Microarray analysis revealed two paternal de novo deletions of 0.7 and 2.5 Mb at two of the breakpoints in 1q24.3 and 6q24.1-q24.2, respectively, which could explain most symptoms of the patient. Subsequent whole-genome mate-pair sequencing confirmed the chromothriptic nature of the rearrangement. The four participating chromosomes were broken into 29 segments longer than 1 kb. Sanger sequencing of all breakpoint junctions revealed additional complexity compatible with the involvement of different repair pathways. We observed translocation of a 33 bp long DNA fragment, which may have implications for the definition of the lower size limit of structural variants. Our observations and literature review indicate that even very small fragments from shattered chromosomes can be detected and handled by the repair machinery during germline chromothriptic chromosome reassembly.

Aggressiveness of non-EMT breast cancer cells relies on FBXO11 activity.

Tumorigenesis is increasingly considered to rely on subclones of cells poised to undergo an epithelial to mesenchymal transition (EMT) program. We and others have provided evidence, however, that the tumorigenesis of human breast cancer is not always restricted to typical EMT cells but is also somewhat paradoxically conveyed by subclones of apparently differentiated, non-EMT cells. Here we characterize such non-EMT-like and EMT-like subclones. Through a loss-of-function screen we found that a member of the E3 ubiquitin ligase complexes, FBXO11, specifically fuels tumor formation of a non-EMT-like clone by restraining the p53/p21 pathway. Interestingly, in the related EMT-like clone, FBXO11 operates through the BCL2 pathway with little or no impact on tumorigenesis. These data command caution in attempts to assess tumorigenesis prospectively based on EMT profiling, and they emphasize the importance of next generation subtyping of tumors, that is at the level of clonal composition.

A functional link between the histone demethylase PHF8 and the transcription factor ZNF711 in X-linked mental retardation.

X-linked mental retardation (XLMR) is an inherited disorder that mostly affects males and is caused by mutations in genes located on the X chromosome. Here, we show that the XLMR protein PHF8 and a C. elegans homolog F29B9.2 catalyze demethylation of di- and monomethylated lysine 9 of histone H3 (H3K9me2/me1). The PHD domain of PHF8 binds to H3K4me3 and colocalizes with H3K4me3 at transcription initiation sites. Furthermore, PHF8 interacts with another XMLR protein, ZNF711, which binds to a subset of PHF8 target genes, including the XLMR gene JARID1C. Of interest, the C. elegans PHF8 homolog is highly expressed in neurons, and mutant animals show impaired locomotion. Taken together, our results functionally link the XLMR gene PHF8 to two other XLMR genes, ZNF711 and JARID1C, indicating that MR genes may be functionally linked in pathways, causing the complex phenotypes observed in patients developing MR.

The small RNA content of human sperm reveals pseudogene-derived piRNAs complementary to protein-coding genes.

At the end of mammalian sperm development, sperm cells expel most of their cytoplasm and dispose of the majority of their RNA. Yet, hundreds of RNA molecules remain in mature sperm. The biological significance of the vast majority of these molecules is unclear. To better understand the processes that generate sperm small RNAs and what roles they may have, we sequenced and characterized the small RNA content of sperm samples from two human fertile individuals. We detected 182 microRNAs, some of which are highly abundant. The most abundant microRNA in sperm is miR-1246 with predicted targets among sperm-specific genes. The most abundant class of small noncoding RNAs in sperm are PIWI-interacting RNAs (piRNAs). Surprisingly, we found that human sperm cells contain piRNAs processed from pseudogenes. Clusters of piRNAs from human testes contain pseudogenes transcribed in the antisense strand and processed into small RNAs. Several human protein-coding genes contain antisense predicted targets of pseudogene-derived piRNAs in the male germline and these piRNAs are still found in mature sperm. Our study provides the most extensive data set and annotation of human sperm small RNAs to date and is a resource for further functional studies on the roles of sperm small RNAs. In addition, we propose that some of the pseudogene-derived human piRNAs may regulate expression of their parent gene in the male germline.

Immature truncated O-glycophenotype of cancer directly induces oncogenic features.

Aberrant expression of immature truncated O-glycans is a characteristic feature observed on virtually all epithelial cancer cells, and a very high frequency is observed in early epithelial premalignant lesions that precede the development of adenocarcinomas. Expression of the truncated O-glycan structures Tn and sialyl-Tn is strongly associated with poor prognosis and overall low survival. The genetic and biosynthetic mechanisms leading to accumulation of truncated O-glycans are not fully understood and include mutation or dysregulation of glycosyltransferases involved in elongation of O-glycans, as well as relocation of glycosyltransferases controlling initiation of O-glycosylation from Golgi to endoplasmic reticulum. Truncated O-glycans have been proposed to play functional roles for cancer-cell invasiveness, but our understanding of the biological functions of aberrant glycosylation in cancer is still highly limited. Here, we used exome sequencing of most glycosyltransferases in a large series of primary and metastatic pancreatic cancers to rule out somatic mutations as a cause of expression of truncated O-glycans. Instead, we found hypermethylation of core 1 ß3-Gal-T-specific molecular chaperone, a key chaperone for O-glycan elongation, as the most prevalent cause. We next used gene editing to produce isogenic cell systems with and without homogenous truncated O-glycans that enabled, to our knowledge, the first polyomic and side-by-side evaluation of the cancer O-glycophenotype in an organotypic tissue model and in xenografts. The results strongly suggest that truncation of O-glycans directly induces oncogenic features of cell growth and invasion. The study provides support for targeting cancer-specific truncated O-glycans with immunotherapeutic measures.

Germline Chromothripsis Driven by L1-Mediated Retrotransposition and Alu/Alu Homologous Recombination.

Chromothripsis (CTH) is a phenomenon where multiple localized double-stranded DNA breaks result in complex genomic rearrangements. Although the DNA-repair mechanisms involved in CTH have been described, the mechanisms driving the localized "shattering" process remain unclear. High-throughput sequence analysis of a familial germline CTH revealed an inserted SVAE retrotransposon associated with a 110-kb deletion displaying hallmarks of L1-mediated retrotransposition. Our analysis suggests that the SVAE insertion did not occur prior to or after, but concurrent with the CTH event. We also observed L1-endonuclease potential target sites in other breakpoints. In addition, we found four Alu elements flanking the 110-kb deletion and associated with an inversion. We suggest that chromatin looping mediated by homologous Alu elements may have brought distal DNA regions into close proximity facilitating DNA cleavage by catalytically active L1-endonuclease. Our data provide the first evidence that active and inactive human retrotransposons can serve as endogenous mutagens driving CTH in the germline.

Epigenetic remodelling and dysregulation of DLGAP4 is linked with early-onset cerebellar ataxia.

Genome instability, epigenetic remodelling and structural chromosomal rearrangements are hallmarks of cancer. However, the coordinated epigenetic effects of constitutional chromosomal rearrangements that disrupt genes associated with congenital neurodevelopmental diseases are poorly understood. To understand the genetic-epigenetic interplay at breakpoints of chromosomal translocations disrupting CG-rich loci, we quantified epigenetic modifications at DLGAP4 (SAPAP4), a key post-synaptic density 95 (PSD95) associated gene, truncated by the chromosome translocation t(8;20)(p12;q11.23), co-segregating with cerebellar ataxia in a five-generation family. We report significant epigenetic remodelling of the DLGAP4 locus triggered by the t(8;20)(p12;q11.23) translocation and leading to dysregulation of DLGAP4 expression in affected carriers. Disruption of DLGAP4 results in monoallelic hypermethylation of the truncated DLGAP4 promoter CpG island. This induced hypermethylation is maintained in somatic cells of carriers across several generations in a t(8;20) dependent-manner however, is erased in the germ cells of the translocation carriers. Subsequently, chromatin remodelling of the locus-perturbed monoallelic expression of DLGAP4 mRNAs and non-coding RNAs in haploid cells having the translocation. Our results provide new mechanistic insight into the way a balanced chromosomal rearrangement associated with a neurodevelopmental disorder perturbs allele-specific epigenetic mechanisms at breakpoints leading to the deregulation of the truncated locus.

Phenotypic subregions within the split-hand/foot malformation 1 locus.

Split-hand/foot malformation 1 (SHFM1) is caused by chromosomal aberrations involving the region 7q21.3, DLX5 mutation, and dysregulation of DLX5/DLX6 expression by long-range position effects. SHFM1 can be isolated or syndromic with incomplete penetrance and a highly variable clinical expression, possibly influenced by sex and imprinting. We report on a new family with five affected individuals with syndromic SHFM1 that includes split-hand/foot malformations, hearing loss, and craniofacial anomalies, and an inv(7)(q21.3q35) present both in the proband and her affected son. The proximal inversion breakpoint, identified by next generation mate-pair sequencing, truncates the SHFM1 locus within the regulatory region of DLX5/6 expression. Through genotype-phenotype correlations of 100 patients with molecularly characterized chromosomal aberrations from 32 SHFM1 families, our findings suggest three phenotypic subregions within the SHFM1 locus associated with (1) isolated SHFM, (2) SHFM and hearing loss, and (3) SHFM, hearing loss, and craniofacial anomalies, respectively (ranked for increasing proximity to DLX5/6), and encompassing previously reported tissue-specific enhancers for DLX5/6. This uniquely well-characterized cohort of SHFM1 patients allowed us to systematically analyze the recently suggested hypothesis of skewed transmission and to confirm a higher penetrance in males vs. females in a subgroup of patients with isolated SHFM.