Myotonic dystrophy (DM), the most common form of muscular dystrophy in adults, is a clinically and genetically heterogeneous neuromuscular disorder. DM is characterized by autosomal dominant inheritance, muscular dystrophy, myotonia, and multisystem involvement. Type 1 DM (DM1) is caused by a (CTG)(n) expansion in the 3' untranslated region of DMPK in 19q13.3. Multiple families, predominantly of German descent and with clinically variable presentation that included proximal myotonic myopathy (PROMM) and type 2 DM (DM2) but without the DM1 mutation, showed linkage to the 3q21 region and were recently shown to segregate a (CCTG)(n) expansion mutation in intron 1 of ZNF9. Here, we present linkage to 3q21 and mutational confirmation in 17 kindreds of European origin with PROMM and proximal myotonic dystrophy, from geographically distinct populations. All patients have the DM2 (CCTG)(n) expansion. To study the evolution of this mutation, we constructed a comprehensive physical map of the DM2 region around ZNF9. High-resolution haplotype analysis of disease chromosomes with five microsatellite and 22 single-nucleotide polymorphism markers around the DM2 mutation identified extensive linkage disequilibrium and a single shared haplotype of at least 132 kb among patients from the different populations. With the exception of the (CCTG)(n) expansion, the available markers indicate that the DM2 haplotype is identical to the most common haplotype in normal individuals. This situation is reminiscent of that seen in DM1. Taken together, these data suggest a single founding mutation in DM2 patients of European origin. We estimate the age of the founding haplotype and of the DM2 (CCTG) expansion mutation to be approximately 200-540 generations.

To date, the contribution of disrupted potentially cis-regulatory conserved non-coding sequences (CNCs) to human disease is most likely underestimated, as no systematic screens for putative deleterious variations in CNCs have been conducted. As a model for monogenic disease we studied the involvement of genetic changes of CNCs in the cis-regulatory domain of FOXL2 in blepharophimosis syndrome (BPES). Fifty-seven molecularly unsolved BPES patients underwent high-resolution copy number screening and targeted sequencing of CNCs. Apart from three larger distant deletions, a de novo deletion as small as 7.4 kb was found at 283 kb 5' to FOXL2. The deletion appeared to be triggered by an H-DNA-induced double-stranded break (DSB). In addition, it disrupts a novel long non-coding RNA (ncRNA) PISRT1 and 8 CNCs. The regulatory potential of the deleted CNCs was substantiated by in vitro luciferase assays. Interestingly, Chromosome Conformation Capture (3C) of a 625 kb region surrounding FOXL2 in expressing cellular systems revealed physical interactions of three upstream fragments and the FOXL2 core promoter. Importantly, one of these contains the 7.4 kb deleted fragment. Overall, this study revealed the smallest distant deletion causing monogenic disease and impacts upon the concept of mutation screening in human disease and developmental disorders in particular.

