Cohen Syndrome (COH1) is a rare autosomal recessive disorder, principally identified by ocular, neural and muscular deficits. We identified three large consanguineous Pakistani families with intellectual disability and in some cases with autistic traits.

Next-generation sequencing (NGS) of exomes and genomes has accelerated the identification of genes involved in Mendelian phenotypes. However, many NGS studies fall short of identifying causal variants, with estimates for success rates as low as 25% for uncovering the pathological variant underlying disease etiology. An important reason for such failures is familial locus heterogeneity, where within a single pedigree causal variants in two or more genes underlie Mendelian trait etiology. As examples of intra- and inter-sibship familial locus heterogeneity, we present 10 consanguineous Pakistani families segregating hearing impairment due to homozygous variants in two different hearing impairment genes and a European-American pedigree in which hearing impairment is caused by four variants in three different genes. We have identified 41 additional pedigrees with syndromic and nonsyndromic hearing impairment for which a single previously reported hearing impairment gene has been identified but only segregates with the phenotype in a subset of affected pedigree members. We estimate that locus heterogeneity occurs in 15.3% (95% confidence interval: 11.9%, 19.9%) of the families in our collection. We demonstrate novel approaches to apply linkage analysis and homozygosity mapping (for autosomal recessive consanguineous pedigrees), which can be used to detect locus heterogeneity using either NGS or SNP array data. Results from linkage analysis and homozygosity mapping can also be used to group sibships or individuals most likely to be segregating the same causal variants and thereby increase the success rate of gene identification.

Ectodermal dysplasias (EDs) are a large heterogeneous group of inherited disorders exhibiting abnormalities in ectodermally derived appendages such as hair, nails, teeth and sweat glands. EDs associated with reticulated pigmentation phenotype are rare entities for which the genetic basis and pathophysiology are not well characterized. The present study describes a five generation consanguineous Pakistani family segregating an autosomal recessive form of a novel type of ectodermal dysplasia. The affected members present with sparse and woolly hair, severe nail dystrophy and reticulate skin pigmentation. After exclusion of known gene loci related with other skin disorders, genome-wide linkage analysis was performed using Illumina HumanOmniExpress beadchip SNP arrays. We linked this form of ED to human chromosome 18p11.32-p11.31 flanked by the SNPs rs9284390 (0.113Mb) and rs4797100 (3.14 Mb). A maximum two-point LOD score of 3.3 was obtained with several markers along the disease interval. The linkage interval of 3.03 Mb encompassed seventeen functional genes. However, sequence analysis of all these genes did not discover any potentially disease causing-variants. The identification of this novel locus provides additional information regarding the mapping of a rare form of ED. Further research, such as the use of whole-genome sequencing, would be expected to reveal any pathogenic mutation within the disease locus.

Discovering the genetic basis of a Mendelian phenotype establishes a causal link between genotype and phenotype, making possible carrier and population screening and direct diagnosis. Such discoveries also contribute to our knowledge of gene function, gene regulation, development, and biological mechanisms that can be used for developing new therapeutics. As of February 2015, 2,937 genes underlying 4,163 Mendelian phenotypes have been discovered, but the genes underlying ∼50% (i.e., 3,152) of all known Mendelian phenotypes are still unknown, and many more Mendelian conditions have yet to be recognized. This is a formidable gap in biomedical knowledge. Accordingly, in December 2011, the NIH established the Centers for Mendelian Genomics (CMGs) to provide the collaborative framework and infrastructure necessary for undertaking large-scale whole-exome sequencing and discovery of the genetic variants responsible for Mendelian phenotypes. In partnership with 529 investigators from 261 institutions in 36 countries, the CMGs assessed 18,863 samples from 8,838 families representing 579 known and 470 novel Mendelian phenotypes as of January 2015. This collaborative effort has identified 956 genes, including 375 not previously associated with human health, that underlie a Mendelian phenotype. These results provide insight into study design and analytical strategies, identify novel mechanisms of disease, and reveal the extensive clinical variability of Mendelian phenotypes. Discovering the gene underlying every Mendelian phenotype will require tackling challenges such as worldwide ascertainment and phenotypic characterization of families affected by Mendelian conditions, improvement in sequencing and analytical techniques, and pervasive sharing of phenotypic and genomic data among researchers, clinicians, and families.

