A system for exact and approximate genetic linkage analysis of SNP data in large pedigrees

Functional variants in a TTTG microsatellite on 15q26.1 cause familial nonautoimmune thyroid abnormalities

Insufficient thyroid hormone production in newborns is referred to as congenital hypothyroidism. Multinodular goiter (MNG), characterized by an enlarged thyroid gland with multiple nodules, is usually seen in adults and is recognized as a separate disorder from congenital hypothyroidism. Here we performed a linkage analysis of a family with both nongoitrous congenital hypothyroidism and MNG and identified a signal at 15q26.1. Follow-up analyses with whole-genome sequencing and genetic screening in congenital hypothyroidism and MNG cohorts showed that changes in a noncoding TTTG microsatellite on 15q26.1 were frequently observed in congenital hypothyroidism (137 in 989) and MNG (3 in 33) compared with controls (3 in 38,722). Characterization of the noncoding variants with epigenomic data and in vitro experiments suggested that the microsatellite is located in a thyroid-specific transcriptional repressor, and its activity is disrupted by the variants. Collectively, we presented genetic evidence linking nongoitrous congenital hypothyroidism and MNG, providing unique insights into thyroid abnormalities.

IHH enhancer variant within neighboring NHEJ1 intron causes microphthalmia anophthalmia and coloboma

Genomic sequences residing within introns of few genes have been shown to act as enhancers affecting expression of neighboring genes. We studied an autosomal recessive phenotypic continuum of microphthalmia, anophthalmia and ocular coloboma, with no apparent coding-region disease-causing mutation. Homozygosity mapping of several affected Jewish Iranian families, combined with whole genome sequence analysis, identified a 0.5 Mb disease-associated chromosome 2q35 locus (maximal LOD score 6.8) harboring an intronic founder variant in NHEJ1, not predicted to affect NHEJ1. The human NHEJ1 intronic variant lies within a known specifically limb-development enhancer of a neighboring gene, Indian hedgehog (Ihh), known to be involved in eye development in mice and chickens. Through mouse and chicken molecular development studies, we demonstrated that this variant is within an Ihh enhancer that drives gene expression in the developing eye and that the identified variant affects this eye-specific enhancer activity. We thus delineate an Ihh enhancer active in mammalian eye development whose variant causes human microphthalmia, anophthalmia and ocular coloboma. The findings highlight disease causation by an intronic variant affecting the expression of a neighboring gene, delineating molecular pathways of eye development.

C-terminal frameshift variant of TDP-43 with pronounced aggregation-propensity causes rimmed vacuole myopathy but not ALS/FTD

Neuronal TDP-43-positive inclusions are neuropathological hallmark lesions in frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS). Pathogenic missense variants in TARDBP, the gene encoding TDP-43, can cause ALS and cluster in the C-terminal prion-like domain (PrLD), where they modulate the liquid condensation and aggregation properties of the protein. TDP-43-positive inclusions are also found in rimmed vacuole myopathies, including sporadic inclusion body myositis, but myopathy-causing TDP-43 variants have not been reported. Using genome-wide linkage analysis and whole exome sequencing in an extended five-generation family with an autosomal dominant rimmed vacuole myopathy, we identified a conclusively linked frameshift mutation in TDP-43 producing a C-terminally altered PrLD (TDP-43p.Trp385IlefsTer10) (maximum multipoint LOD-score 3.61). Patient-derived muscle biopsies showed TDP-43-positive sarcoplasmic inclusions, accumulation of autophagosomes and transcriptomes with abnormally spliced sarcomeric genes (including TTN and NEB) and increased expression of muscle regeneration genes. In vitro phase separation assays demonstrated that TDP-43Trp385IlefsTer10 does not form liquid-like condensates and readily forms solid-like fibrils indicating increased aggregation propensity compared to wild-type TDP-43. In Drosophila TDP-43p.Trp385IlefsTer10 behaved as a partial loss-of-function allele as it was able to rescue the TBPH (fly ortholog of TARDBP) neurodevelopmental lethal null phenotype while showing strongly reduced toxic gain-of-function properties upon overexpression. Accordingly, TDP-43p.Trp385IlefsTer10 showed reduced toxicity in a primary rat neuron disease model. Together, these genetic, pathological, in vitro and in vivo results demonstrate that TDP-43p.Trp385IlefsTer10 is an aggregation-prone partial loss-of-function variant that causes autosomal dominant vacuolar myopathy but not ALS/FTD. Our study genetically links TDP-43 proteinopathy to myodegeneration, and reveals a tissue-specific role of the PrLD in directing pathology.

PRDM10 directs FLCN expression in a novel disorder overlapping with Birt-Hogg-Dubé syndrome and familial lipomatosis

Birt-Hogg-Dubé syndrome (BHD) is an autosomal dominant disorder characterized by fibrofolliculomas, pulmonary cysts, pneumothoraces and renal cell carcinomas. Here, we reveal a novel hereditary disorder in a family with skin and mucosal lesions, extensive lipomatosis and renal cell carcinomas. The proband was initially diagnosed with BHD based on the presence of fibrofolliculomas, but no pathogenic germline variant was detected in FLCN, the gene associated with BHD. By whole exome sequencing we identified a heterozygous missense variant (p.(Cys677Tyr)) in a zinc-finger encoding domain of the PRDM10 gene which co-segregated with the phenotype in the family. We show that PRDM10Cys677Tyr loses affinity for a regulatory binding motif in the FLCN promoter, abrogating cellular FLCN mRNA and protein levels. Overexpressing inducible PRDM10Cys677Tyr in renal epithelial cells altered the transcription of multiple genes, showing overlap but also differences with the effects of knocking out FLCN. We propose that PRDM10 controls an extensive gene program and acts as a critical regulator of FLCN gene transcription in human cells. The germline variant PRDM10Cys677Tyr curtails cellular folliculin expression and underlies a distinguishable syndrome characterized by extensive lipomatosis, fibrofolliculomas and renal cell carcinomas.

Limb girdle muscular disease caused by HMGCR mutation and statin myopathy treatable with mevalonolactone

Myopathy is the main adverse effect of the widely prescribed statin drug class. Statins exert their beneficial effect by inhibiting HMG CoA-reductase, the rate-controlling enzyme of the mevalonate pathway. The mechanism of statin myopathy is yet to be resolved, and its treatment is insufficient. Through homozygosity mapping and whole exome sequencing, followed by functional analysis using confocal microscopy and biochemical and biophysical methods, we demonstrate that a distinct form of human limb girdle muscular disease is caused by a pathogenic homozygous loss-of-function missense mutation in HMG CoA reductase (HMGCR), encoding HMG CoA-reductase. We biochemically synthesized and purified mevalonolactone, never administered to human patients before, and establish the safety of its oral administration in mice. We then show that its oral administration is effective in treating a human patient with no significant adverse effects. Furthermore, we demonstrate that oral mevalonolactone resolved statin-induced myopathy in mice. We conclude that HMGCR mutation causes a late-onset severe progressive muscular disease, which shows similar features to statin-induced myopathy. Our findings indicate that mevalonolactone is effective both in the treatment of hereditary HMGCR myopathy and in a murine model of statin myopathy. Further large clinical trials are in place to enable the clinical use of mevalonolactone both in the rare orphan disease and in the more common statin myopathy.

Familial multiple discoid fibromas is linked to a locus on chromosome 5 including the FNIP1 gene

Previously, we reported a series of families presenting with trichodiscomas, inherited in an autosomal dominant pattern. The phenotype was named familial multiple discoid fibromas (FMDF). The genetic cause of FMDF remained unknown so far. Trichodiscomas are skin lesions previously reported to be part of the same spectrum as the fibrofolliculoma observed in Birt-Hogg-Dubé syndrome (BHD), an inherited disease caused by pathogenic variants in the FLCN gene. Given the clinical and histological differences with BHD and the exclusion of linkage with the FLCN locus, the phenotype was concluded to be distinct from BHD. We performed extensive clinical evaluations and genetic testing in ten families with FMDF. We identified a FNIP1 frameshift variant in nine families and genealogical studies showed common ancestry for eight families. Using whole exome sequencing, we identified six additional rare variants in the haplotype surrounding FNIP1, including a missense variant in the PDGFRB gene that was found to be present in all tested patients with FMDF. Genome-wide linkage analysis showed that the locus on chromosome 5 including FNIP1 was the only region reaching the maximal possible LOD score. We concluded that FMDF is linked to a haplotype on chromosome 5. Additional evaluations in families with FMDF are required to unravel the exact genetic cause underlying the phenotype. When evaluating patients with multiple trichodisomas without a pathogenic variant in the FLCN gene, further genetic testing is warranted and can include analysis of the haplotype on chromosome 5.

Sonic Hedgehog Intron Variant Associated With an Unusual Pediatric Cortical Cataract

Purpose

To identify the genetic basis of an unusual pediatric cortical cataract demonstrating autosomal dominant inheritance in a large European–Australian pedigree.

Methods

DNA from four affected individuals were exome sequenced utilizing a NimbleGen SeqCap EZ Exome V3 kit and HiSeq 2500. DNA from 12 affected and four unaffected individuals were genotyped using Human OmniExpress-24 BeadChips. Multipoint linkage and haplotyping were performed (Superlink-Online SNP). DNA from one affected individual and his unaffected father were whole-genome sequenced on a HiSeq X Ten system. Rare small insertions/deletions and single-nucleotide variants (SNVs) were identified in the disease-linked region (Golden Helix SVS). Combined Annotation Dependent Depletion (CADD) analysis predicted variant deleteriousness. Putative enhancer function and variant effects were determined using the Dual-Glo Luciferase Assay system.

Results

Linkage mapping identified a 6.23-centimorgan support interval at chromosome 7q36. A co-segregating haplotype refined the critical region to 6.03 Mbp containing 21 protein-coding genes. Whole-genome sequencing uncovered 114 noncoding variants from which CADD predicted one was highly deleterious, a novel substitution within intron-1 of the sonic hedgehog signaling molecule (SHH) gene. ENCODE data suggested this site was a putative enhancer, subsequently confirmed by luciferase reporter assays with variant-associated gene overexpression.

Conclusions

In a large pedigree, we have identified a SHH intron variant that co-segregates with an unusual pediatric cortical cataract phenotype. SHH is important for lens formation, and mutations in its receptor (PTCH1) cause syndromic cataract. Our data implicate increased function of an enhancer important for SHH expression primarily within developing eye tissues.

A Novel COCH Mutation Affects the vWFA2 Domain and Leads to a Relatively Mild DFNA9 Phenotype

OBJECTIVE

To study the genotype and phenotype of a Dutch family with autosomal dominantly inherited hearing loss.

STUDY DESIGN

Genotype-phenotype correlation study. Genetic analysis consisted of linkage analysis, variable number of tandem repeats analysis, and Sanger sequencing. Audiovestibular function was examined. Regression analysis was performed on pure tone audiometry and speech recognition scores and correlated with the age and/or level of hearing loss.

SETTING

Tertiary referral center.

PATIENTS

A large Dutch family presenting with sensorineural hearing loss.

MAIN OUTCOME MEASURES

Identification of the underlying genetic defect of the hearing loss in this family. Results of pure tone and speech audiometry, onset age, progression of hearing loss and vestibular (dys)function.

RESULTS

A novel mutation in COCH, c.1312C > T p.(Arg438Cys), cosegregates with hearing loss and a variable degree of vestibular (dys)function in this family. The reported mean age of onset of hearing loss is 33 years (range, 18-49 yr). Hearing loss primarily affects higher frequencies and its progression is relatively mild (0.8 dB/yr). Speech perception is remarkably well preserved in affected family members when compared with other DFNA9 families with different COCH mutations.

CONCLUSION

These findings expand the genotypic and phenotypic spectrum of DFNA9. The c.1312C > T mutation, which affects the vWFA2 domain, causes a relatively mild audiovestibular phenotype when compared with other COCH mutations.

Long-range cis-regulatory elements controlling GDF6 expression are essential for ear development