Ascertaining when and where genes are expressed is of crucial importance to understanding or predicting the physiological role of genes and proteins and how they interact to form the complex networks that underlie organ development and function. It is, therefore, crucial to determine on a genome-wide level, the spatio-temporal gene expression profiles at cellular resolution. This information is provided by colorimetric RNA in situ hybridization that can elucidate expression of genes in their native context and does so at cellular resolution. We generated what is to our knowledge the first genome-wide transcriptome atlas by RNA in situ hybridization of an entire mammalian organism, the developing mouse at embryonic day 14.5. This digital transcriptome atlas, the Eurexpress atlas (http://www.eurexpress.org), consists of a searchable database of annotated images that can be interactively viewed. We generated anatomy-based expression profiles for over 18,000 coding genes and over 400 microRNAs. We identified 1,002 tissue-specific genes that are a source of novel tissue-specific markers for 37 different anatomical structures. The quality and the resolution of the data revealed novel molecular domains for several developing structures, such as the telencephalon, a novel organization for the hypothalamus, and insight on the Wnt network involved in renal epithelial differentiation during kidney development. The digital transcriptome atlas is a powerful resource to determine co-expression of genes, to identify cell populations and lineages, and to identify functional associations between genes relevant to development and disease.

Variation in cellular gene expression levels has been shown to be inherited. Expression is controlled at transcriptional and post-transcriptional levels. Internal ribosome entry sites (IRES) are used by viruses to bypass inhibition of cap-dependent translation, and by eukaryotic cells to control translation under conditions when protein synthesis is inhibited. We aimed at identifying genomic determinants of variability in IRES-mediated translation of viral [Encephalomyocarditis virus (EMCV)] and cellular IRES [X-linked inhibitor-of-apoptosis (XIAP) and c-myc]. Bicistronic lentiviral constructs expressing two fluorescent reporters were used to transduce laboratory and B lymphoblastoid cell lines [15 CEPH pedigrees (n = 205) and 50 unrelated individuals]. IRES efficiency varied according to cell type and among individuals. Control of IRES activity has a significant genetic component (h(2) of 0.47 and 0.36 for EMCV and XIAP, respectively). Quantitative linkage analysis identified a suggestive locus (LOD 2.35) on chromosome 18q21.2, and genome-wide association analysis revealed of a cluster of SNPs on chromosome 3, intronic to the FHIT gene, marginally associated (P = 5.9E-7) with XIAP IRES function. This study illustrates the in vitro generation of intermediate phenotypes by using cell lines for the evaluation of genetic determinants of control of elements such as IRES.

The HSA21 encoded Single-minded 2 (SIM2) transcription factor has key neurological functions and is a good candidate to be involved in the cognitive impairment of Down syndrome. We aimed to explore the functional capacity of SIM2 by mapping its DNA binding sites in mouse embryonic stem cells. ChIP-sequencing revealed 1229 high-confidence SIM2-binding sites. Analysis of the SIM2 target genes confirmed the importance of SIM2 in developmental and neuronal processes and indicated that SIM2 may be a master transcription regulator. Indeed, SIM2 DNA binding sites share sequence specificity and overlapping domains of occupancy with master transcription factors such as SOX2, OCT4 (Pou5f1), NANOG or KLF4. The association between SIM2 and these pioneer factors is supported by co-immunoprecipitation of SIM2 with SOX2, OCT4, NANOG or KLF4. Furthermore, the binding of SIM2 marks a particular sub-category of enhancers known as super-enhancers. These regions are characterized by typical DNA modifications and Mediator co-occupancy (MED1 and MED12). Altogether, we provide evidence that SIM2 binds a specific set of enhancer elements thus explaining how SIM2 can regulate its gene network in neuronal features.

The study of gene expression in mammalian single cells via genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblasts over 133,633 unique heterozygous single-nucleotide variants (hetSNVs). We observed that at the snapshot of analyses, each cell contained mostly transcripts from one allele from the majority of genes; indeed, 76.4% of the hetSNVs displayed stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibited a haplotype-consistent allelic ratio; in contrast, distant sites located in two different genes were independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in the majority of the single cells at a given time point were rare and enriched with highly expressed genes. The relative abundance of each allele in a cell was controlled by some regulatory mechanisms given that we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

DNA double strand breaks (DSBs) are mainly repaired via homologous recombination (HR) or nonhomologous end joining (NHEJ). These breaks pose severe threats to genome integrity but can also be necessary intermediates of normal cellular processes such as immunoglobulin class switch recombination (CSR). During CSR, DSBs are produced in the G1 phase of the cell cycle and are repaired by the classical NHEJ machinery. By studying B lymphocytes derived from patients with Cornelia de Lange Syndrome, we observed a strong correlation between heterozygous loss-of-function mutations in the gene encoding the cohesin loading protein NIPBL and a shift toward the use of an alternative, microhomology-based end joining during CSR. Furthermore, the early recruitment of 53BP1 to DSBs was reduced in the NIPBL-deficient patient cells. Association of NIPBL deficiency and impaired NHEJ was also observed in a plasmid-based end-joining assay and a yeast model system. Our results suggest that NIPBL plays an important and evolutionarily conserved role in NHEJ, in addition to its canonical function in sister chromatid cohesion and its recently suggested function in HR.

Hypophosphatemic rickets (HR) comprises a rare group of inherited diseases. Very recently, mutations in the dentin matrix protein 1 (DMP1) gene were identified in patients with an extremely rare autosomal recessive form of HR (ARHR). To date, very few cases of these mutations were reported.

Down syndrome (DS) is due to increased copy number of human chromosome 21. The contribution of different genetic regions has been tested using mouse models. As shown previously, the Abcg1-U2af1 genetic region contributes to cognitive defects in working and short-term recognition memory in Down syndrome mouse models. Here we analyzed the impact of monosomy of the same genetic interval, using a new mouse model, named Ms2Yah. We used several cognitive paradigms and did not detect defects in the object recognition or the Morris water maze tests. However, surprisingly, Ms2Yah mice displayed increased associative memory in a pure contextual fear-conditioning test and decreased social novelty interaction along with a larger long-term potentiation recorded in the CA1 area following stimulation of Schaffer collaterals. Whole-genome expression studies carried out on hippocampus showed that the transcription of only a small number of genes is affected, mainly from the genetic interval (Cbs, Rsph1, Wdr4), with a few additional ones, including the postsynaptic Gabrr2, Gabbr1, Grid2p, Park2, and Dlg1 and the components of the Ubiquitin-mediated proteolysis (Anapc1, Rnf7, Huwe1, Park2). The Abcg1-U2af1 region is undeniably encompassing dosage-sensitive genes or elements whose change in copy number directly affects learning and memory, synaptic function, and autistic related behavior.

DNA methylation is essential in mammalian development. We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in Down syndrome (DS). In order to determine the methylation differences in T21 without interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. Reduced Representation Bisulfite Sequencing (RRBS) revealed 35 differentially methylated promoter regions (DMRs) (Absolute methylation differences = 25%, FDR < 0.001) in MZ twins discordant for T21 that have also been observed in comparison between unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis (FDR = 1.60 e -03) and include genes of the HOXB and HOXD clusters. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) observed in the iPSC of the T21 versus normal twin. Altogether, the results of this study highlight the epigenetic effects of the extra chromosome 21 in T21 on loci outside of this chromosome that are relevant to DS associated phenotypes.

