The 22q11 low copy repeats are characterized by unprecedented size and structural variability.

Low copy repeats (LCRs) are recognized as a significant source of genomic instability, driving genome variability and evolution. The Chromosome 22 LCRs (LCR22s) mediate nonallelic homologous recombination (NAHR) leading to the 22q11 deletion syndrome (22q11DS). However, LCR22s are among the most complex regions in the genome, and their structure remains unresolved. The difficulty in generating accurate maps of LCR22s has also hindered localization of the deletion end points in 22q11DS patients. Using fiber FISH and Bionano optical mapping, we assembled LCR22 alleles in 187 cell lines. Our analysis uncovered an unprecedented level of variation in LCR22s, including LCR22A alleles ranging in size from 250 to 2000 kb. Further, the incidence of various LCR22 alleles varied within different populations. Additionally, the analysis of LCR22s in 22q11DS patients and their parents enabled further refinement of the rearrangement site within LCR22A and -D, which flank the 22q11 deletion. The NAHR site was localized to a 160-kb paralog shared between the LCR22A and -D in seven 22q11DS patients. Thus, we present the most comprehensive map of LCR22 variation to date. This will greatly facilitate the investigation of the role of LCR variation as a driver of 22q11 rearrangements and the phenotypic variability among 22q11DS patients.

Mutations in THAP11 cause an inborn error of cobalamin metabolism and developmental abnormalities.

CblX (MIM309541) is an X-linked recessive disorder characterized by defects in cobalamin (vitamin B12) metabolism and other developmental defects. Mutations in HCFC1, a transcriptional co-regulator which interacts with multiple transcription factors, have been associated with cblX. HCFC1 regulates cobalamin metabolism via the regulation of MMACHC expression through its interaction with THAP11, a THAP domain-containing transcription factor. The HCFC1/THAP11 complex potentially regulates genes involved in diverse cellular functions including cell cycle, proliferation, and transcription. Thus, it is likely that mutation of THAP11 also results in biochemical and other phenotypes similar to those observed in patients with cblX. We report a patient who presented with clinical and biochemical phenotypic features that overlap cblX, but who does not have any mutations in either MMACHC or HCFC1. We sequenced THAP11 by Sanger sequencing and discovered a potentially pathogenic, homozygous variant, c.240C > G (p.Phe80Leu). Functional analysis in the developing zebrafish embryo demonstrated that both THAP11 and HCFC1 regulate the proliferation and differentiation of neural precursors, suggesting important roles in normal brain development. The loss of THAP11 in zebrafish embryos results in craniofacial abnormalities including the complete loss of Meckel's cartilage, the ceratohyal, and all of the ceratobranchial cartilages. These data are consistent with our previous work that demonstrated a role for HCFC1 in vertebrate craniofacial development. High throughput RNA-sequencing analysis reveals several overlapping gene targets of HCFC1 and THAP11. Thus, both HCFC1 and THAP11 play important roles in the regulation of cobalamin metabolism as well as other pathways involved in early vertebrate development.

Kabuki syndrome genes KMT2D and KDM6A: functional analyses demonstrate critical roles in craniofacial, heart and brain development.

Kabuki syndrome (KS) is a rare multiple congenital anomaly syndrome characterized by distinctive facial features, global developmental delay, intellectual disability and cardiovascular and musculoskeletal abnormalities. While mutations in KMT2D have been identified in a majority of KS patients, a few patients have mutations in KDM6A. We analyzed 40 individuals clinically diagnosed with KS for mutations in KMT2D and KDM6A. Mutations were detected in KMT2D in 12 and KDM6A in 4 cases, respectively. Observed mutations included single-nucleotide variations and indels leading to frame shifts, nonsense, missense or splice-site alterations. In two cases, we discovered overlapping chromosome X microdeletions containing KDM6A. To further elucidate the functional roles of KMT2D and KDM6A, we knocked down the expression of their orthologs in zebrafish. Following knockdown of kmt2d and the two zebrafish paralogs kdm6a and kdm6al, we analyzed morphants for developmental abnormalities in tissues that are affected in individuals with KS, including craniofacial structures, heart and brain. The kmt2d morphants exhibited severe abnormalities in all tissues examined. Although the kdm6a and kdm6al morphants had similar brain abnormalities, kdm6a morphants exhibited craniofacial phenotypes, whereas kdm6al morphants had prominent defects in heart development. Our results provide further support for the similar roles of KMT2D and KDM6A in the etiology of KS by using a vertebrate model organism to provide direct evidence of their roles in the development of organs and tissues affected in KS patients.

