No abstract

Background The concentrations of the highly atherogenic lipoprotein(a) [Lp(a)] are mainly genetically determined by the LPA gene locus. However, up to 70% of the coding sequence is located in the complex so-called kringle IV type 2 (KIV-2) copy number variation, a region hardly accessible by common genotyping and sequencing technologies. Despite its size, little is known about genetic variants in this complex region. The R21X variant is a functional variant located in this region, but it has never been analyzed in large cohorts. Methods We typed R21X in 10,910 individuals from three European populations using a newly developed high-throughput allele-specific qPCR assay. R21X allelic location was determined by separating the LPA alleles using pulsed-field gel electrophoresis (PFGE) and typing them separately. Using GWAS data, we identified a proxy SNP located outside of the KIV-2. Linkage disequilibrium was determined both statistically and by long-range haplotyping using PFGE. Worldwide frequencies were determined by reanalyzing the sequencing data of the 1000 Genomes Project with a dedicated pipeline. Results R21X carriers (frequency 0.016–0.021) showed significantly lower mean Lp(a) concentrations (− 11.7 mg/dL [− 15.5; − 7.82], p = 3.39e−32). The variant is located mostly on medium-sized LPA alleles. In the 1000 Genome data, R21X mostly occurs in Europeans and South Asians, is absent in Africans, and shows varying frequencies in South American populations (0 to 0.022). Of note, the best proxy SNP was another LPA null mutation (rs41272114, D ′ = 0.958, R 2 = 0.281). D ′ was very high in all 1000G populations (0.986–0.996), although rs41272114 frequency varies considerably (0–0.182). Co-localization of both null mutations on the same allele was confirmed by PFGE-based long-range haplotyping. Conclusions We performed the largest epidemiological study on an LPA KIV-2 variant so far, showing that it is possible to assess LPA KIV-2 mutations on a large scale. Surprisingly, in all analyzed populations, R21X was located on the same haplotype as the splice mutation rs41272114, creating “double-null” LPA alleles. Despite being a nonsense variant, the R21X status does not provide additional information beyond the rs41272114 genotype. This has important implications for studies using LPA loss-of-function mutations as genetic instruments and emphasizes the complexity of LPA genetics.

Massive parallel sequencing technologies are promising a highly sensitive detection of low-level mutations, especially in mitochondrial DNA (mtDNA) studies. However, processes from DNA extraction and library construction to bioinformatic analysis include several varying tasks. Further, there is no validated recommendation for the comprehensive procedure. In this study, we examined potential pitfalls on the sequencing results based on two-person mtDNA mixtures. Therefore, we compared three DNA polymerases, six different variant callers in five mixtures between 50% and 0.5% variant allele frequencies generated with two different amplification protocols. In total, 48 samples were sequenced on Illumina MiSeq. Low-level variant calling at the 1% variant level and below was performed by comparing trimming and PCR duplicate removal as well as six different variant callers. The results indicate that sensitivity, specificity, and precision highly depend on the investigated polymerase but also vary based on the analysis tools. Our data highlight the advantage of prior standardization and validation of the individual laboratory setup with a DNA mixture model. Finally, we provide an artificial heteroplasmy benchmark dataset that can help improve somatic variant callers or pipelines, which may be of great interest for research related to cancer and aging.