A hybrid likelihood model for sequence-based disease association studies

Molecular Modeling and In Vitro Functional Analysis of the RGS12 PDZ Domain Variant Associated with High-Penetrance Familial Bipolar Disorder

Bipolar disorder's etiology involves genetics, environmental factors, and gene-environment interactions, underlying its heterogeneous nature and treatment complexity. In 2020, Forstner and colleagues catalogued 378 sequence variants co-segregating with familial bipolar disorder. A notable candidate was an R59Q missense mutation in the PDZ (PSD-95/Dlg1/ZO-1) domain of RGS12. We previously demonstrated that RGS12 loss removes negative regulation on the kappa opioid receptor, disrupting basal ganglia dopamine homeostasis and dampening responses to dopamine-eliciting psychostimulants. Here, we investigated the R59Q variation in the context of potential PDZ domain functional alterations. We first validated a new target for the wildtype RGS12 PDZ domain-the SAPAP3 C-terminus-by molecular docking, surface plasmon resonance (SPR), and co-immunoprecipitation. While initial molecular dynamics (MD) studies predicted negligible effects of the R59Q variation on ligand binding, SPR showed a significant reduction in binding affinity for the three peptide targets tested. AlphaFold2-generated models predicted a modest reduction in protein-peptide interactions, which is consistent with the reduced binding affinity observed by SPR, suggesting that the substituted glutamine side chain may weaken the affinity of RGS12 for its in vivo binding targets, likely through allosteric changes. This difference may adversely affect the CNS signaling related to dynorphin and dopamine in individuals with this R59Q variation, potentially impacting bipolar disorder pathophysiology.

Endoplasmic stress-inducing variants in carboxyl ester lipase and pancreatic cancer risk

BACKGROUND

Endoplasmic reticulum (ER) stress-inducing variants in several pancreatic secretory enzymes have been associated with pancreatic disease. Multiple variants in CEL, encoding carboxyl ester lipase, are known to cause maturity-onset diabetes of the young (MODY8) but have not been implicated in pancreatic cancer risk.

METHODS

The prevalence of ER stress-inducing variants in the CEL gene was compared among pancreatic cancer cases vs. controls. Variants were identified by next-generation sequencing and confirmed by Sanger sequencing. Variants of uncertain significance (VUS) were assessed for their effect on the secretion of CEL protein and variants with reduced protein secretion were evaluated to determine if they induced endoplasmic reticulum stress.

RESULTS

ER stress-inducing CEL variants were found in 34 of 986 cases with sporadic pancreatic ductal adenocarcinoma, and 21 of 1045 controls (P = 0.055). Most of the variants were either the CEL-HYB1 variant, the I488T variant, or the combined CEL-HYB1/I488T variant; one case had a MODY8 variant.

CONCLUSION

This case/control analysis finds ER stress-inducing CEL variants are not associated with an increased likelihood of having pancreatic cancer.

Genomic and neuroimaging approaches to bipolar disorder

BACKGROUND

To date, besides genome-wide association studies, a variety of other genetic analyses (e.g. polygenic risk scores, whole-exome sequencing and whole-genome sequencing) have been conducted, and a large amount of data has been gathered for investigating the involvement of common, rare and very rare types of DNA sequence variants in bipolar disorder. Also, non-invasive neuroimaging methods can be used to quantify changes in brain structure and function in patients with bipolar disorder.

AIMS

To provide a comprehensive assessment of genetic findings associated with bipolar disorder, based on the evaluation of different genomic approaches and neuroimaging studies.

METHOD

We conducted a PubMed search of all relevant literatures from the beginning to the present, by querying related search strings.

RESULTS

ANK3, CACNA1C, SYNE1, ODZ4 and TRANK1 are five genes that have been replicated as key gene candidates in bipolar disorder pathophysiology, through the investigated studies. The percentage of phenotypic variance explained by the identified variants is small (approximately 4.7%). Bipolar disorder polygenic risk scores are associated with other psychiatric phenotypes. The ENIGMA-BD studies show a replicable pattern of lower cortical thickness, altered white matter integrity and smaller subcortical volumes in bipolar disorder.

