The LPA Caller is capable of identifying the LPA Kringle-IV-2 (KIV-2) VNTR unit copy number from whole-genome sequencing (WGS) data. Due to high sequence similarity between the genes, a specialized caller is necessary to resolve the VNTR unit copy number.
The LPA Caller performs the following steps:
Determines total LPA KIV-2 VNTR unit copy number.
Determines the heterozygous LPA KIV-2 VNTR unit copy number if heterozygous KIV-2 markers are present.
Calls small variants in the LPA KIV-2 VNTR region based on the unit copy number along with allele counts from read information.
The LPA Caller requires WGS data aligned to a human reference genome with at least 30x coverage. Reference genome builds must be based on hg19
, GRCh37
, or hg38
.
The first step of LPA calling is to determine the unit copy number of LPA KIV-2. Reads aligned to the LPA KIV-2 are counted. The counts in each region are corrected for GC-bias, and then normalized to a diploid baseline. The GC-bias correction and normalization factors are determined from read counts in 3000 preselected 2 kb regions across the genome. These 3000 normalization regions were randomly selected from the portion of the reference genome having stable coverage across population samples.
The second step of LPA calling is to determine the heterozygous unit copy numbert of LPA KIV-2. Heterozygous unit copy number is determined using two specific linked SNV sites that have been identified as a combined marker allele that is always present in every copy of the repeat unit concordantly. That is, if any copy of the repeat unit in an LPA haplotype contains the ALT alleles at those two SNV sites, then every copy of the repeat unit in that LPA haplotype contains the ALT alleles at those two sites. The relative read coverage for the ALT and REF cases at these sites can therefore be used to determine the proportions of overall copy numbers across the KIV repeat array that belong to each haplotype.
2 small variants are detected from the read alignments. These variants occur in the LPA KIV-2 VNTR region where reads mapping to either of the 6 units in the reference are used for variant calling.
For each variant, reads containing either the variant allele or the nonvariant allele are counted and a binomial model is used to determine the likelihood for each possible variant allele copy number up to the maximum possible as determined from the LPA KIV-2 VNTR unit copy number.
The LPA Caller generates a <output-file-prefix>.targeted.json
file in the output directory. The output file is a JSON formatted file containing the fields below.
The lpa
fields are defined as below.
For the variants
the fields are defined as below.
The LPA Caller also generates a <output-file-prefix>.targeted.vcf[.gz]
file in the output directory. The output file is a VCFv4.2
formatted file possibly compressed.
Examples of the LPA Caller content in the output json file are shown below.
The following are example output files:
Fields in JSON | Explanation | Type and Possible Values |
---|---|---|
Fields in JSON | Explanation | Type and Possible Values |
---|---|---|
Fields in JSON | Explanation | Type and Possible Values |
---|---|---|
sample
The sample name.
string
dragenVersion
The version of DRAGEN.
string
lpa
The LPA targeted caller specific fields.
dictionary
kiv2CopyNumber
Total KIV-2 unit copy number
float
refMarkerAlleleCopyNumber
Null if Homozygous REF/ALT markers call
float, null
Float if Heterozygous markers call and stores the KIV-2 unit copy number of the allele having REF markers
altMarkerAlleleCopyNumber
Null if Homozygous REF/ALT markers call
float, null
Float if Heterozygous markers call and stores the KIV-2 unit copy number of the allele having ALT markers
type
"Heterozygous markers call" if we observe both REF and ALT markers
string, "Heterozygous markers call", "Homozygous REF markers call", "Homozygous ALT markers call"
"Homozygous REF markers call" if we observe only REF markers
"Homozygous ALT markers call" if we observe only ALT markers
variants
List of known variants that were detected in the KIV-2 region.
list of variants
hgvs
HGVS identifier of the variant
string
qual
Phred QUAL score of the variant
double
altCopyNumber
Copy number of the ALT variant
double
altCopyNumberQuality
Phred QUAL copy number of the ALT variant
double