Expansion of GGC Repeat in GIPC1 Is Associated with Oculopharyngodistal Myopathy.

Oculopharyngodistal myopathy (OPDM) is an adult-onset inherited neuromuscular disorder characterized by progressive ptosis, external ophthalmoplegia, and weakness of the masseter, facial, pharyngeal, and distal limb muscles. The myopathological features are presence of rimmed vacuoles (RVs) in the muscle fibers and myopathic changes of differing severity. Inheritance is variable, with either putative autosomal-dominant or autosomal-recessive pattern. Here, using a comprehensive strategy combining whole-genome sequencing (WGS), long-read whole-genome sequencing (LRS), linkage analysis, repeat-primed polymerase chain reaction (RP-PCR), and fluorescence amplicon length analysis polymerase chain reaction (AL-PCR), we identified an abnormal GGC repeat expansion in the 5' UTR of GIPC1 in one out of four families and three sporadic case subjects from a Chinese OPDM cohort. Expanded GGC repeats were further confirmed as the cause of OPDM in an additional 2 out of 4 families and 6 out of 13 sporadic Chinese individuals with OPDM, as well as 7 out of 194 unrelated Japanese individuals with OPDM. Methylation, qRT-PCR, and western blot analysis indicated that GIPC1 mRNA levels were increased while protein levels were unaltered in OPDM-affected individuals. RNA sequencing indicated p53 signaling, vascular smooth muscle contraction, ubiquitin-mediated proteolysis, and ribosome pathways were involved in the pathogenic mechanisms of OPDM-affected individuals with GGC repeat expansion in GIPC1. This study provides further evidence that OPDM is associated with GGC repeat expansions in distinct genes and highly suggests that expanded GGC repeat units are essential in the pathogenesis of OPDM, regardless of the genes in which the expanded repeats are located.

Pinpointing genomic regions associated with root system architecture in rice through an integrative meta-analysis approach.

KEY MESSAGE: Applying an integrated meta-analysis approach led to identification of meta-QTLs/ candidate genes associated with rice root system architecture, which can be used in MQTL-assisted breeding/ genetic engineering of root traits. Root system architecture (RSA) is an important factor for facilitating water and nutrient uptake from deep soils and adaptation to drought stress conditions. In the present research, an integrated meta-analysis approach was employed to find candidate genes and genomic regions involved in rice RSA traits. A whole-genome meta-analysis was performed for 425 initial QTLs reported in 34 independent experiments controlling RSA traits under control and drought stress conditions in the previous twenty years. Sixty-four consensus meta-QTLs (MQTLs) were detected, unevenly distributed on twelve rice chromosomes. The confidence interval (CI) of the identified MQTLs was obtained as 0.11-14.23 cM with an average of 3.79 cM, which was 3.88 times narrower than the mean CI of the original QTLs. Interestingly, 52 MQTLs were co-located with SNP peak positions reported in rice genome-wide association studies (GWAS) for root morphological traits. The genes located in these RSA-related MQTLs were detected and explored to find the drought-responsive genes in the rice root based on the RNA-seq and microarray data. Multiple RSA and drought tolerance-associated genes were found in the MQTLs including the genes involved in auxin biosynthesis or signaling (e.g. YUCCA, WOX, AUX/IAA, ARF), root angle (DRO1-related genes), lateral root development (e.g. DSR, WRKY), root diameter (e.g. OsNAC5), plant cell wall (e.g. EXPA), and lignification (e.g. C4H, PAL, PRX and CAD). The genes located within both the SNP peak positions and the QTL-overview peaks for RSA are suggested as novel candidate genes for further functional analysis. The promising candidate genes and MQTLs can be used as basis for genetic engineering and MQTL-assisted breeding of root phenotypes to improve yield potential, stability and performance in a water-stressed environment.

Methodology and theoretical basis of forward genetic screening for sleep/wakefulness in mice.