CONCLUSIONS

The low amount of explained phenotypic variance highlights the need for further large-scale investigations, especially among non-European populations, to achieve a more complete understanding of the genetic architecture of bipolar disorder and the missing heritability. Combining neuroimaging data with genetic data in large-scale studies might help researchers acquire a better knowledge of the engaged brain regions in bipolar disorder.

The Translational Machine: A novel machine-learning approach to illuminate complex genetic architectures

The Translational Machine (TM) is a machine learning (ML)-based analytic pipeline that translates genotypic/variant call data into biologically contextualized features that richly characterize complex variant architectures and permit greater interpretability and biological replication. It also reduces potentially confounding effects of population substructure on outcome prediction. The TM consists of three main components. First, replicable but flexible feature engineering procedures translate genome-scale data into biologically informative features that appropriately contextualize simple variant calls/genotypes within biological and functional contexts. Second, model-free, nonparametric ML-based feature filtering procedures empirically reduce dimensionality and noise of both original genotype calls and engineered features. Third, a powerful ML algorithm for feature selection is used to differentiate risk variant contributions across variant frequency and functional prediction spectra. The TM simultaneously evaluates potential contributions of variants operative under polygenic and heterogeneous models of genetic architecture. Our TM enables integration of biological information (e.g., genomic annotations) within conceptual frameworks akin to geneset-/pathways-based and collapsing methods, but overcomes some of these methods' limitations. The full TM pipeline is executed in R. Our approach and initial findings from its application to a whole-exome schizophrenia case-control data set are presented. These TM procedures extend the findings of the primary investigation and yield novel results.

An interpretable low-complexity machine learning framework for robust exome-based in-silico diagnosis of Crohn's disease patients

Whole exome sequencing (WES) data are allowing researchers to pinpoint the causes of many Mendelian disorders. In time, sequencing data will be crucial to solve the genome interpretation puzzle, which aims at uncovering the genotype-to-phenotype relationship, but for the moment many conceptual and technical problems need to be addressed. In particular, very few attempts at the in-silico diagnosis of oligo-to-polygenic disorders have been made so far, due to the complexity of the challenge, the relative scarcity of the data and issues such as batch effects and data heterogeneity, which are confounder factors for machine learning (ML) methods. Here, we propose a method for the exome-based in-silico diagnosis of Crohn’s disease (CD) patients which addresses many of the current methodological issues. First, we devise a rational ML-friendly feature representation for WES data based on the gene mutational burden concept, which is suitable for small sample sizes datasets. Second, we propose a Neural Network (NN) with parameter tying and heavy regularization, in order to limit its complexity and thus the risk of over-fitting. We trained and tested our NN on 3 CD case-controls datasets, comparing the performance with the participants of previous CAGI challenges. We show that, notwithstanding the limited NN complexity, it outperforms the previous approaches. Moreover, we interpret the NN predictions by analyzing the learned patterns at the variant and gene level and investigating the decision process leading to each prediction.

Whole-exome sequencing of 81 individuals from 27 multiply affected bipolar disorder families

Bipolar disorder (BD) is a highly heritable neuropsychiatric disease characterized by recurrent episodes of depression and mania. Research suggests that the cumulative impact of common alleles explains 25-38% of phenotypic variance, and that rare variants may contribute to BD susceptibility. To identify rare, high-penetrance susceptibility variants for BD, whole-exome sequencing (WES) was performed in three affected individuals from each of 27 multiply affected families from Spain and Germany. WES identified 378 rare, non-synonymous, and potentially functional variants. These spanned 368 genes, and were carried by all three affected members in at least one family. Eight of the 368 genes harbored rare variants that were implicated in at least two independent families. In an extended segregation analysis involving additional family members, five of these eight genes harbored variants showing full or nearly full cosegregation with BD. These included the brain-expressed genes RGS12 and NCKAP5, which were considered the most promising BD candidates on the basis of independent evidence. Gene enrichment analysis for all 368 genes revealed significant enrichment for four pathways, including genes reported in de novo studies of autism (p_adj < 0.006) and schizophrenia (p_adj = 0.015). These results suggest a possible genetic overlap with BD for autism and schizophrenia at the rare-sequence-variant level. The present study implicates novel candidate genes for BD development, and may contribute to an improved understanding of the biological basis of this common and often devastating disease.

