Position effects at the FGF8 locus are associated with femoral hypoplasia.

Copy-number variations (CNVs) are a common cause of congenital limb malformations and are interpreted primarily on the basis of their effect on gene dosage. However, recent studies show that CNVs also influence the 3D genome chromatin organization. The functional interpretation of whether a phenotype is the result of gene dosage or a regulatory position effect remains challenging. Here, we report on two unrelated families with individuals affected by bilateral hypoplasia of the femoral bones, both harboring de novo duplications on chromosome 10q24.32. The ∼0.5 Mb duplications include FGF8, a key regulator of limb development and several limb enhancer elements. To functionally characterize these variants, we analyzed the local chromatin architecture in the affected individuals' cells and re-engineered the duplications in mice by using CRISPR-Cas9 genome editing. We found that the duplications were associated with ectopic chromatin contacts and increased FGF8 expression. Transgenic mice carrying the heterozygous tandem duplication including Fgf8 exhibited proximal shortening of the limbs, resembling the human phenotype. To evaluate whether the phenotype was a result of gene dosage, we generated another transgenic mice line, carrying the duplication on one allele and a concurrent Fgf8 deletion on the other allele, as a control. Surprisingly, the same malformations were observed. Capture Hi-C experiments revealed ectopic interaction with the duplicated region and Fgf8, indicating a position effect. In summary, we show that duplications at the FGF8 locus are associated with femoral hypoplasia and that the phenotype is most likely the result of position effects altering FGF8 expression rather than gene dosage effects.

Hi-C Identifies Complex Genomic Rearrangements and TAD-Shuffling in Developmental Diseases.

Genome-wide analysis methods, such as array comparative genomic hybridization (CGH) and whole-genome sequencing (WGS), have greatly advanced the identification of structural variants (SVs) in the human genome. However, even with standard high-throughput sequencing techniques, complex rearrangements with multiple breakpoints are often difficult to resolve, and predicting their effects on gene expression and phenotype remains a challenge. Here, we address these problems by using high-throughput chromosome conformation capture (Hi-C) generated from cultured cells of nine individuals with developmental disorders (DDs). Three individuals had previously been identified as harboring duplications at the SOX9 locus and six had been identified with translocations. Hi-C resolved the positions of the duplications and was instructive in interpreting their distinct pathogenic effects, including the formation of new topologically associating domains (neo-TADs). Hi-C was very sensitive in detecting translocations, and it revealed previously unrecognized complex rearrangements at the breakpoints. In several cases, we observed the formation of fused-TADs promoting ectopic enhancer-promoter interactions that were likely to be involved in the disease pathology. In summary, we show that Hi-C is a sensible method for the detection of complex SVs in a clinical setting. The results help interpret the possible pathogenic effects of the SVs in individuals with DDs.

Complete lung agenesis caused by complex genomic rearrangements with neo-TAD formation at the SHH locus.

During human organogenesis, lung development is a timely and tightly regulated developmental process under the control of a large number of signaling molecules. Understanding how genetic variants can disturb normal lung development causing different lung malformations is a major goal for dissecting molecular mechanisms during embryogenesis. Here, through exome sequencing (ES), array CGH, genome sequencing (GS) and Hi-C, we aimed at elucidating the molecular basis of bilateral isolated lung agenesis in three fetuses born to a non-consanguineous family. We detected a complex genomic rearrangement containing duplicated, triplicated and deleted fragments involving the SHH locus in fetuses presenting complete agenesis of both lungs and near-complete agenesis of the trachea, diagnosed by ultrasound screening and confirmed at autopsy following termination. The rearrangement did not include SHH itself, but several regulatory elements for lung development, such as MACS1, a major SHH lung enhancer, and the neighboring genes MNX1 and NOM1. The rearrangement incorporated parts of two topologically associating domains (TADs) including their boundaries. Hi-C of cells from one of the affected fetuses showed the formation of two novel TADs each containing SHH enhancers and the MNX1 and NOM1 genes. Hi-C together with GS indicate that the new 3D conformation is likely causative for this condition by an inappropriate activation of MNX1 included in the neo-TADs by MACS1 enhancer, further highlighting the importance of the 3D chromatin conformation in human disease.