Molecular mechanisms governing the development of the mammalian cochlea, the hearing organ, remain largely unknown. Through genome sequencing in 3 subjects from 2 families with nonsyndromic cochlear aplasia, we identified homozygous 221-kb and 338-kb deletions in a noncoding region on chromosome 8 with an approximately 200-kb overlapping section. Genomic location of the overlapping deleted region started from approximately 350 kb downstream of GDF6, which codes for growth and differentiation factor 6. Otic lineage cells differentiated from induced pluripotent stem cells derived from an affected individual showed reduced expression of GDF6 compared with control cells. Knockout of Gdf6 in a mouse model resulted in cochlear aplasia, closely resembling the human phenotype. We conclude that GDF6 plays a necessary role in early cochlear development controlled by cis-regulatory elements located within an approximately 500-kb region of the genome in humans and that its disruption leads to deafness due to cochlear aplasia.

Pedigree-Based Gene Mapping Supports Previous Loci and Reveals Novel Suggestive Loci in Specific Language Impairment

Purpose Specific language impairment (SLI) is characterized by a delay in language acquisition despite a lack of other developmental delays or hearing loss. Genetics of SLI is poorly understood. The purpose of this study is to identify SLI genetic loci through family-based linkage mapping. Method We performed genome-wide parametric linkage analysis in six families segregating with SLI. An age-appropriate standardized omnibus language measure was used to categorically define the SLI phenotype. Results A suggestive linkage region replicated a previous region of interest with the highest logarithm of odds (LOD) score of 2.40 at 14q11.2-q13.3 in Family 489. A paternal parent-of-origin effect associated with SLI and language phenotypes on a nonsynonymous single nucleotide polymorphism (SNP) in NOP9 (14q12) was reported previously. Linkage analysis identified a new SLI locus at 15q24.3-25.3 with the highest parametric LOD score of 3.06 in Family 315 under a recessive mode of inheritance. Suggestive evidence of linkage was also revealed at 4q31.23-q35.2 in Family 300, with the highest LOD score of 2.41. Genetic linkage was not identified in the other three families included in parametric linkage analysis. Conclusions These results are the first to report genome-wide suggestive linkage with a total language standard score on an age-appropriate omnibus language measure across a wide age range. Our findings confirm previous reports of a language-associated locus on chromosome 14q, report new SLI loci, and validate the pedigree-based parametric linkage analysis approach to mapping genes for SLI. Supplemental Material https://doi.org/10.23641/asha.13203218.

How to Assess Founder Effect in Patients with Congenital Factor XIII Deficiency

Congenital factor XIII (FXIII) deficiency is an extremely rare bleeding disorder (RBD) with estimated prevalence of one per 2 million in the general population. The disorder causes different clinical manifestations such as intracranial hemorrhage (ICH), recurrent miscarriage, umbilical cord bleeding, etc. High incidence of the disorder might be due to founder effect. To assess founder effect, haplotype analysis is an important step. For this purpose, suitable and reliable genetic markers such as microsatellites (Hum FXIIIA01 and HumFXIIIA02) and single nucleotide polymorphisms (SNP) are suggested. In the present study we tried to describe evaluation of founder effect in patients with congenital FXIII deficiency via haplotype analysis using suitable genetic markers.

Germline variants in HEY2 functional domains lead to congenital heart defects and thoracic aortic aneurysms

Purpose

In this study we aimed to establish the genetic cause of a myriad of cardiovascular defects prevalent in individuals from a genetically isolated population, who were found to share a common ancestor in 1728.

Methods

Trio genome sequencing was carried out in an index patient with critical congenital heart disease (CHD), family members had either exome or Sanger sequencing. To confirm enrichment, we performed a gene-based association test and meta-analysis in two independent validation cohorts: one with 2685 CHD cases versus 4370 controls, and the other 326 cases with familial thoracic aortic aneurysms (FTAA) and dissections versus 570 ancestry-matched controls. Functional consequences of identified variants were evaluated using expression studies.

Results

We identified a loss-of-function variant in the Notch target transcription factor-encoding gene HEY2. The homozygous state (n=3) causes life-threatening congenital heart defects, while 80% of heterozygous carriers (n=20) had cardiovascular defects, mainly CHD and FTAA of the ascending aorta. We confirm enrichment of rare risk variants in HEY2 functional domains after meta-analysis (meta-SKAT p=0.018). Furthermore, we show that several identified variants lead to dysregulation of repression by HEY2.

Conclusion

A homozygous germline loss-of-function variant in HEY2 leads to critical CHD. The majority of heterozygotes show a myriad of cardiovascular defects.

Mutations in COMP cause familial carpal tunnel syndrome

Carpal tunnel syndrome (CTS) is the most common peripheral nerve entrapment syndrome, affecting a large proportion of the general population. Genetic susceptibility has been implicated in CTS, but the causative genes remain elusive. Here, we report the identification of two mutations in cartilage oligomeric matrix protein (COMP) that segregate with CTS in two large families with or without multiple epiphyseal dysplasia (MED). Both mutations impair the secretion of COMP by tenocytes, but the mutation associated with MED also perturbs its secretion in chondrocytes. Further functional characterization of the CTS-specific mutation reveals similar histological and molecular changes of tendons/ligaments in patients’ biopsies and the mouse models. The mutant COMP fails to oligomerize properly and is trapped in the ER, resulting in ER stress-induced unfolded protein response and cell death, leading to inflammation, progressive fibrosis and cell composition change in tendons/ligaments. The extracellular matrix (ECM) organization is also altered. Our studies uncover a previously unrecognized mechanism in CTS pathogenesis.

Familial carpal tunnel syndrome (CTS) is common, but causal genes are not characterized. Here the authors report two CTS-related mutations in two large families that impair secretion of COMP in tenocytes, leading to ER stress-induced unfolded protein response, inflammation and fibrosis in patients and mouse models.

Mutations in SPATA13/ASEF2 cause primary angle closure glaucoma

