Zoran Brkanac

Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer's disease

Genome-wide association studies (GWAS) have identified several risk variants for late-onset Alzheimers disease (LOAD). These common variants have replicable but small effects on LOAD risk and generally do not have obvious functional effects. Low-frequency coding variants, not detected by GWAS, are predicted to include functional variants with larger effects on risk. To identify low-frequency coding variants with large effects on LOAD risk, we carried out whole-exome sequencing (WES) in 14 large LOAD families and follow-up analyses of the candidate variants in several large LOAD case–control data sets. A rare variant in PLD3 (phospholipase D3; Val232Met) segregated with disease status in two independent families and doubled risk for Alzheimer’s disease in seven independent case–control series with a total of more than 11,000 cases and controls of European descent. Gene-based burden analyses in 4,387 cases and controls of European descent and 302 African American cases and controls, with complete sequence data for PLD3, reveal that several variants in this gene increase risk for Alzheimer’s disease in both populations. PLD3 is highly expressed in brain regions that are vulnerable to Alzheimer’s disease pathology, including hippocampus and cortex, and is expressed at significantly lower levels in neurons from Alzheimer’s disease brains compared to control brains. Overexpression of PLD3 leads to a significant decrease in intracellular amyloid-β precursor protein (APP) and extracellular Aβ42 and Aβ40 (the 42- and 40-residue isoforms of the amyloid-β peptide), and knockdown of PLD3 leads to a significant increase in extracellular Aβ42 and Aβ40. Together, our genetic and functional data indicate that carriers of PLD3 coding variants have a twofold increased risk for LOAD and that PLD3 influences APP processing. This study provides an example of how densely affected families may help to identify rare variants with large effects on risk for disease or other complex traits.

Missense Mutations in the Regulatory Domain of PKCγ: A New Mechanism for Dominant Nonepisodic Cerebellar Ataxia

We report a nonepisodic autosomal dominant (AD) spinocerebellar ataxia (SCA) not caused by a nucleotide repeat expansion that is, to our knowledge, the first such SCA. The AD SCAs currently comprise a group of > or =16 genetically distinct neurodegenerative conditions, all characterized by progressive incoordination of gait and limbs and by speech and eye-movement disturbances. Six of the nine SCAs for which the genes are known result from CAG expansions that encode polyglutamine tracts. Noncoding CAG, CTG, and ATTCT expansions are responsible for three other SCAs. Approximately 30% of families with SCA do not have linkage to the known loci. We recently mapped the locus for an AD SCA in a family (AT08) to chromosome 19q13.4-qter. A particularly compelling candidate gene, PRKCG, encodes protein kinase C gamma (PKC gamma), a member of a family of serine/threonine kinases. The entire coding region of PRKCG was sequenced in an affected member of family AT08 and in a group of 39 unrelated patients with ataxia not attributable to trinucleotide expansions. Three different nonconservative missense mutations in highly conserved residues in C1, the cysteine-rich region of the protein, were found in family AT08, another familial case, and a sporadic case. The mutations cosegregated with disease in both families. Structural modeling predicts that two of these amino acid substitutions would severely abrogate the zinc-binding or phorbol ester-binding capabilities of the protein. Immunohistochemical studies on cerebellar tissue from an affected member of family AT08 demonstrated reduced staining for both PKC gamma and ataxin 1 in Purkinje cells, whereas staining for calbindin was preserved. These results strongly support a new mechanism for neuronal cell dysfunction and death in hereditary ataxias and suggest that there may be a common pathway for PKC gamma-related and polyglutamine-related neurodegeneration.

Relative Burden of Large CNVs on a Range of Neurodevelopmental Phenotypes