The impact of a fine-scale population stratification on rare variant association test results

Population stratification is a well-known confounding factor in both common and rare variant association analyses. Rare variants tend to be more geographically clustered than common variants, because of their more recent origin. However, it is not yet clear if population stratification at a very fine scale (neighboring administrative regions within a country) would lead to statistical bias in rare variant analyses. As the inclusion of convenience controls from external studies is indeed a common procedure, in order to increase the power to detect genetic associations, this problem is important. We studied through simulation the impact of a fine scale population structure on different rare variant association strategies, assessing type I error and power. We showed that principal component analysis (PCA) based methods of adjustment for population stratification adequately corrected type I error inflation at the largest geographical scales, but not at finest scales. We also showed in our simulations that adding controls obviously increased power, but at a considerably lower level when controls were drawn from another population.

Exome sequencing in large, multiplex bipolar disorder families from Cuba

Bipolar disorder (BD) is a major psychiatric illness affecting around 1% of the global population. BD is characterized by recurrent manic and depressive episodes, and has an estimated heritability of around 70%. Research has identified the first BD susceptibility genes. However, the underlying pathways and regulatory networks remain largely unknown. Research suggests that the cumulative impact of common alleles with small effects explains only around 25-38% of the phenotypic variance for BD. A plausible hypothesis therefore is that rare, high penetrance variants may contribute to BD risk. The present study investigated the role of rare, nonsynonymous, and potentially functional variants via whole exome sequencing in 15 BD cases from two large, multiply affected families from Cuba. The high prevalence of BD in these pedigrees renders them promising in terms of the identification of genetic risk variants with large effect sizes. In addition, SNP array data were used to calculate polygenic risk scores for affected and unaffected family members. After correction for multiple testing, no significant increase in polygenic risk scores for common, BD-associated genetic variants was found in BD cases compared to healthy relatives. Exome sequencing identified a total of 17 rare and potentially damaging variants in 17 genes. The identified variants were shared by all investigated BD cases in the respective pedigree. The most promising variant was located in the gene SERPING1 (p.L349F), which has been reported previously as a genome-wide significant risk gene for schizophrenia. The present data suggest novel candidate genes for BD susceptibility, and may facilitate the discovery of disease-relevant pathways and regulatory networks.

Mutations in the pancreatic secretory enzymes CPA1 and CPB1 are associated with pancreatic cancer

To evaluate whether germline variants in genes encoding pancreatic secretory enzymes contribute to pancreatic cancer susceptibility, we sequenced the coding regions of CPB1 and other genes encoding pancreatic secretory enzymes and known pancreatitis susceptibility genes (PRSS1, CPA1, CTRC, and SPINK1) in a hospital series of pancreatic cancer cases and controls. Variants in CPB1, CPA1 (encoding carboxypeptidase B1 and A1), and CTRC were evaluated in a second set of cases with familial pancreatic cancer and controls. More deleterious CPB1 variants, defined as having impaired protein secretion and induction of endoplasmic reticulum (ER) stress in transfected HEK 293T cells, were found in the hospital series of pancreatic cancer cases (5/986, 0.5%) than in controls (0/1,045, P = 0.027). Among familial pancreatic cancer cases, ER stress-inducing CPB1 variants were found in 4 of 593 (0.67%) vs. 0 of 967 additional controls (P = 0.020), with a combined prevalence in pancreatic cancer cases of 9/1,579 vs. 0/2,012 controls (P < 0.01). More ER stress-inducing CPA1 variants were also found in the combined set of hospital and familial cases with pancreatic cancer than in controls [7/1,546 vs. 1/2,012; P = 0.025; odds ratio, 9.36 (95% CI, 1.15-76.02)]. Overall, 16 (1%) of 1,579 pancreatic cancer cases had an ER stress-inducing CPA1 or CPB1 variant, compared with 1 of 2,068 controls (P < 0.00001). No other candidate genes had statistically significant differences in variant prevalence between cases and controls. Our study indicates ER stress-inducing variants in CPB1 and CPA1 are associated with pancreatic cancer susceptibility and implicate ER stress in pancreatic acinar cells in pancreatic cancer development.