Current estimates suggest 50% of glaucoma blindness worldwide is caused by primary angle-closure glaucoma (PACG) but the causative gene is not known. We used genetic linkage and whole genome sequencing to identify Spermatogenesis Associated Protein 13, SPATA13 (NM_001166271; NP_001159743, SPATA13 isoform I), also known as ASEF2 (Adenomatous polyposis coli-stimulated guanine nucleotide exchange factor 2), as the causal gene for PACG in a large seven-generation white British family showing variable expression and incomplete penetrance. The 9 bp deletion, c.1432_1440del; p.478_480del was present in all affected individuals with angle-closure disease. We show ubiquitous expression of this transcript in cell lines derived from human tissues and in iris, retina, retinal pigment and ciliary epithelia, cornea and lens. We also identified eight additional mutations in SPATA13 in a cohort of 189 unrelated PACS/PAC/PACG samples. This gene encodes a 1277 residue protein which localises to the nucleus with partial co-localisation with nuclear speckles. In cells undergoing mitosis SPATA13 isoform I becomes part of the kinetochore complex co-localising with two kinetochore markers, polo like kinase 1 (PLK-1) and centrosome-associated protein E (CENP-E). The 9 bp deletion reported in this study increases the RAC1-dependent guanine nucleotide exchange factors (GEF) activity. The increase in GEF activity was also observed in three other variants identified in this study. Taken together, our data suggest that SPATA13 is involved in the regulation of mitosis and the mutations dysregulate GEF activity affecting homeostasis in tissues where it is highly expressed, influencing PACG pathogenesis.

Identification of ADHD risk genes in extended pedigrees by combining linkage analysis and whole-exome sequencing

Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder with a complex genetic background, hampering identification of underlying genetic risk factors. We hypothesized that combining linkage analysis and whole-exome sequencing (WES) in multi-generation pedigrees with multiple affected individuals can point toward novel ADHD genes. Three families with multiple ADHD-affected members (N_total = 70) and apparent dominant inheritance pattern were included in this study. Genotyping was performed in 37 family members, and WES was additionally carried out in 10 of those. Linkage analysis was performed using multi-point analysis in Superlink Online SNP 1.1. From prioritized linkage regions with a LOD score ≥ 2, a total of 24 genes harboring rare variants were selected. Those genes were taken forward and were jointly analyzed in gene-set analyses of exome-chip data using the MAGMA software in an independent sample of patients with persistent ADHD and healthy controls (N = 9365). The gene-set including all 24 genes together, and particularly the gene-set from one of the three families (12 genes), were significantly associated with persistent ADHD in this sample. Among the latter, gene-wide analysis for the AAED1 gene reached significance. A rare variant (rs151326868) within AAED1 segregated with ADHD in one of the families. The analytic strategy followed here is an effective approach for identifying novel ADHD risk genes. Additionally, this study suggests that both rare and more frequent variants in multiple genes act together in contributing to ADHD risk, even in individual multi-case families.

Loss of ADAMTS19 causes progressive non-syndromic heart valve disease

Valvular heart disease is observed in approximately 2% of the general population¹. Although the initial observation is often localized (for example, to the aortic or mitral valve), disease manifestations are regularly observed in the other valves and patients frequently require surgery. Despite the high frequency of heart valve disease, only a handful of genes have so far been identified as the monogenic causes of disease^2-7. Here we identify two consanguineous families, each with two affected family members presenting with progressive heart valve disease early in life. Whole-exome sequencing revealed homozygous, truncating nonsense alleles in ADAMTS19 in all four affected individuals. Homozygous knockout mice for Adamts19 show aortic valve dysfunction, recapitulating aspects of the human phenotype. Expression analysis using a lacZ reporter and single-cell RNA sequencing highlight Adamts19 as a novel marker for valvular interstitial cells; inference of gene regulatory networks in valvular interstitial cells positions Adamts19 in a highly discriminatory network driven by the transcription factor lymphoid enhancer-binding factor 1 downstream of the Wnt signaling pathway. Upregulation of endocardial Krüppel-like factor 2 in Adamts19 knockout mice precedes hemodynamic perturbation, showing that a tight balance in the Wnt-Adamts19-Klf2 axis is required for proper valve maturation and maintenance.

A genome-wide analysis in consanguineous families reveals new chromosomal loci in specific language impairment (SLI)

Language is a uniquely human ability, and failure to attain this ability can have a life-long impact on the affected individuals. This is particularly true for individuals with specific language impairment (SLI), which is defined as an impairment in normal language development in the absence of any other developmental disability. Although SLI displays high heritability, family-based linkage studies have been hampered by an unclear mode of Mendelian segregation, variable disease penetrance, and heterogeneity of diagnostic criteria. We performed genome-wide parametric linkage analysis and homozygosity mapping in 14 consanguineous families from Pakistan segregating SLI. Linkage analysis revealed a multipoint LOD score of 4.18 at chromosome 2q in family PKSLI05 under a recessive mode of inheritance. A second linkage score of 3.85 was observed in family PKSLI12 at a non-overlapping locus on chromosome 2q. Two other suggestive linkage loci were found in family PKSLI05 on 14q and 22q with LOD scores of 2.37 and 2.23, respectively, that were also identified in homozygosity mapping. Reduction to homozygosity was observed on chromosomes 2q, 5p, 8q, 14q, 17q, and 22q. Each homozygosity region occurred in multiple PKSLI families. We report new SLI loci on chromosomes 2 and 8 and confirm suggestive SLI linkage loci on chromosomes 5, 14, 17, and 22 reported previously in the population of Robinson Crusoe Island. These findings indicate that linkage and homozygosity mapping in consanguineous families can improve genetic analyses in SLI and suggest the involvement of additional genes in the causation of this disorder.

Hypomorphic caspase activation and recruitment domain 11 (CARD11) mutations associated with diverse immunologic phenotypes with or without atopic disease

BACKGROUND

Caspase activation and recruitment domain 11 (CARD11) encodes a scaffold protein in lymphocytes that links antigen receptor engagement with downstream signaling to nuclear factor κB, c-Jun N-terminal kinase, and mechanistic target of rapamycin complex 1. Germline CARD11 mutations cause several distinct primary immune disorders in human subjects, including severe combined immune deficiency (biallelic null mutations), B-cell expansion with nuclear factor κB and T-cell anergy (heterozygous, gain-of-function mutations), and severe atopic disease (loss-of-function, heterozygous, dominant interfering mutations), which has focused attention on CARD11 mutations discovered by using whole-exome sequencing.

OBJECTIVES

We sought to determine the molecular actions of an extended allelic series of CARD11 and to characterize the expanding range of clinical phenotypes associated with heterozygous CARD11 loss-of-function alleles.

METHODS

Cell transfections and primary T-cell assays were used to evaluate signaling and function of CARD11 variants.

RESULTS

