Investigation of Rare Single-Nucleotide PCDH15 Variants in Schizophrenia and Autism Spectrum Disorders

Both schizophrenia (SCZ) and autism spectrum disorders (ASD) are neuropsychiatric disorders with overlapping genetic etiology. Protocadherin 15 (PCDH15), which encodes a member of the cadherin super family that contributes to neural development and function, has been cited as a risk gene for neuropsychiatric disorders. Recently, rare variants of large effect have been paid attention to understand the etiopathology of these complex disorders. Thus, we evaluated the impacts of rare, single-nucleotide variants (SNVs) in PCDH15 on SCZ or ASD. First, we conducted coding exon-targeted resequencing of PCDH15 with next-generation sequencing technology in 562 Japanese patients (370 SCZ and 192 ASD) and detected 16 heterozygous SNVs. We then performed association analyses on 2,096 cases (1,714 SCZ and 382 ASD) and 1,917 controls with six novel variants of these 16 SNVs. Of these six variants, four (p.R219K, p.T281A, p.D642N, c.3010-1G>C) were ultra-rare variants (minor allele frequency < 0.0005) that may increase disease susceptibility. Finally, no statistically significant association between any of these rare, heterozygous PCDH15 point variants and SCZ or ASD was found. Our results suggest that a larger sample size of resequencing subjects is necessary to detect associations between rare PCDH15 variants and neuropsychiatric disorders.     

DCLK1 variants are associated across schizophrenia and attention deficit/hyperactivity disorder

Doublecortin and calmodulin like kinase 1 (DCLK1) is implicated in synaptic plasticity and neurodevelopment. Genetic variants in DCLK1 are associated with cognitive traits, specifically verbal memory and general cognition. We investigated the role of DCLK1 variants in three psychiatric disorders that have neuro-cognitive dysfunctions: schizophrenia (SCZ), bipolar affective disorder (BP) and attention deficit/hyperactivity disorder (ADHD). We mined six genome wide association studies (GWASs) that were available publically or through collaboration; three for BP, two for SCZ and one for ADHD. We also genotyped the DCLK1 region in additional samples of cases with SCZ, BP or ADHD and controls that had not been whole-genome typed. In total, 9895 subjects were analysed, including 5308 normal controls and 4,587 patients (1,125 with SCZ, 2,496 with BP and 966 with ADHD). Several DCLK1 variants were associated with disease phenotypes in the different samples. The main effect was observed for rs7989807 in intron 3, which was strongly associated with SCZ alone and even more so when cases with SCZ and ADHD were combined (P-value = 4 × 10(-5) and 4 × 10(-6), respectively). Associations were also observed with additional markers in intron 3 (combination of SCZ, ADHD and BP), intron 19 (SCZ+BP) and the 3'UTR (SCZ+BP). Our results suggest that genetic variants in DCLK1 are associated with SCZ and, to a lesser extent, with ADHD and BP. Interestingly the association is strongest when SCZ and ADHD are considered together, suggesting common genetic susceptibility. Given that DCLK1 variants were previously found to be associated with cognitive traits, these results are consistent with the role of DCLK1 in neurodevelopment and synaptic plasticity.     

Integrated Model of De Novo and Inherited Genetic Variants Yields Greater Power to Identify Risk Genes

De novo mutations affect risk for many diseases and disorders, especially those with early-onset. An example is autism spectrum disorders (ASD). Four recent whole-exome sequencing (WES) studies of ASD families revealed a handful of novel risk genes, based on independent de novo loss-of-function (LoF) mutations falling in the same gene, and found that de novo LoF mutations occurred at a twofold higher rate than expected by chance. However successful these studies were, they used only a small fraction of the data, excluding other types of de novo mutations and inherited rare variants. Moreover, such analyses cannot readily incorporate data from case-control studies. An important research challenge in gene discovery, therefore, is to develop statistical methods that accommodate a broader class of rare variation. We develop methods that can incorporate WES data regarding de novo mutations, inherited variants present, and variants identified within cases and controls. TADA, for Transmission And De novo Association, integrates these data by a gene-based likelihood model involving parameters for allele frequencies and gene-specific penetrances. Inference is based on a Hierarchical Bayes strategy that borrows information across all genes to infer parameters that would be difficult to estimate for individual genes. In addition to theoretical development we validated TADA using realistic simulations mimicking rare, large-effect mutations affecting risk for ASD and show it has dramatically better power than other common methods of analysis. Thus TADA''s integration of various kinds of WES data can be a highly effective means of identifying novel risk genes. Indeed, application of TADA to WES data from subjects with ASD and their families, as well as from a study of ASD subjects and controls, revealed several novel and promising ASD candidate genes with strong statistical support.     