The regulatory network of genes and molecules in sleep/wakefulness remains to be elucidated. Here we describe the methodology and workflow of the dominant screening of randomly mutagenized mice and discuss theoretical basis of forward genetics research for sleep in mice. Our high-throughput screening employs electroencephalogram (EEG) and electromyogram (EMG) to stage vigilance states into a wake, rapid eye movement sleep (REMS) and non-REM sleep (NREMS). Based on their near-identical sleep/wake behavior, C57BL/6J (B6J) and C57BL/6N (B6N) are chosen as mutagenized and counter strains, respectively. The total time spent in the wake and NREMS, as well as the REMS episode duration, shows sufficient reproducibility with small coefficients of variance, indicating that these parameters are most suitable for quantitative phenotype-driven screening. Coarse linkage analysis of the quantitative trait, combined with whole-exome sequencing, can identify the gene mutation associated with sleep abnormality. Our simulations calculate the achievable LOD score as a function of the phenotype strength and the numbers of mice examined. A pedigree showing a mild decrease in total wake time resulting from a heterozygous point mutation in the Cacna1a gene is described as an example.

Whole-genome sequencing identifies functional noncoding variation in SEMA3C that cosegregates with dyslexia in a multigenerational family.

Dyslexia is a common heritable developmental disorder involving impaired reading abilities. Its genetic underpinnings are thought to be complex and heterogeneous, involving common and rare genetic variation. Multigenerational families segregating apparent monogenic forms of language-related disorders can provide useful entrypoints into biological pathways. In the present study, we performed a genome-wide linkage scan in a three-generational family in which dyslexia affects 14 of its 30 members and seems to be transmitted with an autosomal dominant pattern of inheritance. We identified a locus on chromosome 7q21.11 which cosegregated with dyslexia status, with the exception of two cases of phenocopy (LOD = 2.83). Whole-genome sequencing of key individuals enabled the assessment of coding and noncoding variation in the family. Two rare single-nucleotide variants (rs144517871 and rs143835534) within the first intron of the SEMA3C gene cosegregated with the 7q21.11 risk haplotype. In silico characterization of these two variants predicted effects on gene regulation, which we functionally validated for rs144517871 in human cell lines using luciferase reporter assays. SEMA3C encodes a secreted protein that acts as a guidance cue in several processes, including cortical neuronal migration and cellular polarization. We hypothesize that these intronic variants could have a cis-regulatory effect on SEMA3C expression, making a contribution to dyslexia susceptibility in this family.

Bi-allelic Loss-of-Function Variants in DNMBP Cause Infantile Cataracts.

Infantile and childhood-onset cataracts form a heterogeneous group of disorders; among the many genetic causes, numerous pathogenic variants in additional genes associated with autosomal-recessive infantile cataracts remain to be discovered. We identified three consanguineous families affected by bilateral infantile cataracts. Using exome sequencing, we found homozygous loss-of-function variants in DNMBP: nonsense variant c.811C>T (p.Arg271∗) in large family F385 (nine affected individuals; LOD score = 5.18 at θ = 0), frameshift deletion c.2947_2948del (p.Asp983∗) in family F372 (two affected individuals), and frameshift variant c.2852_2855del (p.Thr951Metfs∗41) in family F3 (one affected individual). The phenotypes of all affected individuals include infantile-onset cataracts. RNAi-mediated knockdown of the Drosophila ortholog still life (sif), enriched in lens-secreting cells, affects the development of these cells as well as the localization of E-cadherin, alters the distribution of septate junctions in adjacent cone cells, and leads to a ∼50% reduction in electroretinography amplitudes in young flies. DNMBP regulates the shape of tight junctions, which correspond to the septate junctions in invertebrates, as well as the assembly pattern of E-cadherin in human epithelial cells. E-cadherin has an important role in lens vesicle separation and lens epithelial cell survival in humans. We therefore conclude that DNMBP loss-of-function variants cause infantile-onset cataracts in humans.

