PSEUDOMARKER 2.0: efficient computation of likelihoods using NOMAD

Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.

Joint Linkage and Association Analysis Using GENEHUNTER-MODSCORE with an Application to Familial Pancreatic Cancer

INTRODUCTION

Joint linkage and association (JLA) analysis combines two disease gene mapping strategies: linkage information contained in families and association information contained in populations. Such a JLA analysis can increase mapping power, especially when the evidence for both linkage and association is low to moderate. Similarly, an association analysis based on haplotypes instead of single markers can increase mapping power when the association pattern is complex.

METHODS

In this paper, we present an extension to the GENEHUNTER-MODSCORE software package that enables a JLA analysis based on haplotypes and uses information from arbitrary pedigree types and unrelated individuals. Our new JLA method is an extension of the MOD score approach for linkage analysis, which allows the estimation of trait-model and linkage disequilibrium (LD) parameters, i.e., penetrance, disease-allele frequency, and haplotype frequencies. LD is modeled between alleles at a single diallelic disease locus and up to three diallelic test markers. Linkage information is contributed by additional multi-allelic flanking markers. We investigated the statistical properties of our JLA implementation using extensive simulations, and we compared our approach to another commonly used single-marker JLA test. To demonstrate the applicability of our new method in practice, we analyzed pedigree data from the German National Case Collection for Familial Pancreatic Cancer (FaPaCa).

RESULTS

Based on the simulated data, we demonstrated the validity of our JLA-MOD score analysis implementation and identified scenarios in which haplotype-based tests outperformed the single-marker test. The estimated trait-model and LD parameters were in good accordance with the simulated values. Our method outperformed another commonly used JLA single-marker test when the LD pattern was complex. The exploratory analysis of the FaPaCa families led to the identification of a promising genetic region on chromosome 22q13.33, which can serve as a starting point for future mutation analysis and molecular research in pancreatic cancer.

CONCLUSION

Our newly proposed JLA-MOD score method proves to be a valuable gene mapping and characterization tool, especially when either linkage or association information alone provide insufficient power to identify the disease-causing genetic variants.

Linkage analysis using whole exome sequencing data implicates SLC17A1, SLC17A3, TATDN2 and TMEM131L in type 1 diabetes in Kuwaiti families

Type 1 diabetes (T1D) is characterized by the progressive destruction of pancreatic β-cells, leading to insulin deficiency and lifelong dependency on exogenous insulin. Higher estimates of heritability rates in monozygotic twins, followed by dizygotic twins and sib-pairs, indicate the role of genetics in the pathogenesis of T1D. The incidence and prevalence of T1D are alarmingly high in Kuwait. Consanguineous marriages account for 50-70% of all marriages in Kuwait, leading to an excessive burden of recessive allele enrichment and clustering of familial disorders. Thus, genetic studies from this Arab region are expected to lead to the identification of novel gene loci for T1D. In this study, we performed linkage analyses to identify the recurrent genetic variants segregating in high-risk Kuwaiti families with T1D. We studied 18 unrelated Kuwaiti native T1D families using whole exome sequencing data from 86 individuals, of whom 37 were diagnosed with T1D. The study identified three potential loci with a LOD score of ≥ 3, spanning across four candidate genes, namely SLC17A1 (rs1165196:pT269I), SLC17A3 (rs942379: p.S370S), TATDN2 (rs394558:p.V256I), and TMEM131L (rs6848033:p.R190R). Upon examination of missense variants from these genes in the familial T1D dataset, we observed a significantly increased enrichment of the genotype homozygous for the minor allele at SLC17A3 rs56027330_p.G279R accounting for 16.2% in affected children from 6 unrelated Kuwaiti T1D families compared to 1000 genomes Phase 3 data (0.9%). Data from the NephQTL database revealed that the rs1165196, rs942379, rs394558, and rs56027330 SNPs exhibited genotype-based differential expression in either glomerular or tubular tissues. Data from the GTEx database revealed rs942379 and rs394558 as QTL variants altering the expression of TRIM38 and IRAK2 respectively. Global genome-wide association studies indicated that SLC17A1 rs1165196 and other variants from SLC17A3 are associated with uric acid concentrations and gout. Further evidence from the T1D Knowledge portal supported the role of shortlisted variants in T1D pathogenesis and urate metabolism. Our study suggests the involvement of SLC17A1, SLC17A3, TATDN2, and TMEM131L genes in familial T1D in Kuwait. An enrichment selection of genotype homozygous for the minor allele is observed at SLC17A3 rs56027330_p.G279R variant in affected members of Kuwaiti T1D families. Future studies may focus on replicating the findings in a larger T1D cohort and delineate the mechanistic details of the impact of these novel candidate genes on the pathophysiology of T1D.

