eQTL colocalization analysis highlights novel susceptibility genes in Autism Spectrum Disorders (ASD)

Autism Spectrum Disorders (ASD) are a group of neurodevelopmental disorders (NDDs) characterized by difficulties in social interaction and communication, repetitive behavior, and restricted interests. ASD has proven to have a strong genetic component. However, defining causal genes is still one of the main challenges in GWAS, since the vast majority (>90%) of detected signals lie within the non-coding genome. Expression quantitative trait locus (eQTL) colocalization analysis determines whether a specific variant is responsible for both a local eQTL and GWAS association and has helped leverage data and rendering gene discovery for a wide array of diseases. Here we further mine the largest ASD GWAS performed to date (18,381 cases and 27,969 controls) altogether with GWAS summary statistics from the main PGC studies (Schizophrenia, MD (Major Depression) and ADHD (Attention Deficit/Hyperactivity Disorder)), by using eQTpLot, a newly developed tool that illustrates the colocalization of GWAS and eQTL signals in a locus, and the enrichment of and correlation between the candidate gene eQTLs and trait-significant variants. This analysis points up 8 genes with a significant eQTL colocalization signal in ASD (CRHR1, KANSL1, MANBA, MAPT, MMP12, NKX2-2, PTPRE and WNT3) and one gene (SRPK2) with a marginally significant colocalization signal (r = 0.69, p < 1 × 10-6), and specifically highlights the potentially causal role of MAPT (r = 0.76, p < 1 × 10-6), NKX2-2 (r = 0.71, p-value = 2.26-02) and PTPRE (r = 0.97, p-value = 2.63-04) when restricting the analysis to brain tissue.

The non-coding genome in Autism Spectrum Disorders

Autism Spectrum Disorders (ASD) are a group of neurodevelopmental disorders (NDDs) characterized by difficulties in social interaction and communication, repetitive behavior, and restricted interests. While ASD have been proven to have a strong genetic component, current research largely focuses on coding regions of the genome. However, non-coding DNA, which makes up for ∼99% of the human genome, has recently been recognized as an important contributor to the high heritability of ASD, and novel sequencing technologies have been a milestone in opening up new directions for the study of the gene regulatory networks embedded within the non-coding regions. Here, we summarize current progress on the contribution of non-coding alterations to the pathogenesis of ASD and provide an overview of existing methods allowing for the study of their functional relevance, discussing potential ways of unraveling ASD's "missing heritability".

Deep sequencing of GDF5 reveals the absence of rare variants at this important osteoarthritis susceptibility locus

OBJECTIVE

The common single nucleotide polymorphism (SNP) rs143383 in the 5' untranslated region (5'UTR) of growth and differentiation factor 5 (GDF5) is strongly associated with osteoarthritis (OA) and influences GDF5 allelic expression in vitro and in the joint tissues of OA patients. This effect is modulated in cis by another common SNP, also located within the 5'UTR, whilst a common SNP in the 3'UTR influences allelic expression independent of rs143383. DNA variants can be common, rare or extremely rare/unique. To therefore enhance our understanding of the allelic architecture of this very important OA susceptibility locus we sequenced the gene for potentially functional and novel rare variants.

METHOD

Using the Sanger method we sequenced GDF5 in 992 OA patients and 944 controls, with DNA changes identified by sequencing software. We encompassed the protein-coding region of the two GDF5 exons, both untranslated regions and approximately 100 bp of the proximal promoter of the gene.

RESULTS

We detected 13 variants. Six were extremely rare with minor allele frequencies (MAFs) of ≤ 0.0006. One is in a predicted transcription factor binding site in the GDF5 promoter whilst two substitute conserved amino acids. The remaining seven variants were common and are previously known variants, with MAFs ranging from 0.025 to 0.39. There was a complete absence of variants with frequencies in-between the extremely rare (n=6) and the common (n=7).

CONCLUSIONS

This is the first report of the deep sequencing of an OA susceptibility locus. The absence of rare variants informs us that within the regions of the gene that we have sequenced GDF5 does not harbour any novel variants that are able to contribute, at a population level, to the OA association signal mediated by rs143383 nor does it harbour, at a population level, any novel variants that can influence OA susceptibility independent of rs143383.