Genetic linkage analysis in the age of whole-genome sequencing.

For many years, linkage analysis was the primary tool used for the genetic mapping of Mendelian and complex traits with familial aggregation. Linkage analysis was largely supplanted by the wide adoption of genome-wide association studies (GWASs). However, with the recent increased use of whole-genome sequencing (WGS), linkage analysis is again emerging as an important and powerful analysis method for the identification of genes involved in disease aetiology, often in conjunction with WGS filtering approaches. Here, we review the principles of linkage analysis and provide practical guidelines for carrying out linkage studies using WGS data.

Genome-Wide Linkage Study Meta-Analysis of Male Sexual Orientation.

Male sexual orientation is a scientifically and socially important trait shown by family and twin studies to be influenced by environmental and complex genetic factors. Individual genome-wide linkage studies (GWLS) have been conducted, but not jointly analyzed. Two main datasets account for > 90% of the published GWLS concordant sibling pairs on the trait and are jointly analyzed here: MGSOSO (Molecular Genetic Study of Sexual Orientation; 409 concordant sibling pairs in 384 families, Sanders et al. (2015)) and Hamer (155 concordant sibling pairs in 145 families, Mustanski et al. (2005)). We conducted multipoint linkage analyses with Merlin on the datasets separately since they were genotyped differently, integrated genetic marker positions, and combined the resultant LOD (logarithm of the odds) scores at each 1 cM grid position. We continue to find the strongest linkage support at pericentromeric chromosome 8 and chromosome Xq28. We also incorporated the remaining published GWLS dataset (on 55 families) by using meta-analytic approaches on published summary statistics. The meta-analysis has maximized the positional information from GWLS of currently available family resources and can help prioritize findings from genome-wide association studies (GWAS) and other approaches. Although increasing evidence highlights genetic contributions to male sexual orientation, our current understanding of contributory loci is still limited, consistent with the complexity of the trait. Further increasing genetic knowledge about male sexual orientation, especially via large GWAS, should help advance our understanding of the biology of this important trait.

Introgression and QTL mapping conferring resistance for Alternaria brassicae in the backcross progeny of Sinapis alba + Brassica juncea somatic hybrids.

KEY MESSAGE: A total of three QTLs, responsible for A. brassicae resistance were introgressed into S. alba - B. juncea introgression lines from S. alba and mapped through donor genome-specific SSR markers. Alternaria brassicae is a key pathogen of the Brassicaceae family causing severe blight disease to oilseed crops that leads to heavy yield losses due to lack of resistance source within cultivated Brassicas. However, the host resistance present in the Sinapis alba, an allied member of the Brassicaceae family is still unexplored precisely due to the unavailability of segregating population for Alternaria blight resistance and scarcity of donor genome-specific genetic markers. The present study was undertaken to identify quantitative trait loci governing resistance to Alternaria blight which was introgressed from S. alba to the backcross population of stable S. alba + B. juncea somatic hybrids (2n = 60; AABBSS). The second backcross population showed significant phenotypic variations for Alternaria blight ranging from immune to highly susceptible phenotype, thus suggesting quantitative nature of resistance for the disease. A subset of 154 BC2F3-4 lines was used for disease screening and genotyping with 234 S. alba genome-specific microsatellite markers. As a result of the study, twelve linkage groups were developed corresponding to 12 chromosomes of S. alba (n = 12) covering a length of 1694.02 cM. The chromosomes 5 and 11 harbored a total of 1 (Abr-01), and 2 (Abr-02, and Abr-03) QTLs detected by ICIM-ADD mapping method at LOD score values 3.7, 5.12, and 6.74, respectively. The QTLs identified during the study have a range of 5.51-10.87 percent phenotypic variations for disease resistance. To the best of our knowledge, this is the first report of QTLs introgression for A. brassicae resistance in cultivated Brassica from an allied member of Brassicaceae.

