Visitar URL original
dbvar/Structural_Variant_Sets/Nonredundant_Structural_Variants at master · ncbi/dbvar · GitHub
Skip to content

Latest commit

 

History

History

Nonredundant_Structural_Variants

Folders and files

NameName
Last commit message
Last commit date

parent directory

..
 
 
 
 
 
 
 
 
 
 
 
 

README.md

dbVar Human Nonredundant Structural Variants (NR SVs)

Work in progress - data subject to change

Documentation updated: 04/23/2020

Data Summary

See also https://ftp.ncbi.nlm.nih.gov/pub/dbVar/sandbox/sv_datasets/nonredundant/release_notes/NR_stats.latest.txt

FTP Directory: nonredundant
Last modified: Oct 30, 2023
File types: bed, bedpe, tsv

Deletions

Type and FTP Directory GRCh37 GRCh38
All 2,176,754 2,178,660
Common 211,888 211,465
Pathogenic 20,697 20,503
Somatic 23,385 23,344
All-ACMG 6,471 6,482
Common-ACMG 392 403
Pathogenic-ACMG 3,199 3,201
Somatic-ACMG 966 965

Duplications

Type and FTP Directory GRCh37 GRCh38
All 784,471 787,474
Common 59,576 59,150
Pathogenic 5,290 5,153
Somatic 15,105 15,079
All-ACMG 5,262 5,340
Common-ACMG 125 129
Pathogenic-ACMG 1,021 1,019
Somatic-ACMG 811 810

Insertions

Type and FTP Directory GRCh37 GRCh38
All 1,640,496 1,651,014
Common 121,240 121,187
Pathogenic 250 250
Somatic 4 5
All-ACMG 3,370 3,445
Common-ACMG 255 261
Pathogenic-ACMG 65 65

All files are available in bed, bedpe, and tsv formats:

https://ftp.ncbi.nlm.nih.gov/pub/dbVar/sandbox/sv_datasets/nonredundant/deletions/ https://ftp.ncbi.nlm.nih.gov/pub/dbVar/sandbox/sv_datasets/nonredundant/insertions/ https://ftp.ncbi.nlm.nih.gov/pub/dbVar/sandbox/sv_datasets/nonredundant/duplications/

Description of NR SV data files:

  • Sets of "non-redundant structural variations" (NR SVs) derived from dbVar are available via FTP as tab delimited files by assembly, GRCh37 & GRCh38, and type of variant.

  • Non-redundant refers to variant coordinates, i.e. chr, outermost start, and outermost stop. Please note: the non-redundant coordinates are based strictly on exact overlap of coordinates, not on partial overlaps.

  • Other features of NR SV files:

    • variant calls are from germline samples only (no somatic)
    • placements are "BestAvailable" on the assembly (guarantees no duplicate placements for a variant)
    • placements are on finished chromosomes only (not on NT_ or NW_ contigs)
    • placements are 1-based in the .tsv files
    • placements are zero-based start and 1-based stop in .bed and .bedpe files
    • insertion_length is set to sequence length if the sequence was submitted to dbVar without a specific insertion_length
    • insertions submitted to dbVar without insertion_length or submitted sequence are not included in the NR files
  • Other files based on NR SV files:

    • NR SV files annotated with overlapping ACMG genes
    • NR SV files in .bed format
    • NR SV files in .bedpe format

File format:

Column NR SV TSV File BED File BEDPE File
1 chr chr chr
2 outermost_start (1-based) outermost_start (0-based) outermost_start (0-based)
3 outermost_stop (1-based) outermost_stop (1-based) outermost_stop (1-based)
4 variant_count NR_SV_id .  (chrom2)
5 variant_type   -1  (start2)
6 method   -1  (end2)
7 analysis   NR_SV_id
8 platform   .  (score)
9 study   .  (strand1)
10 variant   .  (stramd2)
11 clinical_assertion   variant_count
12 clinvar_accession   variant_type
13 bin_size   method
14 min_insertion_length*   analysis
15 max_insertion_length*   platform
16     study
17     variant
18     clinical_assertion
19     clinvar_accession
20     bin_size
21     min_insertion_length
22     max_insertion_length

Please note:

  • * = NR_SV TSV fields 14 and 15 are in nr_insertion.tsv files only
  • NR_SV_id = chr_outermost_start_outermost_stop_type where type is del, dup, or ins
  • bin_size = small (length < 50 bp), medium (< 1000000), large (>= 1000000). Length = outermost_stop - outermost_start + 1.
  • In all cases, bedpe columns 4 through 6, and 8 through 10, are populated with default values per the bedpe specification
  • The bed and bedpe specifications are found here: https://bedtools.readthedocs.io/en/latest/content/general-usage.html

Some fields may have multiple values:

  • The fields type, method, analysis, platform, variant, study, clinical_significance, clinvar_accession, and gene may contain multiple values.
  • Each of the values is associated with one or more calls found in the variant field.
  • The values in the variant field are "dbVar call accessions".

Records in the NR SV files:

Records in the deletions or duplications NR SV files, e.g.:

