Find research datasets worth reusing
Search datasets from major research repositories and use ShareScore to quickly assess how well each record supports discovery, access, and reuse.
176
datasets available to search
ShareScore release 0.9.0
Dataset results
176 results for “QTL”
QTL summary statistics from the DIRECT consortium
<p>These are the complete summary statistics for DIRECT genotype-phenotypes associations (QTLs). The project performed genotype-phenotype associations for gene expression (RNAseq), targeted proteins (Olink), targeted metabolites (Biocrates) and untargeted metabolites (Metabolon) derived from 3,029 blood and plasma samples from the DIRECT cohort. This submission includes supplementary files and nominal pvalues (as uncorrected pvalues) for all associations included in the manuscript. Trans associations included are typically limited to pvalues <1e-04. Network tables are also included, with information to load and use Cytoscape to visualize them. This is version 2, some files were missing on version 1.</p>
QTL Mapping for Resistance to Cankers Induced by Pseudomonas syringae pv. actinidiae (Psa) in a Tetraploid Actinidia chinensis Kiwifruit Population
<p>Raw Illumina R1 sequence reads for individual plants genotyped for the study entitled "QTL Mapping for Resistance to Cankers Induced by <em>Pseudomonas syringae</em> pv. <em>actinidiae</em> (Psa) in a Tetraploid <em>Actinidia chinensis</em> Kiwifruit Population" accepted in MDPI Pathogen journals, Special issue <em>"</em><em>Pseudomonas syringae</em> Species Complex"</p>
QTL Mapping for Resistance to Cankers Induced by Pseudomonas syringae pv. actinidiae (Psa) in a Tetraploid Actinidia chinensis Kiwifruit Population
<p>Raw Illumina R2 sequence reads for individual plants genotyped for the study entitled "QTL Mapping for Resistance to Cankers Induced by <em>Pseudomonas syringae</em> pv. <em>actinidiae</em> (Psa) in a Tetraploid <em>Actinidia chinensis</em> Kiwifruit Population" accepted in MDPI Pathogen journals, Special issue <em>"</em><em>Pseudomonas syringae</em> Species Complex"</p>
Strawberry volatile organic compounds metabolomic data and QTL study
<p>This dataset contains the supplementary materials of the publication "Multivariate QTL approach reveals a major regulator of terpenoid production and other volatiles in strawberry" of the same authors. In this study we extracted volatile organic compounds from several strawberry samples and analysed their identity and abundance. We used this volatile data to perform an extensive multivariate QTL study, the results of which can be found in this dataset.</p> <p>All analysis, results and figures can be reproduced using the folder included in the <strong>supplementary data 1</strong>. If you want to reproduce part or all of our analysis, only download sup data 1. </p> <p>If you only need one of our results or data table you can download them individually:</p> <ul> <li>Sup data 2: abundances of volatile organic compounds from a biparental and diverse panel (GWAS) population.</li> <li>Sup data 3 and 4: p-value tables for all metabolites as well as multivariate traits (see publication for more information).</li> <li>Sup table 1: Summary of metabolite abundances and heritabilities across both populations.</li> <li>Sup tables 2 and 3: significant QTL signals for each trait individually and summarised per QTL locus.</li> <li>Sup table 4: previously reported VOC QTLs in strawberry, with positions imputed in the Royal Royce genome.</li> <li>Sup table 5: metadata about all the identified compounds on this and previous studies.</li> <li>Sup table 6: SNP array positions imputed in the "Royal Royce" genome assembly.</li> <li>Sup table 7: Number of markers per chromosome in this analysis.</li> </ul> <h3>Update 2025</h3> <p>We updated the underlying code and datasets to reflect several revisions made to this work. Most notably, the QTL results have been reworked. They now do not include Blink or FarmCPU results (only mixed model results, obtained through statgenGWAS). Additionally, heritability estimations, QQ-plots and other figures have been added to the reproducible results code.</p>
Mapping of promoter usage QTL using RNA-seq data reveals their contributions to complex traits
<p>Supporting data associated with the manuscript "Mapping of promoter usage QTL using RNA-seq data reveals their contributions to complex traits". These include:</p> <ul> <li><strong>StringTie2.438.GEUVADIS.gtf.gz</strong> - Assembled transcriptome using RNA-seq from 438 samples from the GEUVADIS project by StringTie2</li> <li><strong>puQTL_nominal_output.txt.gz</strong> - puQTL results from QTLtools nominal pass</li> <li><strong>puQTL_conditional_output.txt.gz</strong> - puQTL results from QTLtools conditional pass</li> <li><strong>eQTL_nominal_output.txt.gz</strong> - eQTL results from QTLtools nominal pass</li> <li><strong>eQTL_conditional_output.txt.gz</strong> - eQTL results from QTLtools conditional pass</li> </ul> <p>See <a href="https://qtltools.github.io/qtltools/">https://qtltools.github.io/qtltools/</a> for column descriptions of QTL mapping results.</p>
Gene expression QTL mapping in stimulated iPSC-derived macrophages provides insights into common complex diseases.
<p>Many disease-associated variants are thought to be regulatory but are not present in existing catalogues of expression quantitative trait loci (eQTL). We hypothesise that these variants may regulate expression in specific biological contexts, such as stimulated immune cells. Here, we used human iPSC-derived macrophages to map eQTLs across 24 cellular conditions. We found that 76% of eQTLs detected in at least one stimulated condition were also found in naive cells. The percentage of response eQTLs (reQTLs) varied widely across conditions (3.7% - 28.4%), with reQTLs specific to a single condition being rare (1.11%). Despite their relative rarity, reQTLs were overrepresented (p=0.05, Fisher's exact test) among disease-colocalizing eQTLs. We nominated an additional 21.7% of disease effector genes at GWAS loci via colocalization of reQTLs, with 38.6% of these not found in the Genotype–Tissue Expression (GTEx) catalogue. Our study highlights the diversity of genetic effects on expression and demonstrates how condition-specific regulatory variation can enhance our understanding of common disease risk alleles.</p>
GWAS and GTEx QTL integration
<p># Data usage policy</p> <p>When using this data, you must acknowledge the source by citing the publication "Widespread dose-dependent effects of RNA expression and splicing on complex diseases and traits" (https://doi.org/10.1101/814350).</p> <p> </p> <pre><em># GTEx GWAS integration </em> This package contains the application of several GWAS-QTL integration methods. The results were analyzed in [this preprint](<em>https://www.biorxiv.org/content/10.1101/814350v1</em>) about GTEx v8 application to several GWAS traits. <em>``` </em><em>. </em><em>|-- colocalization </em><em>| |-- coloc </em><em>| | `-- coloc_enloc_priors_eqtl.tar.gz </em><em>| |-- enloc </em><em>| | |-- enloc_eqtl_eur.tar.gz </em><em>| | `-- enloc_sqtl_eur.tar.gz </em><em>| `-- eur_ld.bed.gz </em><em>|-- prediction_models </em><em>| |-- gtex_v8_expression_mashr_snp_smultixcan_covariance.txt.gz </em><em>| |-- gtex_v8_splicing_mashr_snp_smultixcan_covariance.txt.gz </em><em>| |-- mashr_eqtl.tar </em><em>| `-- mashr_sqtl.tar </em><em>|-- smr </em><em>| |-- SMR_gtex_v8_README.txt </em><em>| `-- SMRresults_GTEx_v8_peQTL5e-08.tar.gz </em><em>|-- smultixcan </em><em>| |-- smultixcan_eqtl.tar.gz </em><em>| `-- smultixcan_sqtl.tar.gz </em><em>`-- spredixcan </em><em> |-- spredixcan_eqtl.tar.gz </em><em> `-- spredixcan_sqtl.tar.gz </em> <em> ``` </em><em> </em>You can uncompress gzipped tarball packages <em>`*.tar.gz` </em>in a UNIX command line with an instruction such as: <em>```bash </em><em>tar -xzvpf smultixcan_eqtl.tar.gz </em><em>``` </em>, and the tar packages (<em>`*.tar`</em>) with an analogous instruction: <em>```bash </em><em>tar -xvpf mashr_eqtl.tar </em><em>``` </em> <em>## Preliminaries </em> <strong>**</strong>Finemapping<strong>** </strong>results are contained in a separate release due to size constraints. GWAS summary statistics for 114 traits were harmonized and imputed to GTEx v8 variants with MAF>0.01 using only european samples. (summary imputation software [here](<em>https://github.com/hakyimlab/summary-gwas-imputation</em>)). Some of the following analyses used the full set of 114 traits, while some focused only on 87 traits whose imputed associations showed no deflation (the imputation algorithm is conservative, and studies with too few available variants have a depleted distribution of association p-values after imputation). The harmonized and imputed GWAS summary statistics are contained in a separate release due to size constraints. For completeness' sake, the imputed summary statistics look like: <em>``` </em><em>variant_id panel_variant_id chromosome position effect_allele non_effect_allele current_build frequency sample_size zscore pvalue effect_size standard_error imputation_status n_cases </em><em>rs554008981 chr1_13550_G_A_b38 chr1 13550 A G hg38 0.017316017316017316 336474 -2.2919929353647097 0.021906050841240293 NA NA imputed NA </em><em>rs201055865 chr1_14671_G_C_b38 chr1 14671 C G hg38 0.012987012987012988 336474 -0.9559192804440632 0.33911301727494103 NA NA imputed NA </em><em>... </em><em>``` </em> The GWAS were split in approximately independent LD regions (Berisa-Pickrell)/ GWAS regions are defined in <em>`eur_ld.bed.gz` </em>(note that a few of them are ill-defined in hg38 and where ignored; only completely defined regions were used). <em>## Colocalization </em> <em>### Enloc </em> ENLOC ([see fotware here](<em>https://github.com/xqwen/integrative</em>)) was run for sQTLs and eQTLs using individuals of european ancestry and DAP-G QTL enrichment results on 87 traits. Result files are included in <em>`enloc_eqtl_eur.tar.gz` </em>and <em>`enloc_sqtl_eur.tar.gz` </em>Each file contains a particular tissue-trait combination. Each row details colocalization between a GWAS region (Berisa-Pickrell) and gene's or intron's cis-window. A region might overlap multiple genes/introns or viceversa. Each ENLOC file contains the following columns: <strong>* </strong>gwas_locus: GWAS LD region <strong>* </strong>molecular_qtl_trait: gene or intron <strong>* </strong>locus_gwas_pip: posterior inclusion probability of variants in the GWAS LD region <strong>* </strong>locus_rcp: regional colocalization probability (main colocalization measure) <strong>* </strong>lead_coloc_SNP: snp with highest RCP <strong>* </strong>lead_snp_rcp: rcp of the lead coloc snp <em>### Coloc </em> Coloc ([see software here](<em>https://cran.r-project.org/web/packages/coloc/index.html</em>)) was run using prior probabilities estimated from QTL enrichment of GWAS variants (computed via ENLOC). Results for eQTL are available in <em>`coloc_enloc_priors_eqtl.tar.gz`</em>. Each file contains results for a trait-tissue combination. Columns are: <strong>* </strong>gene_id: gene or intron id <strong>* </strong>p0: probability that neither QTL nor GWAS contain a causal variant <strong>* </strong>p1: probability that only GWAS contains a causal variant <strong>* </strong>p2: probability that only QTL has a causal variant <strong>* </strong>p3: probability that GWAS and QTL have a causal variant and it's distinct <strong>* </strong>p4: probability that GWAS and QTL have a causal variant and it's the same (main colocalization measure) <em>## PrediXcan </em> <em>`mashr_eqtl.tar` </em>and <em>`mashr_sqtl.tar` </em>contain prediction models (trained on expression or splicing data respectively, for 49 GTEx tissues) and LD compilations to be used with PrediXcan, S-PrediXcan, MultiXcan and S-MultiXcan. For every tissue, the <em>`mashr_{tissue}.db` </em>file is a SQLite file with the prediction model definitions. <em>`mashr_{tissue}.txt.gz` </em>is a gzipped-text file with the upper triangular matrices of covariance between snps within a gene/intron prediction model. Many variants in these models don't have an rsid. To fully leverage the information in these models, it is advised to at least harmonize to GTEx variants, and if possible impute as we did [here](<em>https://github.com/hakyimlab/summary-gwas-imputation</em>). <em>### S-PrediXcan </em> S-PrediXcan was run for the 114 harmonized and imputed traits, on eQTL and sQTL mashr prediction models. All of the GWAS traits had the same format, so that the following format parameters were used with S-PrediXcan: <em>``` </em><em>--snp_column panel_variant_id --effect_allele_column effect_allele --non_effect_allele_column non_effect_allele --zscore_column zscore \ </em><em>--keep_non_rsid --additional_output --model_db_snp_key varID \ </em><em>``` </em> Each file is a CSV, with each row containing a gene/intron association at a given trait-tissue combination: <strong>* </strong>gene: ENSEMBLE ID or intron id <strong>* </strong>gene_name: HUGO name or intron id <strong>* </strong>zscore: predicted association z-score <strong>* </strong>effect_size: estimated effect size <strong>* </strong>pvalue: association p-value <strong>* </strong>var_g: estimated variance of predicted expression or splicing <strong>* </strong>pred_perf_r2: prediction model cross-validated performance <strong>* </strong>pred_perf_pval: prediction model cross-validated performance <strong>* </strong>pred_perf_qval: deprecated, empty field left for compatibility <strong>* </strong>n_snps_used: number of snps in the intersection of GWAS and model <strong>* </strong>n_snps_in_cov: number of snps in the LD compilation <strong>* </strong>n_snps_in_model: number of snps in the model <strong>* </strong>best_gwas_p: smallest p-value acros GWAS snps used in this model <strong>* </strong>largest_weight: largest prediction model weight <em>### S-Multixcan </em> S-MultiXcan results were generated from the above S-PrediXcan results. Each fiel contains multi-tissue associations for a given trait: <strong>* </strong>gene: ENSEMBLE ID or intron id <strong>* </strong>gene_name: HUGO name or intron id <strong>* </strong>pvalue: multi-tissue association p-value <strong>* </strong>n: number of models avialble for this gene/intron <strong>* </strong>n_indep: number of independent components of variation in predicted expression/splicing (surviving principal components) <strong>* </strong>p_i_best: highest single-tissue p-value (S-PrediXcan) <strong>* </strong>t_i_best: tissue of highest p-value <strong>* </strong>p_i_worst: lowest single-tissue p-value (S-PrediXcan) <strong>* </strong>t_i_worst: tissue of lowest p-value <strong>* </strong>eigen_max: maximum eigenvalue of SVD <strong>* </strong>eigen_min: minimum eigenvalue of SVD <strong>* </strong>eigen_min_kept: smallest eigenvalue retained after discarding smallest variations <strong>* </strong>z_min: minimum single-tissue z-score <strong>* </strong>z_max: maximum single-tissue z-score <strong>* </strong>z_mean: mean single-tissue zscre <strong>* </strong>z_sd: standard deviation of the single-tissue z-scores <strong>* </strong>tmi: trace of M * M_i where M is predicted expression/splicing covariance across tissues for a gene, and M_i is its SVD pseudo-inverse <strong>* </strong>status: computation status, 0 if no errors <em>## SMR </em> See <em>`SMR_gtex_v8_README.txt` </em>for details.</pre> <p> </p> <p> </p> <p># Disclaimer</p> <p>The data is provided "as is", and the authors assume no responsibility for errors or omissions. <br> The User assumes the entire risk associated with its use of these data. <br> The authors shall not be held liable for any use or misuse of the data described and/or contained herein. <br> The User bears all responsibility in determining whether these data are fit for the User's intended use. </p> <p>The information contained in these data is not better than the original sources from which they were derived,<br> and both scale and accuracy may vary across the data set. <br> These data may not have the accuracy, resolution, completeness, timeliness, or other characteristics<br> appropriate for applications that potential users of the data may contemplate. <br> <br> The user is responsible to comply with any data usage policy from the original GWAS studies;<br> refer to the list of traits described [here](https://www.biorxiv.org/content/10.1101/814350v1)<br> to identify their respective Consortia's requirements.</p> <p><br> THE DATA IS PROVIDED WITHOUT WARRANTY OF ANY KIND,<br> EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,<br> FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.<br> IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,<br> WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,<br> OUT OF OR IN CONNECTION WITH THE DATA OR THE USE OR OTHER DEALINGS IN THE DATA.</p>
Full expression and splicing QTL summary statistics for three human spinal cord segments
<p>This dataset is part of the manuscript: <strong>"Integrative genetic analysis of the amyotrophic lateral sclerosis spinal cord implicates glial activation and suggests new risk genes</strong>" by Jack Humphrey, et al.</p> <p>The files below contain nominal and permuted quantitative trait loci (QTL) associations between common genetic variants derived from whole genome sequencing and either gene expression or splicing phenotypes generated from RNA-seq of post-mortem spinal cord sections. All QTLs were mapped with TensorQTL.</p> <p>Nominal association files are coordinate sorted, bgzip-compressed and tabix-indexed tab-separated variable files. Top association files are gzip-compressed tab-separated variable files.</p> <p>Sample sizes for each dataset are as follows:</p> <p>Cervical spinal cord n = 216.</p> <p>Lumbar spinal cord n = 197.</p> <p>Thoracic spinal cord n = 68.</p> <p>Description of files:<br> CervicalSpinalCord_expression_peer30_gene.cis_qtl.txt.gz - <strong>Top per-gene </strong>eQTL summary statistics from cervical spinal cord <br> CervicalSpinalCord_expression_peer30_gene.cis_qtl_nominal_tabixed.tsv.gz - <strong>Full nominal</strong> eQTL summary statistics from cervical spinal cord <br> CervicalSpinalCord_splicing_peer15_cluster.cis_qtl.txt.gz - <strong>Top per-cluster </strong>sQTL summary statistics from cervical spinal cord<br> CervicalSpinalCord_splicing_peer15_gene.cis_qtl_nominal_tabixed.tsv.gz - <strong>Full nominal</strong> sQTL summary statistics from cervical spinal cord<br> LumbarSpinalCord_expression_peer30_gene.cis_qtl.txt.gz - <strong>Top per-gene </strong>eQTL summary statistics from lumbar spinal cord<br> LumbarSpinalCord_expression_peer30_gene.cis_qtl_nominal_tabixed.tsv.gz -<strong> Full nominal</strong> eQTL summary statistics from lumbar spinal cord<br> LumbarSpinalCord_splicing_peer10_cluster.cis_qtl.txt.gz - <strong>Top per-cluster </strong>sQTL summary statistics from lumbar spinal cord<br> LumbarSpinalCord_splicing_peer10_gene.cis_qtl_nominal_tabixed.tsv.gz - <strong>Full nominal</strong> sQTL summary statistics from lumbar spinal cord<br> ThoracicSpinalCord_expression_peer10_gene.cis_qtl.txt.gz - <strong>Top per-gene </strong>eQTL summary statistics from thoracic spinal cord<br> ThoracicSpinalCord_expression_peer10_gene.cis_qtl_nominal_tabixed.tsv.gz - <strong>Full nominal</strong> eQTL summary statistics from thoracic spinal cord<br> ThoracicSpinalCord_splicing_peer5_cluster.cis_qtl.txt.gz - <strong>Top per-cluster </strong>sQTL summary statistics from thoracic spinal cord<br> ThoracicSpinalCord_splicing_peer5_gene.cis_qtl_nominal_tabixed.tsv.gz - <strong>Full nominal</strong> sQTL summary statistics from thoracic spinal cord</p> <p>als_snps_maf0.01_alleles.tsv.gz<strong> - Allele information</strong> for all SNPs tested in the eQTL and sQTL analysis</p> <p>Table columns are formatted as follows:</p> <p>Nominal QTL results include all SNP-gene pairs tested using either a 1Mb window from each side of the transcription start site (TSS) of the gene (eQTLs) or 100kb either side of the middle of the splicing cluster (sQTLs). Table columns are formatted as follows:</p> <ul> <li>phenotype_id - ensembl ID of the gene tested (GENCODE v30) for eQTL, or the ID of the splice junction for sQTLs, splicing junction position, splicing cluster id, and gene Ensembl id (GENCODE v30),separated by a colon</li> <li>variant_id - SNP tested for association (rsid or chr:position:ref:alt)</li> <li>tss_distance - distance of the SNP to the gene transcription start site (TSS)</li> <li>maf - minor allele frequency in cohort</li> <li>ma_samples - number of samples carrying the minor allele</li> <li>ma_count - total number of minor alleles across individuals</li> <li>pval_nominal - nominal <em>P</em>-value from linear regression</li> <li>slope - slope of the linear regression</li> <li>slope_se - standard error of the slope</li> </ul> <p>Permuted QTL results include only the top SNP-gene association for each gene (eQTL) or cluster (sQTL). Table columns are formatted as follows:</p> <ul> <li>phenotype_id - ensembl ID of the gene tested (GENCODE v30), or the coordinates of the cluster, formatted as coordinate, splicing cluster id, and gene Ensembl id (GENCODE v30), separated by a colon</li> <li>num_var - total number of variants tested in <em>cis</em></li> <li>beta_shape1 - first parameter value of the fitted beta distribution</li> <li>beta_shape2 - second parameter value of the fitted beta distribution</li> <li>true_df - effective degrees of freedom the beta distribution approximation</li> <li>pval_true_df - empirical <em>P</em>-value for the beta distribution approximation</li> <li>variant_id - ID of the top variant (rsid or chr:position:ref:alt)</li> <li>tss_distance - distance of the SNP to the gene transcription start site (TSS)</li> <li>ma_samples - number of samples carrying the minor allele</li> <li>ma_count - total number of minor alleles across individuals</li> <li>maf -minor allele frequency in MiGA cohort</li> <li>ref_factor - flag indicating if the alternative allele is the minor allele in the cohort (1 if AF <= 0.5, -1 if not)</li> <li>pval_nominal - nominal <em>P</em>-value from linear regression</li> <li>slope - slope of the linear regression</li> <li>slope_se - standard error of the slope</li> <li>pval_perm - first permutation <em>P</em>-value directly obtained from the permutations with the direct method</li> <li>pval_beta - second permutation <em>P</em>-value obtained via beta approximation. This is the one to use for downstream analysis</li> <li>qval - Storey q-value derived from pval_beta (FDR adjusted)</li> <li>pval_nominal_threshold - nominal <em>P</em>-value threshold for calling a variant-gene pair significant for the gene</li> </ul> <p>Allele Information for each variant:</p> <ul> <li>CHROM - chromosome position of the variant</li> <li>POS - position of the variant in the chromosome</li> <li>REF - reference allele (GRCh38)</li> <li>ALT - alternative allele (this is the <strong>effect</strong> allele in the eQTL and sQTL analysis)</li> <li>ID - variant id (rsid or chr:position:ref:alt)</li> </ul>
QTL mapping and transcriptome analysis of Sclerotinia-resistance in the wild cabbage species Brassica oleracea var. villosa [Main code]
<p>This is the main code supplement for my computational analysis for the manuscript: "QTL mapping and transcriptome analysis of Sclerotinia-resistance in the wild cabbage species <em>Brassica oleracea </em>var<em>. villosa".</em> The main code is availabe in separate html-files. DOI will be added if available.</p>
Pooled DNA sequencing to identify SNPs associated with a major QTL for bacterial wilt resistance in Italian ryegrass (Lolium multiflorum Lam.)
<p>We used pooled DNA sequencing to characterize a major QTL for bacterial wilt resistance of Italian ryegrass and to develop inexpensive sequence-based markers to efficiently target resistance alleles for marker-assisted recurrent selection. From the mapping population segregating for the QTL, DNA of 44 of the most resistant and 44 of the most susceptible F<sub>1</sub> individuals were pooled and sequenced using the Illumina HiSeq2000 platform. Allele frequencies of 18 x 10<sup>6</sup> single nucleotide polymorphisms (SNP) were determined in the resistant and susceptible pool. A total of 271 SNPs on 140 scaffold sequences of the reference parental genome showed significantly different allele frequencies in both pools. We converted 44 selected SNPs to KASP markers, genetically mapped these proximal to the major QTL and thus validated their association with bacterial wilt resistance.</p>
Epistatic QTL for yield heterosis in tomato
<p>Controlled population development and genome-wide association studies have proven powerful in uncovering genes and alleles underlying complex traits. An underexplored dimension of such studies is the phenotypic contribution of non-additive interactions between quantitative trait loci (QTL). Capturing such epistasis in a genome-wide manner requires very large populations to represent replicated combinations of loci whose interactions determine phenotypic outcomes. Here, we dissect epistasis using a densely genotyped population of 1400 backcross-inbred lines (BILs) between a modern processing tomato inbred (<em>Solanum lycopersicum</em>) and the Lost Accession of a distant, green-fruited, drought-tolerant wild species, <em>Solanum pennellii</em>. The homozygous BILs, each representing an average of 11 introgressions, and their hybrids with the recurrent parents were phenotyped for tomato yield components. Population-wide mean yield of the BILs was less than 50% of that of their hybrids (BILH), and correspondingly, homozygous introgressions across the genome reduced yield relative to recurrent parent while several QTL of the BILHs independently improved productivity. Analysis of two QTL scans showed 32 cases of less than additive and 48 cases of more than additive interactions. Strikingly, one epistatic interaction involving <em>S. pennellii</em> QTL on chromosomes 1 and 7, which independently did not affect yield, increased fruit yield by 20-50% in the double introgression hybrid grown in irrigated and dry fields over a period of three years and varied genetic backgrounds. Our work demonstrates the power of large, interspecific controlled population development to uncover hidden QTL phenotypes, and how rare epistatic interactions can improve crop productivity via heterosis.</p>
Molecular mapping of quantitative trait loci (QTL) for resistance to early blight in tomato
<p>Molecular mapping of quantitative trait loci (QTL) for resistance to early blight in tomato (1135884)</p>
isoMIGA: Expression and Splicing QTL Summary Statistics in Human Microglia
<p>https://github.com/RajLabMSSM/isoMiGA </p> <p>Full associations and top association summary statistics for QTLs mapped in a multi-ethnic meta-analysis of human microglia, as part of the isoform-centric Microglia Genomic Atlas (isoMiGA) project.</p> <p>Sample size = 555 samples from 391 unique donors</p> <p>6 cohorts (Raj MFG, Raj STG, Raj, SVZ, Raj, THA, Roussos, Gaffney) meta-analysed using the linear mixed model random-effects meta-analysis software mmQTL (PMID: 35058635).</p> <p>Each file has the following naming convention:</p> <p>{REFERENCE}_{PHENOTYPE}_{ASSOC}.tsv.gz</p> <p><strong>References</strong></p> <p>The following two transcript references were used for mapping phenotypes:</p> <p>1. GENCODE - GENCODE v38 comprehensive transcripts (https://www.gencodegenes.org/human/release_38.html)</p> <p>2. Union - a union of all GENCODE v38 transcripts with an additional 35,879 novel transcripts identified from long-read RNA-seq in 30 human microglia samples. GTF to be uploaded in Zenodo.</p> <p><strong>Phenotypes</strong></p> <p>The following phenotypes were tested for genetic association:</p> <ul> <li>expression: total gene expression of each gene following voom normalization of read counts.</li> <li>transcript: transcript usage - each transcript expression normalized by TPM divided by the total TPM for that gene</li> <li>leafcutter: junction usage using the Leafcutter framework to count and cluster intron-splicing junction reads. Each phenotype is the relative usage of the intron against the total count of the introns within that cluster.</li> <li>SUPPA_A3: alternate 3' splice site usage, as identified by SUPPA2 from transcript TPMs</li> <li>SUPPA_A5: alternate 5' splice site usage</li> <li>SUPPA_AF: alternate first exon usage</li> <li>SUPPA_AL: alternate last exon usage</li> <li>SUPPA_RI: intron retention usage</li> <li>SUPPA_SE: exon skipping usage</li> </ul> <p>All phenotype matrices were scaled and centred and then quantile normalized.</p> <p><strong>Associations</strong></p> <p>Top associations (top_assoc.tsv.gz) list the SNP-feature pair with the lowest adjusted P-value (qval) for that feature.</p> <p>Full associations (full_assoc.tsv.gz) list all tested SNP-feature pairs.</p> <p><strong>Data dictionary</strong></p> <p>The columns of the two association files only differ by the presence of the qval column in the top associations.</p> <p>feature: the phenotype being tested</p> <p>variant_id: the genetic variant being tested</p> <p>chr: chromosome</p> <p>pos: position (hg38)</p> <p>ref: reference allele</p> <p>alt: alternate allele</p> <p>Allele: the effect allele that the beta is relative to</p> <p>beta_tissue_0: Gaffney cohort beta</p> <p>sd_tissue_0: Gaffney cohort standard error</p> <p>z_tissue_0: Gafney cohort Z-score</p> <p>beta_tissue_1: Roussos cohort beta</p> <p>sd_tissue_1: Roussos cohort standard error</p> <p>z_tissue_1: Roussos cohort Z-core</p> <p>beta_tissue_2: Raj MFG cohort beta</p> <p>sd_tissue_2: Raj MFG cohort standard error</p> <p>z_tissue_2: Raj MFG cohort Z-score</p> <p>beta_tissue_3: Raj STG cohort beta</p> <p>sd_tissue_3: Raj STG cohort standard error</p> <p>z_tissue_3: Raj STG cohort Z-score</p> <p>beta_tissue_4: Raj SVZ cohort beta</p> <p>sd_tissue_4: Raj SVZ cohort standard error</p> <p>z_tissue_4: Raj SVZ cohort Z-score</p> <p>beta_tissue_5: Raj THA cohort beta</p> <p>sd_tissue_5: Raj THA cohort standard error</p> <p>z_tissue_5: Raj THA cohort Z-score</p> <p>fixed_beta: Fixed effect meta-analysis estimate of the beta</p> <p>fixed_sd: Fixed effect meta-analysis standard error of the beta</p> <p>fixed_z: Fixed effect meta-analysis Z-score</p> <p>Random_Z: Random effect meta-analysis Z-score</p> <p>Fixed_P: Fixed effect meta-analysis P-value</p> <p>Random_P: Random effect meta-analysis P-value</p> <p>Fixed_bonf: Fixed effect meta-analysis P-value adjusted for the number of variants tested in that feature (Bonferroni)</p> <p>Random_bonf: Random effect meta-analysis P-value adjusted for the number of variants tested in that feature (Bonferroni)</p> <p>Fixed_FDR: Fixed effect meta-analysis P-value adjusted for the number of variants tested in that feature (FDR)</p> <p>Random_FDR: Random effect meta-analysis P-value adjusted for the number of variants tested in that feature (FDR)</p> <p>qval: Random effect meta-analysis P-value adjusted for the number of variants tested in that feature (FDR) and for the number of features tested in the dataset (Storey's q value)</p> <p> </p> <p><br> </p> <p> </p> <p> </p> <p> </p> <p> </p> <p> </p> <p> </p>
Integrative QTL mapping and selection signatures in Groningen White Headed cattle inferred from whole-genome sequences
<p>Here, we aimed to identify and characterize genomic regions that differ between Groningen White Headed (GWH) breed and other cattle, and in particular to identify candidate genes associated with coat color and/or eye-protective phenotypes. Firstly, whole genome sequences of 170 animals from eight breeds were used to evaluate the genetic structure of the GWH in relation to other cattle breeds by carrying out principal components and model-based clustering analyses. Secondly, the candidate genomic regions were identified by integrating the findings from: a) a genome-wide association study using GWH, other white headed breeds (Hereford and Simmental), and breeds with a non-white headed phenotype (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel, Dutch Belted, and Holstein Friesian); b) scans for specific signatures of selection in GWH cattle by comparison with four other Dutch traditional breeds (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel and Dutch Belted) and the commercial Holstein Friesian; and c) detection of candidate genes identified via these approaches. The alignment of the filtered reads to the reference genome (ARS-UCD1.2) resulted in a mean depth of coverage of 8.7X. After variant calling, the lowest number of breed-specific variants was detected in Holstein Friesian (148,213), and the largest in Deep Red (558,909). By integrating the results, we identified five genomic regions under selection on BTA4 (70.2–71.3 Mb), BTA5 (10.0–19.7 Mb), BTA20 (10.0–19.9 and 20.0–22.7 Mb), and BTA25 (0.5–9.2 Mb). These regions contain positional and functional candidate genes associated with retinal degeneration (e.g., <em>CWC27</em> and <em>CLUAP1</em>), ultraviole<em>t</em> protection (e.g., <em>ERCC8</em>), and pigmentation (e.g. <em>PDE4D</em>) which are probably associated with the GWH specific pigmentation and/or eye-protective phenotypes, e.g. Ambilateral Circumocular Pigmentation (ACOP). Our results will assist in characterizing the molecular basis of GWH phenotypes and the biological implications of its adaptation.</p>
Epistatic QTL for yield heterosis in tomato
Open the record for dataset details and reuse information.
Data from: Major QTL controls adaptation to serpentine soils in Mimulus guttatus
Spatially varying selection is a critical driver of adaptive differentiation. Yet, there are few examples where the fitness effects of naturally segregating variants that contribute to local adaptation have been measured in the field. Plant adaptation to harsh soil habitats provides an ideal study system for investigating the genetic basis of local adaptation. The work presented here identifies a major locus underlying adaptation to serpentine soils in Mimulus guttatus and estimates the strength of selection on this locus in native field sites. Reciprocal transplant and common garden studies show that serpentine and non-serpentine populations of M. guttatus differ in their ability to survive on serpentine soils. We directly mapped these field survival differences by performing a bulk segregant analysis with F2 survivors from a field transplant study and identify a single QTL where individuals that are homozygous for the non-serpentine allele do not survive on serpentine soils. Genotyping the survivors from an independent mapping population reveals that this same QTL controls serpentine tolerance in a second, geographically distant population. Finally, we show that this QTL controls tolerance to soil properties, as opposed to some other aspect of the field sites that may differ, by performing a lab-based common-garden experiment in native serpentine soils that replicates the survival differences observed in the field. These results indicate that despite the myriad chemical and physical challenges that plants face in serpentine habitats, adaptation to these soils in M. guttatus has a simple genetic basis.
Dissecting the genetic basis of variation in Drosophila sleep using a multiparental QTL mapping resource
There is considerable variation in sleep duration, timing and quality in human populations, and sleep dysregulation has been implicated as a risk factor for a range of health problems. Human sleep traits are known to be regulated by genetic factors, but also by an array of environmental and social factors. These uncontrolled, non-genetic effects complicate powerful identification of the loci contributing to sleep directly in humans. The model system, Drosophila melanogaster, exhibits a behavior that shows the hallmarks of mammalian sleep, and here we use a multitiered approach, encompassing high-resolution QTL mapping, expression QTL data, and functional validation with RNAi to investigate the genetic basis of sleep under highly controlled environmental conditions. We measured a battery of sleep phenotypes in >750 genotypes derived from a multiparental mapping panel and identified several, modest-effect QTL contributing to natural variation for sleep. Merging sleep QTL data with a large head transcriptome eQTL mapping dataset from the same population allowed us to refine the list of plausible candidate causative sleep loci. This set includes genes with previously characterized effects on sleep and circadian rhythms, in addition to novel candidates. Finally, we employed adult, nervous system-specific RNAi on the Dopa decarboxylase, dyschronic, and timeless genes, finding significant effects on sleep phenotypes for all three. The genes we resolve are strong candidates to harbor causative, regulatory variation contributing to sleep.
Data from: Validation of grain yield QTL from soft winter wheat using a CIMMYT spring wheat panel
Validation of quantitative trait loci (QTLs) is an essential step in marker-assisted breeding. The objectives of this study were to validate grain yield (GY) QTLs previously identified in soft red winter wheat (Triticum aestivum L.) through biparental and association mapping using the spring wheat association mapping initiative (WAMI) panel from CIMMYT, Mexico, and to identify allele combinations of the validated QTLs that resulted to the highest GY. Linked single-nucleotide polymorphisms for IWA3560 (3A), IWA1818 (4B), and IWA755 (6B) were significantly associated (P < 0.001) with GY, grain number, and thousand-grain weight in the WAMI. Lines possessing the favorable allele for the QTL at the 3A, 4B, and 6B loci (ACG allele combination) validated on the WAMI had the highest mean GY at 4.55 t ha−1, but three other haplotypes (ACA, GCA, and GCG) differing by one or two alleles in the validated QTL regions were not significantly different. These results validate GY QTLs across winter and spring wheat through genome-wide association analysis and further demonstrate the potential for pyramiding favorable alleles for the genetic improvement of wheat breeding populations.
QTL script and VCF files - GENETIC ARCHITECTURE OF FLOWERING-TIME VARIATION IN BRACHYPODIUM DISTACHYON
<p>Supplemental data for the article:</p> <p>GENETIC ARCHITECTURE OF FLOWERING-TIME VARIATION IN BRACHYPODIUM DISTACHYON (Woods et al., 2016).</p> <p>Supplemental data includes:</p> <ul> <li> <p>A folder named “Data”, which contains:</p> <ol> <li> <p><strong>Raw genotypic data</strong>. Data showing the parental genotype in a RIL population (F7; Bd21 X Bd1-1)</p> </li> <li> <p><strong>Raw phenotypic data</strong>. Flowering time (leaves and days to flowering) of the RIL population in different growth environments.</p> </li> <li> <p><strong>The genetic map </strong>(see material and methods section of Woods et al., 2016 for additional information).</p> </li> <li> <p><strong>Gene positions.</strong> A file containing the position of annotated genes on the Brachypodium distachyon genome V2.1. Data from Phytozome (https://phytozome.jgi.doe.gov/pz/portal.html).</p> </li> </ol> </li> <li> <p>The<strong> R script</strong> used for the QTL analysis (Final_script.R).</p> </li> <li> <p><strong>VCF files of the loci of interest</strong>. A folder called “VCF” includes the VCF files of the genes presented in Fig. 6 (VRN1, PHYC, VRN2, and FD). </p> </li> </ul> <p> </p>
QTL mapping for seedling and adult plant resistance to stripe and leaf rust in two winter wheat populations
<p><span>The two recombinant inbred lines (RIL) populations developed by crossing Almaly × Avocet S (206 RILs) and Almaly × Anza (162 RILs) were used to detect the novel genomic regions associated with adult plant resistance (APR) and seedling or all-stage resistance (ASR) to yellow rust (YR) and leaf rust (LR). Both the populations were evaluated for YR APR in two environments (2018 and 2019) and LR APR in three environments (2018, 2019, and 2020) in the Anza population and two environments (2018 and 2019) in the Avocet population; both the populations were phenotyped for one environment during 2020 for LR and YR ASR and genotyped using high throughput DArTseq technology. A set of 51 QTLs including 22 for YR APR, nine for LR APR, nine for YR ASR, and 11 for LR ASR were identified. Also, a set of 13 stable QTLs including nine QTLs (<em>QYR-APR-2A.1, QYR-APR-2A.2, QYR-APR-4D.2, QYR-APR-1B, QYR-APR-2B.1, QYR-APR-2B.2, QYR-APR-3D, QYR-APR-4D.1, </em>and<em> QYR-APR-4D.2</em>) for YR APR and four QTLs (<em>QLR-APR-4A, QLR-APR-2B, QLR-APR-3B, </em>and<em> </em></span><em>QLR-APR-5A.2</em>) <span>for LR APR were identified. </span><span>In silico analysis revealed that the key putative candidate genes such as <em>Cytochrome P450</em></span><em><span>, Protein kinase-like domain superfamily</span><span>, Zinc-binding ribosomal protein</span><span>, SANT/Myb domain</span><span>, WRKY transcription factor</span><span>, Nucleotide-sugar transporter,</span></em><span> and <em>NAC</em> </span><em><span>domain superfamily</span></em><span> were in the QTL regions and involved in the regulation of host response towards the pathogen infection. </span><span>The stable QTLs identified in this study are useful for developing rust-resistant varieties through marker-assisted selection (MAS).</span></p>
ScienceDex guides
Understand access before you commit
These curated guides explain access requirements, typical timelines, costs, and reuse considerations for widely used research datasets.
Allen Brain Atlas
Allen Brain Atlas is an Allen Institute collection of brain map atlases, datasets, APIs, and analysis tools covering mouse, human, and non-human primate brain resources.
Annotated Behaviour and Observability Dataset (ABODe)
ABODe is a University of Edinburgh DataShare dataset for behavior classification in group-housed mice using home-cage video, identities, bounding boxes, ground-plate positions, and annotator labels.
DANDI Archive for NWB datasets
DANDI is a BRAIN Initiative archive for publishing and sharing neurophysiology data, including electrophysiology, optophysiology, and behavioral data packaged as NWB and related standards.
International Brain Laboratory public data
The International Brain Laboratory public data releases expose standardized mouse decision-making experiments, including Neuropixels recordings, widefield calcium imaging, behavior, and session metadata accessed through the ONE API.
OpenNeuro
OpenNeuro is a free, open platform for sharing neuroimaging datasets, with public search, dataset pages, and download paths for web, S3, DataLad, and the OpenNeuro CLI.