An X-linked cobalamin disorder caused by mutations in transcriptional coregulator HCFC1.

Derivatives of vitamin B12 (cobalamin) are essential cofactors for enzymes required in intermediary metabolism. Defects in cobalamin metabolism lead to disorders characterized by the accumulation of methylmalonic acid and/or homocysteine in blood and urine. The most common inborn error of cobalamin metabolism, combined methylmalonic acidemia and hyperhomocysteinemia, cblC type, is caused by mutations in MMACHC. However, several individuals with presumed cblC based on cellular and biochemical analysis do not have mutations in MMACHC. We used exome sequencing to identify the genetic basis of an X-linked form of combined methylmalonic acidemia and hyperhomocysteinemia, designated cblX. A missense mutation in a global transcriptional coregulator, HCFC1, was identified in the index case. Additional male subjects were ascertained through two international diagnostic laboratories, and 13/17 had one of five distinct missense mutations affecting three highly conserved amino acids within the HCFC1 kelch domain. A common phenotype of severe neurological symptoms including intractable epilepsy and profound neurocognitive impairment, along with variable biochemical manifestations, was observed in all affected subjects compared to individuals with early-onset cblC. The severe reduction in MMACHC mRNA and protein within subject fibroblast lines suggested a role for HCFC1 in transcriptional regulation of MMACHC, which was further supported by the identification of consensus HCFC1 binding sites in MMACHC. Furthermore, siRNA-mediated knockdown of HCFC1 expression resulted in the coordinate downregulation of MMACHC mRNA. This X-linked disorder demonstrates a distinct disease mechanism by which transcriptional dysregulation leads to an inborn error of metabolism with a complex clinical phenotype.

Mutations in the mitochondrial cysteinyl-tRNA synthase gene, CARS2, lead to a severe epileptic encephalopathy and complex movement disorder.

BACKGROUND: Mitochondrial disease is often suspected in cases of severe epileptic encephalopathy especially when a complex movement disorder, liver involvement and progressive developmental regression are present. Although mutations in either mitochondrial DNA or POLG are often present, other nuclear defects in mitochondrial DNA replication and protein translation have been associated with a severe epileptic encephalopathy. METHODS AND RESULTS: We identified a proband with an epileptic encephalopathy, complex movement disorder and a combined mitochondrial respiratory chain enzyme deficiency. The child presented with neurological regression, complex movement disorder and intractable seizures. A combined deficiency of mitochondrial complexes I, III and IV was noted in liver tissue, along with increased mitochondrial DNA content in skeletal muscle. Incomplete assembly of complex V, using blue native polyacrylamide gel electrophoretic analysis and complex I, using western blotting, suggested a disorder of mitochondrial transcription or translation. Exome sequencing identified compound heterozygous mutations in CARS2, a mitochondrial aminoacyl-tRNA synthetase. Both mutations affect highly conserved amino acids located within the functional ligase domain of the cysteinyl-tRNA synthase. A specific decrease in the amount of charged mt-tRNA(Cys) was detected in patient fibroblasts compared with controls. Retroviral transfection of the wild-type CARS2 into patient skin fibroblasts led to the correction of the incomplete assembly of complex V, providing functional evidence for the role of CARS2 mutations in disease aetiology. CONCLUSIONS: Our findings indicate that mutations in CARS2 result in a mitochondrial translational defect as seen in individuals with mitochondrial epileptic encephalopathy.