Here we report on an expanded cohort of patients harboring novel heterozygous CARD11 mutations that extend beyond atopy to include other immunologic phenotypes not previously associated with CARD11 mutations. In addition to (and sometimes excluding) severe atopy, heterozygous missense and indel mutations in CARD11 presented with immunologic phenotypes similar to those observed in signal transducer and activator of transcription 3 loss of function, dedicator of cytokinesis 8 deficiency, common variable immunodeficiency, neutropenia, and immune dysregulation, polyendocrinopathy, enteropathy, X-linked-like syndrome. Pathogenic variants exhibited dominant negative activity and were largely confined to the CARD or coiled-coil domains of the CARD11 protein.

CONCLUSION

These results illuminate a broader phenotypic spectrum associated with CARD11 mutations in human subjects and underscore the need for functional studies to demonstrate that rare gene variants encountered in expected and unexpected phenotypes must nonetheless be validated for pathogenic activity.

Genetic linkage analysis of a large family identifies FIGN as a candidate modulator of reduced penetrance in heritable pulmonary arterial hypertension

BACKGROUND

Mapping the genetic component of molecular mechanisms responsible for the reduced penetrance (RP) of rare disorders constitutes one of the most challenging problems in human genetics. Heritable pulmonary arterial hypertension (PAH) is one such disorder characterised by rare mutations mostly occurring in the bone morphogenetic protein receptor type 2 (BMPR2) gene and a wide heterogeneity of penetrance modifier mechanisms. Here, we analyse 32 genotyped individuals from a large Iberian family of 65 members, including 22 carriers of the pathogenic BMPR2 mutation c.1472G>A (p.Arg491Gln), 8 of them diagnosed with PAH by right-heart catheterisation, leading to an RP rate of 36.4%.

METHODS

We performed a linkage analysis on the genotyping data to search for genetic modifiers of penetrance. Using functional genomics data, we characterised the candidate region identified by linkage analysis. We also predicted the haplotype segregation within the family.

RESULTS

We identified a candidate chromosome region in 2q24.3, 38 Mb upstream from BMPR2, with significant linkage (LOD=4.09) under a PAH susceptibility model. This region contains common variants associated with vascular aetiology and shows functional evidence that the putative genetic modifier is located in the upstream distal promoter of the fidgetin (FIGN) gene.

CONCLUSION

Our results suggest that the genetic modifier acts through FIGN transcriptional regulation, whose expression variability would contribute to modulating heritable PAH. This finding may help to advance our understanding of RP in PAH across families sharing the p.Arg491Gln pathogenic mutation in BMPR2.

Homozygous frameshift mutations in FAT1 cause a syndrome characterized by colobomatous-microphthalmia, ptosis, nephropathy and syndactyly

A failure in optic fissure fusion during development can lead to blinding malformations of the eye. Here, we report a syndrome characterized by facial dysmorphism, colobomatous microphthalmia, ptosis and syndactyly with or without nephropathy, associated with homozygous frameshift mutations in FAT1. We show that Fat1 knockout mice and zebrafish embryos homozygous for truncating fat1a mutations exhibit completely penetrant coloboma, recapitulating the most consistent developmental defect observed in affected individuals. In human retinal pigment epithelium (RPE) cells, the primary site for the fusion of optic fissure margins, FAT1 is localized at earliest cell-cell junctions, consistent with a role in facilitating optic fissure fusion during vertebrate eye development. Our findings establish FAT1 as a gene with pleiotropic effects in human, in that frameshift mutations cause a severe multi-system disorder whereas recessive missense mutations had been previously associated with isolated glomerulotubular nephropathy.

Loss of the cadherin FAT1 has been associated with nephropathy and epithelial cell adhesion defects. Here, the authors report five families with a syndromic form of coloboma associated with homozygous frameshift variants in FAT1 and recapitulate the phenotype in mutant mice and zebrafish.

SCAPER localizes to primary cilia and its mutation affects cilia length, causing Bardet-Biedl syndrome

Studies of ciliopathies have served in elucidating much of our knowledge of structure and function of primary cilia. We report humans with Bardet-Biedl syndrome who display intellectual disability, retinitis pigmentosa, obesity, short stature and brachydactyly, stemming from a homozyogous truncation mutation in SCAPER, a gene previously associated with mitotic progression. Our findings, based on linkage analysis and exome sequencing studies of two remotely related large consanguineous families, are in line with recent reports of SCAPER variants associated with intellectual disability and retinitis pigmentosa. Using immuno-fluorescence and live cell imaging in NIH/3T3 fibroblasts and SH-SY5Y neuroblastoma cell lines over-expressing SCAPER, we demonstrate that both wild type and mutant SCAPER are expressed in primary cilia and co-localize with tubulin, forming bundles of microtubules. While wild type SCAPER was rarely localized along the ciliary axoneme and basal body, the aberrant protein remained sequestered to the cilia, mostly at the ciliary tip. Notably, longer cilia were demonstrated both in human affected fibroblasts compared to controls, as well as in NIH/3T3 cells transfected with mutant versus wildtype SCAPER. As SCAPER expression is known to peak at late G1 and S phase, overlapping the timing of ciliary resorption, our data suggest a possible role of SCAPER in ciliary dynamics and disassembly, also affecting microtubule-related mitotic progression. Thus, we outline a human ciliopathy syndrome and demonstrate that it is caused by a mutation in SCAPER, affecting primary cilia.

Combined CNV, haplotyping and whole exome sequencing enable identification of two distinct novel EYS mutations causing RP in a single inbred tribe

Whole exome sequencing (WES) has become routine in clinical practice, especially in studies of recessive hereditary diseases in inbred consanguineous families, where homozygosity of a founder mutation is assumed. Multiple members of two consanguineous families of a single Bedouin tribe were diagnosed with apparently autosomal recessive/pseudo-dominant retinitis pigmentosa (RP). Affected individuals exhibited severe visual impairment with nyctalopia, marked constriction of visual fields, markedly reduced and delayed responses on electro-retinography (ERG) and eventual loss of central vision. Combined copy-number variant (CNV) analysis, haplotype reconstruction and WES of the kindred identified two distinct novel mutations in EYS (RP25): a p.(W1817*) nonsense mutation (identified through WES) and a large deletion encompassing 9 of the 43 exons, that was missed by WES and was identified through microarray CNV analysis. Segregation analysis of both mutations demonstrated that all affected individuals were either homozygous for one of the mutations, or compound heterozygous for both. The two mutations are predicted to cause loss of function of the encoded protein and were not present in screening of 200 ethnically-matched controls. Our findings of two distinct mutations in the same gene in a single inbred kindred, identified only through combined WES and microarray CNV analysis, highlight the limitations of either CNV or WES alone, as the heterozygous deletion had normal WES read-depth values. Moreover, they demonstrate pitfalls in homozygosity mapping for disease-causing variant identification in inbred communities.

