A challenge in clinical genomics is to predict whether copy number variation (CNV) affecting a gene or multiple genes will manifest as disease. Increasing recognition of gene dosage effects in neurodevelopmental disorders prompted us to develop a computational approach based on critical-exon (highly expressed in brain, highly conserved) examination for potential etiologic effects. Using a large CNV dataset, our updated analyses revealed significant (P < 1.64 × 10(-15)) enrichment of critical-exons within rare CNVs in cases compared to controls. Separately, we used a weighted gene co-expression network analysis (WGCNA) to construct an unbiased protein module from prenatal and adult tissues and found it significantly enriched for critical exons in prenatal (P < 1.15 × 10(-50), OR = 2.11) and adult (P < 6.03 × 10(-18), OR = 1.55) tissues. WGCNA yielded 1,206 proteins for which we prioritized the corresponding genes as likely to have a role in neurodevelopmental disorders. We compared the gene lists obtained from critical-exon and WGCNA analysis and found 438 candidate genes associated with CNVs annotated as pathogenic, or as variants of uncertain significance (VOUS), from among 10,619 developmental delay cases. We identified genes containing CNVs previously considered to be VOUS to be new candidate genes for neurodevelopmental disorders (GIT1, MVB12B and PPP1R9A) demonstrating the utility of this strategy to index the clinical effects of CNVs.

Genetic studies of autism spectrum disorder (ASD) have established that de novo duplications and deletions contribute to risk. However, ascertainment of structural variants (SVs) has been restricted by the coarse resolution of current approaches. By applying a custom pipeline for SV discovery, genotyping, and de novo assembly to genome sequencing of 235 subjects (71 affected individuals, 26 healthy siblings, and their parents), we compiled an atlas of 29,719 SV loci (5,213/genome), comprising 11 different classes. We found a high diversity of de novo mutations, the majority of which were undetectable by previous methods. In addition, we observed complex mutation clusters where combinations of de novo SVs, nucleotide substitutions, and indels occurred as a single event. We estimate a high rate of structural mutation in humans (20%) and propose that genetic risk for ASD is attributable to an elevated frequency of gene-disrupting de novo SVs, but not an elevated rate of genome rearrangement.

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.

To assess the relative impact of inherited and de novo variants on autism risk, we generated a comprehensive set of exonic single-nucleotide variants (SNVs) and copy number variants (CNVs) from 2,377 families with autism. We find that private, inherited truncating SNVs in conserved genes are enriched in probands (odds ratio = 1.14, P = 0.0002) in comparison to unaffected siblings, an effect involving significant maternal transmission bias to sons. We also observe a bias for inherited CNVs, specifically for small (<100 kb), maternally inherited events (P = 0.01) that are enriched in CHD8 target genes (P = 7.4 × 10(-3)). Using a logistic regression model, we show that private truncating SNVs and rare, inherited CNVs are statistically independent risk factors for autism, with odds ratios of 1.11 (P = 0.0002) and 1.23 (P = 0.01), respectively. This analysis identifies a second class of candidate genes (for example, RIMS1, CUL7 and LZTR1) where transmitted mutations may create a sensitized background but are unlikely to be completely penetrant.

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.

Next-generation sequencing has made possible the detection of rare variant (RV) associations with quantitative traits (QT). Due to high sequencing cost, many studies can only sequence a modest number of selected samples with extreme QT. Therefore association testing in individual studies can be underpowered. Besides the primary trait, many clinically important secondary traits are often measured. It is highly beneficial if multiple studies can be jointly analyzed for detecting associations with commonly measured traits. However, analyzing secondary traits in selected samples can be biased if sample ascertainment is not properly modeled. Some methods exist for analyzing secondary traits in selected samples, where some burden tests can be implemented. However p-values can only be evaluated analytically via asymptotic approximations, which may not be accurate. Additionally, potentially more powerful sequence kernel association tests, variable selection-based methods, and burden tests that require permutations cannot be incorporated. To overcome these limitations, we developed a unified method for analyzing secondary trait associations with RVs (STAR) in selected samples, incorporating all RV tests. Statistical significance can be evaluated either through permutations or analytically. STAR makes it possible to apply more powerful RV tests to analyze secondary trait associations. It also enables jointly analyzing multiple cohorts ascertained under different study designs, which greatly boosts power. The performance of STAR and commonly used RV association tests were comprehensively evaluated using simulation studies. STAR was also implemented to analyze a dataset from the SardiNIA project where samples with extreme low-density lipoprotein levels were sequenced. A significant association between LDLR and systolic blood pressure was identified, which is supported by pharmacogenetic studies. In summary, for sequencing studies, STAR is an important tool for detecting secondary-trait RV associations.

Inherited deafness is clinically and genetically heterogeneous. We recently mapped DFNB86, a locus associated with nonsyndromic deafness, to chromosome 16p. In this study, whole-exome sequencing was performed with genomic DNA from affected individuals from three large consanguineous families in which markers linked to DFNB86 segregate with profound deafness. Analyses of these data revealed homozygous mutation c.208G>T (p.Asp70Tyr) or c.878G>C (p.Arg293Pro) in TBC1D24 as the underlying cause of deafness in the three families. Sanger sequence analysis of TBC1D24 in an additional large family in which deafness segregates with DFNB86 identified the c.208G>T (p.Asp70Tyr) substitution. These mutations affect TBC1D24 amino acid residues that are conserved in orthologs ranging from fruit fly to human. Neither variant was observed in databases of single-nucleotide variants or in 634 chromosomes from ethnically matched control subjects. TBC1D24 in the mouse inner ear was immunolocalized predominantly to spiral ganglion neurons, indicating that DFNB86 deafness might be an auditory neuropathy spectrum disorder. Previously, six recessive mutations in TBC1D24 were reported to cause seizures (hearing loss was not reported) ranging in severity from epilepsy with otherwise normal development to epileptic encephalopathy resulting in childhood death. Two of our four families in which deafness segregates with mutant alleles of TBC1D24 were available for neurological examination. Cosegregation of epilepsy and deafness was not observed in these two families. Although the causal relationship between genotype and phenotype is not presently understood, our findings, combined with published data, indicate that recessive alleles of TBC1D24 can cause either epilepsy or nonsyndromic deafness.

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.

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.

Synesthesia is a perceptual condition in which sensory stimulation triggers anomalous sensory experiences. In colored sequence synesthesia (CSS), color experiences are triggered by sequences such as letters or numbers. We performed a family based linkage analysis to identify genetic loci responsible for the increased neural crosstalk underlying CSS. Our results implicate a 23 MB region at 16q12.2-23.1, providing the first step in understanding the molecular basis of CSS.

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.