Clinical and genetic analysis of familial neuromyelitis optica spectrum disorder in Chinese: associated with ubiquitin-specific peptidase USP18 gene variants

BACKGROUND

Familial clustering of neuromyelitis optica spectrum disorder (NMOSD) was present in Chinese. This study was to investigate the clinical characteristics and genetic background of familial NMOSD.

METHODS

Through questionnaires in four medical centres in 2016-2020, we identified 10 families with NMOSD aggregation. The statistical differences of clinical characteristics between familial and sporadic NMOSD (22 cases and 459 cases) were summarised. The whole-exome sequencing (WES) for seven families (13 cases and 13 controls) was analysed, compared with our previous WES data for sporadic NMOSD (228 cases and 1 400 controls). The family-based and population-based association and linkage analysis were conducted to identify the pathogenetic genes, the variant impacts were predicted.

RESULTS

The familial occurrence was 0.87% in Chinese. Familial patients had higher expanded disability status scale score than sporadic patients (p=0.03). The single-nucleotide polymorphism (SNP) rs2252257 in the promoter and enhancer of ubiquitin-specific peptidase USP18 was linked to familial NMOSD (p=7.8E-05, logarithm of the odds (LOD)=3.1), SNPs rs361553, rs2252257 and rs5746523 were related to sporadic NMOSD (p=1.29E-10, 3.45E-07 and 2.01E-09, respectively). Patients with the SNP rs361553 T/T genotype had higher recurrence rate than C/T or C/C genotype (1.22±0.85 vs 0.69±0.57 and 0.81±0.65, p=0.003 and 0.001, respectively). SNPs rs361553 and rs2252257 altered USP18 expression in brain and nerve tissues.

CONCLUSION

Most clinical characteristics of familial NMOSD were indistinguishable from sporadic NMOSD except for the worst episodes severity. USP18 with impaired intronic regulatory function contributed to the pathogenesis of NMOSD.

Whole genome sequencing identifies candidate genes for familial essential tremor and reveals biological pathways implicated in essential tremor aetiology

Background

Essential tremor (ET), one of the most common neurological disorders, has a phenotypically heterogeneous presentation characterized by bilateral kinetic tremor of the arms and, in some patients, tremor involving other body regions (e.g., head, voice). Genetic studies suggest that ET is genetically heterogeneous.

Methods

We analyzed whole genome sequence data (WGS) generated on 104 multi-generational white families with European ancestry affected by ET. Genome-wide parametric linkage and association scans were analyzed using adjusted logistic regression models through the application of the Pseudomarker software. To investigate the additional contribution of rare variants in familial ET, we also performed an aggregate variant non-parametric linkage (NPL) analysis using the collapsed haplotype method implemented in CHP-NPL software.

Findings

Parametric linkage analysis of common variants identified several loci with significant evidence of linkage (HLOD ≥3.6). Among the gene regions within the strongest ET linkage peaks were BTC (4q13.3, HLOD=4.53), N6AMT1 (21q21.3, HLOD=4.31), PCDH9 (13q21.32, HLOD=4.21), EYA1 (8q13.3, HLOD=4.04), RBFOX1 (16p13.3, HLOD=4.02), MAPT (17q21.31, HLOD=3.99) and SCARB2 (4q21.1, HLOD=3.65). CHP-NPL analysis identified fifteen additional genes with evidence of significant linkage (LOD ≥3.8). These genes include TUBB2A, VPS33B, STEAP1B, SPINK5, ZRANB1, TBC1D3C, PDPR, NPY4R, ETS2, ZNF736, SPATA21, ARL17A, PZP, BLK and CCDC94. In one ET family contributing to the linkage peak on chromosome 16p13.3, we identified a likely pathogenic heterozygous canonical splice acceptor variant in exon 2 of RBFOX1 (ENST00000547372; c.4-2A>G), that co-segregated with the ET phenotype in the family.