DoEstRare: A statistical test to identify local enrichments in rare genomic variants associated with disease

Next-generation sequencing technologies made it possible to assay the effect of rare variants on complex diseases. As an extension of the "common disease-common variant" paradigm, rare variant studies are necessary to get a more complete insight into the genetic architecture of human traits. Association studies of these rare variations show new challenges in terms of statistical analysis. Due to their low frequency, rare variants must be tested by groups. This approach is then hindered by the fact that an unknown proportion of the variants could be neutral. The risk level of a rare variation may be determined by its impact but also by its position in the protein sequence. More generally, the molecular mechanisms underlying the disease architecture may involve specific protein domains or inter-genic regulatory regions. While a large variety of methods are optimizing functionality weights for each single marker, few evaluate variant position differences between cases and controls. Here, we propose a test called DoEstRare, which aims to simultaneously detect clusters of disease risk variants and global allele frequency differences in genomic regions. This test estimates, for cases and controls, variant position densities in the genetic region by a kernel method, weighted by a function of allele frequencies. We compared DoEstRare with previously published strategies through simulation studies as well as re-analysis of real datasets. Based on simulation under various scenarios, DoEstRare was the sole to consistently show highest performance, in terms of type I error and power both when variants were clustered or not. DoEstRare was also applied to Brugada syndrome and early-onset Alzheimer's disease data and provided complementary results to other existing tests. DoEstRare, by integrating variant position information, gives new opportunities to explain disease susceptibility. DoEstRare is implemented in a user-friendly R package.

Incorporation of protein binding effects into likelihood ratio test for exome sequencing data

Statistical association studies are an important tool in detecting novel disease genes. However, for sequencing data, association studies confront the challenge of low power because of relatively small data sample size and rare variants. Incorporating biological information that reflects disease mechanism is likely to strengthen the association evidence of disease genes, and thus increase the power of association studies. In this paper, we annotate non-synonymous single-nucleotide variants according to protein binding sites (BSs) by using a more accurate BS prediction method. We then incorporate this information into association study through a statistical framework of likelihood ratio test (LRT) based on weighted burden score of single-nucleotide variants (SNVs). The strategy is applied to Genetic Analysis Workshop 19 exome-sequencing data for detecting novel genes associated to hypotension. The SNV-weighting LRT idea is empirically verified by the simulated phenotypes (336 cases and 1607 controls), and the weights based on BS annotation are applied to the real phenotypes (394 cases and 1457 controls). Such strategy of weighting the prior information on protein functional sites is shown to be superior to the unweighted LRT and serves as a good complement to the existing association tests. Several putative genes are reported; some of them are functionally related to hypertension according to the previous evidence in the literature.

Filtering genetic variants and placing informative priors based on putative biological function

