A simple method for analyzing exome sequencing data shows distinct levels of nonsynonymous variation for human immune and nervous system genes

Functional analysis of an R311C variant of Ca2+ -calmodulin-dependent protein kinase kinase-2 (CaMKK2) found as a de novo mutation in a patient with bipolar disorder

OBJECTIVES

Loss-of-function mutations in the gene encoding the calcium-calmodulin (Ca²⁺ -CaM)-dependent protein kinase kinase-2 (CaMKK2) enzyme are linked to bipolar disorder. Recently, a de novo arginine to cysteine (R311C) mutation in CaMKK2 was identified from a whole exome sequencing study of bipolar patients and their unaffected parents. The aim of the present study was to determine the functional consequences of the R311C mutation on CaMKK2 activity and regulation by Ca²⁺ -CaM.

METHODS

The effects of the R311C mutation on CaMKK2 activity and Ca²⁺ -CaM activation were examined using a radiolabeled adenosine triphosphate (ATP) kinase assay. We performed immunoblot analysis to determine whether the R311C mutation impacts threonine-85 (T85) autophosphorylation, an activating phosphorylation site on CaMKK2 that has also been implicated in bipolar disorder. We also expressed the R311C mutant in CaMKK2 knockout HAP1 cells and used immunoblot analysis and an MTS reduction assay to study its effects on Ca²⁺ -dependent downstream signaling and cell viability, respectively.

RESULTS

The R311C mutation maps to the conserved HRD motif within the catalytic loop of CaMKK2 and caused a marked reduction in kinase activity and Ca²⁺ -CaM activation. The R311C mutation virtually abolished T85 autophosphorylation in response to Ca²⁺ -CaM and exerted a dominant-negative effect in cells as it impaired the ability of wild-type CaMKK2 to initiate downstream signaling and maintain cell viability.

CONCLUSIONS

The highly disruptive, loss-of-function impact of the de novo R311C mutation in human CaMKK2 provides a compelling functional rationale for being considered a potential rare monogenic cause of bipolar disorder.

The population genomics of rhesus macaques (Macaca mulatta) based on whole-genome sequences

Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primate in biomedical research, have the largest natural geographic distribution of any nonhuman primate, and have been the focus of much evolutionary and behavioral investigation. Consequently, rhesus macaques are one of the most thoroughly studied nonhuman primate species. However, little is known about genome-wide genetic variation in this species. A detailed understanding of extant genomic variation among rhesus macaques has implications for the use of this species as a model for studies of human health and disease, as well as for evolutionary population genomics. Whole-genome sequencing analysis of 133 rhesus macaques revealed more than 43.7 million single-nucleotide variants, including thousands predicted to alter protein sequences, transcript splicing, and transcription factor binding sites. Rhesus macaques exhibit 2.5-fold higher overall nucleotide diversity and slightly elevated putative functional variation compared with humans. This functional variation in macaques provides opportunities for analyses of coding and noncoding variation, and its cellular consequences. Despite modestly higher levels of nonsynonymous variation in the macaques, the estimated distribution of fitness effects and the ratio of nonsynonymous to synonymous variants suggest that purifying selection has had stronger effects in rhesus macaques than in humans. Demographic reconstructions indicate this species has experienced a consistently large but fluctuating population size. Overall, the results presented here provide new insights into the population genomics of nonhuman primates and expand genomic information directly relevant to primate models of human disease.

Improved variant calling accuracy by merging replicates in whole-exome sequencing studies

In large scale population-based whole-exome sequencing (WES) studies, there are some samples occasionally sequenced two or more times due to a variety of reasons. To investigate how to efficiently utilize these duplicated sequencing data, we conducted comprehensive evaluation of variant calling strategies. 92 samples subjected to WES twice were selected from a large population study. These 92 duplicated samples were divided into two groups: group H consisting of the higher sequencing depth for each subject and group L consisting of the lower depth for each subject. The merged samples for each subject were put in a third group M. Using the GATK multisample toolkit, we compared variant calling accuracy among three strategies. Hierarchical clustering analysis indicated that the two replicates for each subject showed high homogeneity. The comparative analyses on the basis of heterozygous-homozygous ratio (Hete/Homo), transition-transversion ratio (Ti/Tv), and overlapping rate with the 1000 Genomes Project consistently showed that the data quality of the SNPs detected from the M group was more accurate than that of SNPs detected from the H and L groups. These results suggested that merging homogeneous duplicated exomes instead of using one of them could improve variant calling accuracy.