Biallelic UBE4A loss-of-function variants cause intellectual disability and global developmental delay.

PURPOSE: To identify novel genes associated with intellectual disability (ID) in four unrelated families. METHODS: Here, through exome sequencing and international collaboration, we report eight individuals from four unrelated families of diverse geographic origin with biallelic loss-of-function variants in UBE4A. RESULTS: Eight evaluated individuals presented with syndromic intellectual disability and global developmental delay. Other clinical features included hypotonia, short stature, seizures, and behavior disorder. Characteristic features were appreciated in some individuals but not all; in some cases, features became more apparent with age. We demonstrated that UBE4A loss-of-function variants reduced RNA expression and protein levels in clinical samples. Mice generated to mimic patient-specific Ube4a loss-of-function variant exhibited muscular and neurological/behavioral abnormalities, some of which are suggestive of the clinical abnormalities seen in the affected individuals. CONCLUSION: These data indicate that biallelic loss-of-function variants in UBE4A cause a novel intellectual disability syndrome, suggesting that UBE4A enzyme activity is required for normal development and neurological function.

A novel complex neurological phenotype due to a homozygous mutation in FDX2.

Defects in iron-sulphur [Fe-S] cluster biogenesis are increasingly recognized as causing neurological disease. Mutations in a number of genes that encode proteins involved in mitochondrial [Fe-S] protein assembly lead to complex neurological phenotypes. One class of proteins essential in the early cluster assembly are ferredoxins. FDX2 is ubiquitously expressed and is essential in the de novo formation of [2Fe-2S] clusters in humans. We describe and genetically define a novel complex neurological syndrome identified in two Brazilian families, with a novel homozygous mutation in FDX2. Patients were clinically evaluated, underwent MRI, nerve conduction studies, EMG and muscle biopsy. To define the genetic aetiology, a combination of homozygosity mapping and whole exome sequencing was performed. We identified six patients from two apparently unrelated families with autosomal recessive inheritance of a complex neurological phenotype involving optic atrophy and nystagmus developing by age 3, followed by myopathy and recurrent episodes of cramps, myalgia and muscle weakness in the first or second decade of life. Sensory-motor axonal neuropathy led to progressive distal weakness. MRI disclosed a reversible or partially reversible leukoencephalopathy. Muscle biopsy demonstrated an unusual pattern of regional succinate dehydrogenase and cytochrome c oxidase deficiency with iron accumulation. The phenotype was mapped in both families to the same homozygous missense mutation in FDX2 (c.431C > T, p.P144L). The deleterious effect of the mutation was validated by real-time reverse transcription polymerase chain reaction and western blot analysis, which demonstrated normal expression of FDX2 mRNA but severely reduced expression of FDX2 protein in muscle tissue. This study describes a novel complex neurological phenotype with unusual MRI and muscle biopsy features, conclusively mapped to a mutation in FDX2, which encodes a ubiquitously expressed mitochondrial ferredoxin essential for early [Fe-S] cluster biogenesis.

Homozygous missense mutation in MED25 segregates with syndromic intellectual disability in a large consanguineous family.

BACKGROUND: Intellectual disability (ID) is a highly heterogeneous condition affecting 2% of the population worldwide. In a field study conducted in a highly inbred area of Northeastern Brazil, we investigated a consanguineous family in which seven adults presented syndromic ID. METHODS: Genome-Wide Human SNP Array 6.0 (Affymetrix) microarray was used to determine regions of homozygosity-by-descent and whole exome sequencing (WES) was performed in one affected individual using Extended Nextera Rapid-Capture Exome and Illumina HiSeq2500. RESULTS: We found two regions with an logarithm of the odds (LOD) score of 3.234: a region spanning 4.0 Mb in 19q13.32-q13.33 and a pericentromeric 20 Mb area in chromosome 2 (2p12-q11.2). WES disclosed in the critical region of chromosome 19 a homozygous variant (c.418C>T, p.Arg140Trp) in Mediator complex subunit 25 (MED25), predicted as deleterious by PolyPhen-2, Provean, Mutation Taster and Sorting Intolerant From Tolerant (SIFT). MED25 is a component of the Mediator complex, involved in regulation of transcription of nearly all RNA polymerase II-dependent genes. Deleterious mutations in MED12, MED17 and MED23 have already been associated with ID. CONCLUSIONS: These findings demonstrate that the combination of field investigation of families in highly inbred regions with modern technologies is an effective way for identifying new genes associated with ID.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects.

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