Characterization of genetic determinants of the resistance to phylloxera, Daktulosphaira vitifoliae, and the dagger nematode Xiphinema index from muscadine background.

BACKGROUND: Muscadine (Muscadinia rotundifolia) is known as a resistance source to many pests and diseases in grapevine. The genetics of its resistance to two major grapevine pests, the phylloxera D. vitifoliae and the dagger nematode X. index, vector of the Grapevine fanleaf virus (GFLV), was investigated in a backcross progeny between the F1 resistant hybrid material VRH8771 (Vitis-Muscadinia) derived from the muscadine R source 'NC184-4' and V. vinifera cv. 'Cabernet-Sauvignon' (CS). RESULTS: In this pseudo-testcross, parental maps were constructed using simple-sequence repeats markers and single nucleotide polymorphism markers from a GBS approach. For the VRH8771 map, 2271 SNP and 135 SSR markers were assembled, resulting in 19 linkage groups (LG) and an average distance between markers of 0.98 cM. Phylloxera resistance was assessed by monitoring root nodosity number in an in planta experiment and larval development in a root in vitro assay. Nematode resistance was studied using 10-12 month long tests for the selection of durable resistance and rating criteria based on nematode reproduction factor and gall index. A major QTL for phylloxera larval development, explaining more than 70% of the total variance and co-localizing with a QTL for nodosity number, was identified on LG 7 and designated RDV6. Additional QTLs were detected on LG 3 (RDV7) and LG 10 (RDV8), depending on the in planta or in vitro experiments, suggesting that various loci may influence or modulate nodosity formation and larval development. Using a Bulked Segregant Analysis approach and a proportion test, markers clustered in three regions on LG 9, LG 10 and LG 18 were shown to be associated to the nematode resistant phenotype. QTL analysis confirmed the results and QTLs were thus designated respectively XiR2, XiR3 and XiR4, although a LOD-score below the significant threshold value was obtained for the QTL on LG 18. CONCLUSIONS: Based on a high-resolution linkage map and a segregating grapevine backcross progeny, the first QTLs for resistance to D. vitifoliae and to X. index were identified from a muscadine source. All together these results open the way to the development of marker-assisted selection in grapevine rootstock breeding programs based on muscadine derived resistance to phylloxera and to X. index in order to delay GFLV transmission.

Grain protein content and thousand kernel weight QTLs identified in a durum × wild emmer wheat mapping population tested in five environments.

KEY MESSAGE: Genetic dissection of GPC and TKW in tetraploid durum × WEW RIL population, based on high-density SNP genetic map, revealed 12 GPC QTLs and 11 TKW QTLs, with favorable alleles for 11 and 5 QTLs, respectively, derived from WEW. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW) was shown to exhibit high grain protein content (GPC) and therefore possess a great potential for improvement of cultivated wheat nutritional value. Genetic dissection of thousand kernel weight (TKW) and grain protein content (GPC) was performed using a high-density genetic map constructed based on a recombinant inbred line (RIL) population derived from a cross between T. durum var. Svevo and WEW acc. Y12-3. Genotyping of 208 F6 RILs with a 15 K wheat single nucleotide polymorphism (SNP) array yielded 4166 polymorphic SNP markers, of which 1510 were designated as skeleton markers. A total map length of 2169 cM was obtained with an average distance of 1.5 cM between SNPs. A total of 12 GPC QTLs and 11 TKW QTLs were found under five different environments. No significant correlations were found between GPC and TKW across all environments. Four major GPC QTLs with favorable alleles from WEW were found on chromosomes 4BS, 5AS, 6BS and 7BL. The 6BS GPC QTL coincided with the physical position of the NAC transcription factor TtNAM-B1, underlying the cloned QTL, Gpc-B1. Comparisons of the physical intervals of the GPC QTLs described here with the results previously reported in other durum × WEW RIL population led to the discovery of seven novel GPC QTLs. Therefore, our research emphasizes the importance of GPC QTL dissection in diverse WEW accessions as a source of novel alleles for improvement of GPC in cultivated wheat.