While numerous studies have implicated copy number variants (CNVs) in a range of neurological phenotypes, the impact relative to disease severity has been difficult to ascertain due to small sample sizes, lack of phenotypic details, and heterogeneity in platforms used for discovery. Using a customized microarray enriched for genomic hotspots, we assayed for large CNVs among 1,227 individuals with various neurological deficits including dyslexia (376), sporadic autism (350), and intellectual disability (ID) (501), as well as 337 controls. We show that the frequency of large CNVs (>1 Mbp) is significantly greater for ID–associated phenotypes compared to autism (p = 9.58×10−11, odds ratio = 4.59), dyslexia (p = 3.81×10−18, odds ratio = 14.45), or controls (p = 2.75×10−17, odds ratio = 13.71). There is a striking difference in the frequency of rare CNVs (>50 kbp) in autism (10%, p = 2.4×10−6, odds ratio = 6) or ID (16%, p = 3.55×10−12, odds ratio = 10) compared to dyslexia (2%) with essentially no difference in large CNV burden among dyslexia patients compared to controls. Rare CNVs were more likely to arise de novo (64%) in ID when compared to autism (40%) or dyslexia (0%). We observed a significantly increased large CNV burden in individuals with ID and multiple congenital anomalies (MCA) compared to ID alone (p = 0.001, odds ratio = 2.54). Our data suggest that large CNV burden positively correlates with the severity of childhood disability: ID with MCA being most severely affected and dyslexics being indistinguishable from controls. When autism without ID was considered separately, the increase in CNV burden was modest compared to controls (p = 0.07, odds ratio = 2.33).

Linkage analyses of four regions previously implicated in dyslexia: Confirmation of a locus on chromosome 15q

Dyslexia is a common, complex disorder, which is thought to have a genetic component. There have been numerous reports of linkage to several regions of the genome for dyslexia and continuous dyslexia‐related phenotypes. We attempted to confirm linkage of continuous measures of (1) accuracy and efficiency of phonological decoding; and (2) accuracy of single word reading (WID) to regions on chromosomes 2p, 6p, 15q, and 18p, using 111 families with a total of 898 members. We used both single‐marker and multipoint variance components linkage analysis and Markov Chain Monte Carlo (MCMC) joint segregation and linkage analysis for initial inspection of these regions. Positive results were followed with traditional parametric lod score analysis using a model estimated by MCMC segregation analysis. No positive linkage signals were found on chromosomes 2p, 6p, or 18p. Evidence of linkage of WID to chromosome 15q was found with both methods of analysis. The maximum single‐marker parametric lod score of 2.34 was obtained at a distance of 3 cM from D15S143. Multipoint analyses localized the putative susceptibility gene to the interval between markers GATA50C03 and D15S143, which falls between a region implicated in a recent genome screen for attention‐deficit/hyperactivity disorder, and DYX1C1, a candidate gene for dyslexia. This apparent multiplicity of linkage signals in the region for developmental disorders may be the result of errors in map and/or model specification obscuring the pleiotropic effect of a single gene on different phenotypes, or it may reflect the presence of multiple genes. This article contains supplementary material, which may be viewed at the American Journal of Medical Genetics website at http://www.interscience.wiley.com/jpages/0148‐7299:1/suppmat/index.html.

Replication of CNTNAP2 association with nonword repetition and support for FOXP2 association with timed reading and motor activities in a dyslexia family sample

Two functionally related genes, FOXP2 and CNTNAP2, influence language abilities in families with rare syndromic and common nonsyndromic forms of impaired language, respectively. We investigated whether these genes are associated with component phenotypes of dyslexia and measures of sequential motor ability. Quantitative transmission disequilibrium testing (QTDT) and linear association modeling were used to evaluate associations with measures of phonological memory (nonword repetition, NWR), expressive language (sentence repetition), reading (real word reading efficiency, RWRE; word attack, WATT), and timed sequential motor activities (rapid alternating place of articulation, RAPA; finger succession in the dominant hand, FS-D) in 188 family trios with a child with dyslexia. Consistent with a prior study of language impairment, QTDT in dyslexia showed evidence of CNTNAP2 single nucleotide polymorphism (SNP) association with NWR. For FOXP2, we provide the first evidence for SNP association with component phenotypes of dyslexia, specifically NWR and RWRE but not WATT. In addition, FOXP2 SNP associations with both RAPA and FS-D were observed. Our results confirm the role of CNTNAP2 in NWR in a dyslexia sample and motivate new questions about the effects of FOXP2 in neurodevelopmental disorders.

Evaluation of candidate genes for DYX1 and DYX2 in families with dyslexia.