Non-secretor status due to homozygosity for the common FUT2 variant c.461G>A (p.Trp154∗) is associated with either risk for autoimmune diseases or protection against viral diarrhea and HIV. We determined the role of FUT2 in otitis media susceptibility by obtaining DNA samples from 609 multi-ethnic families and simplex case subjects with otitis media. Exome and Sanger sequencing, linkage analysis, and Fisher exact and transmission disequilibrium tests (TDT) were performed. The common FUT2 c.604C>T (p.Arg202∗) variant co-segregates with otitis media in a Filipino pedigree (LOD = 4.0). Additionally, a rare variant, c.412C>T (p.Arg138Cys), is associated with recurrent/chronic otitis media in European-American children (p = 1.2 × 10-5) and US trios (TDT p = 0.01). The c.461G>A (p.Trp154∗) variant was also over-transmitted in US trios (TDT p = 0.01) and was associated with shifts in middle ear microbiota composition (PERMANOVA p < 10-7) and increased biodiversity. When all missense and nonsense variants identified in multi-ethnic US trios with CADD > 20 were combined, FUT2 variants were over-transmitted in trios (TDT p = 0.001). Fut2 is transiently upregulated in mouse middle ear after inoculation with non-typeable Haemophilus influenzae. Four FUT2 variants-namely p.Ala104Val, p.Arg138Cys, p.Trp154∗, and p.Arg202∗-reduced A antigen in mutant-transfected COS-7 cells, while the nonsense variants also reduced FUT2 protein levels. Common and rare FUT2 variants confer susceptibility to otitis media, likely by modifying the middle ear microbiome through regulation of A antigen levels in epithelial cells. Our families demonstrate marked intra-familial genetic heterogeneity, suggesting that multiple combinations of common and rare variants plus environmental factors influence the individual otitis media phenotype as a complex trait.

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.