Hcfc1b, a zebrafish ortholog of HCFC1, regulates craniofacial development by modulating mmachc expression.

Mutations in HCFC1 (MIM300019), have been recently associated with cblX (MIM309541), an X-linked, recessive disorder characterized by multiple congenital anomalies including craniofacial abnormalities. HCFC1 is a transcriptional co-regulator that modulates the expression of numerous downstream target genes including MMACHC, but it is not clear how these HCFC1 targets play a role in the clinical manifestations of cblX. To begin to elucidate the mechanism by which HCFC1 modulates disease phenotypes, we have carried out loss of function analyses in the developing zebrafish. Of the two HCFC1 orthologs in zebrafish, hcfc1a and hcfc1b, the loss of hcfc1b specifically results in defects in craniofacial development. Subsequent analysis revealed that hcfc1b regulates cranial neural crest cell differentiation and proliferation within the posterior pharyngeal arches. Further, the hcfc1b-mediated craniofacial abnormalities were rescued by expression of human MMACHC, a downstream target of HCFC1 that is aberrantly expressed in cblX. Furthermore, we tested distinct human HCFC1 mutations for their role in craniofacial development and demonstrated variable effects on MMACHC expression in humans and craniofacial development in zebrafish. Notably, several individuals with mutations in either HCFC1 or MMACHC have been reported to have mild to moderate facial dysmorphia. Thus, our data demonstrates that HCFC1 plays a role in craniofacial development, which is in part mediated through the regulation of MMACHC expression.

Copy number variation at 1q21.1 associated with neuroblastoma.

Common copy number variations (CNVs) represent a significant source of genetic diversity, yet their influence on phenotypic variability, including disease susceptibility, remains poorly understood. To address this problem in human cancer, we performed a genome-wide association study of CNVs in the childhood cancer neuroblastoma, a disease in which single nucleotide polymorphism variations are known to influence susceptibility. We first genotyped 846 Caucasian neuroblastoma patients and 803 healthy Caucasian controls at approximately 550,000 single nucleotide polymorphisms, and performed a CNV-based test for association. We then replicated significant observations in two independent sample sets comprised of a total of 595 cases and 3,357 controls. Here we describe the identification of a common CNV at chromosome 1q21.1 associated with neuroblastoma in the discovery set, which was confirmed in both replication sets. This CNV was validated by quantitative polymerase chain reaction, fluorescent in situ hybridization and analysis of matched tumour specimens, and was shown to be heritable in an independent set of 713 cancer-free parent-offspring trios. We identified a previously unknown transcript within the CNV that showed high sequence similarity to several neuroblastoma breakpoint family (NBPF) genes and represents a new member of this gene family (NBPF23). This transcript was preferentially expressed in fetal brain and fetal sympathetic nervous tissues, and the expression level was strictly correlated with CNV state in neuroblastoma cells. These data demonstrate that inherited copy number variation at 1q21.1 is associated with neuroblastoma and implicate a previously unknown neuroblastoma breakpoint family gene in early tumorigenesis of this childhood cancer.

An individual with blepharophimosis-ptosis-epicanthus inversus syndrome (BPES) and additional features expands the phenotype associated with mutations in KAT6B.

Blepharophimosis-ptosis-epicanthus inversus syndrome (BPES) is an autosomal dominant disorder caused by mutations in FOXL2. We identified an individual with BPES and additional phenotypic features who did not have a FOXL2 mutation. We used whole exome sequencing to identify a de novo mutation in KAT6B (lysine acetyltransferase 6B) in this individual. The mutation was a 2-bp insertion leading to a frameshift which resulted in a premature stop codon. The resulting truncated protein does not have the C-terminal serine/methionine transcription activation domain necessary for interaction with other transcriptional and epigenetic regulators. This mutation likely has a dominant-negative or gain-of-function effect, similar to those observed in other genetic disorders resulting from KAT6B mutations, including Say-Barber-Biesecker-Young-Simpson (SBBYSS) and genitopatellar syndrome (GTPTS). Thus, our subject's phenotype broadens the spectrum of clinical findings associated with mutations in KAT6B. Furthermore, our results suggest that individuals with BPES without a FOXL2 mutation should be tested for KAT6B mutations. The transcriptional and epigenetic regulation mediated by KAT6B appears crucial to early developmental processes, which when perturbed can lead to a wide spectrum of phenotypic outcomes.

