Permalink
Branch: master
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
274 lines (197 sloc) 10.5 KB

CyVerse logo

Home_Icon Learning Center Home

The Annotation Tutorial

This tutorial is a step-by-step guide for using SciApps to perform MAKER based annotation.


Example Data

Input Description Example
Assembled genome A scaled-down genome that is comprised of the first 300kb of three chromosomes of rice test_genome.fasta.gz
Annotated gene models MAKER output in GFF3 format maker_out.gff.gz

Apps:

App name Version Description App link Notes/other links
MAKER 2.32 A portable and easily configurable genome annotation pipeline MAKER-2.32 MAKER documentation
SNAP 0.0.1 Semi-HMM-based Nucleic Acid Parser SNAP-0.0.1 SNAP documentation

Step 1: Requiring access to SciApps

This is one-time operation. Please login to SciApps directly if you have completed this step before.

  1. Log into CyVerse User portal at https://user.cyverse.org.

  2. By default, you will be under the 'Services' page, click on 'AVAILABLE', then 'REQUEST ACCESS' to SciApps.

    cyverse_user

  3. Click on 'MY SERVICES', then click on 'LAUNCH' for Discovery Environment.

  4. Once in Discovery Environment, click to open the 'Data' window. You should see the sci_data folder under your root folder:/iplant/home/YOUR_USER_NAME.

    de_data


Step 2: Uploading data for SciApps

This step will demo how to upload data to the sci_data folder for accessing from SciApps.

  1. Click sci_data folder to open it.

  2. Click 'Upload', then 'Import from URL' to import this URL: https://data.sciapps.org/example_data/maker/my.all.gff.gz

    url_window

    Note

    This may take a few minutes. You can check the status by clicking the 'Bell' on the top right corner of DE. Once importing completed, 'Refresh' the window to see the file. This is a GFF3 formatted file from MAKER.

  3. Alternatively, use Cyberduck or iCommands for bulk data transfer to the sci_data folder.

Step 3: HMM parameters estimation with SNAP

  1. Login to SciApps at https://www.SciApps.org.

  2. Click Prediciton category (left panel) or search for SNAP, then click SNAP to load SNAP-0.0.1.

  3. Under “GFF file”, click Browse DataStore, then navigate to the maker folder (example > maker); select maker_out.gff.gz and click 'Select and Close'.

    data_window2

    Tip

    Click 'Refresh' if you can not see any newly uploaded files.

  4. Leave other parameters as default, and click Submit Job. You will be asked to confirm; click "Submit". You will be prompted to check the job status in the right panel.

Step 4: Running MAKER with SNAP output

This step will demo how to use SNAP output with MAKER to do a second around annotation.

  1. Click Annotation category (left panel) or search for MAKER, then click to load MAKER-2.32.

  2. Under “Genome sequence file” click Browse DataStore, then navigate to the maker folder (example > maker); select test_genome.fasta.gz and click 'Select and Close'.

  3. Click SNAP-0.0.1 in the History panel to expand its outputs, then drag and drop snap_out.hmm into the SNAP HMM file field.

    build_workflow2

  4. Under “Maker annotations” click Browse DataStore, then navigate to the maker folder (example > maker); select maker_out.gff.gz and click 'Select and Close'.

  5. Leave others as defaults, then click the "Submit Job" button.

  6. Once COMPLETED, click the Visualization icon for MAKER-2.32 in the History panel to bring up its outputs. Select jbrowse_out.view.tgz from the list of outputs, then click Visualize, you will be directed to a genome browser to visualize your annotation results.

    jbrowse_out

Step 5: Creating a Workflow

This step will demo how to build a two-step workflow with previously completed MAKER and SNAP jobs.

  1. Check the checkboxes for step 1 (SNAP), and step 2 (MAKER) in the History panel, then click the ‘build a workflow’ link to load the Workflow building page.

    build_workflow3

    Tip

    History panel Checkboxes and the workflow building page are interactive. Use the 'Select All' or 'Reset' button to simplify the selection process.

  2. Modify Workflow Name and Workflow Description, then click the ‘Build Workflow’ button to visualize the workflow.

    annotation_workflow2

    Note

    The connection between SNAP-0.0.1 and MAKER-2.32 (via my.all.hmm) is recorded through dragging and dropping, which feeds the output of SNAP as an input for MAKER.

  3. On the 'Workflow Diagram', you can save the workflow. Your saved workflows will appear in 'My workflows' (under the 'Workflow' menu from top navigation panel).

    Note

    You can select the workfrom from 'My workflows' and click share to get a direct link, which can be passed to others for sharing the entire analysis. workflow_URL


Step 6: Running a Workflow

This step will demo how to run a workflow you created or someone shared with you.

  1. Navigate to ‘Workflow’, then ‘My workflows’, to load the workflow you created and saved (in Step 2).

    myworkflows_window

    Note

    Alternatively, you can load the app forms and job histories directly if you have the direct link for a workflow.

  2. Scroll down the main panel, then click Submit Workflow. You will be asked to confirm and prompted to check the job status in the right panel. Then a live workflow diagram will be displayed with real-time analysis status updates.


Step 7: Using Apollo for Community Annotation

In reality, annotated genes from MAKER will be further filtered or even manually annotated before being released (for example, by Gramene/Plant Ensembl). For manual annotation of the MAKER results with Apollo, we set up a demo at http://data.maizecode.org/apollo. You can login with username: demo@demo.com, and password: demo.

Summary

This tutorial covers how to use SciApps for your annotaition work, including accessing data in CyVerse Data Store, launching jobs, building workflows, running workflows, visualizing results, and importing workflows to re-run.

More help and additional information

Search for an answer:
CyVerse Learning Center or CyVerse Wiki
Post your question to the user forum:
Ask CyVerse

More resources:


Fix or improve this documentation:


Home_Icon Learning Center Home