Interpretation

Linkage and association analyses of WGS identified several novel ET candidate genes, which are implicated in four major pathways that include 1) the epidermal growth factor receptor-phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha-AKT serine/threonine kinase 1 (EGFR-PI3K-AKT) and Mitogen-activated protein Kinase 1 (ERK) pathways, 2) Reactive oxygen species (ROS) and DNA repair, 3) gamma-aminobutyric acid-ergic (GABAergic) system and 4) RNA binding and regulation of RNA processes. Our study provides evidence for a possible overlap in the genetic architecture of ET, neurological disease, cancer and aging. The genes and pathways identified can be prioritized in future genetic and functional studies.

Funding

National Institutes of Health, NINDS, NS073872 (USA) and NIA AG058131(USA).

Novel Linkage Peaks Discovered for Diabetic Nephropathy in Individuals With Type 1 Diabetes

Genome-wide association studies (GWAS) and linkage studies have had limited success in identifying genome-wide significantly linked regions or risk loci for diabetic nephropathy (DN) in individuals with type 1 diabetes (T1D). As GWAS cohorts have grown, they have also included more documented and undocumented familial relationships. Here we computationally inferred and manually curated pedigrees in a study cohort of >6,000 individuals with T1D and their relatives without diabetes. We performed a linkage study for 177 pedigrees consisting of 452 individuals with T1D and their relatives using a genome-wide genotyping array with >300,000 single nucleotide polymorphisms and PSEUDOMARKER software. Analysis resulted in genome-wide significant linkage peaks on eight chromosomal regions from five chromosomes (logarithm of odds score >3.3). The highest peak was localized at the HLA region on chromosome 6p, but whether the peak originated from T1D or DN remained ambiguous. Of other significant peaks, the chromosome 4p22 region was localized on top of ARHGAP24, a gene associated with focal segmental glomerulosclerosis, suggesting this gene may play a role in DN as well. Furthermore, rare variants have been associated with DN and chronic kidney disease near the 4q25 peak, localized on top of CCSER1.

Variants in regulatory elements of PDE4D associate with major mental illness in the Finnish population

We have previously reported a replicable association between variants at the PDE4D gene and familial schizophrenia in a Finnish cohort. In order to identify the potential functional mutations underlying these previous findings, we sequenced 1.5 Mb of the PDE4D genomic locus in 20 families (consisting of 96 individuals and 79 independent chromosomes), followed by two stages of genotyping across 6668 individuals from multiple Finnish cohorts for major mental illnesses. We identified 4570 SNPs across the PDE4D gene, with 380 associated to schizophrenia (p ≤ 0.05). Importantly, two of these variants, rs35278 and rs165940, are located at transcription factor-binding sites, and displayed replicable association in the two-stage enlargement of the familial schizophrenia cohort (combined statistics for rs35278 p = 0.0012; OR = 1.18, 95% CI: 1.06-1.32; and rs165940 p = 0.0016; OR = 1.27, 95% CI: 1.13-1.41). Further analysis using additional cohorts and endophenotypes revealed that rs165940 principally associates within the psychosis (p = 0.025, OR = 1.18, 95% CI: 1.07-1.30) and cognitive domains of major mental illnesses (g-score p = 0.044, β = -0.033). Specifically, the cognitive domains represented verbal learning and memory (p = 0.0091, β = -0.044) and verbal working memory (p = 0.0062, β = -0.036). Moreover, expression data from the GTEx database demonstrated that rs165940 significantly correlates with the mRNA expression levels of PDE4D in the cerebellum (p-value = 0.04; m-value = 0.9), demonstrating a potential functional consequence for this variant. Thus, rs165940 represents the most likely functional variant for major mental illness at the PDE4D locus in the Finnish population, increasing risk broadly to psychotic disorders.