chr outermost_start outermost_stop variant_count variant_type method analysis platform study variant clincical_assertion clinvar_accession bin_size
15 98085101 101843270 1 copy_number_loss Oligo_aCGH Probe_signal_intensity Agilent ISCA 44K ClinGen_Laboratory-Submitted nssv14082018 Pathogenic SCV000586438

Records in the insertions NR SV files, e.g.:

chr outermost_start outermost_stop variant_count variant_type method analysis platform study variant clincical_assertion clinvar_accession bin_size min_insertion_length max_insertion_length
1 1889055 1889055 1 alu_insertion Sequencing Split_read_and_paired-end_mapping HiSeq 2000 Gardner2017 nssv14051747 small 258 258
  • only insertion SVs have minimum_insertion_length and maximum_insertion_length fields

Records in the NR SV .bed files, e.g.

chr outermost_start outermost_stop NR_SV_id
chr1 0 10000 chr1_0_10000_del
  • Placements in bed files are zero-based start and one-based stop
  • name is comprised of chromosome, outermost_start, outermost_stop, and type (del, dup, or ins)
  • NR SV .bed files may be used with a variety of tools as shown in the tutorial:

https://github.com/ncbi/dbvar/blob/master/Structural_Variant_Sets/Nonredundant_Structural_Variants/ToolGuide.md

Records in the NR SV .bedpe files, e.g.

chr outermost_start outermost_stop chrom2 start2 end2 NR_SV_id score strand1 strand2 variant_count variant_type method analysis platform study variant clinical_assertion clinvar_accession
chr1 14873 7527302 . -1 -1 chr1_14873_7527302_del . . . 1 copy_number_loss Oligo_aCGH Probe_signal_intensity NA ClinGen_Laboratory-Submitted nssv13638713 Pathogenic SCV000495999
  • Placements in bedpe files are zero-based start and one-based stop
  • bedpe files are normally used for disjointed genomic sequences
  • the NR bedpe files do not contain disjointed genomic sequences, and instead use default values for chrom2, start2 and end2
  • the NR bedpe files use the optional fields starting at field 11 to hold: chr, outermost_start, outermost_stop, name, and all the additional fields which are in the NR SV files

ACMG files as "proof of principle" for a "use case"

Records in the NR SV ACMG files contain a field for the ACMG gene that overlaps the variant, e.g.:

chr outermost_start outermost_stop variant_count variant_type method analysis platform study variant clinical_significance clinvar_accession bin_size min_insertion_length max_insertion_length gene
3 30621876 30621876 2 alu_insertion Merging;Sequencing Merging;Split_read_and_paired-end_mapping See merged experiments;HiSeq 2000 1000_Genomes_Consortium_Phase_3_SV_Submission;Gardner2017 essv18243203;nssv14059593 small 279 280 TGFBR2

Caveats for the ACMG files:

  • ACMG files are provided as a "proof of principle" for a "use case" for the NR SV files.
  • ACMG files are based on region/gene overlaps, and are missing a few call/gene overlaps in the cases where the parent variant region does not include all of the variant call placement and the region does not overlap the gene or is upstream or downstream of the gene.
  • Placements in the ACMG files do not account for confidence intervals, even though the overlaps reported in the file were determined using confidence intervals. This results in a few of the overlaps that are not supported by the placements as reported in the file.

For information on ACMG genes please see: https://www.ncbi.nlm.nih.gov/clinvar/docs/acmg/

Methods and Analyses

example values

Methods include, e.g. Analyses include, e.g.
BAC_aCGH Curated
Curated Genotyping
MLPA Local_sequence_assembly
Merging Merging
Multiple Multiple
Not_provided Not_provided
Oligo_aCGH Optical_mapping
Optical_mapping Other
ROMA Paired-end_mapping
SNP_array Probe_signal_intensity
Sequencing Read_depth
qPCR SNP_genotyping_analysis
ROMA Sequence_alignment
SNP_array Split_read_mapping
Sequencing de_novo_sequence_assembly

README files for deletions, insertions, and duplications

Please see README files for deletions, insertions, and duplications for example records and additional details.

Brief Outline of algorithm used to generate NR-SVs.

The algorithm makes use of previously existing scripts.

Input files are generated from the dbVar database with tab separated values and contain SVs by assembly, type, and other relevant fields.

Selected type files are grouped into "aggregated type files" as specified above, by chr.

The "aggregated type files" are converted into XML records containing all the neccessary fields required by the nr process.

The XML is then parsed to generate SV records with coordinates, type, method, analysis, platform, insertion_length, SV accession and study.

The SV records are then processed to generate NR SV tab-separated value (tsv) files by assembly and type, as described above, e.g.

  • GRCh38.nr_deletions.tsv.gz

Tutorials

Examples of using the NR files with various tools and browsers can be found in: https://github.com/ncbi/dbvar/blob/master/Structural_Variant_Sets/Nonredundant_Structural_Variants/ToolGuide.md

Questions or feedback

Thanks!

Thanks for your interest in the dbVar human "non-redundant structural variations" (NR SVs) data files from NCBI.

Please check back soon for further updates.