Bioinformatics. Группа авторов

Чтение книги онлайн.

Читать онлайн книгу Bioinformatics - Группа авторов страница 47

Bioinformatics - Группа авторов

Скачать книгу

presents genomic annotation in the form of tracks. Each track provides a different type of feature, from genes to SNPs to predicted gene regulatory regions to expression data. Each organism has its own set of tracks, some created by the UCSC Genome Bioinformatics team and others provided by members of the bioinformatics community. Over 200 tracks are available for the GRCh37 version of the human genome assembly. The newer human genome assembly, GRCh38, has fewer tracks, as not all the data have been remapped from the older assembly. Other genomes are not as well annotated as human; for example, fewer than 20 tracks are available for the sea hare. Some tracks, such as those created from NCBI transcript data, are updated weekly, while others, such as the SNP tracks created from NCBI variant data (Sayers et al. 2019), are updated less frequently, depending on the release schedule of the underlying data. For ease of use, tracks are organized into subsections. For example, depending on the organism, the Genes and Gene Predictions section may include evidence-based gene predictions, ab initio gene predictions, and/or alignment of protein sequences from other species.

      The home page of the UCSC Genome Browser provides a stepping-off point for many of the resources developed by the Genome Bioinformatics group at UCSC, including the Genome Browser, BLAT, and the Table Browser, which will be described in detail later in this chapter. The Tools menu provides a link to liftOver, a widely used tool that converts genomic coordinates from one assembly to another. Using this tool, it is possible to update annotation files so that old data can be integrated into a new genome assembly. The Download menu provides an option to download all the sequence and annotation data for each genome assembly hosted by UCSC, as well as some of the source code. The What's New section provides updates on new genome assemblies, as well as new tools and features. Finally, there is an extensive Help menu, with detailed documentation as well as videos. Users may also submit questions to a mailing list, and most queries are answered within a day.

      Box 4.1 Common File Types for Genomic Data

      Both the UCSC and Ensembl Genome Browsers allow users to upload their own data so that they can be viewed in context with other genome-scale data. User data must be formatted in a commonly used data structure in order to be interpreted correctly by the browser.

       Browser Extensible Data (BED) format is a tab-delimited format that is flexible enough to display many types of data. It can be used to display fairly simple features like the location of transcription binding factor sites, as well more complex ones like transcripts and their exons.

       Binary Alignment/Map (BAM) format is the compressed binary version of the Sequence Alignment/Map (SAM) format. It is a compact format designed for use with very large files of nucleotide sequence alignments. Because it can be indexed, only the portion of the file that is needed for display is transferred to the browser. Many tools for next generation sequence analysis use BAM format as output or input.

       Variant Call Format (VCF) is a flexible format for large files of variation data including single-nucleotide variants, insertions/deletions, copy number variants, and structural variants. Like BAM format, it is compressed and indexed, and only the portion of the file that is needed for display is transferred to the browser. Many tools for variant analysis use VCF format as output or input.

      To perform a search, enter text into the Position/Search Term box. If the query maps to a unique position in the genome, such as a search for a particular chromosome and position, the Go button links directly to the Genome Browser. However, if there is more than one hit for the query, such as a search for the term metalloprotease, the resulting page will contain a list of results that all contain that term. For some species, the terms have been indexed, and typing a gene symbol into the search box will bring up a list of possible matches. In this example, we will search for the human hypoxia inducible factor 1 alpha subunit (HIF1A) gene (Figure 4.1), which produces a single hit on GRCh38.

Snapshot depicts the home page of the UCSC Genome Browser, showing a query for the gene HIF1A on the human GRCh38 genome assembly in which the organism can be selected by clicking on its name in the phylogenetic tree. For many organisms, more than one genome assembly is available.

      Below the browser window illustrated in Figure 4.2, one would find a list of tracks that are available for display on the assembly. The tracks are separated

Скачать книгу