ABSTRACT
Genetic variation within intron 3 of CACNA1C surrounding a variable-number tandem repeat (VNTR) is associated with schizophrenia and bipolar disorder, but the causal variant(s) and their effects remain unclear. We fine-mapped the association at CACNA1C including this VNTR using sequences from 155 long-read genome assemblies. Across global populations, we found 7 alleles (called Types) of the CACNA1C VNTR where sequence differences within repeat units revealed distinct VNTR structure. Within the most common VNTR Type, a previously identified polymorphism of repeat unit composition (termed Variable Region 2) was in complete linkage disequilibrium with fine-mapped schizophrenia SNPs. Applying a genotyping strategy to GTEx data that capitalizes on sequence differences between repeat units, we found that the eQTL in brain tissues at Variable Region 2 had a similar effect size and significance as SNP eQTLs. We show that the risk allele of Variable Region 2 decreases CACNA1C gene expression. Our long-read-informed approach to genotype structurally complex VNTR alleles in large cohorts permits investigation of other variants missed by short-read sequencing. Our work suggests an effect on gene expression arising from sequence variation within a VNTR and provides a detailed characterization of new alleles at a flagship psychiatric GWAS locus.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This study was funded in part by grants from the National Institutes of Health
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
All data were publicly available at the start of the study: Long-read haplotype assemblies were downloaded from the HGSVC2 FTP site, the HPRC S3 bucket, and the T2T GitHub site. Summary statistics from the schizophrenia GWAS are available at the PGC site. Variant calls for HGSVC2 (Ebert et al. 2021) and 1000 Genomes individuals (Byrska-Bishop et al. 2022) are available through the 1000 Genomes FTP site. Ancestry information is available on IGSR. WGS alignment files for archaic human individuals can be found at the public FTP sites.
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Data Availability
All data produced in the present study are available upon reasonable request to the authors
ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data_collections/HGSVC2/release/v1.0/assemblies
https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html?prefix=working/
http://walters.psycm.cf.ac.uk/clozuk_pgc2.meta.sumstats.txt.gz
https://www.internationalgenome.org/data-portal/sample
http://cdna.eva.mpg.de/neandertal/altai/AltaiNeandertal/bam/
http://ftp.eva.mpg.de/neandertal/Chagyrskaya/BAM/
http://cdna.eva.mpg.de/denisova/alignments/
http://ftp.eva.mpg.de/neandertal/Vindija/bam/Pruefer_etal_2017/