Aphid resistance in Capsicum maps to a locus containing LRR-RLK gene analogues.

KEY MESSAGE: A QTL for aphid resistance on pepper chromosome 2 was identified and validated. This QTL affects aphid survival and reproduction, and was fine mapped to a locus containing LRR-RLK analogues. Myzus persicae is one of the most threatening insect pests that adversely affects pepper (Capsicum) cultivation. Resistance to aphids was previously identified in Capsicum baccatum. This study aimed at elucidating the genetics of aphid resistance in C. baccatum. A QTL analysis was carried out for M. persicae resistance in an F2 population derived from an intraspecific cross between a highly resistant plant and a susceptible plant. Survival and reproduction were used as resistance parameters. Interval mapping detected two QTLs affecting aphid survival (Rmpas-1) and reproduction (Rmprp-1), respectively, both localized in the same area and sharing the same top marker on chromosome 2. Use of this marker as co-factor in multiple-QTL mapping analysis revealed a second, minor QTL (Rmprp-2) only affecting aphid reproduction, on chromosome 4. Fine mapping confirmed the effects of Rmpas-1 and Rmprp-1 and narrowed the major QTL Rmprp-1 down to a genomic region of 96 kb which is predicted to encode four analogues of resistance genes of the receptor-like kinase family containing a leucine-rich repeat domain (LRR-RLKs). This work provides not only initial information for breeding aphid-resistant pepper varieties, but also forms the basis for future molecular analysis of gene(s) involved in aphid resistance.

Properties and Evaluation of the MOBIT - a novel Linkage-based Test Statistic and Quantification Method for Imprinting.

Genomic imprinting is a parent-of-origin effect apparent in an appreciable number of human diseases. We have proposed the new imprinting test statistic MOBIT, which is based on MOD score analysis. We were interested in the properties of the MOBIT concerning its distribution under three hypotheses: (1) H0,a: no linkage, no imprinting; (2) H0,b: linkage, no imprinting; (3) H1: linkage and imprinting. More specifically, we assessed the confounding between imprinting and sex-specific recombination frequencies, which presents a major difficulty in linkage-based testing for imprinting, and evaluated the power of the test. To this end, we have performed a linkage simulation study of affected sib-pairs and a three-generation pedigree with two trait models, many two- and multipoint marker scenarios, three genetic map ratios, two sample sizes, and five imprinting degrees. We also investigated the ability of the MOBIT to quantify the degree of imprinting and applied the MOBIT using a real data example on house dust mite allergy. We further proposed and evaluated two approaches to obtain empiric p values for the MOBIT. Our results showed that twopoint analyses assuming a sex-averaged marker map led to an inflated type I error due to confounding, especially for a larger marker-trait locus distance. When the correct sex-specific marker map was assumed, twopoint analyses have a reduced power to detect imprinting, compared to sex-averaged analyses with an appropriate correction for the inflation of the test statistic. However, confounding was not an issue in multipoint analysis unless the map ratio was extreme and marker spacing was sparse. With multipoint analysis, power as well as the ability to quantify the imprinting degree were almost equally high when a sex-averaged or the correct sex-specific map was used in the analysis. We recommend to obtain empiric p values for the MOBIT using genotype simulations based on the best-fitting nonimprinting model of the real dataset analysis. In addition, an implementation of a method based on the permutation of parental sexes is also available. In summary, we propose to perform multipoint analyses using densely spaced markers to efficiently discover new imprinted loci and to reliably quantify the degree of imprinting.

NAA10 polyadenylation signal variants cause syndromic microphthalmia.