Genes and biological processes commonly disrupted in rare and heterogeneous developmental delay syndromes.

Rare copy number variations (CNVs) are a recognized cause of common human disease. Predicting the genetic element(s) within a small CNV whose copy number loss or gain underlies a specific phenotype might be achieved reasonably rapidly for single patients. Identifying the biological processes that are commonly disrupted within a large patient cohort which possess larger CNVs, however, requires a more objective approach that exploits genomic resources. In this study, we first identified 98 large, rare CNVs within patients exhibiting multiple congenital anomalies. All patients presented with global developmental delay (DD), while other secondary symptoms such as cardiac defects, craniofacial features and seizures were varyingly presented. By applying a robust statistical procedure that matches patients' clinical phenotypes to laboratory mouse gene knockouts, we were able to strongly implicate anomalies in brain morphology and, separately, in long-term potentiation as manifestations of these DD patients' disorders. These and other significantly enriched model phenotypes provide insights into the pathoetiology of human DD and behavioral and anatomical secondary symptoms that are specific to DD patients. These enrichments set apart 103 genes, from among thousands overlapped by these CNVs, as strong candidates whose copy number change causally underlies approximately 46% of the cohort's DD syndromes and between 59 and 80% of the cohort's secondary symptoms. We also identified significantly enriched model phenotypes among genes overlapped by CNVs in both DD and learning disability cohorts, indicating a congruent etiology. These results demonstrate the high predictive potential of model organism phenotypes when implicating candidate genes for rare genomic disorders.

High-resolution mapping and analysis of copy number variations in the human genome: a data resource for clinical and research applications.

We present a database of copy number variations (CNVs) detected in 2026 disease-free individuals, using high-density, SNP-based oligonucleotide microarrays. This large cohort, comprised mainly of Caucasians (65.2%) and African-Americans (34.2%), was analyzed for CNVs in a single study using a uniform array platform and computational process. We have catalogued and characterized 54,462 individual CNVs, 77.8% of which were identified in multiple unrelated individuals. These nonunique CNVs mapped to 3272 distinct regions of genomic variation spanning 5.9% of the genome; 51.5% of these were previously unreported, and >85% are rare. Our annotation and analysis confirmed and extended previously reported correlations between CNVs and several genomic features such as repetitive DNA elements, segmental duplications, and genes. We demonstrate the utility of this data set in distinguishing CNVs with pathologic significance from normal variants. Together, this analysis and annotation provides a useful resource to assist with the assessment of CNVs in the contexts of human variation, disease susceptibility, and clinical molecular diagnostics.

Genomic regions associated with microdeletion/microduplication syndromes exhibit extreme diversity of structural variation.

Segmental duplications (SDs) are a class of long, repetitive DNA elements whose paralogs share a high level of sequence similarity with each other. SDs mediate chromosomal rearrangements that lead to structural variation in the general population as well as genomic disorders associated with multiple congenital anomalies, including the 7q11.23 (Williams-Beuren Syndrome, WBS), 15q13.3, and 16p12.2 microdeletion syndromes. Population-level characterization of SDs has generally been lacking because most techniques used for analyzing these complex regions are both labor and cost intensive. In this study, we have used a high-throughput technique to genotype complex structural variation with a single molecule, long-range optical mapping approach. We characterized SDs and identified novel structural variants (SVs) at 7q11.23, 15q13.3, and 16p12.2 using optical mapping data from 154 phenotypically normal individuals from 26 populations comprising five super-populations. We detected several novel SVs for each locus, some of which had significantly different prevalence between populations. Additionally, we localized the microdeletion breakpoints to specific paralogous duplicons located within complex SDs in two patients with WBS, one patient with 15q13.3, and one patient with 16p12.2 microdeletion syndromes. The population-level data presented here highlights the extreme diversity of large and complex SVs within SD-containing regions. The approach we outline will greatly facilitate the investigation of the role of inter-SD structural variation as a driver of chromosomal rearrangements and genomic disorders.