Dyslexia is a common heterogeneous disorder with a significant genetic component. Multiple studies have replicated the evidence for linkage between variously defined phenotypes of dyslexia and chromosomal regions on 15q21 (DYX1) and 6p22.2 (DYX2). Based on association studies and the possibility for functional significance of several polymorphisms, candidate genes responsible for the observed linkage signal have been proposed—DYX1C1 for 15q21, and KIAA0319 and DCDC2 for 6p22.2. We investigated the evidence for contribution of these candidate genes to dyslexia in our sample of multigenerational families. Our previous quantitative linkage analyses in this dataset provided supportive evidence for linkage of dyslexia to the locus on chromosome 15, but not to the locus on chromosome 6. In the current study, we used probands from 191 families for a case control analysis, and proband‐parent trios for family‐based TDT analyses. The observation of weak evidence for transmission disequilibrium for one of the two studied polymorphisms in DYX1C1 suggests involvement of this gene in dyslexia in our dataset. We did not find evidence for the association of KIAA0319 or DCDC2 alleles to dyslexia in our sample. We observed a slight tendency for an intronic deletion in DCDC2 to be associated with worse performance on some quantitative measures of dyslexia in the probands in our sample, but not in their parents.

A genome scan in multigenerational families with dyslexia: Identification of a novel locus on chromosome 2q that contributes to phonological decoding efficiency.

Dyslexia is a common and complex developmental disorder manifested by unexpected difficulty in learning to read. Multiple different measures are used for diagnosis, and may reflect different biological pathways related to the disorder. Impaired phonological decoding (translation of written words without meaning cues into spoken words) is thought to be a core deficit. We present a genome scan of two continuous measures of phonological decoding ability: phonemic decoding efficiency (PDE) and word attack (WA). PDE measures both accuracy and speed of phonological decoding, whereas WA measures accuracy alone. Multipoint variance component linkage analyses (VC) and Markov chain Monte-Carlo (MCMC) multipoint joint linkage and segregation analyses were performed on 108 families. A strong signal was observed on chromosome 2 for PDE using both VC (LOD=2.65) and MCMC methods (intensity ratio (IR)=32.1). The IR is an estimate of the ratio of the posterior to prior probability of linkage in MCMC analysis. The chromosome 2 signal was not seen for WA. More detailed mapping with additional markers provided statistically significant evidence for linkage of PDE to chromosome 2, with VC-LOD=3.0 and IR=59.6 at D2S1399. Parametric analyses of PDE, using a model obtained by complex segregation analysis, provided a multipoint maximum LOD=2.89. The consistency of results from three analytic approaches provides strong evidence for a locus on chromosome 2 that influences speed but not accuracy of phonological decoding.

Genomewide scan for real-word reading subphenotypes of dyslexia: Novel chromosome 13 locus and genetic complexity

Dyslexia is a common learning disability exhibited as a delay in acquiring reading skills despite adequate intelligence and instruction. Reading single real words (real‐word reading, RWR) is especially impaired in many dyslexics. We performed a genome scan, using variance components (VC) linkage analysis and Bayesian Markov chain Monte Carlo (MCMC) joint segregation and linkage analysis, for three quantitative measures of RWR in 108 multigenerational families, with follow up of the strongest signals with parametric LOD score analyses. We used single‐word reading efficiency (SWE) to assess speed and accuracy of RWR, and word identification (WID) to assess accuracy alone. Adjusting SWE for WID provided a third measure of RWR efficiency. All three methods of analysis identified a strong linkage signal for SWE on chromosome 13q. Based on multipoint analysis with 13 markers we obtained a MCMC intensity ratio (IR) of 53.2 (chromosome‐wide P < 0.004), a VC LOD score of 2.29, and a parametric LOD score of 2.94, based on a quantitative‐trait model from MCMC segregation analysis (SA). A weaker signal for SWE on chromosome 2q occurred in the same location as a significant linkage peak seen previously in a scan for phonological decoding. MCMC oligogenic SA identified three models of transmission for WID, which could be assigned to two distinct linkage peaks on chromosomes 12 and 15. Taken together, these results indicate a locus for efficiency and accuracy of RWR on chromosome 13, and a complex model for inheritance of RWR accuracy with loci on chromosomes 12 and 15.

IFRD1 Is a Candidate Gene for SMNA on Chromosome 7q22-q23