The power comparison of the haplotype-based collapsing tests and the variant-based collapsing tests for detecting rare variants in pedigrees

Background

Both common and rare genetic variants have been shown to contribute to the etiology of complex diseases. Recent genome-wide association studies (GWAS) have successfully investigated how common variants contribute to the genetic factors associated with common human diseases. However, understanding the impact of rare variants, which are abundant in the human population (one in every 17 bases), remains challenging. A number of statistical tests have been developed to analyze collapsed rare variants identified by association tests. Here, we propose a haplotype-based approach. This work inspired by an existing statistical framework of the pedigree disequilibrium test (PDT), which uses genetic data to assess the effects of variants in general pedigrees. We aim to compare the performance between the haplotype-based approach and the rare variant-based approach for detecting rare causal variants in pedigrees.

Results

Extensive simulations in the sequencing setting were carried out to evaluate and compare the haplotype-based approach with the rare variant methods that drew on a more conventional collapsing strategy. As assessed through a variety of scenarios, the haplotype-based pedigree tests had enhanced statistical power compared with the rare variants based pedigree tests when the disease of interest was mainly caused by rare haplotypes (with multiple rare alleles), and vice versa when disease was caused by rare variants acting independently. For most of other situations when disease was caused both by haplotypes with multiple rare alleles and by rare variants with similar effects, these two approaches provided similar power in testing for association.

Conclusions

The haplotype-based approach was designed to assess the role of rare and potentially causal haplotypes. The proposed rare variants-based pedigree tests were designed to assess the role of rare and potentially causal variants. This study clearly documented the situations under which either method performs better than the other. All tests have been implemented in a software, which was submitted to the Comprehensive R Archive Network (CRAN) for general use as a computer program named rvHPDT.     

Variants in GATA4 are a rare cause of familial and sporadic congenital diaphragmatic hernia

Congenital diaphragmatic hernia (CDH) is characterized by incomplete formation of the diaphragm occurring as either an isolated defect or in association with other anomalies. Genetic factors including aneuploidies and copy number variants are important in the pathogenesis of many cases of CDH, but few single genes have been definitively implicated in human CDH. In this study, we used whole exome sequencing (WES) to identify a paternally inherited novel missense GATA4 variant (c.754C>T; p.R252W) in a familial case of CDH with incomplete penetrance. Phenotypic characterization of the family included magnetic resonance imaging of the chest and abdomen demonstrating asymptomatic defects in the diaphragm in the two “unaffected” missense variant carriers. Screening 96 additional CDH patients identified a de novo heterozygous GATA4 variant (c.848G>A; p.R283H) in a non-isolated CDH patient. In summary, GATA4 is implicated in both familial and sporadic CDH, and our data suggests that WES may be a powerful tool to discover rare variants for CDH.     

Association and linkage analysis of RGS4 polymorphisms with schizophrenia and bipolar disorder in Brazil

Linkage and association studies in five independently ascertained samples have suggested that polymorphisms of the regulator of G-protein signaling 4 (RGS4) may confer risk for schizophrenia (SCZ). Suggestive evidence for association with bipolar disorder (BD) has also been presented. However, the associated alleles and haplotypes have differed among the samples. Data from other independent samples may clarify the putative associations. Hence, we investigated an independent, ethnically diverse Brazilian population comprising patients with SCZ (n=271) or BD1 (n=306), who were contrasted with 576 community-based controls. Parents of 49 SCZ cases and 44 BD cases were available for transmission disequilibrium tests (TDTs). Four RGS4 single-nucleotide polymorphisms (SNPs) 1, 4, 7 and 18 putatively associated with SCZ were investigated. In the SCZ samples, significant case-control differences were not observed for individual SNPs or haplotypes, though the TDT suggested transmission distortion similar to that observed in the initial report. For the BD sample, case-control comparisons revealed no significant differences for individual SNPs, but an omnibus test suggested differences in the overall distribution of haplotypes bearing all four SNPs (SNP-EM Omnibus likelihood ratio test; P=0.003). The TDT revealed over-transmission of allele A at SNP7 (P=0.016), as well as haplotypes incorporating this allele. However, global tests incorporating all haplotypes yielded only suggestive trends for association (P=0.19). In conclusion, association with SCZ was not detected in the present analyses. The failure to detect an association may be related to inadequate power or to confounds related to ethnic admixture. Suggestive associations with BD detected here require further investigation in a larger sample.     