Gordon syndrome (GS), or distal arthrogryposis type 3, is a rare, autosomal-dominant disorder characterized by cleft palate and congenital contractures of the hands and feet. Exome sequencing of five GS-affected families identified mutations in piezo-type mechanosensitive ion channel component 2 (PIEZO2) in each family. Sanger sequencing revealed PIEZO2 mutations in five of seven additional families studied (for a total of 10/12 [83%] individuals), and nine families had an identical c.8057G>A (p.Arg2686His) mutation. The phenotype of GS overlaps with distal arthrogryposis type 5 (DA5) and Marden-Walker syndrome (MWS). Using molecular inversion probes for targeted sequencing to screen PIEZO2, we found mutations in 24/29 (82%) DA5-affected families and one of two MWS-affected families. The presence of cleft palate was significantly associated with c.8057G>A (Fisher's exact test, adjusted p value < 0.0001). Collectively, although GS, DA5, and MWS have traditionally been considered separate disorders, our findings indicate that they are etiologically related and perhaps represent variable expressivity of the same condition.

Massively parallel whole-genome sequencing (WGS) data have ushered in a new era in human genetics. These data are now being used to understand the role of rare variants in complex traits and to advance the goals of precision medicine. The technological and computing advances that have enabled us to generate WGS data on thousands of individuals have also outpaced our ability to perform analyses in scientifically and statistically rigorous and thoughtful ways. The past several years have witnessed the application of whole-exome sequencing (WES) to complex traits and diseases. From our analysis of NHLBI Exome Sequencing Project (ESP) data, not only have a number of important disease and complex trait association findings emerged, but our collective experience offers some valuable lessons for WGS initiatives. These include caveats associated with generating automated pipelines for quality control and analysis of rare variants; the importance of studying minority populations; sample size requirements and efficient study designs for identifying rare-variant associations; and the significance of incidental findings in population-based genetic research. With the ESP as an example, we offer guidance and a framework on how to conduct a large-scale association study in the era of WGS.

Joubert syndrome (JS) is a recessive neurodevelopmental disorder characterized by hypotonia, ataxia, abnormal eye movements, and variable cognitive impairment. It is defined by a distinctive brain malformation known as the "molar tooth sign" on axial MRI. Subsets of affected individuals have malformations such as coloboma, polydactyly, and encephalocele, as well as progressive retinal dystrophy, fibrocystic kidney disease, and liver fibrosis. More than 35 genes have been associated with JS, but in a subset of families the genetic cause remains unknown. All of the gene products localize in and around the primary cilium, making JS a canonical ciliopathy. Ciliopathies are unified by their overlapping clinical features and underlying mechanisms involving ciliary dysfunction. In this work, we identify biallelic rare, predicted-deleterious ARMC9 variants (stop-gain, missense, splice-site, and single-exon deletion) in 11 individuals with JS from 8 families, accounting for approximately 1% of the disorder. The associated phenotypes range from isolated neurological involvement to JS with retinal dystrophy, additional brain abnormalities (e.g., heterotopia, Dandy-Walker malformation), pituitary insufficiency, and/or synpolydactyly. We show that ARMC9 localizes to the basal body of the cilium and is upregulated during ciliogenesis. Typical ciliopathy phenotypes (curved body shape, retinal dystrophy, coloboma, and decreased cilia) in a CRISPR/Cas9-engineered zebrafish mutant model provide additional support for ARMC9 as a ciliopathy-associated gene. Identifying ARMC9 mutations as a cause of JS takes us one step closer to a full genetic understanding of this important disorder and enables future functional work to define the central biological mechanisms underlying JS and other ciliopathies.

Autism spectrum disorder (ASD) is a highly heritable disorder of complex and heterogeneous aetiology. It is primarily characterized by altered cognitive ability including impaired language and communication skills and fundamental deficits in social reciprocity. Despite some notable successes in neuropsychiatric genetics, overall, the high heritability of ASD (~90%) remains poorly explained by common genetic risk variants. However, recent studies suggest that rare genomic variation, in particular copy number variation, may account for a significant proportion of the genetic basis of ASD. We present a large scale analysis to identify candidate genes which may contain low-frequency recessive variation contributing to ASD while taking into account the potential contribution of population differences to the genetic heterogeneity of ASD. Our strategy, homozygous haplotype (HH) mapping, aims to detect homozygous segments of identical haplotype structure that are shared at a higher frequency amongst ASD patients compared to parental controls. The analysis was performed on 1,402 Autism Genome Project trios genotyped for 1 million single nucleotide polymorphisms (SNPs). We identified 25 known and 1,218 novel ASD candidate genes in the discovery analysis including CADM2, ABHD14A, CHRFAM7A, GRIK2, GRM3, EPHA3, FGF10, KCND2, PDZK1, IMMP2L and FOXP2. Furthermore, 10 of the previously reported ASD genes and 300 of the novel candidates identified in the discovery analysis were replicated in an independent sample of 1,182 trios. Our results demonstrate that regions of HH are significantly enriched for previously reported ASD candidate genes and the observed association is independent of gene size (odds ratio 2.10). Our findings highlight the applicability of HH mapping in complex disorders such as ASD and offer an alternative approach to the analysis of genome-wide association data.

Previously, we have described the clinical and molecular characterization of a de novo 14q13.1-q21.1 microdeletion, less than 3.5 Mb in size, in a patient with severe microcephaly, psychomotor retardation, and other clinical anomalies. Here we report the characterization of the genomic structure of the human tuberin-like protein gene 1 (TULIP1; approved gene symbol GARNL1), a CpGisland-associated, brain-expressed candidate gene for the neurological findings in our patient, and its murine homologue. The human TULIP1 gene was mapped to chromosome band 14q13.2 by fluorescence in situ hybridization of BAC clone RP11-355C3 (GenBank Accession No. AL160231), containing the 3' region of the gene. TULIP1 spans about 271 kb of human genomic DNA and is divided into 41 exons. An untranscribed, processed pseudogene of TULIP1 was found on human chromosome band 9q31.1. The active locus TULIP1, encoding a predicted protein of 2036 amino acids, is expressed ubiquitously in pre- and postnatal human tissues. The murine homologue Tulip1 spans about 220 kb of mouse genomic DNA and is also divided into 41 exons, encoding a predicted protein of 2035 amino acids. No pseudogene could be found in the available mouse sequence data. Several splicing variants were found. Considering the location, expression profile, and predicted function, TULIP1 is a strong candidate for several neurological features seen in 14q deletion patients. Additionally we searched for mutations in the coding region of TULIP1 in subjects from a family with idiopathic basal ganglia calcification (IBGC; Fahr disease), previously linked to chromosome 14q. We identified two novel SNPs in the intron-exon boundaries; however, they did not segregate only with affected subjects in the predicted model of an autosomal dominant disease such as IBGC.

We identified overlapping homozygous regions within the DFNB25 locus in two Dutch and ten Pakistani families with sensorineural autosomal-recessive nonsyndromic hearing impairment (arNSHI). Only one of the families, W98-053, was not consanguineous, and its sibship pointed toward a reduced critical region of 0.9 Mb. This region contained the GRXCR1 gene, and the orthologous mouse gene was described to be mutated in the pirouette (pi) mutant with resulting hearing loss and circling behavior. Sequence analysis of the GRXCR1 gene in hearing-impaired family members revealed splice-site mutations in two Dutch families and a missense and nonsense mutation, respectively, in two Pakistani families. The splice-site mutations are predicted to cause frameshifts and premature stop codons. In family W98-053, this could be confirmed by cDNA analysis. GRXCR1 is predicted to contain a GRX-like domain. GRX domains are involved in reversible S-glutathionylation of proteins and thereby in the modulation of activity and/or localization of these proteins. The missense mutation is located in this domain, whereas the nonsense and splice-site mutations may result in complete or partial absence of the GRX-like domain or of the complete protein. Hearing loss in patients with GRXCR1 mutations is congenital and is moderate to profound. Progression of the hearing loss was observed in family W98-053. Vestibular dysfunction was observed in some but not all affected individuals. Quantitative analysis of GRXCR1 transcripts in fetal and adult human tissues revealed a preferential expression of the gene in fetal cochlea, which may explain the nonsyndromic nature of the hearing impairment.

Joubert syndrome and related disorders (JSRD) are primarily autosomal-recessive conditions characterized by hypotonia, ataxia, abnormal eye movements, and intellectual disability with a distinctive mid-hindbrain malformation. Variable features include retinal dystrophy, cystic kidney disease, and liver fibrosis. JSRD are included in the rapidly expanding group of disorders called ciliopathies, because all six gene products implicated in JSRD (NPHP1, AHI1, CEP290, RPGRIP1L, TMEM67, and ARL13B) function in the primary cilium/basal body organelle. By using homozygosity mapping in consanguineous families, we identify loss-of-function mutations in CC2D2A in JSRD patients with and without retinal, kidney, and liver disease. CC2D2A is expressed in all fetal and adult tissues tested. In ciliated cells, we observe localization of recombinant CC2D2A at the basal body and colocalization with CEP290, whose cognate gene is mutated in multiple hereditary ciliopathies. In addition, the proteins can physically interact in vitro, as shown by yeast two-hybrid and GST pull-down experiments. A nonsense mutation in the zebrafish CC2D2A ortholog (sentinel) results in pronephric cysts, a hallmark of ciliary dysfunction analogous to human cystic kidney disease. Knockdown of cep290 function in sentinel fish results in a synergistic pronephric cyst phenotype, revealing a genetic interaction between CC2D2A and CEP290 and implicating CC2D2A in cilium/basal body function. These observations extend the genetic spectrum of JSRD and provide a model system for studying extragenic modifiers in JSRD and other ciliopathies.

We have identified a large multigenerational Austrian family displaying a novel form of X-linked recessive myopathy. Affected individuals develop an adult-onset scapulo-axio-peroneal myopathy with bent-spine syndrome characterized by specific atrophy of postural muscles along with pseudoathleticism or hypertrophy and cardiac involvement. Known X-linked myopathies were excluded by simple-tandem-repeat polymorphism (STRP) and single-nucleotide polymorphism (SNP) analysis, direct gene sequencing, and immunohistochemical analysis. STRP analysis revealed significant linkage at Xq25-q27.1. Haplotype analysis based on SNP microarray data from selected family members confirmed this linkage region on the distal arm of the X chromosome, thereby narrowing down the critical interval to 12 Mb. Sequencing of functional candidate genes led to the identification of a missense mutation within the four and a half LIM domain 1 gene (FHL1), which putatively disrupts the fourth LIM domain of the protein. Mutation screening of FHL1 in a myopathy family from the UK exhibiting an almost identical phenotype revealed a 3 bp insertion mutation within the second LIM domain. FHL1 on Xq26.3 is highly expressed in skeletal and cardiac muscles. Western-blot analysis of muscle biopsies showed a marked decrease in protein expression of FHL1 in patients, in concordance with the genetic data. In summary, we have to our knowledge characterized a new disorder, X-linked myopathy with postural muscle atrophy (XMPMA), and identified FHL1 as the causative gene. This is the first FHL protein to be identified in conjunction with a human genetic disorder and further supports the role of FHL proteins in the development and maintenance of muscle tissue. Mutation screening of FHL1 should be considered for patients with uncharacterized myopathies and cardiomyopathies.

Lapachone is a natural naphthoquinone-derived compound found in Tabebuia avellanedae. It is well-known for its analgesic, anti-inflammatory, anti-microbial, diuretic, and anti-cancerous effects. However, the wound-healing effects of this compound are not known yet. The aim of this study was to investigate the wound healing activity of naphthoquinones (α-lapachone and β-lapachone) from Handroanthus impetiginosus.

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.

Germline pathogenic variants in chromatin-modifying enzymes are a common cause of pediatric developmental disorders. These enzymes catalyze reactions that regulate epigenetic inheritance via histone post-translational modifications and DNA methylation. Cytosine methylation (5-methylcytosine [5mC]) of DNA is the quintessential epigenetic mark, yet no human Mendelian disorder of DNA demethylation has yet been delineated. Here, we describe in detail a Mendelian disorder caused by the disruption of DNA demethylation. TET3 is a methylcytosine dioxygenase that initiates DNA demethylation during early zygote formation, embryogenesis, and neuronal differentiation and is intolerant to haploinsufficiency in mice and humans. We identify and characterize 11 cases of human TET3 deficiency in eight families with the common phenotypic features of intellectual disability and/or global developmental delay; hypotonia; autistic traits; movement disorders; growth abnormalities; and facial dysmorphism. Mono-allelic frameshift and nonsense variants in TET3 occur throughout the coding region. Mono-allelic and bi-allelic missense variants localize to conserved residues; all but one such variant occur within the catalytic domain, and most display hypomorphic function in an assay of catalytic activity. TET3 deficiency and other Mendelian disorders of the epigenetic machinery show substantial phenotypic overlap, including features of intellectual disability and abnormal growth, underscoring shared disease mechanisms.

Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.

Nucleolin (NCL/C23; OMIM: 164035) is a major nucleolar protein that plays a critical role in multiple processes, including ribosome assembly and maturation, chromatin decondensation, and pre-rRNA transcription. Due to its diverse functions, nucleolin has frequently been implicated in pathological processes, including cancer and viral infection. We recently identified a de novo frameshifting indel mutation of NCL, p.Gly664Glufs*70, through whole-exome sequencing of autism spectrum disorder trios. Through the transfection of constructs encoding either a wild-type human nucleolin or a mutant nucleolin with the same C-terminal sequence predicted for the autism proband, and by using co-localization with the nucleophosmin (NPM; B23) protein, we have shown that the nucleolin mutation leads to mislocalization of the NCL protein from the nucleolus to the nucleoplasm. Moreover, a construct with a nonsense mutation at the same residue, p.Gly664*, shows a very similar effect on the location of the NCL protein, thus confirming the presence of a predicted nucleolar location signal in this region of the NCL protein. Real-time fluorescence recovery experiments show significant changes in the kinetics and mobility of mutant NCL protein in the nucleoplasm of HEK293Tcells. Several other studies also report de novoNCL mutations in ASD or neurodevelopmental disorders. The altered mislocalization and dynamics of mutant NCL (p.G664Glufs*70/p.G664*) may have relevance to the etiopathlogy of NCL-related ASD and other neurodevelopmental phenotypes.

TET3 encodes an essential dioxygenase involved in epigenetic regulation through DNA demethylation. TET3 deficiency, or Beck-Fahrner syndrome (BEFAHRS; MIM: 618798), is a recently described neurodevelopmental disorder of the DNA demethylation machinery with a nonspecific phenotype resembling other chromatin-modifying disorders, but inconsistent variant types and inheritance patterns pose diagnostic challenges. Given TET3's direct role in regulating 5-methylcytosine and recent identification of syndrome-specific DNA methylation profiles, we analyzed genome-wide DNA methylation in whole blood of TET3-deficient individuals and identified an episignature that distinguishes affected and unaffected individuals and those with mono-allelic and bi-allelic pathogenic variants. Validation and testing of the episignature correctly categorized known TET3 variants and determined pathogenicity of variants of uncertain significance. Clinical utility was demonstrated when the episignature alone identified an affected individual from over 1000 undiagnosed cases and was confirmed upon distinguishing TET3-deficient individuals from those with 46 other disorders. The TET3-deficient signature - and the signature resulting from activating mutations in DNMT1 which normally opposes TET3 - are characterized by hypermethylation, which for BEFAHRS involves CpG sites that may be biologically relevant. This work expands the role of epi-phenotyping in molecular diagnosis and reveals genome-wide DNA methylation profiling as a quantitative, functional readout for characterization of this new biochemical category of disease.