PedMiner: a tool for linkage analysis-based identification of disease-associated variants using family based whole-exome sequencing data

With the advances of next-generation sequencing technology, the field of disease research has been revolutionized. However, pinpointing the disease-causing variants from millions of revealed variants is still a tough task. Here, we have reviewed the existing linkage analysis tools and presented PedMiner, a web-based application designed to narrow down candidate variants from family based whole-exome sequencing (WES) data through linkage analysis. PedMiner integrates linkage analysis, variant annotation and prioritization in one automated pipeline. It provides graphical visualization of the linked regions along with comprehensive annotation of variants and genes within these linked regions. This efficient and comprehensive application will be helpful for the scientific community working on Mendelian inherited disorders using family based WES data.

A rare TTC30B variant is identified as a candidate for synpolydactyly in a Chinese pedigree

BACKGROUND

Syndactyly type II (synpolydactyly, SPD) is a rare autosomal dominant inherited disease with higher incomplete penetrance. Currently, several variants in HOXD13 and one deletion in FBLN1 have been associated with SPD. However, the causative variants in several SPD families and their etiological mechanism are still largely unknown.

METHODS

Whole exome and PCR-sanger sequencing followed by two-point linkage analysis were performed to identify the pathogenic variant in a six-generation Chinese pedigree. Homology modeling in combination with the RNAi and qRT-PCR experiments was used for revealing the pathogenic mechanism of the TTC30B variant.

RESULTS

A six-generation SPD family was reported. The affected subjects in this family had no other clinical malformation beyond SPD. A rare missense variant c.1157C>T [p.Ala375Val] (chr2:178416368, hg19) in TTC30B was demonstrated to be responsible for this SPD family. The modeling structure indicated that the Ala375 was evolutionarily and structurally conserved. The variant p.Ala375Val was predicted to be deleterious for protein structure and/or stability. Two-point linkage analysis resulted in a maximum LOD score of 3.1444 (P = 0.000071). Furthermore, we found that TTC30B was regulated by the Shh signaling pathway and the abnormal expression of TTC30B will affect the activation of the Shh signaling pathway in human retinal pigment epithelial cells.

CONCLUSIONS

This study demonstrates for the first time that an IFT (intraflagellar transport) - related gene TTC30B is implicated with SPD.

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.

A novel association of rs13334070 in the RPGRIP1L gene with adiposity factors discovered by joint linkage and linkage disequilibrium analysis in Iranian pedigrees: Tehran Cardiometabolic Genetic Study (TCGS)

Understanding the genetic and metabolic bases of obesity is helpful in planning and developing health strategies. Therefore, the first family-based joint linkage and linkage disequilibrium study was conducted in Iranian pedigrees to assess the relationship between obesity and single-nucleotide polymorphisms (SNPs) located in the 16q12.2 region. In the present study, a total of 13,344 individuals were included, of whom 12,502 individuals were within 3,109 pedigrees and 842 were unrelated singletons. To investigate the relationship between obesity and genetic variants, a joint model of linkage and linkage disequilibrium was applied. Moreover, a sequence kernel association test (SKAT) was used to evaluate the association of the SNP set with body size and lipid profile measurements. The joint model showed that rs13334070, in the intron 4 of the RPGRIP1L gene, has a significant association with obesity. According to the 4-gamete rule, which is a procedure for constructing SNP sets by considering recombination occurrence between SNPs, this polymorphism has a high correlation with six nearby SNPs that make an SNP set. SKAT showed that this SNP set has a significant association with body size factors, but almost no association with most of the lipid profile measurements. In conclusion, from the result of this study, it might be reasonable to consider RPGRIP1L as an important gene whose variations could be associated with obesity risk factors.

Genetic variants of TREML2 are associated with HLA-B27-positive ankylosing spondylitis