MPZL2 is a novel gene associated with autosomal recessive nonsyndromic moderate hearing loss

While recent studies have revealed a substantial portion of the genes underlying human hearing loss, the extensive genetic landscape has not been completely explored. Here, we report a loss-of-function variant (c.72delA) in MPZL2 in three unrelated multiplex families from Turkey and Iran with autosomal recessive nonsyndromic hearing loss. The variant co-segregates with moderate sensorineural hearing loss in all three families. We show a shared haplotype flanking the variant in our families implicating a single founder. While rare in other populations, the allele frequency of the variant is ~ 0.004 in Ashkenazi Jews, suggesting that it may be an important cause of moderate hearing loss in that population. We show that Mpzl2 is expressed in mouse inner ear, and the protein localizes in the auditory inner and outer hair cells, with an asymmetric subcellular localization. We thus present MPZL2 as a novel gene associated with sensorineural hearing loss.

Heterozygous missense variants of LMX1A lead to nonsyndromic hearing impairment and vestibular dysfunction

Unraveling the causes and pathomechanisms of progressive disorders is essential for the development of therapeutic strategies. Here, we identified heterozygous pathogenic missense variants of LMX1A in two families of Dutch origin with progressive nonsyndromic hearing impairment (HI), using whole exome sequencing. One variant, c.721G > C (p.Val241Leu), occurred de novo and is predicted to affect the homeodomain of LMX1A, which is essential for DNA binding. The second variant, c.290G > C (p.Cys97Ser), predicted to affect a zinc-binding residue of the second LIM domain that is involved in protein–protein interactions. Bi-allelic deleterious variants of Lmx1a are associated with a complex phenotype in mice, including deafness and vestibular defects, due to arrest of inner ear development. Although Lmx1a mouse mutants demonstrate neurological, skeletal, pigmentation and reproductive system abnormalities, no syndromic features were present in the participating subjects of either family. LMX1A has previously been suggested as a candidate gene for intellectual disability, but our data do not support this, as affected subjects displayed normal cognition. Large variability was observed in the age of onset (a)symmetry, severity and progression rate of HI. About half of the affected individuals displayed vestibular dysfunction and experienced symptoms thereof. The late-onset progressive phenotype and the absence of cochleovestibular malformations on computed tomography scans indicate that heterozygous defects of LMX1A do not result in severe developmental abnormalities in humans. We propose that a single LMX1A wild-type copy is sufficient for normal development but insufficient for maintenance of cochleovestibular function. Alternatively, minor cochleovestibular developmental abnormalities could eventually lead to the progressive phenotype seen in the families.

Fatal thoracic aortic aneurysm and dissection in a large family with a novel MYLK gene mutation: delineation of the clinical phenotype

Background

Thoracic and abdominal aortic aneurysms and dissection often develop in hypertensive elderly patients. At higher risk are smokers and those who have a family history of aortic aneurysms. In most affected families, the aortic aneurysms and dissection is inherited in an autosomal dominant manner with decreased penetrance and variable expressivity. Mutations at two chromosomal loci, TAA1 at 11q23 and the TAA2 at 5q13–14, and eight genes, MYLK, MYH11, TGFBR2, TGFBR1, ACTA2, SMAD3, TGFB2, and MAT2A, have been identified as being responsible for the disease in 23% of affected families.

Results

Herein, we inform on the clinical, genetic and pathological characteristics of nine living and deceased members of a large consanguineous Arab family with thoracic aortic aneurysm and dissection who carry a missense mutation c.4471G > T (Ala1491Ser), in exon 27 of MYLK gene. We show a reduced kinase activity of the Ala1491Ser protein compared to wildtype protein. This mutation is expressed as aortic aneurysm and dissection in one of two distinct phenotypes. A severe fatal and early onset symptom in homozygous or mild late onset in heterozygous genotypes.

Conclusions

We found that MYLK gene Ala1491Ser mutation affect the kinase activity and clinically, it presents with vascular aneurysms and dissection. We describe a distinct genotype phenotype correlation where; heterozygous patients have mild late onset and incomplete penetrance disease compared with the early onset severe and generally fatal outcome in homozygous patients.

Progressive hereditary spastic paraplegia caused by a homozygous KY mutation

Twelve individuals of consanguineous Bedouin kindred presented with autosomal recessive progressive spastic paraplegia evident as of age 0-24 months, with spasticity of lower limbs, hyperreflexia, toe walking and equinus deformity. Kyphoscolisois was evident in older patients. Most had atrophy of the lateral aspects of the tongue and few had intellectual disability. Nerve conduction velocity, electromyography and head and spinal cord magnetic resonance imaging were normal in tested subjects. Muscle biopsy showed occasional central nuclei and fiber size variability with small angular fibers. Genome-wide linkage analysis identified a 6.7Mbp disease-associated locus on chromosome 3q21.3-3q22.2 (LOD score 9.02; D3S1290). Whole-exome sequencing identified a single homozygous variant within this locus, c.51_52ins(28); p.(V18fs56*) in KY, segregating in the family as expected and not found in 190 Bedouin controls. High KY transcript levels were demonstrated in muscular organs with lower expression in the CNS. The phenotype is reminiscent of kyphoscoliosis seen in Ky null mice. Two recent studies done independently and parallel to ours describe somewhat similar phenotypes in one and two patients with KY mutations. KY encodes a tranglutaminase-like peptidase, which interacts with muscle cytoskeletal proteins and is part of a Z-band protein complex, suggesting the disease mechanism may resemble myofibrillar myopathy. However, the mixed myopathic-neurologic features caused by human and mouse Ky mutations are difficult to explain by loss of KY sarcomere stabilizing function alone. KY transcription in CNS tissues may imply that it also has a role in neuromotor function, in line with the irregularity of neuromuscular junction in Ky null mutant mice.

Family Based Whole Exome Sequencing Reveals the Multifaceted Role of Notch Signaling in Congenital Heart Disease

