Project description:We used microarrays and a previously established linkage map to localize the genetic determinants of brain gene expression for a backcross family of lake whitefish species pairs (Coregonus sp.). Our goals were to elucidate the genomic distribution and sex-specificity of brain expression QTL (eQTL) and to determine the extent to which genes controlling transcriptional variation may underlie adaptive divergence in the recently evolved dwarf (limnetic) and normal (benthic) whitefish. We observed a sex-bias in transcriptional genetic architecture, with more eQTL observed in males, as well as divergence in genome location of eQTL between sexes. Hotspots of nonrandom aggregations of up to 32 eQTL in one location were observed. We identified candidate genes for species pair divergence involved with energetic metabolism, protein synthesis, and neural development based on co-localization of eQTL for these genes with eight previously identified adaptive phenotypic QTL and four previously identified outlier loci from a genome scan in natural populations. 88% of eQTL-phenotypic QTL co-localization involved growth rate and condition factor QTL, two traits central to adaptive divergence between whitefish species pairs. Hotspots co-localized with phenotypic QTL in several cases, revealing possible locations where master regulatory genes, such as a zinc finger protein in one case, control gene expression directly related to adaptive phenotypic divergence. We observed little evidence of co-localization of brain eQTL with behavioral QTL, which provides insight on the genes identified by behavioral QTL studies. These results extend to the transcriptome level previous work illustrating that selection has shaped recent parallel divergence between dwarf and normal lake whitefish species pairs and that metabolic, more than morphological differences appear to play a key role in this divergence. Keywords: eQTL mapping, gene expression, linkage mapping, adaptive radiation, Coregonus, microarrays

Project description:Genome-wide association studies (GWAS) have successfully identified 145 genomic regions that contribute to schizophrenia risk, but linkage disequilibrium (LD) makes it challenging to discern causal variants. Computational finemapping prioritized thousands of credible variants, ~98% of which lie within poorly characterized non-coding regions. To functionally validate their regulatory effects, we performed a massively parallel reporter assay (MPRA) on 5,173 finemapped schizophrenia GWAS variants in primary human neural progenitors (HNPs). We identified 440 variants with allelic regulatory effects (MPRA-positive variants), with 72% of GWAS loci containing at least one MPRA-positive variant. Transcription factor binding had modest predictive power for predicting the allelic activity of MPRA-positive variants, while GWAS association, finemap posterior probability, enhancer overlap, and evolutionary conservation failed to predict MPRA-positive variants. Furthermore, 64% of MPRA-positive variants did not exhibit eQTL signature, suggesting that MPRA could identify yet unexplored variants with regulatory potentials. MPRA-positive variants differed from eQTLs, as they were more frequently located in distal neuronal enhancers. Therefore, we leveraged neuronal 3D chromatin architecture to identify 273 genes that physically interact with MPRA-positive variants. These genes annotated by chromatin interactome displayed higher mutational constraints and regulatory complexity than genes annotated by eQTLs, recapitulating a recent finding that eQTL- and GWAS-detected variants map to genes with different properties. Finally, we propose a model in which allelic activity of multiple variants within a GWAS locus can be aggregated to predict gene expression by taking chromatin contact frequency and accessibility into account. In conclusion, we demonstrate that MPRA can effectively identify functional regulatory variants and delineate previously unknown regulatory principles of schizophrenia.

Project description:We used microarrays and a previously established linkage map to localize the genetic determinants of brain gene expression for a backcross family of lake whitefish species pairs (Coregonus sp.). Our goals were to elucidate the genomic distribution and sex-specificity of brain expression QTL (eQTL) and to determine the extent to which genes controlling transcriptional variation may underlie adaptive divergence in the recently evolved dwarf (limnetic) and normal (benthic) whitefish. We observed a sex-bias in transcriptional genetic architecture, with more eQTL observed in males, as well as divergence in genome location of eQTL between sexes. Hotspots of nonrandom aggregations of up to 32 eQTL in one location were observed. We identified candidate genes for species pair divergence involved with energetic metabolism, protein synthesis, and neural development based on co-localization of eQTL for these genes with eight previously identified adaptive phenotypic QTL and four previously identified outlier loci from a genome scan in natural populations. 88% of eQTL-phenotypic QTL co-localization involved growth rate and condition factor QTL, two traits central to adaptive divergence between whitefish species pairs. Hotspots co-localized with phenotypic QTL in several cases, revealing possible locations where master regulatory genes, such as a zinc finger protein in one case, control gene expression directly related to adaptive phenotypic divergence. We observed little evidence of co-localization of brain eQTL with behavioral QTL, which provides insight on the genes identified by behavioral QTL studies. These results extend to the transcriptome level previous work illustrating that selection has shaped recent parallel divergence between dwarf and normal lake whitefish species pairs and that metabolic, more than morphological differences appear to play a key role in this divergence. Keywords: eQTL mapping, gene expression, linkage mapping, adaptive radiation, Coregonus, microarrays The objective of this study was to elucidate the genomic distribution and sex-specificity of brain eQTL in dwarf and normal lake whitefish. Dissected brain tissue (250-350 mg) was sampled for 55 individuals from a hybrid x dwarf backcross mapping family. We used a loop design (YANG and SPEED 2002; CHURCHILL 2002) to maximize the number of sampled meioses. Each of 55 samples was technically replicated on two distinct slides, while performing dye swapping (Cy3 and Alexa) to estimate the dye intensity variation bias. After correcting for local background, raw intensity values were both log2 transformed and normalized using the regional LOWESS method implemented in the R/MANOVA software (KERR et al. 2000). We used a previously generated linkage map based on the same backcross individuals for which gene expression was measured. eQTL mapping was performed with QTL Cartographer.

Dataset Information

Integrative QTL analysis of gene expression and chromatin accessibility identifies multi-tissue patterns of genetic regulation (RNA-seq)

Dataset's files

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets