Expression analysis and genotyping of DGKZ: a GWAS-derived risk gene for schizophrenia

Whole-genome sequencing identifies novel genes for autism in Chinese trios

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with high genetic heritability but heterogeneity. Fully understanding its genetics requires whole-genome sequencing (WGS), but the ASD studies utilizing WGS data in Chinese population are limited. In this study, we present a WGS study for 334 individuals, including 112 ASD patients and their non-ASD parents. We identified 146 de novo variants in coding regions in 85 cases and 60 inherited variants in coding regions. By integrating these variants with an association model, we identified 33 potential risk genes (P<0.001) enriched in neuron and regulation related biological process. Besides the well-known ASD genes (SCN2A, NF1, SHANK3, CHD8 etc.), several high confidence genes were highlighted by a series of functional analyses, including CTNND1, DGKZ, LRP1, DDN, ZNF483, NR4A2, SMAD6, INTS1, and MRPL12, with more supported evidence from GO enrichment, expression and network analysis. We also integrated RNA-seq data to analyze the effect of the variants on the gene expression and found 12 genes in the individuals with the related variants had relatively biased expression. We further presented the clinical phenotypes of the proband carrying the risk genes in both our samples and Caucasian samples to show the effect of the risk genes on phenotype. Regarding variants in non-coding regions, a total of 74 de novo variants and 30 inherited variants were predicted as pathogenic with high confidence, which were mapped to specific genes or regulatory features. The number of de novo variants found in patient was significantly associated with the parents' ages at the birth of the child, and gender with trend. We also identified small de novo structural variants in ASD trios. The results in this study provided important evidence for understanding the genetic mechanism of ASD.

Comparative analysis of PDZ-binding motifs in the diacylglycerol kinase family

Diacylglycerol kinases (DGKs) control local and temporal amounts of diacylglycerol (DAG) and phosphatidic acid (PA) by converting DAG to PA through phosphorylation in cells. Certain DGK enzymes possess C-terminal sequences that encode potential PDZ-binding motifs (PBMs), which could be involved in their recruitment into supramolecular signaling complexes. In this study, we used two different interactomic approaches, quantitative native holdup (nHU) and qualitative affinity purification (AP), both coupled to mass spectrometry (MS) to investigate the PDZ partners associated with the potential PBMs of DGKs. Complementing these results with site-specific affinity interactomic data measured on isolated PDZ domain fragments and PBM motifs, as well as evolutionary conservation analysis of the PBMs of DGKs, we explored functional differences within different DGK groups. All our results indicate that putative PBM sequences of type II enzymes, namely DGKδ, DGKη, and DGKκ, are likely to be nonfunctional. In contrast, type IV enzymes, namely DGKζ and DGKι, possess highly promiscuous PBMs that interact with a set of PDZ proteins with very similar affinity interactomes. The combination of various interactomic assays and evolutionary analyses provides a useful strategy for identifying functional domains and motifs within diverse enzyme families.

Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder

Three-dimensional chromatin interactions regulate gene expressions. The significance of de novo mutations (DNMs) in chromatin interactions remains poorly understood for autism spectrum disorder (ASD). We generated 813 whole-genome sequences from 242 Korean simplex families to detect DNMs, and identified target genes which were putatively affected by non-coding DNMs in chromatin interactions. Non-coding DNMs in chromatin interactions were significantly involved in transcriptional dysregulations related to ASD risk. Correspondingly, target genes showed spatiotemporal expressions relevant to ASD in developing brains and enrichment in biological pathways implicated in ASD, such as histone modification. Regarding clinical features of ASD, non-coding DNMs in chromatin interactions particularly contributed to low intelligence quotient levels in ASD probands. We further validated our findings using two replication cohorts, Simons Simplex Collection (SSC) and MSSNG, and showed the consistent enrichment of non-coding DNM-disrupted chromatin interactions in ASD probands. Generating human induced pluripotent stem cells in two ASD families, we were able to demonstrate that non-coding DNMs in chromatin interactions alter the expression of target genes at the stage of early neural development. Taken together, our findings indicate that non-coding DNMs in ASD probands lead to early neurodevelopmental disruption implicated in ASD risk via chromatin interactions.

Genetic control of RNA splicing and its distinct role in complex trait variation

Most genetic variants identified from genome-wide association studies (GWAS) in humans are noncoding, indicating their role in gene regulation. Previous studies have shown considerable links of GWAS signals to expression quantitative trait loci (eQTLs) but the links to other genetic regulatory mechanisms, such as splicing QTLs (sQTLs), are underexplored. Here, we introduce an sQTL mapping method, testing for heterogeneity between isoform-eQTLeffects (THISTLE), with improved power over competing methods. Applying THISTLE together with a complementary sQTL mapping strategy to brain transcriptomic (n = 2,865) and genotype data, we identified 12,794 genes with cis-sQTLs at P < 5 × 10⁻⁸, approximately 61% of which were distinct from eQTLs. Integrating the sQTL data into GWAS for 12 brain-related complex traits (including diseases), we identified 244 genes associated with the traits through cis-sQTLs, approximately 61% of which could not be discovered using the corresponding eQTL data. Our study demonstrates the distinct role of most sQTLs in the genetic regulation of transcription and complex trait variation.

A powerful method for splicing quantitative trait loci (sQTL) mapping, THISTLE, is presented and applied to a collection of 2,865 brain samples. Integration with GWAS identifies 244 genes associated via cis-sQTLs, of which 61% were not identified using expression QTLs.

Finding associations in a heterogeneous setting: statistical test for aberration enrichment

Most two-group statistical tests find broad patterns such as overall shifts in mean, median, or variance. These tests may not have enough power to detect effects in a small subset of samples, e.g., a drug that works well only on a few patients. We developed a novel statistical test targeting such effects relevant for clinical trials, biomarker discovery, feature selection, etc. We focused on finding meaningful associations in complex genetic diseases in gene expression, miRNA expression, and DNA methylation. Our test outperforms traditional statistical tests in simulated and experimental data and detects potentially disease-relevant genes with heterogeneous effects.