Left-ventricular outflow tract obstructions (LVOTO) encompass a wide spectrum of phenotypically heterogeneous heart malformations which frequently cluster in families. We performed family based whole-exome and targeted re-sequencing on 182 individuals from 51 families with multiple affected members. Central to our approach is the family unit which serves as a reference to identify causal genotype-phenotype correlations. Screening a multitude of 10 overlapping phenotypes revealed disease associated and co-segregating variants in 12 families. These rare or novel protein altering mutations cluster predominantly in genes (NOTCH1, ARHGAP31, MAML1, SMARCA4, JARID2, JAG1) along the Notch signaling cascade. This is in line with a significant enrichment (Wilcoxon, p< 0.05) of variants with a higher pathogenicity in the Notch signaling pathway in patients compared to controls. The significant enrichment of novel protein truncating and missense mutations in NOTCH1 highlights the allelic and phenotypic heterogeneity in our pediatric cohort. We identified novel co-segregating pathogenic mutations in NOTCH1 associated with left and right-sided cardiac malformations in three independent families with a total of 15 affected individuals. In summary, our results suggest that a small but highly pathogenic fraction of family specific mutations along the Notch cascade are a common cause of LVOTO.

Left-ventricular outflow tract obstructions comprise a group of developmental heart disorders that are genetically and phenotypically heterogeneous, with no single gene accounting for the majority of cases. In order to identify mutations contributing to disease, we selected patients from 51 families with a history of congenital cardiac malformations. We interrogated the entire coding sequences of 106 patients and identified a small but highly pathogenic fraction of mutations that are likely to contribute to disease in 12 families (23.5%). Furthermore, we present a strategy for identifying candidate mutations based on familial segregation in a genetically heterogeneous disorder.

Unexpected identification of a recurrent mutation in the DLX3 gene causing amelogenesis imperfecta

OBJECTIVE

To identify the molecular genetic aetiology of a family with autosomal dominant amelogenesis imperfecta (AI).

SUBJECTS AND METHODS

DNA samples were collected from a six-generation family, and the candidate gene approach was used to screen for the enamelin (ENAM) gene. Whole-exome sequencing and linkage analysis with SNP array data identified linked regions, and candidate gene screening was performed.

RESULTS

Mutational analysis revealed a mutation (c.561_562delCT and p.Tyr188Glnfs*13) in the DLX3 gene. After finding a recurrent DLX3 mutation, the clinical phenotype of the family members was re-examined. The proband's mother had pulp elongation in the third molars. The proband had not hair phenotype, but her cousin had curly hair at birth.

CONCLUSIONS

In this study, we identified a recurrent 2-bp deletional DLX3 mutation in a new family. The clinical phenotype was the mildest one associated with the DLX3 mutations. These results will advance the understanding of the functional role of DLX3 in developmental processes.

Increased Probability of Co-Occurrence of Two Rare Diseases in Consanguineous Families and Resolution of a Complex Phenotype by Next Generation Sequencing

Massively parallel sequencing of whole genomes and exomes has facilitated a direct assessment of causative genetic variation, now enabling the identification of genetic factors involved in rare diseases (RD) with Mendelian inheritance patterns on an almost routine basis. Here, we describe the illustrative case of a single consanguineous family where this strategy suffered from the difficulty to distinguish between two etiologically distinct disorders, namely the co-occurrence of hereditary hypophosphatemic rickets (HRR) and congenital myopathies (CM), by their phenotypic manifestation alone. We used parametric linkage analysis, homozygosity mapping and whole exome-sequencing to identify mutations underlying HRR and CM. We also present an approximate approach for assessing the probability of co-occurrence of two unlinked recessive RD in a single family as a function of the degree of consanguinity and the frequency of the disease-causing alleles. Linkage analysis and homozygosity mapping yielded elusive results when assuming a single RD, but whole-exome sequencing helped to identify two mutations in two genes, namely SLC34A3 and SEPN1, that segregated independently in this family and that have previously been linked to two etiologically different diseases. We assess the increase in chance co-occurrence of rare diseases due to consanguinity, i.e. under circumstances that generally favor linkage mapping of recessive disease, and show that this probability can increase by several orders of magnitudes. We conclude that such potential co-occurrence represents an underestimated risk when analyzing rare or undefined diseases in consanguineous families and should be given more consideration in the clinical and genetic evaluation.

PBAP: a pipeline for file processing and quality control of pedigree data with dense genetic markers

MOTIVATION

Huge genetic datasets with dense marker panels are now common. With the availability of sequence data and recognition of importance of rare variants, smaller studies based on pedigrees are again also common. Pedigree-based samples often start with a dense marker panel, a subset of which may be used for linkage analysis to reduce computational burden and to limit linkage disequilibrium between single-nucleotide polymorphisms (SNPs). Programs attempting to select markers for linkage panels exist but lack flexibility.

RESULTS

We developed a pedigree-based analysis pipeline (PBAP) suite of programs geared towards SNPs and sequence data. PBAP performs quality control, marker selection and file preparation. PBAP sets up files for MORGAN, which can handle analyses for small and large pedigrees, typically human, and results can be used with other programs and for downstream analyses. We evaluate and illustrate its features with two real datasets.

AVAILABILITY AND IMPLEMENTATION

PBAP scripts may be downloaded from http://faculty.washington.edu/wijsman/software.shtml.

CONTACT

wijsman@uw.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

UNC80 mutation causes a syndrome of hypotonia, severe intellectual disability, dyskinesia and dysmorphism, similar to that caused by mutations in its interacting cation channel NALCN

BACKGROUND

A syndrome of profound hypotonia, intellectual disability, intrauterine growth retardation with subsequent failure to thrive, dyskinesia and epilepsy was diagnosed in Bedouin Israeli families. Mild dysmorphism was evident: plagiocephaly, broad forehead with prominent nose, smooth philtrum and congenital esotropia. We set out to decipher the molecular basis of this syndrome.

METHODS

Genome-wide linkage analysis and fine mapping were done. Whole exome sequencing data were filtered for candidate variants within locus. Validation and segregation of the mutation was assayed via Sanger sequencing. UNC80 expression pattern was analysed through reverse transcription PCR.

RESULTS