Autosomal recessive, complete TYK2 deficiency was previously described in a patient (P1) with intracellular bacterial and viral infections and features of hyper-IgE syndrome (HIES), including atopic dermatitis, high serum IgE levels, and staphylococcal abscesses. We identified seven other TYK2-deficient patients from five families and four different ethnic groups. These patients were homozygous for one of five null mutations, different from that seen in P1. They displayed mycobacterial and/or viral infections, but no HIES. All eight TYK2-deficient patients displayed impaired but not abolished cellular responses to (a) IL-12 and IFN-α/β, accounting for mycobacterial and viral infections, respectively; (b) IL-23, with normal proportions of circulating IL-17(+) T cells, accounting for their apparent lack of mucocutaneous candidiasis; and (c) IL-10, with no overt clinical consequences, including a lack of inflammatory bowel disease. Cellular responses to IL-21, IL-27, IFN-γ, IL-28/29 (IFN-λ), and leukemia inhibitory factor (LIF) were normal. The leukocytes and fibroblasts of all seven newly identified TYK2-deficient patients, unlike those of P1, responded normally to IL-6, possibly accounting for the lack of HIES in these patients. The expression of exogenous wild-type TYK2 or the silencing of endogenous TYK2 did not rescue IL-6 hyporesponsiveness, suggesting that this phenotype was not a consequence of the TYK2 genotype. The core clinical phenotype of TYK2 deficiency is mycobacterial and/or viral infections, caused by impaired responses to IL-12 and IFN-α/β. Moreover, impaired IL-6 responses and HIES do not appear to be intrinsic features of TYK2 deficiency in humans.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

