I want to know more about pathogenicity of different variant types on a large scale¶
The instructions in each section include links to the relevant pages in the documentation. Links are tagged as:
- Tools - descriptive
- Data - descriptive
- Pre-made workflows
- Reference lists/tables
You can use LabKey or Participant Explorer to create a cohort of interest. You will need to make yourself familiar with the clinical data we have available and the LabKey API which you will use to access it. You can use both of these to create a list of gVCF or BAM filepaths. We also have a tutorial on cohort building to work through:
- 100kGP clinical and phenotype data
- Participant Explorer
- Search for participants
- Labkey API
- Building cohorts
From either of these, you can create a list of filepaths to the gVCFs or BAM files for each participant:
Work with the HPC¶
You can use the gVCF files directly, or you can call the variants yourself from the BAM files, using some of the tools installed on the HPC. You can learn more about the HPC and how to work with it:
- High Performance Cluster (HPC)
- Accessing the HPC
- Using software on the HPC
- How to submit jobs to LSF
There are folders on the HPC for your GECIP domain or Discovery forum. You should use your relevant folder as your working directory. These are also accessible from the desktop:
For functional annotation of VCFs, you can use the VEP or the functional annotation pipeline:
If you have other functional annotation software you wish to use, you can bring this in using containers.
Alternatively, functionality annotated versions of our aggreagate VCFs are available. Code-books are available with examples of how to query them:
- Functional annotation of germline aggregated variant calls
- AggV2 code book functional annotation queries
- Somatic aggregated variant calls
- somAgg code book functional annotation queries
You can now filter the VCFs to find variant types of interest using your preferred programming language. We provide support for coding in Python or R in the RE. You can use interactive coding tools such as RStudio and Jupyter notebooks, which you can use on the HPC:
Combine with publicly available data¶
A number of publicly datasets have been made available in the RE, such as gnomAD and ClinVar. You can include these in your analyses.
Compile text and figures¶
Use your preferred programming language and statistical tools to compare, correlate, verify and model your results. You can use LO Calc to create figures and tables. You can also write any notes in LO Writer.
The only way to get the results of your analysis out is using Airlock. You should include any notes you may have made by hand. It is your responsibility to ensure your data conforms to the Airlock rules and does not contain any identifying data.