AML with complex karyotype: extreme genomic complexity revealed by combined long-read sequencing and Hi-C technology.

Acute myeloid leukemia with complex karyotype (CK-AML) is associated with poor prognosis, which is only in part explained by underlying TP53 mutations. Especially in the presence of complex chromosomal rearrangements, such as chromothripsis, the outcome of CK-AML is dismal. However, this degree of complexity of genomic rearrangements contributes to the leukemogenic phenotype and treatment resistance of CK-AML remains largely unknown. Applying an integrative workflow for the detection of structural variants (SVs) based on Oxford Nanopore (ONT) genomic DNA long-read sequencing (gDNA-LRS) and high-throughput chromosome confirmation capture (Hi-C) in a well-defined cohort of CK-AML identified regions with an extreme density of SVs. These rearrangements consisted to a large degree of focal amplifications enriched in the proximity of mammalian-wide interspersed repeat elements, which often result in oncogenic fusion transcripts, such as USP7::MVD, or the deregulation of oncogenic driver genes as confirmed by RNA-seq and ONT direct complementary DNA sequencing. We termed this novel phenomenon chromocataclysm. Thus, our integrative SV detection workflow combing gDNA-LRS and Hi-C enables to unravel complex genomic rearrangements at a very high resolution in regions hard to analyze by conventional sequencing technology, thereby providing an important tool to identify novel important drivers underlying cancer with complex karyotypic changes.

Enhancer hijacking at the ARHGAP36 locus is associated with connective tissue to bone transformation.

Heterotopic ossification is a disorder caused by abnormal mineralization of soft tissues in which signaling pathways such as BMP, TGFß and WNT are known key players in driving ectopic bone formation. Identifying novel genes and pathways related to the mineralization process are important steps for future gene therapy in bone disorders. In this study, we detect an inter-chromosomal insertional duplication in a female proband disrupting a topologically associating domain and causing an ultra-rare progressive form of heterotopic ossification. This structural variant lead to enhancer hijacking and misexpression of ARHGAP36 in fibroblasts, validated here by orthogonal in vitro studies. In addition, ARHGAP36 overexpression inhibits TGFß, and activates hedgehog signaling and genes/proteins related to extracellular matrix production. Our work on the genetic cause of this heterotopic ossification case has revealed that ARHGAP36 plays a role in bone formation and metabolism, outlining first details of this gene contributing to bone-formation and -disease.

Integration of Hi-C with short and long-read genome sequencing reveals the structure of germline rearranged genomes.

Structural variants are a common cause of disease and contribute to a large extent to inter-individual variability, but their detection and interpretation remain a challenge. Here, we investigate 11 individuals with complex genomic rearrangements including germline chromothripsis by combining short- and long-read genome sequencing (GS) with Hi-C. Large-scale genomic rearrangements are identified in Hi-C interaction maps, allowing for an independent assessment of breakpoint calls derived from the GS methods, resulting in >300 genomic junctions. Based on a comprehensive breakpoint detection and Hi-C, we achieve a reconstruction of whole rearranged chromosomes. Integrating information on the three-dimensional organization of chromatin, we observe that breakpoints occur more frequently than expected in lamina-associated domains (LADs) and that a majority reshuffle topologically associating domains (TADs). By applying phased RNA-seq, we observe an enrichment of genes showing allelic imbalanced expression (AIG) within 100 kb around the breakpoints. Interestingly, the AIGs hit by a breakpoint (19/22) display both up- and downregulation, thereby suggesting different mechanisms at play, such as gene disruption and rearrangements of regulatory information. However, the majority of interpretable genes located 200 kb around a breakpoint do not show significant expression changes. Thus, there is an overall robustness in the genome towards large-scale chromosome rearrangements.