A de novo 8.8-Mb deletion of 21q21.1-q21.3 in an autistic male with a complex rearrangement involving chromosomes 6, 10, and 21.

We report here on a normal-appearing male with pervasive developmental disorder who was found to have a de novo, apparently balanced complex rearrangement involving chromosomes 6, 10, and 21: 46,XY,ins(21;10)(q11.2;p11.2p13)t(6;21)(p23;q11.2). Further analysis by high-density oligonucleotide microarray was performed, showing an 8.8-Mb heterozygous deletion at 21q21.1-q21.3. Interestingly, the deletion is distal to the translocation breakpoint on chromosome 21. The deletion involves 19 genes, including NCAM2 and GRIK1, both of which are associated with normal brain development and function, and have been considered as possible candidate genes in autism and other neurobehavioral disorders. This case underscores the utility of genomewide microarray analysis for the detection of copy number alterations in patients with apparently balanced complex rearrangements and abnormal phenotypes.

Discovery of a potentially deleterious variant in TMEM87B in a patient with a hemizygous 2q13 microdeletion suggests a recessive condition characterized by congenital heart disease and restrictive cardiomyopathy.

Restrictive cardiomyopathy (RCM) is a rare cause of heart muscle disease with the highest mortality rate among cardiomyopathy types. The etiology of RCM is poorly understood, although genetic causes have been implicated, and syndromic associations have been described. Here, we describe a patient with an atrial septal defect and restrictive cardiomyopathy along with craniofacial anomalies and intellectual disabilities. Initial screening using chromosomal microarray analysis (CMA) identified a maternally inherited 2q13 microdeletion. The patient had many of the features reported in previous cases with the recurrent 2q13 microdeletion syndrome. However, the inheritance of the microdeletion from an unaffected mother combined with the low incidence (10%) and milder forms of cardiac defects in previously reported cases made the clinical significance of the CMA results unclear. Whole-exome sequencing (WES) with trio-based analysis was performed and identified a paternally inherited TMEM87B mutation (c.1366A>G, p.Asn456Asp) in the patient. TMEM87B, a highly conserved, transmembrane protein of currently unknown function, lies within the critical region of the recurrent 2q13 microdeletion syndrome. Furthermore, a recent study had demonstrated that depletion of TMEM87B in zebrafish embryos affected cardiac development and led to cardiac hypoplasia. Thus, by combining CMA and WES, we potentially uncover an autosomal-recessive disorder characterized by a severe cardiac phenotype caused by mutations in TMEM87B. This study expands the spectrum of phenotypes associated with the recurrent 2q13 microdeletion syndrome and also further suggests the role of TMEM87B in its etiology, especially the cardiac pathology.

A 3.1-Mb microdeletion of 3p21.31 associated with cortical blindness, cleft lip, CNS abnormalities, and developmental delay.

We report a 3.1-Mb de novo deletion of 3p21.31 in a 3.5-year-old female with cortical blindness, cleft lip, CNS abnormalities, and gross developmental delays. Examination of the region showed approximately 80 genes to be involved in the deletion. Functional analysis of the deleted genes suggests that several of them may be important in normal neuronal maturation and function. Thus, haploinsufficiency of one or more of these genes could potentially contribute to the observed phenotype. Our patient does not have clinical features that overlap completely with either proximal or distal 3p deletions, suggesting that the deletion seen in our patient leads to a distinct clinical phenotype not described previously.