Although ankylosing spondylitis (AS) is a common, highly heritable arthropathy, the precise genetic mechanism underlying the disease remains elusive. Here, we investigate the disease-causing mutations in a large AS family with distinguished complexity, consisting of 23 patients covering four generations and exhibiting a mixed HLA-B27 (+) and (-) status. Linkage analysis with 32 members using three methods and whole-exome sequencing analysis with three HLA-B27 (+) patients, one HLA-B27 (-) patient, and one healthy individual did not identify a mutation common to all of the patients, strongly suggesting the existence of genetic heterogeneity in this large pedigree. However, if only B27-positive patients were analyzed, the linkage analysis located a 22-Mb region harboring the HLA gene cluster in chromosome 6 (LOD = 4.2), and the subsequent exome analysis identified two non-synonymous mutations in the TREML2 and IP6K3 genes. These genes were resequenced among 370 sporadic AS patients and 487 healthy individuals. A significantly higher mutation frequency of TREML2 was observed in AS patients (1.51% versus 0.21%). The results obtained for the AS pedigree and sporadic patients suggest that mutation of TREML2 is a major factor leading to AS for HLA-B27 (+) members in this large family and that TREML2 is also a susceptibility gene promoting the development of ankylosing spondylitis in HLA-B27 (+) individuals.

MGST2 and WNT2 are candidate genes for comitant strabismus susceptibility in Japanese patients

BACKGROUND/AIM

Strabismus is a common condition with misalignment between two eyes that may lead to decrease of visual acuity, lack of binocularity, and diplopia. It is caused by heterogeneous environmental and genetic risk factors. Our previous research has identified new chromosomal susceptibility loci in 4q28.3 and 7q31.2 regions for comitant strabismus in Japanese families. We conducted a verification study by linkage analysis to narrow the chromosomal loci down to a single gene.

METHODS

From Japanese and U.S. databases, 24 rsSNPs and 233 rsSNPs were chosen from the 4q28.3 and 7q31.2 region, respectively, and were typed in 108 affected subjects and 96 unaffected subjects of 58 families with primary and non-syndromic comitant strabismus. Three major analytical methods were used: transmission disequilibrium test (TDT), TDT allowing for errors (TDTae), and linkage analysis under dominant and recessive inheritance.

RESULTS

The SNPs with significant P values in TDT and TDTae were located solely at the gene, microsomal glutathione S-transferase 2 (MGST2), on chromosome 4q28.3 locus. In contrast, significant SNPs were dispersed in a few genes, containing wingless-type MMTV integration site family member 2 (WNT2), on chromosome 7q31.2 locus. The distribution of significant SNPs on the 7q31.2 locus showed that only the ST7 to WNT2 region in the same big haplotype block contained significant SNPs for all three methods of linkage analysis.

CONCLUSIONS

This study suggests that MGST2 and WNT2 are potential candidates for comitant strabismus in Japanese population.

Neuregulin signaling pathway in smoking behavior

Understanding molecular processes that link comorbid traits such as addictions and mental disorders can provide novel therapeutic targets. Neuregulin signaling pathway (NSP) has previously been implicated in schizophrenia, a neurodevelopmental disorder with high comorbidity to smoking. Using a Finnish twin family sample, we have previously detected association between nicotine dependence and ERBB4 (a neuregulin receptor), and linkage for smoking initiation at the ERBB4 locus on 2q33. Further, Neuregulin3 has recently been shown to associate with nicotine withdrawal in a behavioral mouse model. In this study, we scrutinized association and linkage between 15 036 common, low frequency and rare genetic variants in 10 NSP genes and phenotypes encompassing smoking and alcohol use. Using the Finnish twin family sample (N=1998 from 740 families), we detected 66 variants (representing 23 LD blocks) significantly associated (false discovery rate P<0.05) with smoking initiation, nicotine dependence and nicotine withdrawal. We comprehensively annotated the associated variants using expression (eQTL) and methylation quantitative trait loci (meQTL) analyses in a Finnish population sample. Among the 66 variants, we identified 25 eQTLs (in NRG1 and ERBB4), 22 meQTLs (in NRG3, ERBB4 and PSENEN), a missense variant in NRG1 (rs113317778) and a splicing disruption variant in ERBB4 (rs13385826). Majority of the QTLs in blood were replicated in silico using publicly available databases, with additional QTLs observed in brain. In conclusion, our results support the involvement of NSP in smoking behavior but not in alcohol use and abuse, and disclose functional potential for 56 of the 66 associated single-nucleotide polymorphism.