Diagnostic Yield of Whole Genome Sequencing After Nondiagnostic Exome Sequencing or Gene Panel in Developmental and Epileptic Encephalopathies.

OBJECTIVE: To assess the benefits and limitations of whole genome sequencing (WGS) compared to exome sequencing (ES) or multigene panel (MGP) in the molecular diagnosis of developmental and epileptic encephalopathies (DEE). METHODS: We performed WGS of 30 comprehensively phenotyped DEE patient trios that were undiagnosed after first-tier testing, including chromosomal microarray and either research ES (n = 15) or diagnostic MGP (n = 15). RESULTS: Eight diagnoses were made in the 15 individuals who received prior ES (53%): 3 individuals had complex structural variants; 5 had ES-detectable variants, which now had additional evidence for pathogenicity. Eleven diagnoses were made in the 15 MGP-negative individuals (68%); the majority (n = 10) involved genes not included in the panel, particularly in individuals with postneonatal onset of seizures and those with more complex presentations including movement disorders, dysmorphic features, or multiorgan involvement. A total of 42% of diagnoses were autosomal recessive or X-chromosome linked. CONCLUSION: WGS was able to improve diagnostic yield over ES primarily through the detection of complex structural variants (n = 3). The higher diagnostic yield was otherwise better attributed to the power of re-analysis rather than inherent advantages of the WGS platform. Additional research is required to assist in the assessment of pathogenicity of novel noncoding and complex structural variants and further improve diagnostic yield for patients with DEE and other neurogenetic disorders.

A fast method to reprogram and CRISPR/Cas9 gene editing from erythroblasts.

An efficient one-step procedure to reprogram fibroblasts into human induced pluripotent stem cells (hiPSC) and perform CRISPR/Cas9 gene editing simultaneously was recently reported. Here we show that such simultaneous reprogramming and gene editing can be efficiently done from erythroblasts. We successfully obtained human induced pluripotent stem cells colonies together with in frame and out of frame CAPN1 mutations in one or both alleles. We did not identify off-targets in edited cell lines. The entire process, from blood collection to mutated hiPSC took approximately 5 weeks, a much shorter period than standard multi-step methodologies using fibroblasts. Noteworthy, blood drawing is a less invasive procedure than a skin biopsy.

Origin and age of the causative mutations in KLC2, IMPA1, MED25 and WNT7A unravelled through Brazilian admixed populations.

The mutation age and local ancestry of chromosomal segments harbouring mutations associated with autosomal recessive (AR) disorders in Brazilian admixed populations remain unknown; additionally, inbreeding levels for these affected individuals continue to be estimated based on genealogical information. Here, we calculated inbreeding levels using a runs of homozygosity approach, mutation age and local ancestry to infer the origin of each chromosomal segments containing disorder-causing mutations in KLC2, IMPA1, MED25 and WNT7A. Genotyped data were generated from 18 patients affected by AR diseases and combined to the 1000 genome project (1KGP) and Simons genome diversity project (SGDP) databases to infer local ancestry. We found a major European contribution for mutated haplotypes with recent mutation age and inbreeding values found only in Native American and Middle East individuals. These results contribute to identifying the origin of and to understanding how these diseases are maintained and spread in Brazilian and world populations.

Waardenburg syndrome: Novel mutations in a large Brazilian sample.

This paper deals with the molecular investigation of Waardenburg syndrome (WS) in a sample of 49 clinically diagnosed probands (most from southeastern Brazil), 24 of them having the type 1 (WS1) variant (10 familial and 14 isolated cases) and 25 being affected by the type 2 (WS2) variant (five familial and 20 isolated cases). Sequential Sanger sequencing of all coding exons of PAX3, MITF, EDN3, EDNRB, SOX10 and SNAI2 genes, followed by CNV detection by MLPA of PAX3, MITF and SOX10 genes in selected cases revealed many novel pathogenic variants. Molecular screening, performed in all patients, revealed 19 causative variants (19/49 = 38.8%), six of them being large whole-exon deletions detected by MLPA, seven (four missense and three nonsense substitutions) resulting from single nucleotide substitutions (SNV), and six representing small indels. A pair of dizygotic affected female twins presented the c.430delC variant in SOX10, but the mutation, imputed to gonadal mosaicism, was not found in their unaffected parents. At least 10 novel causative mutations, described in this paper, were found in this Brazilian sample. Copy-number-variation detected by MLPA identified the causative mutation in 12.2% of our cases, corresponding to 31.6% of all causative mutations. In the majority of cases, the deletions were sporadic, since they were not present in the parents of isolated cases. Our results, as a whole, reinforce the fact that the screening of copy-number-variants by MLPA is a powerful tool to identify the molecular cause in WS patients.