BACKGROUND: A single variant in NAA10 (c.471+2T>A), the gene encoding N-acetyltransferase 10, has been associated with Lenz microphthalmia syndrome. In this study, we aimed to identify causative variants in families with syndromic X-linked microphthalmia. METHODS: Three families, including 15 affected individuals with syndromic X-linked microphthalmia, underwent analyses including linkage analysis, exome sequencing and targeted gene sequencing. The consequences of two identified variants in NAA10 were evaluated using quantitative PCR and RNAseq. RESULTS: Genetic linkage analysis in family 1 supported a candidate region on Xq27-q28, which included NAA10. Exome sequencing identified a hemizygous NAA10 polyadenylation signal (PAS) variant, chrX:153,195,397T>C, c.*43A>G, which segregated with the disease. Targeted sequencing of affected males from families 2 and 3 identified distinct NAA10 PAS variants, chrX:g.153,195,401T>C, c.*39A>G and chrX:g.153,195,400T>C, c.*40A>G. All three variants were absent from gnomAD. Quantitative PCR and RNAseq showed reduced NAA10 mRNA levels and abnormal 3' UTRs in affected individuals. Targeted sequencing of NAA10 in 376 additional affected individuals failed to identify variants in the PAS. CONCLUSION: These data show that PAS variants are the most common variant type in NAA10-associated syndromic microphthalmia, suggesting reduced RNA is the molecular mechanism by which these alterations cause microphthalmia/anophthalmia. We reviewed recognised variants in PAS associated with Mendelian disorders and identified only 23 others, indicating that NAA10 harbours more than 10% of all known PAS variants. We hypothesise that PAS in other genes harbour unrecognised pathogenic variants associated with Mendelian disorders. The systematic interrogation of PAS could improve genetic testing yields.

Genetic Dissection of Phomopsis Stem Canker Resistance in Cultivated Sunflower Using High Density SNP Linkage Map.

Phomopsis stem canker (PSC) caused by Diaporthe helianthi is increasingly becoming a global threat for sunflower production. In this study, the genetic basis of PSC resistance was investigated in a recombinant inbred line (RIL) population developed from a cross between HA 89 (susceptible) and HA-R3 (resistant). The RIL population was evaluated for PSC disease incidence (DI) in seven screening trials at multiple locations during 2016-2018. The distribution of PSC DI in the RIL population was continuous, confirming a polygenic inheritance of the trait. A moderately high broad-sense heritability (H2, 0.76) was estimated for the trait across environments. In the combined analysis, both the genotype and the genotype × environment interactions were highly significant. A linkage map spanning 1505.33 cM was constructed using genotyping-by-sequencing derived markers. Marker-trait association analysis identified a total of 15 quantitative trait loci (QTL) associated with PSC resistance on 11 sunflower chromosomes, each explaining between 5.24 and 17.39% of the phenotypic variation. PSC resistance QTL were detected in two genomic regions each on chromosomes 3, 5, 13, and 17, while one QTL each was detected in the remaining seven chromosomes. Tightly linked single nucleotide polymorphism (SNP) markers flanking the PSC resistance QTL will facilitate marker-assisted selection in PSC resistance sunflower breeding.

A Novel Locus and Candidate Gene for Familial Developmental Dyslexia on Chromosome 4q.