We have established strong linkage evidence that supports mapping autosomal-dominant sensory/motor neuropathy with ataxia (SMNA) to chromosome 7q22-q32. SMNA is a rare neurological disorder whose phenotype encompasses both the central and the peripheral nervous system. In order to identify a gene responsible for SMNA, we have undertaken a comprehensive genomic evaluation of the region of linkage, including evaluation for repeat expansion and small deletions or duplications, capillary sequencing of candidate genes, and massively parallel sequencing of all coding exons. We excluded repeat expansion and small deletions or duplications as causative, and through microarray-based hybrid capture and massively parallel short-read sequencing, we identified a nonsynonymous variant in the human interferon-related developmental regulator gene 1 (IFRD1) as a disease-causing candidate. Sequence conservation, animal models, and protein structure evaluation support the involvement of IFRD1 in SMNA. Mutation analysis of IFRD1 in additional patients with similar phenotypes is needed for demonstration of causality and further evaluation of its importance in neurological diseases.

R47H Variant of TREM2 Associated With Alzheimer Disease in a Large Late-Onset Family Clinical, Genetic, and Neuropathological Study

IMPORTANCE The R47H variant in the triggering receptor expressed on myeloid cells 2 gene (TREM2), a modulator of the immune response of microglia, is a strong genetic risk factor for Alzheimer disease (AD) and possibly other neurodegenerative disorders. OBJECTIVE To investigate a large family with late-onset AD (LOAD), in which R47H cosegregated with 75% of cases. DESIGN, SETTING, AND PARTICIPANTS This study includes genetic and pathologic studies of families with LOAD from 1985 to 2014. A total of 131 families with LOAD (751 individuals) were included from the University of Washington Alzheimer Disease Research Center. To identify LOAD genes/risk factors in the LOAD123 family with 21 affected members and 12 autopsies, we sequenced 4 exomes. Candidate variants were tested for cosegregation with the disease. TREM2 R47H was genotyped in an additional 130 families with LOAD. We performed clinical and neuropathological assessments of patients with and without R47H and evaluated the variants effect on brain pathology, cellular morphology, and expression of microglial markers. MAIN OUTCOMES AND MEASURES We assessed the effect of TREM2 genotype on age at onset and disease duration. We compared Braak and Consortium to Establish a Registry for Alzheimers Disease scores, presence of α-synuclein and TAR DNA-binding protein 43 aggregates, and additional vascular or Parkinson pathology in TREM2 R47H carriers vs noncarriers. Microglial activation was assessed by quantitative immunohistochemistry and morphometry. RESULTS Twelve of 16 patients with AD in the LOAD123 family carried R47H. Eleven patients with dementia had apolipoprotein E 4 (ApoE4) and R47H genotypes. We also found a rare missense variant, D353N, in a nominated AD risk gene, unc-5 homolog C (UNC5C), in 5 affected individuals in the LOAD123 family. R47H carriers demonstrated a shortened disease duration (mean [SD], 6.7 [2.8] vs 11.1 [6.6] years; 2-tailed t test; P = .04) and more frequent α-synucleinopathy. The panmicroglial marker ionized calcium-binding adapter molecule 1 was decreased in all AD cases and the decrease was most pronounced in R47H carriers (mean [SD], in the hilus: 0.114 [0.13] for R47H_AD vs 0.574 [0.26] for control individuals; 2-tailed t test; P = .005 and vs 0.465 [0.32] for AD; P = .02; in frontal cortex gray matter: 0.006 [0.004] for R47H_AD vs 0.016 [0.01] for AD; P = .04 and vs 0.033 [0.013] for control individuals; P < .001). Major histocompatibility complex class II, a marker of microglial activation, was increased in all patients with AD (AD: 2.5, R47H_AD: 2.7, and control: 1.0; P < .01). CONCLUSIONS AND RELEVANCE Our results demonstrate a complex genetic landscape of LOAD, even in a single pedigree with an apparent autosomal dominant pattern of inheritance. ApoE4, TREM2 R47H, and rare variants in other genes, such as UNC5C D353N, are likely responsible for the notable occurrence of AD in this family. Our findings support the role of the TREM2 receptor in microglial clearance of aggregation-prone proteins that is compromised in R47H carriers and may accelerate the course of disease.