A genome-wide screen for acrophobia susceptibility loci in a Finnish isolate

Acrophobia, an abnormal fear of heights, is a specific phobia characterized as apprehension cued by the occurrence or anticipation of elevated spaces. It is considered a complex trait with onset influenced by both genetic and environmental factors. Identification of genetic risk variants would provide novel insight into the genetic basis of the fear of heights phenotype and contribute to the molecular-level understanding of its aetiology. Genetic isolates may facilitate identification of susceptibility alleles due to reduced genetic heterogeneity. We took advantage of an internal genetic isolate in Finland in which a distinct acrophobia phenotype appears to be segregating in pedigrees originally ascertained for schizophrenia. We conducted parametric, nonparametric, joint linkage and linkage disequilibrium analyses using a microsatellite marker panel, genotyped in families to search for chromosomal regions correlated with acrophobia. Our results implicated a few regions with suggestive evidence for linkage on chromosomes 4q28 (LOD = 2.17), 8q24 (LOD = 2.09) and 13q21-q22 (LOD = 2.22). We observed no risk haplotypes shared between different families. These results suggest that genetic predisposition to acrophobia in this genetic isolate is unlikely to be mediated by a small number of shared high-risk alleles, but rather has a complex genetic architecture.

Creative Activities in Music--A Genome-Wide Linkage Analysis

Creative activities in music represent a complex cognitive function of the human brain, whose biological basis is largely unknown. In order to elucidate the biological background of creative activities in music we performed genome-wide linkage and linkage disequilibrium (LD) scans in musically experienced individuals characterised for self-reported composing, arranging and non-music related creativity. The participants consisted of 474 individuals from 79 families, and 103 sporadic individuals. We found promising evidence for linkage at 16p12.1-q12.1 for arranging (LOD 2.75, 120 cases), 4q22.1 for composing (LOD 2.15, 103 cases) and Xp11.23 for non-music related creativity (LOD 2.50, 259 cases). Surprisingly, statistically significant evidence for linkage was found for the opposite phenotype of creative activity in music (neither composing nor arranging; NCNA) at 18q21 (LOD 3.09, 149 cases), which contains cadherin genes like CDH7 and CDH19. The locus at 4q22.1 overlaps the previously identified region of musical aptitude, music perception and performance giving further support for this region as a candidate region for broad range of music-related traits. The other regions at 18q21 and 16p12.1-q12.1 are also adjacent to the previously identified loci with musical aptitude. Pathway analysis of the genes suggestively associated with composing suggested an overrepresentation of the cerebellar long-term depression pathway (LTD), which is a cellular model for synaptic plasticity. The LTD also includes cadherins and AMPA receptors, whose component GSG1L was linked to arranging. These results suggest that molecular pathways linked to memory and learning via LTD affect music-related creative behaviour. Musical creativity is a complex phenotype where a common background with musicality and intelligence has been proposed. Here, we implicate genetic regions affecting music-related creative behaviour, which also include genes with neuropsychiatric associations. We also propose a common genetic background for music-related creative behaviour and musical abilities at chromosome 4.

AMIGO-Kv2.1 Potassium Channel Complex Is Associated With Schizophrenia-Related Phenotypes

The enormous variability in electrical properties of neurons is largely affected by a multitude of potassium channel subunits. Kv2.1 is a widely expressed voltage-dependent potassium channel and an important regulator of neuronal excitability. The Kv2.1 auxiliary subunit AMIGO constitutes an integral part of the Kv2.1 channel complex in brain and regulates the activity of the channel. AMIGO and Kv2.1 localize to the distinct somatodendritic clusters at the neuronal plasma membrane. Here we have created and characterized a mouse line lacking the AMIGO gene. Absence of AMIGO clearly reduced the amount of the Kv2.1 channel protein in mouse brain and altered the electrophysiological properties of neurons. These changes were accompanied by behavioral and pharmacological abnormalities reminiscent of those identified in schizophrenia. Concomitantly, we have detected an association of a rare, population-specific polymorphism of KV2.1 (KCNB1) with human schizophrenia in a genetic isolate enriched with schizophrenia. Our study demonstrates the involvement of AMIGO-Kv2.1 channel complex in schizophrenia-related behavioral domains in mice and identifies KV2.1 (KCNB1) as a strong susceptibility gene for schizophrenia spectrum disorders in humans.