Homozygosity mapping followed by fine mapping identified a 7.5 Mb disease-associated locus (logarithm of odds score 3.5) on chromosome 2. Whole exome and Sanger sequencing identified a single homozygous nonsense mutation within this locus, segregating within the families as expected for recessive heredity and not found in a homozygous state in 150 Bedouin controls: c.151C>T, p.(R51*) in UNC80.

CONCLUSIONS

The syndrome described is caused by a mutation in UNC80, truncating most of the 3258 amino acids highly conserved encoded protein, that has no known motifs. UNC80 bridges between UNC79 and the cation channel NALCN, enabling NALCN's role in basal Na(+) leak conductance in neurons, essential for neuronal function. The phenotype caused by the UNC80 mutation resembles that previously described for homozygous NALCN mutations.

Allelic Mutations of KITLG, Encoding KIT Ligand, Cause Asymmetric and Unilateral Hearing Loss and Waardenburg Syndrome Type 2

Linkage analysis combined with whole-exome sequencing in a large family with congenital and stable non-syndromic unilateral and asymmetric hearing loss (NS-UHL/AHL) revealed a heterozygous truncating mutation, c.286_303delinsT (p.Ser96Ter), in KITLG. This mutation co-segregated with NS-UHL/AHL as a dominant trait with reduced penetrance. By screening a panel of probands with NS-UHL/AHL, we found an additional mutation, c.200_202del (p.His67_Cys68delinsArg). In vitro studies revealed that the p.His67_Cys68delinsArg transmembrane isoform of KITLG is not detectable at the cell membrane, supporting pathogenicity. KITLG encodes a ligand for the KIT receptor. Also, KITLG-KIT signaling and MITF are suggested to mutually interact in melanocyte development. Because mutations in MITF are causative of Waardenburg syndrome type 2 (WS2), we screened KITLG in suspected WS2-affected probands. A heterozygous missense mutation, c.310C>G (p.Leu104Val), that segregated with WS2 was identified in a small family. In vitro studies revealed that the p.Leu104Val transmembrane isoform of KITLG is located at the cell membrane, as is wild-type KITLG. However, in culture media of transfected cells, the p.Leu104Val soluble isoform of KITLG was reduced, and no soluble p.His67_Cys68delinsArg and p.Ser96Ter KITLG could be detected. These data suggest that mutations in KITLG associated with NS-UHL/AHL have a loss-of-function effect. We speculate that the mechanism of the mutation underlying WS2 and leading to membrane incorporation and reduced secretion of KITLG occurs via a dominant-negative or gain-of-function effect. Our study unveils different phenotypes associated with KITLG, previously associated with pigmentation abnormalities, and will thereby improve the genetic counseling given to individuals with KITLG variants.

Variants in Nebulin (NEB) Are Linked to the Development of Familial Primary Angle Closure Glaucoma in Basset Hounds

Several dog breeds are susceptible to developing primary angle closure glaucoma (PACG), which suggests a genetic basis for the disease. We have identified a four-generation Basset Hound pedigree with characteristic autosomal recessive PACG that closely recapitulates PACG in humans. Our aim is to utilize gene mapping and whole exome sequencing approaches to identify PACG-causing sequence variants in the Basset. Extensive clinical phenotyping of all pedigree members was conducted. SNP-chip genotyping was carried out in 9 affected and 15 unaffected pedigree members. Two-point and multipoint linkage analyses of genome-wide SNP data were performed using Superlink-Online SNP-1.1 and a locus was mapped to chromosome 19q with a maximum LOD score of 3.24. The locus contains 12 Ensemble predicted canine genes and is syntenic to a region on chromosome 2 in the human genome. Using exome-sequencing analysis, a possibly damaging, non-synonymous variant in the gene Nebulin (NEB) was found to segregate with PACG which alters a phylogenetically conserved Lysine residue. The association of this variants with PACG was confirmed in a secondary cohort of unrelated Basset Hounds (p = 3.4 × 10^-4, OR = 15.3 for homozygosity). Nebulin, a protein that promotes the contractile function of sarcomeres, was found to be prominently expressed in the ciliary muscles of the anterior segment. Our findings may provide insight into the molecular mechanisms that underlie PACG. The phenotypic similarities of disease presentation in dogs and humans may enable the translation of findings made in this study to patients with PACG.

A novel variant in CDKN1C is associated with intrauterine growth restriction, short stature, and early-adulthood-onset diabetes

CONTEXT

CDKN1C, a cyclin-dependent kinase inhibitor and negative regulator of cellular proliferation, is paternally imprinted and has been shown to regulate β-cell proliferation. CDKN1C mutations are associated with growth disorders, including Beckwith-Wiedemann syndrome and IMAGe syndrome.

OBJECTIVE

To investigate the genetic basis for a familial disorder characterized by intrauterine growth restriction, short stature, and early-adulthood-onset diabetes.

DESIGN, SETTING, AND PARTICIPANTS

Genomic DNA samples (15 affected and 26 unaffected from a six-generation pedigree) were analyzed by genome-wide single nucleotide polymorphism arrays, whole exome and Sanger sequencing, and multiplex ligation-dependent probe amplification.

MAIN OUTCOME MEASURE(S)

Subjects were assessed for height, weight, adrenal gland size, ACTH, diabetes status, and testis volume. Linkage and sequence analyses were performed, and the identified genetic variant was functionally evaluated in reconstitution studies.

RESULTS

The pedigree followed a paternally imprinted pattern of inheritance, and genetic linkage analysis identified a single significant 2.6-megabase locus on chromosome 11p15, within the imprinting center region 2. Multiplex ligation-dependent probe amplification did not detect copy number variants or methylation abnormalities. Whole exome sequencing revealed a single novel variant in the proliferating cell nuclear antigen-binding region of CDKN1C (c.842G>T, p.R281I) that co-segregated with affected status and, unlike variants found in IMAGe, did not entirely abrogate proliferating cell nuclear antigen binding. Clinical assessments revealed that affected individuals had low testicular volume but normal adrenal function.

CONCLUSIONS

We report a novel CDKN1C mutation associated with features of IMAGe syndrome, but without adrenal insufficiency or metaphyseal dysplasia, and characterized by early-adulthood-onset diabetes. Our data expand the range of phenotypes observed with CDKN1C defects and suggest that CDKN1C mutations may represent a novel monogenic form of diabetes.