High-density genetic marker data, especially sequence data, imply an immense multiple testing burden. This can be ameliorated by filtering genetic variants, exploiting or accounting for correlations between variants, jointly testing variants, and by incorporating informative priors. Priors can be based on biological knowledge or predicted variant function, or even be used to integrate gene expression or other omics data. Based on Genetic Analysis Workshop (GAW) 19 data, this article discusses diversity and usefulness of functional variant scores provided, for example, by PolyPhen2, SIFT, or RegulomeDB annotations. Incorporating functional scores into variant filters or weights and adjusting the significance level for correlations between variants yielded significant associations with blood pressure traits in a large family study of Mexican Americans (GAW19 data set). Marker rs218966 in gene PHF14 and rs9836027 in MAP4 significantly associated with hypertension; additionally, rare variants in SNUPN significantly associated with systolic blood pressure. Variant weights strongly influenced the power of kernel methods and burden tests. Apart from variant weights in test statistics, prior weights may also be used when combining test statistics or to informatively weight p values while controlling false discovery rate (FDR). Indeed, power improved when gene expression data for FDR-controlled informative weighting of association test p values of genes was used. Finally, approaches exploiting variant correlations included identity-by-descent mapping and the optimal strategy for joint testing rare and common variants, which was observed to depend on linkage disequilibrium structure.

Whole Genome Sequencing Defines the Genetic Heterogeneity of Familial Pancreatic Cancer

Pancreatic cancer is projected to become the second leading cause of cancer-related death in the United States by 2020. A familial aggregation of pancreatic cancer has been established, but the cause of this aggregation in most families is unknown. To determine the genetic basis of susceptibility in these families, we sequenced the germline genomes of 638 patients with familial pancreatic cancer and the tumor exomes of 39 familial pancreatic adenocarcinomas. Our analyses support the role of previously identified familial pancreatic cancer susceptibility genes such as BRCA2, CDKN2A, and ATM, and identify novel candidate genes harboring rare, deleterious germline variants for further characterization. We also show how somatic point mutations that occur during hematopoiesis can affect the interpretation of genome-wide studies of hereditary traits. Our observations have important implications for the etiology of pancreatic cancer and for the identification of susceptibility genes in other common cancer types.

SIGNIFICANCE

The genetic basis of disease susceptibility in the majority of patients with familial pancreatic cancer is unknown. We whole genome sequenced 638 patients with familial pancreatic cancer and demonstrate that the genetic underpinning of inherited pancreatic cancer is highly heterogeneous. This has significant implications for the management of patients with familial pancreatic cancer.

Whole genome/exome sequencing in mood and psychotic disorders

Recent developments in DNA sequencing technologies have allowed for genetic studies using whole genome or exome analysis, and these have been applied in the study of mood and psychotic disorders, including bipolar disorder, depression, schizophrenia, and schizoaffective disorder. In this review, the current situation, recent findings, methodological problems, and future directions of whole genome/exome analysis studies of these disorders are summarized. Whole genome/exome studies of bipolar disorder have included pedigree analysis and case-control studies, demonstrating the role of previously implicated pathways, such as calcium signaling, cyclic adenosine monophosphate response element binding protein (CREB) signaling, and potassium channels. Extensive analysis of trio families and case-control studies showed that de novo mutations play a role in the genetic architecture of schizophrenia and indicated that mutations in several molecular pathways, including chromatin regulation, activity-regulated cytoskeleton, post-synaptic density, N-methyl-D-aspartate receptor, and targets of fragile X mental retardation protein, are associated with this disorder. Depression is a heterogeneous group of diseases and studies using exome analysis have been conducted to identify rare mutations causing Mendelian diseases that accompany depression. In the near future, clarification of the genetic architecture of bipolar disorder and schizophrenia is expected. Identification of causative mutations using these new technologies will facilitate neurobiological studies of these disorders.

New developments in the genetics of bipolar disorder