The recent availability of the chicken genome sequence poses the question of whether there are human protein-coding genes conserved in chicken that are currently not included in the human gene catalog. Here, we show, using comparative gene finding followed by experimental verification of exon pairs by RT-PCR, that the addition to the multi-exonic subset of this catalog could be as little as 0.2%, suggesting that we may be closing in on the human gene set. Our protocol, however, has two shortcomings: (i) the bioinformatic screening of the predicted genes, applied to filter out false positives, cannot handle intronless genes; and (ii) the experimental verification could fail to identify expression at a specific developmental time. This highlights the importance of developing methods that could provide a reliable estimate of the number of these two types of genes.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Down syndrome (DS) is characterized by extensive phenotypic variability, with most traits occurring in only a fraction of affected individuals. Substantial gene-expression variation is present among unaffected individuals, and this variation has a strong genetic component. Since DS is caused by genomic-dosage imbalance, we hypothesize that gene-expression variation of human chromosome 21 (HSA21) genes in individuals with DS has an impact on the phenotypic variability among affected individuals. We studied gene-expression variation in 14 lymphoblastoid and 17 fibroblast cell lines from individuals with DS and an equal number of controls. Gene expression was assayed using quantitative real-time polymerase chain reaction on 100 and 106 HSA21 genes and 23 and 26 non-HSA21 genes in lymphoblastoid and fibroblast cell lines, respectively. Surprisingly, only 39% and 62% of HSA21 genes in lymphoblastoid and fibroblast cells, respectively, showed a statistically significant difference between DS and normal samples, although the average up-regulation of HSA21 genes was close to the expected 1.5-fold in both cell types. Gene-expression variation in DS and normal samples was evaluated using the Kolmogorov-Smirnov test. According to the degree of overlap in expression levels, we classified all genes into 3 groups: (A) nonoverlapping, (B) partially overlapping, and (C) extensively overlapping expression distributions between normal and DS samples. We hypothesize that, in each cell type, group A genes are the most dosage sensitive and are most likely involved in the constant DS traits, group B genes might be involved in variable DS traits, and group C genes are not dosage sensitive and are least likely to participate in DS pathological phenotypes. This study provides the first extensive data set on HSA21 gene-expression variation in DS and underscores its role in modulating the outcome of gene-dosage imbalance.

Split-hand/foot malformation with long-bone deficiency (SHFLD) is a rare, severe limb deformity characterized by tibia aplasia with or without split-hand/split-foot deformity. Identification of genetic susceptibility loci for SHFLD has been unsuccessful because of its rare incidence, variable phenotypic expression and associated anomalies, and uncertain inheritance pattern. SHFLD is usually inherited as an autosomal dominant trait with reduced penetrance, although recessive inheritance has also been postulated. We conducted a genomewide linkage analysis, using a 10K SNP array in a large consanguineous family (UR078) from the United Arab Emirates (UAE) who had disease transmission consistent with an autosomal dominant inheritance pattern. The study identified two novel SHFLD susceptibility loci at 1q42.2-q43 (nonparametric linkage [NPL] 9.8, P=.000065) and 6q14.1 (NPL 7.12, P=.000897). These results were also supported by multipoint parametric linkage analysis. Maximum multipoint LOD scores of 3.20 and 3.78 were detected for genomic locations 1q42.2-43 and 6q14.1, respectively, with the use of an autosomal dominant mode of inheritance with reduced penetrance. Haplotype analysis with informative crossovers enabled mapping of the SHFLD loci to a region of approximately 18.38 cM (8.4 Mb) between single-nucleotide polymorphisms rs1124110 and rs535043 on 1q42.2-q43 and to a region of approximately 1.96 cM (4.1 Mb) between rs623155 and rs1547251 on 6q14.1. The study identified two novel loci for the SHFLD phenotype in this UAE family.

The importance of natural gene expression variation for human behavior is undisputed, but its impact on circadian physiology remains mostly unexplored. Using umbilical cord fibroblasts, we have determined by genome-wide association how common genetic variation impacts upon cellular circadian function. Gene set enrichment points to differences in protein catabolism as one major source of clock variation in humans. The two most significant alleles regulated expression of COPS7B, a subunit of the COP9 signalosome. We further show that the signalosome complex is imported into the nucleus in timed fashion to stabilize the essential circadian protein BMAL1, a novel mechanism to oppose its proteasome-mediated degradation. Thus, circadian clock properties depend in part upon a genetically-encoded competition between stabilizing and destabilizing forces, and genetic alterations in these mechanisms provide one explanation for human chronotype.