Objective: Developmental dyslexia is a highly heritable specific reading and writing disability. To identify a possible new locus and candidate gene for this disability, we investigated a four-generation pedigree where transmission of dyslexia is consistent with an autosomal dominant inheritance pattern. Methods: We performed genome wide array-based SNP genotyping and parametric linkage analysis and sequencing analysis of protein-coding exons, exon-intron boundaries and conserved extragenic regions within the haplotype cosegregating with dyslexia in DNA from one affected and one unaffected family member. Cosegregation was confirmed by sequencing all available family members. Additionally, we analyzed 96 dyslexic individuals who had previously shown positive LOD scores on chromosome 4q28 as well as an even larger sample (n = 2591). Results: We found a single prominent linkage interval on chromosome 4q, where sequence analysis revealed a nucleotide variant in the 3' UTR of brain expressed SPRY1 in the dyslexic family member that cosegregated with dyslexia. This sequence alteration might affect the binding efficiency of the IGF2BP1 RNA-binding protein and thus influence the expression level of the SPRY1 gene product. An analysis of 96 individuals from a cohort of dyslexic individuals revealed a second heterozygous variant in this gene, which was absent in the unaffected sister of the proband. An investigation of the region in a much larger sample further found a nominal p-value of 0.0016 for verbal short-term memory (digit span) in 2,591 individuals for a neighboring SNV. After correcting for the local number of analyzed SNVs, and after taking into account linkage disequilibrium, we found this corresponds to a p-value of 0.0678 for this phenotype. Conclusions: We describe a new locus for familial dyslexia and discuss the possibility that SPRY1 might play a role in the etiology of a monogenic form of dyslexia.

Analysis of pedigree data in populations with multiple ancestries: Strategies for dealing with admixture in Caribbean Hispanic families from the ADSP.

Multipoint linkage analysis is an important approach for localizing disease-associated loci in pedigrees. Linkage analysis, however, is sensitive to misspecification of marker allele frequencies. Pedigrees from recently admixed populations are particularly susceptible to this problem because of the challenge of accurately accounting for population structure. Therefore, increasing emphasis on use of multiethnic samples in genetic studies requires reevaluation of best practices, given data currently available. Typical strategies have been to compute allele frequencies from the sample, or to use marker allele frequencies determined by admixture proportions averaged over the entire sample. However, admixture proportions vary among pedigrees and throughout the genome in a family-specific manner. Here, we evaluate several approaches to model admixture in linkage analysis, providing different levels of detail about ancestral origin. To perform our evaluations, for specification of marker allele frequencies, we used data on 67 Caribbean Hispanic admixed families from the Alzheimer's Disease Sequencing Project. Our results show that choice of admixture model has an effect on the linkage analysis results. Variant-specific admixture proportions, computed for individual families, provide the most detailed regional admixture estimates, and, as such, are the most appropriate allele frequencies for linkage analysis. This likely decreases the number of false-positive results, and is straightforward to implement.

Linkage analysis and whole exome sequencing reveals AHNAK2 as a novel genetic cause for autosomal recessive CMT in a Malaysian family.

Charcot-Marie-Tooth (CMT) disease is a form of inherited peripheral neuropathy that affects motor and sensory neurons. To identify the causative gene in a consanguineous family with autosomal recessive CMT (AR-CMT), we employed a combination of linkage analysis and whole exome sequencing. After excluding known AR-CMT genes, genome-wide linkage analysis mapped the disease locus to a 7.48-Mb interval on chromosome 14q32.11-q32.33, flanked by the markers rs2124843 and rs4983409. Whole exome sequencing identified two non-synonymous variants (p.T40P and p.H915Y) in the AHNAK2 gene that segregated with the disease in the family. Pathogenic predictions indicated that p.T40P is the likely causative allele. Analysis of AHNAK2 expression in the AR-CMT patient fibroblasts showed significantly reduced mRNA and protein levels. AHNAK2 binds directly to periaxin which is encoded by the PRX gene, and PRX mutations are associated with another form of AR-CMT (CMT4F). The altered expression of mutant AHNAK2 may disrupt the AHNAK2-PRX interaction in which one of its known functions is to regulate myelination.

Novel truncation mutations in MYRF cause autosomal dominant high hyperopia mapped to 11p12-q13.3.

