https://github.com/broadinstitute/eqtl_annotations
Science Score: 26.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
○CITATION.cff file
-
✓codemeta.json file
Found codemeta.json file -
✓.zenodo.json file
Found .zenodo.json file -
○DOI references
-
○Academic publication links
-
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (5.8%) to scientific vocabulary
Repository
Basic Info
- Host: GitHub
- Owner: broadinstitute
- Language: Python
- Default Branch: main
- Size: 84 KB
Statistics
- Stars: 1
- Watchers: 4
- Forks: 1
- Open Issues: 0
- Releases: 0
Metadata Files
README.md
eQTL_annotations
Village finemapped eQTL annotation pipeline.
Main pipeline
Overall: Annotates variants from our finemapped sets (eQTL pipeline) with information about if they are in peaks (e.g. ATAC or CHIP, etc.) and GTEx annotated regions (e.g. promoter, missense, splice site, etc.) * Peak overlaps - annotates all variants in the VCF with their closest peak for each peakfile * Also, takes in ABC specific input, for gene-specific peak pairs (peaks are associated with specific genes in this peakfile) * Merges all peak variant annotations with finemapped variants * Merges GTEx annotations (enhancer, promoter, etc.) with finemapped variants * Plots enrichment for each annotation category by pip bins: PIP < 0.01, 0.01 < PIP <0.1, 0.1 < PIP < 0.5, 0.5 < PIP < 0.9, 0.9 < PIP, for each day (each fiinemapped file) * Plot all finemapped files together, with enrichment & proportion of variants in each annotation category across all the variants
Peak Predictions
Optional section of the pipeline, dependent on the peak predictor file input to the pipeline (which for now, should be an ATAC peak file in bed format.) Model: a logistic regression model that takes in each peak-gene pair in the format outlined below.
- Variants each have a set of annotations, a distance, a gene link, and a PIP
- Label each variant-annot-gene-pip with a peakname
- For each peak-gene pair:
- take "or" of all annots (if any variant linked to that peak-gene is an enhancer, promoter, etc.)
- take average of distance to gene TSS
- take max of PIP
- Output = peak-maxpip-annots-distance-gene
Train on high and low PIP peak-gene pairs (PIP > 0.1, PIP < 0.03>) in all other chromosomes (leave-one-out chr) to get AUC. Score all peak-gene pairs w/ leave-one-out chromosome training.
Similar to ABC eQTL models, we remove any peaks that have variants in a promoter region, 3' or 5' UTR, frameshift, missense, stop-gain and splice variants. We also remove variants < 250 bp from the TSS.
Each peak-gene gets assigned a predicted probability of being a 'regulatory region', AKA containing a high-PIP variant.
Owner
- Name: Broad Institute
- Login: broadinstitute
- Kind: organization
- Location: Cambridge, MA
- Website: http://www.broadinstitute.org/
- Twitter: broadinstitute
- Repositories: 1,083
- Profile: https://github.com/broadinstitute
Broad Institute of MIT and Harvard
GitHub Events
Total
Last Year
Issues and Pull Requests
Last synced: over 1 year ago
All Time
- Total issues: 0
- Total pull requests: 2
- Average time to close issues: N/A
- Average time to close pull requests: about 8 hours
- Total issue authors: 0
- Total pull request authors: 1
- Average comments per issue: 0
- Average comments per pull request: 0.0
- Merged pull requests: 2
- Bot issues: 0
- Bot pull requests: 0
Past Year
- Issues: 0
- Pull requests: 2
- Average time to close issues: N/A
- Average time to close pull requests: about 8 hours
- Issue authors: 0
- Pull request authors: 1
- Average comments per issue: 0
- Average comments per pull request: 0.0
- Merged pull requests: 2
- Bot issues: 0
- Bot pull requests: 0
Top Authors
Issue Authors
Pull Request Authors
- thouis (4)