Developmental eye defects often severely reduce vision. Despite extensive efforts, for a substantial fraction of these cases the molecular causes are unknown. Recessive eye disorders are frequent in consanguineous populations and such large families with multiple affected individuals provide an opportunity to identify recessive causative genes. We studied a Pakistani consanguineous family with three affected individuals with congenital vision loss and progressive eye degeneration. The family was analyzed by exome sequencing of one affected individual and genotyping of all family members. We have identified a non-synonymous homozygous variant (NM_001128918.2: c.1708C > G: p.Arg570Gly) in the MARK3 gene as the likely cause of the phenotype. Given that MARK3 is highly conserved in flies (I: 55%; S: 67%) we knocked down the MARK3 homologue, par-1, in the eye during development. This leads to a significant reduction in eye size, a severe loss of photoreceptors and loss of vision based on electroretinogram (ERG) recordings. Expression of the par-1 p.Arg792Gly mutation (equivalent to the MARK3 variant found in patients) in developing fly eyes also induces loss of eye tissue and reduces the ERG signals. The data in flies and human indicate that the MARK3 variant corresponds to a loss of function. We conclude that the identified mutation in MARK3 establishes a new gene-disease link, since it likely causes structural abnormalities during eye development and visual impairment in humans, and that the function of MARK3/par-1 is evolutionarily conserved in eye development.

Cargo transport along the cytoplasmic microtubular network is essential for neuronal function, and cytoplasmic dynein-1 is an established molecular motor that is critical for neurogenesis and homeostasis. We performed whole-exome sequencing, homozygosity mapping, and chromosomal microarray studies in five individuals from three independent pedigrees and identified likely-pathogenic variants in DYNC1I2 (Dynein Cytoplasmic 1 Intermediate Chain 2), encoding a component of the cytoplasmic dynein 1 complex. In a consanguineous Pakistani family with three affected individuals presenting with microcephaly, severe intellectual disability, simplification of cerebral gyration, corpus callosum hypoplasia, and dysmorphic facial features, we identified a homozygous splice donor site variant (GenBank: NM_001378.2:c.607+1G>A). We report two additional individuals who have similar neurodevelopmental deficits and craniofacial features and harbor deleterious variants; one individual bears a c.740A>G (p.Tyr247Cys) change in trans with a 374 kb deletion encompassing DYNC1I2, and an unrelated individual harbors the compound-heterozygous variants c.868C>T (p.Gln290∗) and c.740A>G (p.Tyr247Cys). Zebrafish larvae subjected to CRISPR-Cas9 gene disruption or transient suppression of dync1i2a displayed significantly altered craniofacial patterning with concomitant reduction in head size. We monitored cell death and cell cycle progression in dync1i2a zebrafish models and observed significantly increased apoptosis, likely due to prolonged mitosis caused by abnormal spindle morphology, and this finding offers initial insights into the cellular basis of microcephaly. Additionally, complementation studies in zebrafish demonstrate that p.Tyr247Cys attenuates gene function, consistent with protein structural analysis. Our genetic and functional data indicate that DYNC1I2 dysfunction probably causes an autosomal-recessive microcephaly syndrome and highlight further the critical roles of the dynein-1 complex in neurodevelopment.

Spontaneously arising mouse mutations have served as the foundation for understanding gene function for more than 100 years. We have used exome sequencing in an effort to identify the causative mutations for 172 distinct, spontaneously arising mouse models of Mendelian disorders, including a broad range of clinically relevant phenotypes. To analyze the resulting data, we developed an analytics pipeline that is optimized for mouse exome data and a variation database that allows for reproducible, user-defined data mining as well as nomination of mutation candidates through knowledge-based integration of sample and variant data. Using these new tools, putative pathogenic mutations were identified for 91 (53%) of the strains in our study. Despite the increased power offered by potentially unlimited pedigrees and controlled breeding, about half of our exome cases remained unsolved. Using a combination of manual analyses of exome alignments and whole-genome sequencing, we provide evidence that a large fraction of unsolved exome cases have underlying structural mutations. This result directly informs efforts to investigate the similar proportion of apparently Mendelian human phenotypes that are recalcitrant to exome sequencing.