A Novel Relationship for Schizophrenia,Bipolar, and Major Depressive Disorder. Part 8: a Hint from Chromosome 8 High Density Association Screen

Convergent evidence from genetics, symptomatology, and psychopharmacology implies that there are intrinsic connections between schizophrenia (SCZ), bipolar disorder (BPD), and major depressive disorder (MDD); for example, any two or even three of these disorders could co-exist in some families. A total of 48,753 single nucleotide polymorphisms (SNPs) on chromosome 8 were genotyped by Affymetrix Genome-Wide Human SNP array 6.0 on 119 SCZ, 253 BPD (type I), 177 MDD patients, and 1000 controls. Associated SNP loci were comprehensively revealed, and outstanding susceptibility genes were identified including CSMD1, NRG1, PXDNL, SGCZ, and TMEM66. Unexpectedly, flanking genes for up to 95.9 % of the associated SNPs were replicated (P ≤ 9.9E?8) in the enlarged cohort of 986 SCZ patients. Considering convergent evidence, our results implicate that bipolar disorder and major depressive disorder might be subtypes of schizophrenia.     

Whole Exome Sequencing of Distant Relatives in Multiplex Families Implicates Rare Variants in Candidate Genes for Oral Clefts

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.     

Alternative-Splicing in the Exon-10 Region of GABAA Receptor β2 Subunit Gene: Relationships between Novel Isoforms and Psychotic Disorders

Background

Non-coding single nucleotide polymorphisms (SNPs) in GABRB2, the gene for β₂-subunit of gamma-aminobutyric acid type A (GABA_A) receptor, have been associated with schizophrenia (SCZ) and quantitatively correlated to mRNA expression and alternative splicing.

Methods and Findings

Expression of the Exon 10 region of GABRB2 from minigene constructs revealed this region to be an “alternative splicing hotspot” that readily gave rise to differently spliced isoforms depending on intron sequences. This led to a search in human brain cDNA libraries, and the discovery of two novel isoforms, β_2S1 and β_2S2, bearing variations in the neighborhood of Exon-10. Quantitative real-time PCR analysis of postmortem brain samples showed increased β_2S1 expression and decreased β_2S2 expression in both SCZ and bipolar disorder (BPD) compared to controls. Disease-control differences were significantly correlated with SNP rs187269 in BPD males for both β_2S1 and β_2S2 expressions, and significantly correlated with SNPs rs2546620 and rs187269 in SCZ males for β_2S2 expression. Moreover, site-directed mutagenesis indicated that Thr³⁶⁵, a potential phosphorylation site in Exon-10, played a key role in determining the time profile of the ATP-dependent electrophysiological current run-down.

Conclusion

This study therefore provided experimental evidence for the importance of non-coding sequences in the Exon-10 region in GABRB2 with respect to β₂-subunit splicing diversity and the etiologies of SCZ and BPD.     

Pooled Sequencing of 531 Genes in Inflammatory Bowel Disease Identifies an Associated Rare Variant in BTNL2 and Implicates Other Immune Related Genes

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn’s disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10⁻¹⁰, OR = 2.3[95% CI = 1.75–3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1–5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.     

Relationships between substance abuse/dependence and psychiatric disorders based on polygenic scores

Genetic susceptibility to substance use disorders (SUDs) is partially shared between substances. Heritability of any substance dependence, estimated as 54%, is partly explained by additive effects of common variants. Comorbidity between SUDs and other psychiatric disorders is frequent. The present study aims to analyze the additive role of common variants in this comorbidity using polygenic scores (PGSs) based on genome‐wide association study discovery samples of schizophrenia (SCZ), bipolar disorder, attention‐deficit/hyperactivity disorder, autism spectrum disorder, major depressive disorder and anxiety disorders, available from large consortia. PGSs were calculated for 534 patients meeting DSM‐IV criteria for dependence of a substance and abuse/dependence of another substance between alcohol, tobacco, cannabis, cocaine, opiates, hypnotics, stimulants, hallucinogens and solvents; and 587 blood donors from the same population, Iberians from Galicia, as controls. Significance of the PGS and percentage of variance explained were calculated by logistic regression. Using discovery samples of similar size, significant associations with SUDs were detected for SCZ PGS. SCZ PGS explained more variance in SUDs than in most psychiatric disorders. Cross‐disorder PGS based on five psychiatric disorders was significant after adjustment for the effect of SCZ PGS. SCZ PGS was significantly higher in women than in men abusing alcohol. Our findings indicate that SUDs share genetic susceptibility with SCZ to a greater extent than with other psychiatric disorders, including externalizing disorders such as attention‐deficit/hyperactivity disorder. Women have lower probability to develop substance abuse/dependence than men at similar PGS probably because of a higher social pressure against excessive drug use in women.     

Negative linkage disequilibrium between amino acid changing variants reveals interference among deleterious mutations in the human genome

Evolutionary forces like Hill-Robertson interference and negative epistasis can lead to deleterious mutations being found on distinct haplotypes. However, the extent to which these forces depend on the selection and dominance coefficients of deleterious mutations and shape genome-wide patterns of linkage disequilibrium (LD) in natural populations with complex demographic histories has not been tested. In this study, we first used forward-in-time simulations to predict how negative selection impacts LD. Under models where deleterious mutations have additive effects on fitness, deleterious variants less than 10 kb apart tend to be carried on different haplotypes relative to pairs of synonymous SNPs. In contrast, for recessive mutations, there is no consistent ordering of how selection coefficients affect LD decay, due to the complex interplay of different evolutionary effects. We then examined empirical data of modern humans from the 1000 Genomes Project. LD between derived alleles at nonsynonymous SNPs is lower compared to pairs of derived synonymous variants, suggesting that nonsynonymous derived alleles tend to occur on different haplotypes more than synonymous variants. This result holds when controlling for potential confounding factors by matching SNPs for frequency in the sample (allele count), physical distance, magnitude of background selection, and genetic distance between pairs of variants. Lastly, we introduce a new statistic H_R^(j) which allows us to detect interference using unphased genotypes. Application of this approach to high-coverage human genome sequences confirms our finding that nonsynonymous derived alleles tend to be located on different haplotypes more often than are synonymous derived alleles. Our findings suggest that interference may play a pervasive role in shaping patterns of LD between deleterious variants in the human genome, and consequently influences genome-wide patterns of LD.     

OLA-DRB1 microsatellite variants are associated with ovine growth and reproduction traits

The DRB1 intron 2 (GT)_n (GA)_m microsatellite was genotyped in experimental flocks of seven Merinoland rams and 249 ewes as well as their offspring (381 lambs) from consecutive lambings. A total of 16 DRB1 alleles were detected, ranging between 353 and 857 bp. In comparison with carriers of other alleles, the ewes carrying the predominant 411 bp allele had higher values of all the recorded fertility traits. For ewes carrying the 394 and 857 bp alleles, the birth weight of lambs was about 400 g higher as compared to the residual group of ewes. The observed associations could be due to differences in disease resistance, cell recognition or tissue differentiation between carriers of various MHC haplotypes which can in turn affect individual fertility and growth performance.     

Association of PDE4B polymorphisms and schizophrenia in Northwestern Han Chinese

The phosphodiesterase 4B (PDE4B) is a candidate susceptibility gene for schizophrenia (SCZ), interacting with DISC1, a known genetic risk factor for SCZ. To examine if variants within PDE4B gene are associated with SCZ in Northwestern Han Chinese, and if these effects vary in gender-specific subgroup, we analyzed 20 SNPs, selected from previous studies and preliminary HapMap data analyses with minor allele frequency (MAF) ≥ 20%, in a cohort of 428 cases and 572 controls from genetically independent Northwestern Han Chinese. Single SNP association, haplotype association and sex-specific association analysis were performed. We found that rs472952 is significantly associated with SCZ and rs7537440 is associated with SCZ in females. Further analysis indicated that a haplotype block spanning PDE4B2 splice site is highly associated with SCZ and several haplotypes in this block have about twofold to threefold increase in cases. Our results provide further evidence that PDE4B may play important roles in the etiology of SCZ.     

HapFABIA: Identification of very short segments of identity by descent characterized by rare variants in large sequencing data

Identity by descent (IBD) can be reliably detected for long shared DNA segments, which are found in related individuals. However, many studies contain cohorts of unrelated individuals that share only short IBD segments. New sequencing technologies facilitate identification of short IBD segments through rare variants, which convey more information on IBD than common variants. Current IBD detection methods, however, are not designed to use rare variants for the detection of short IBD segments. Short IBD segments reveal genetic structures at high resolution. Therefore, they can help to improve imputation and phasing, to increase genotyping accuracy for low-coverage sequencing and to increase the power of association studies. Since short IBD segments are further assumed to be old, they can shed light on the evolutionary history of humans. We propose HapFABIA, a computational method that applies biclustering to identify very short IBD segments characterized by rare variants. HapFABIA is designed to detect short IBD segments in genotype data that were obtained from next-generation sequencing, but can also be applied to DNA microarray data. Especially in next-generation sequencing data, HapFABIA exploits rare variants for IBD detection. HapFABIA significantly outperformed competing algorithms at detecting short IBD segments on artificial and simulated data with rare variants. HapFABIA identified 160 588 different short IBD segments characterized by rare variants with a median length of 23 kb (mean 24 kb) in data for chromosome 1 of the 1000 Genomes Project. These short IBD segments contain 752 000 single nucleotide variants (SNVs), which account for 39% of the rare variants and 23.5% of all variants. The vast majority—152 000 IBD segments—are shared by Africans, while only 19 000 and 11 000 are shared by Europeans and Asians, respectively. IBD segments that match the Denisova or the Neandertal genome are found significantly more often in Asians and Europeans but also, in some cases exclusively, in Africans. The lengths of IBD segments and their sharing between continental populations indicate that many short IBD segments from chromosome 1 existed before humans migrated out of Africa. Thus, rare variants that tag these short IBD segments predate human migration from Africa. The software package HapFABIA is available from Bioconductor. All data sets, result files and programs for data simulation, preprocessing and evaluation are supplied at http://www.bioinf.jku.at/research/short-IBD.     

Haplotype Diversity in “Source-Sink” Dynamics of <Emphasis Type="Italic">Escherichia coli</Emphasis> Urovirulence

FimH, the mannose-specific, type 1 fimbrial adhesin of Escherichia coli, acquires amino acid replacements adaptive in extraintestinal niches (the genitourinary tract) but detrimental in the main habitat (the large intestine). This microevolutionary dynamics is reminiscent of an ecological “source-sink” model of continuous species spread from a stable primary habitat (source) into transient secondary niches (sink), with eventual extinction of the sink-evolved populations. Here, we have adapted two ecological analytical tools—diversity indexes D _S and α—to compare size and frequency distributions of fimH haplotypes between evolutionarily conserved FimH variants (“source” haplotypes) and FimH variants with adaptive mutations (putative “sink” haplotypes). Both indexes show two- to threefold increased diversity of the sink fimH haplotypes relative to the source haplotypes, a pattern that ran opposite to those seen with nonstructural fimbrial genes (fimC and fimI) and housekeeping loci (adk and fumC) but similar to that seen with another fimbrial adhesin of E. coli, papG-II, also implicated in extraintestinal infections. The increased diversity of the sink pool of adhesin genes is due to the increased richness of the haplotypes (the number of unique haplotypes), rather than their evenness (the extent of similarity in relative abundances). Taken together, this pattern supports a continuous emergence and extinction of the gene alleles adaptive to virulence sink habitats of E. coli, rather than a one-time change in the habitat conditions. Thus, ecological methods of species diversity analysis can be successfully adapted to characterize the emergence of microbial virulence in bacterial pathogens subject to source-sink dynamics. [Reviewing Editor: Dr. Margaret Riley]     

Exome Sequencing in 53 Sporadic Cases of Schizophrenia Identifies 18 Putative Candidate Genes

Schizophrenia (SCZ) is a severe, debilitating mental illness which has a significant genetic component. The identification of genetic factors related to SCZ has been challenging and these factors remain largely unknown. To evaluate the contribution of de novo variants (DNVs) to SCZ, we sequenced the exomes of 53 individuals with sporadic SCZ and of their non-affected parents. We identified 49 DNVs, 18 of which were predicted to alter gene function, including 13 damaging missense mutations, 2 conserved splice site mutations, 2 nonsense mutations, and 1 frameshift deletion. The average number of exonic DNV per proband was 0.88, which corresponds to an exonic point mutation rate of 1.7×10⁻⁸ per nucleotide per generation. The non-synonymous-to-synonymous mutation ratio of 2.06 did not differ from neutral expectations. Overall, this study provides a list of 18 putative candidate genes for sporadic SCZ, and when combined with the results of similar reports, identifies a second proband carrying a non-synonymous DNV in the RGS12 gene.