The last several years have been breakthrough ones in bipolar disorder (BPD) genetics, as the field has identified robust risk variants for the first time. Leading the way have been genome-wide association studies (GWAS) that have assessed common genetic markers across very large groups of patients and controls. These have resulted in findings in genes including ANK3, CACNA1C, SYNE1, ODZ4, and TRANK1. Additional studies have begun to examine the biology of these genes and how risk variants influence aspects of brain and behavior that underlie BPD. For example, carriers of the CACNA1C risk variant have been found to exhibit hippocampal and anterior cingulate dysfunction during episodic memory recall. This work has shed additional light on the relationship of bipolar susceptibility variants to other disorders, particularly schizophrenia. Even larger BPD GWAS are expected with samples now amassed of 21,035 cases and 28,758 controls. Studies have examined the pharmacogenomics of BPD with studies of lithium response, yielding high profile results that remain to be confirmed. The next frontier in the field is the identification of rare bipolar susceptibility variants through large-scale DNA sequencing. While only a couple of papers have been published to date, many studies are underway. The Bipolar Sequencing Consortium has been formed to bring together all of the groups working in this area, and to perform meta-analyses of the data generated. The consortium, with 13 member groups, now has exome data on ~3,500 cases and ~5,000 controls, and on ~162 families. The focus will likely shift within several years from exome data to whole genome data as costs of obtaining such data continue to drop. Gene-mapping studies are now providing clear results that provide insights into the pathophysiology of the disorder. Sequencing studies should extend this process further. Findings could eventually set the stage for rational therapeutic development.

Identifying rare and common disease associated variants in genomic data using Parkinson's disease as a model

BACKGROUND

Genome-wide association studies have been successful in identifying common genetic variants for human diseases. However, much of the heritable variation associated with diseases such as Parkinson's disease remains unknown suggesting that many more risk loci are yet to be identified. Rare variants have become important in disease association studies for explaining missing heritability. Methods for detecting this type of association require prior knowledge on candidate genes and combining variants within the region. These methods may suffer from power loss in situations with many neutral variants or causal variants with opposite effects.

RESULTS

We propose a method capable of scanning genetic variants to identify the region most likely harbouring disease gene with rare and/or common causal variants. Our method assigns a score at each individual variant based on our scoring system. It uses aggregate scores to identify the region with disease association. We evaluate performance by simulation based on 1000 Genomes sequencing data and compare with three commonly used methods. We use a Parkinson's disease case-control dataset as a model to demonstrate the application of our method. Our method has better power than CMC and WSS and similar power to SKAT-O with well-controlled type I error under simulation based on 1000 Genomes sequencing data. In real data analysis, we confirm the association of α-synuclein gene (SNCA) with Parkinson's disease (p = 0.005). We further identify association with hyaluronan synthase 2 (HAS2, p = 0.028) and kringle containing transmembrane protein 1 (KREMEN1, p = 0.006). KREMEN1 is associated with Wnt signalling pathway which has been shown to play an important role for neurodegeneration in Parkinson's disease.

CONCLUSIONS

Our method is time efficient and less sensitive to inclusion of neutral variants and direction effect of causal variants. It can narrow down a genomic region or a chromosome to a disease associated region. Using Parkinson's disease as a model, our method not only confirms association for a known gene but also identifies two genes previously found by other studies. In spite of many existing methods, we conclude that our method serves as an efficient alternative for exploring genomic data containing both rare and common variants.

The Growing Importance of CNVs: New Insights for Detection and Clinical Interpretation

Differences between genomes can be due to single nucleotide variants, translocations, inversions, and copy number variants (CNVs, gain or loss of DNA). The latter can range from sub-microscopic events to complete chromosomal aneuploidies. Small CNVs are often benign but those larger than 500 kb are strongly associated with morbid consequences such as developmental disorders and cancer. Detecting CNVs within and between populations is essential to better understand the plasticity of our genome and to elucidate its possible contribution to disease. Hence there is a need for better-tailored and more robust tools for the detection and genome-wide analyses of CNVs. While a link between a given CNV and a disease may have often been established, the relative CNV contribution to disease progression and impact on drug response is not necessarily understood. In this review we discuss the progress, challenges, and limitations that occur at different stages of CNV analysis from the detection (using DNA microarrays and next-generation sequencing) and identification of recurrent CNVs to the association with phenotypes. We emphasize the importance of germline CNVs and propose strategies to aid clinicians to better interpret structural variations and assess their clinical implications.