High hyperopia is a common and severe form of refractive error. Genetic factors play important roles in the development of high hyperopia but the exact gene responsible for this condition is mostly unknown. We identified a large Chinese family with autosomal dominant high hyperopia. A genome-wide linkage scan mapped the high hyperopia to chromosome 11p12-q13.3, with maximum log of the odds scores of 4.68 at theta = 0 for D11S987. Parallel whole-exome sequencing detected a novel c.3377delG (p.Gly1126Valfs*31) heterozygous mutation in the MYRF gene within the linkage interval. Whole-exome sequencing in other 121 probands with high hyperopia identified additional novel mutations in MYRF within two other families: a de novo c.3274_3275delAG (p.Leu1093Profs*22) heterozygous mutation and a c.3194+2T>C heterozygous mutation. All three mutations are located in the C-terminal region of MYRF and are predicted to result in truncation of that portion. Two patients from two of the three families developed angle-closure glaucoma. These three mutations were present in neither the ExAC database nor our in-house whole-exome sequencing data from 3280 individuals. No other truncation mutations in MYRF were detected in the 3280 individuals. Knockdown of myrf resulted in small eye size in zebrafish. These evidence all support that truncation mutations in the C-terminal region of MYRF are responsible for autosomal dominant high hyperopia in these families. Our results may provide useful clues for further understanding the functional role of the C-terminal region of this critical myelin regulatory factor, as well as the molecular pathogenesis of high hyperopia and its associated angle-closure glaucoma.

Autosomal-Recessive Hearing Impairment Due to Rare Missense Variants within S1PR2.

The sphingosine-1-phosphate receptors (S1PRs) are a well-studied class of transmembrane G protein-coupled sphingolipid receptors that mediate multiple cellular processes. However, S1PRs have not been previously reported to be involved in the genetic etiology of human traits. S1PR2 lies within the autosomal-recessive nonsyndromic hearing impairment (ARNSHI) locus DFNB68 on 19p13.2. From exome sequence data we identified two pathogenic S1PR2 variants, c.323G>C (p.Arg108Pro) and c.419A>G (p.Tyr140Cys). Each of these variants co-segregates with congenital profound hearing impairment in consanguineous Pakistani families with maximum LOD scores of 6.4 for family DEM4154 and 3.3 for family PKDF1400. Neither S1PR2 missense variant was reported among ∼120,000 chromosomes in the Exome Aggregation Consortium database, in 76 unrelated Pakistani exomes, or in 720 Pakistani control chromosomes. Both DNA variants affect highly conserved residues of S1PR2 and are predicted to be damaging by multiple bioinformatics tools. Molecular modeling predicts that these variants affect binding of sphingosine-1-phosphate (p.Arg108Pro) and G protein docking (p.Tyr140Cys). In the previously reported S1pr2(-/-) mice, stria vascularis abnormalities, organ of Corti degeneration, and profound hearing loss were observed. Additionally, hair cell defects were seen in both knockout mice and morphant zebrafish. Family PKDF1400 presents with ARNSHI, which is consistent with the lack of gross malformations in S1pr2(-/-) mice, whereas family DEM4154 has lower limb malformations in addition to hearing loss. Our findings suggest the possibility of developing therapies against hair cell damage (e.g., from ototoxic drugs) through targeted stimulation of S1PR2.

Exome sequencing and CRISPR/Cas genome editing identify mutations of ZAK as a cause of limb defects in humans and mice.

The CRISPR/Cas technology enables targeted genome editing and the rapid generation of transgenic animal models for the study of human genetic disorders. Here we describe an autosomal recessive human disease in two unrelated families characterized by a split-foot defect, nail abnormalities of the hands, and hearing loss, due to mutations disrupting the SAM domain of the protein kinase ZAK. ZAK is a member of the MAPKKK family with no known role in limb development. We show that Zak is expressed in the developing limbs and that a CRISPR/Cas-mediated knockout of the two Zak isoforms is embryonically lethal in mice. In contrast, a deletion of the SAM domain induces a complex hindlimb defect associated with down-regulation of Trp63, a known split-hand/split-foot malformation disease gene. Our results identify ZAK as a key player in mammalian limb patterning and demonstrate the rapid utility of CRISPR/Cas genome editing to assign causality to human mutations in the mouse in <10 wk.