Publisher Correction: Discordant congenital Zika syndrome twins show differential in vitro viral susceptibility of neural progenitor cells.

The original PDF version of this Article contained errors in the spelling of Luiz Carlos Caires-Júnior, Uirá Souto Melo, Bruno Henrique Silva Araujo, Alessandra Soares-Schanoski, Murilo Sena Amaral, Kayque Alves Telles-Silva, Vanessa van der Linden, Helio van der Linden, João Ricardo Mendes de Oliveira, Nivia Maria Rodrigues Arrais, Joanna Goes Castro Meira, Ana Jovina Barreto Bispo, Esper Abrão Cavalheiro, and Robert Andreata-Santos, which were incorrectly given as Luiz Carlos de Caires Jr., UiráSouto Melo, Bruno Silva Henrique Araujo, Alessandra Soares Schanoski, MuriloSena Amaral, Kayque Telles Alves Silva, Vanessa Van der Linden, Helio Van der Linden, João Mendes Ricardo de Oliveira, Nivia Rodrigues Maria Arrais, Joanna Castro Goes Meira, Ana JovinaBarreto Bispo, EsperAbrão Cavalheiro, and Robert Andreata Santos. Furthermore, in both the PDF and HTML versions of the Article, the top panel of Fig. 3e was incorrectly labeled '10608-1' and should have been '10608-4', and financial support from CAPES and DECIT-MS was inadvertently omitted from the Acknowledgements section. These errors have now been corrected in both the PDF and HTML versions of the Article.

Discordant congenital Zika syndrome twins show differential in vitro viral susceptibility of neural progenitor cells.

Congenital Zika syndrome (CZS) causes early brain development impairment by affecting neural progenitor cells (NPCs). Here, we analyze NPCs from three pairs of dizygotic twins discordant for CZS. We compare by RNA-Seq the NPCs derived from CZS-affected and CZS-unaffected twins. Prior to Zika virus (ZIKV) infection the NPCs from CZS babies show a significantly different gene expression signature of mTOR and Wnt pathway regulators, key to a neurodevelopmental program. Following ZIKV in vitro infection, cells from affected individuals have significantly higher ZIKV replication and reduced cell growth. Whole-exome analysis in 18 affected CZS babies as compared to 5 unaffected twins and 609 controls excludes a monogenic model to explain resistance or increased susceptibility to CZS development. Overall, our results indicate that CZS is not a stochastic event and depends on NPC intrinsic susceptibility, possibly related to oligogenic and/or epigenetic mechanisms.

[Could endogamy explain the higher prevalence of disabilities in the population of the Brazilian Northeast?]. / A endogamia explicaria a elevada prevalência de deficiências em populações do Nordeste brasileiro?

Despite conceptual inaccuracies and methods of measurement, the WHO estimates that about 10% of the world population has some form of disability. In order to investigate the prevalence and etiology of disabilities and to evaluate if they could be associated with inbreeding, a cross-sectional epidemiological study was conducted using the informant method in five communities of the state Rio Grande do Norte in the Brazilian Northeast, in which consanguineous marriages frequencies varied between 9 and 32%. The average prevalence of disabilities in the five sampled communities was 4.53%, obtained by interviews that involved 37.87% of a population of 39,054 inhabitants. On average, 25% of consanguineous and 12% of non-consanguineous couples had one or more children with disabilities. The increased rate of individuals with disabilities in the Brazilian Northeast could be associated with the maintenance of the tradition of consanguineous marriages in these populations and some of these disabilities may be caused by genetic disorders.