Linkage analyses in Caribbean Hispanic families identify novel loci associated with familial late-onset Alzheimer's disease

INTRODUCTION

We performed linkage analyses in Caribbean Hispanic families with multiple late-onset Alzheimer's disease (LOAD) cases to identify regions that may contain disease causative variants.

METHODS

We selected 67 LOAD families to perform genome-wide linkage scan. Analysis of the linked regions was repeated using the entire sample of 282 families. Validated chromosomal regions were analyzed using joint linkage and association.

RESULTS

We identified 26 regions linked to LOAD (HLOD ≥3.6). We validated 13 of the regions (HLOD ≥2.5) using the entire family sample. The strongest signal was at 11q12.3 (rs2232932: HLODmax = 4.7, Pjoint = 6.6 × 10(-6)), a locus located ∼2 Mb upstream of the membrane-spanning 4A gene cluster. We additionally identified a locus at 7p14.3 (rs10255835: HLODmax = 4.9, Pjoint = 1.2 × 10(-5)), a region harboring genes associated with the nervous system (GARS, GHRHR, and NEUROD6).

DISCUSSION

Future sequencing efforts should focus on these regions because they may harbor familial LOAD causative mutations.

Whole exome sequencing for handedness in a large and highly consanguineous family

Pinpointing genes involved in non-right-handedness has the potential to clarify developmental contributions to human brain lateralization. Major-gene models have been considered for human handedness which allow for phenocopy and reduced penetrance, i.e. an imperfect correspondence between genotype and phenotype. However, a recent genome-wide association scan did not detect any common polymorphisms with substantial genetic effects. Previous linkage studies in families have also not yielded significant findings. Genetic heterogeneity and/or polygenicity are therefore indicated, but it remains possible that relatively rare, or even unique, major-genetic effects may be detectable in certain extended families with many non-right-handed members. Here we applied whole exome sequencing to 17 members from a single, large consanguineous family from Pakistan. Multipoint linkage analysis across all autosomes did not yield clear candidate genomic regions for involvement in the trait and single-point analysis of exomic variation did not yield clear candidate mutations/genes. Any genetic contribution to handedness in this unusual family is therefore likely to have a complex etiology, as at the population level.

Genetic variants associated with susceptibility to psychosis in late-onset Alzheimer's disease families

Psychotic symptoms are frequent in late-onset Alzheimer's disease (LOAD) patients. Although the risk for psychosis in LOAD is genetically mediated, no genes have been identified. To identify loci potentially containing genetic variants associated with risk of psychosis in LOAD, a total of 263 families from the National Institute of Aging-LOAD cohort were classified into psychotic (LOAD+P, n = 215) and nonpsychotic (LOAD-P, n = 48) families based on the presence/absence of psychosis during the course of LOAD. The LOAD+P families yielded strong evidence of linkage on chromosome 19q13 (two-point [2-pt] ​logarithm of odds [LOD] = 3.8, rs2285513 and multipoint LOD = 2.7, rs541169). Joint linkage and association in 19q13 region detected strong association with rs2945988 (p = 8.7 × 10(-7)). Linkage results for the LOAD-P families yielded nonsignificant 19q13 LOD scores. Several 19q13 single-nucleotide polymorphisms generalized the association of LOAD+P in a Caribbean Hispanic (CH) cohort, and the strongest signal was rs10410711 (pmeta = 5.1 × 10(-5)). A variant located 24 kb upstream of rs10410711 and rs10421862 was strongly associated with LOAD+P (pmeta = 1.0 × 10(-5)) in a meta-analysis of the CH cohort and an additional non-Hispanic Caucasian dataset. Identified variants rs2945988 and rs10421862 affect brain gene expression levels. Our results suggest that genetic variants in genes on 19q13, some of which are involved in brain development and neurodegeneration, may influence the susceptibility to psychosis in LOAD patients.

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.