The kringle IV type 2 domain variant 4925G>A causes the elusive association signal of the LPA pentanucleotide repeat.

Lipoprotein(a) [Lp(a)] concentrations are regulated by the LPA gene mainly via the large kringle IV-type 2 (KIV-2) copy number variation and multiple causal variants. Early studies suggested an effect of long pentanucleotide repeat (PNR) alleles (10 and 11 repeats, PNR10 and PNR11) in the LPA promoter on gene transcription and found an association with lower Lp(a). Subsequent in vitro studies showed no effects on mRNA transcription, but the association with strongly decreased Lp(a) remained consistent. We investigated the isolated and combined effect of PNR10, PNR11, and the frequent splice site variant KIV-2 4925G>A on Lp(a) concentrations in the Cooperative Health Research in the Region of Augsburg F4 study by multiple quantile regression in single-SNP and joint models. Data on Lp(a), apolipoprotein(a) Western blot isoforms, and variant genotypes were available for 2,858 individuals. We found a considerable linkage disequilibrium between KIV-2 4925G>A and the alleles PNR10 and PNR11. In single-variant analysis adjusted for age, sex, and the shorter apo(a) isoform, we determined that both PNR alleles were associated with a highly significant Lp(a) decrease (PNR10: ß = -14.43 mg/dl, 95% CI: -15.84, -13.02, P = 3.33e-84; PNR11: ß = -17.21 mg/dl, 95% CI: -20.19, -14.23, P = 4.01e-29). However, a joint model, adjusting the PNR alleles additionally for 4925G>A, abolished the effect on Lp(a) (PNR10: ß = +0.44 mg/dl, 95% CI: -1.73, 2.60, P = 0.69; PNR11: ß = -1.52 mg/dl, 95% CI: -6.05, 3.00, P = 0.51). Collectively, we conclude that the previously reported Lp(a) decrease observed in pentanucleotide alleles PNR10 or PNR11 carriers results from a linkage disequilibrium with the frequent splicing mutation KIV-2 4925G>A.

Lysis reagents, cell numbers, and calculation method influence high-throughput measurement of HDL-mediated cholesterol efflux capacity.

HDL-mediated cholesterol efflux capacity (CEC) may protect against cardiovascular disease. However, CEC assays are not standardized, hampering their application in large cohorts and comparison between studies. To improve standardization, we systematically investigated technical differences between existing protocols that influence assay performance that have not been previously addressed. CEC was measured in 96-well plates using J774A.1 macrophages labeled with BODIPY-cholesterol and incubated for 4 h with 2% apolipoprotein B-depleted human serum. The time zero method, which calculates CEC using control wells, and the per-well method, which calculates CEC based on the actual content of BODIPY-cholesterol in each well, were compared in 506 samples. We showed that the per-well method had a considerably lower sample rejection rate (4.74% vs. 13.44%) and intra-assay (4.48% vs. 5.28%) and interassay coefficients of variation (two controls: 7.85%, 9.86% vs. 13.58%, 15.29%) compared with the time zero method. Correction for plate-to-plate differences using four controls on each plate also improved assay performance of both methods. In addition, we observed that the lysis reagent used had a significant effect. Compared with cholic acid, lysis with sodium hydroxide results in higher (P = 0.0082) and Triton X-100 in lower (P = 0.0028) CEC values. Furthermore, large cell seeding errors (30% variation) greatly biased CEC for both referencing methods (P < 0.0001) as measured by a resazurin assay. In conclusion, lysis reagents, cell numbers, and assay setup greatly impact the quality and reliability of CEC quantification and should be considered when this method is newly established in a laboratory.

Analyzing Low-Level mtDNA Heteroplasmy-Pitfalls and Challenges from Bench to Benchmarking.

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.

Resolving intra-repeat variation in medically relevant VNTRs from short-read sequencing data using the cardiovascular risk gene LPA as a model.

BACKGROUND: Variable number tandem repeats (VNTRs) are highly polymorphic DNA regions harboring many potentially disease-causing variants. However, VNTRs often appear unresolved ("dark") in variation databases due to their repetitive nature. One particularly complex and medically relevant VNTR is the KIV-2 VNTR located in the cardiovascular disease gene LPA which encompasses up to 70% of the coding sequence. RESULTS: Using the highly complex LPA gene as a model, we develop a computational approach to resolve intra-repeat variation in VNTRs from largely available short-read sequencing data. We apply the approach to six protein-coding VNTRs in 2504 samples from the 1000 Genomes Project and developed an optimized method for the LPA KIV-2 VNTR that discriminates the confounding KIV-2 subtypes upfront. This results in an F1-score improvement of up to 2.1-fold compared to previously published strategies. Finally, we analyze the LPA VNTR in > 199,000 UK Biobank samples, detecting > 700 KIV-2 mutations. This approach successfully reveals new strong Lp(a)-lowering effects for KIV-2 variants, with protective effect against coronary artery disease, and also validated previous findings based on tagging SNPs. CONCLUSIONS: Our approach paves the way for reliable variant detection in VNTRs at scale and we show that it is transferable to other dark regions, which will help unlock medical information hidden in VNTRs.

KLKB1 and CLSTN2 are associated with HDL-mediated cholesterol efflux capacity in a genome-wide association study.

BACKGROUND AND AIMS: HDL-mediated cholesterol efflux capacity (CEC) may protect from cardiovascular disease. Thus, we aimed to identify its genetic and non-genetic determinants. METHODS: We measured CEC to 2% apolipoprotein B-depleted serum using BODIPY-cholesterol and cAMP-stimulated J774A.1 macrophages using serum samples from 4,981 participants in the German Chronic Kidney Disease (GCKD) study. Variance of CEC explained by clinical and biochemical parameters in a multivariable linear regression model was calculated by proportional marginal variance decomposition. A genome-wide association study with 7,746,917 variants was performed based on an additive genetic model. The main model was adjusted for age, sex and principal components 1-10. Further models were selected for sensitivity analysis and to reduce residual variance by known CEC pathways. RESULTS: Variables that explained 1% and more of the variance of CEC were concentrations of triglycerides (12.9%), HDL-cholesterol (11.8%), LDL-cholesterol (3.0%), apolipoprotein A-IV (2.8%), PCSK9 (1.0%), and eGFR (1.0%). The KLKB1 (chr4) and APOE/C1 (chr19) loci were genome-wide significantly (p < 5x10-8) associated with CEC in our main model (p = 8.8x10-10 and p = 3.3x10-10, respectively). KLKB1 remained significantly associated after additional adjustment for either kidney parameters, HDL-cholesterol, triglycerides or apolipoprotein A-IV concentrations, while the APOE/C1 locus was not significantly associated anymore after adjustment for triglycerides. Adjustment for triglycerides also revealed an association with the CLSTN2 locus (chr3; p = 6.0x10-9). CONCLUSIONS: We identified HDL-cholesterol and triglycerides as the main determinants of CEC. Furthermore, we newly found a significant association of CEC with the KLKB1 and the CLSTN2 locus and confirmed the association with the APOE/C1 locus, likely mediated by triglycerides.

The effect of LPA Thr3888Pro on lipoprotein(a) and coronary artery disease is modified by the LPA KIV-2 variant 4925G>A.

BACKGROUND AND AIMS: High lipoprotein(a) [Lp(a)] concentrations are associated with increased coronary artery disease (CAD) risk. Lp(a) is regulated mainly genetically by the LPA gene but involved genetic variants have not been fully elucidated. Improved understanding of the entanglements of genetic Lp(a) regulation may enhance genetic prediction of Lp(a) and CAD risk. We investigated an interaction between the well-known LPA missense SNP rs41272110 (known as Thr3888Pro) and the frequent LPA splicing mutation KIV-2 4925G>A. METHODS: Effects on Lp(a) concentrations were investigated by multiple quantile regression in the German Chronic Kidney Disease (GCKD) study, KORA-F3 and KORA-F4 (ntotal = 10,405) as well as in the UK Biobank (UKB) 200k exome dataset (n = 173,878). The impact of the interaction on CAD risk was assessed by survival analysis in UKB. RESULTS: We observed a significant SNP-SNP interaction in all studies (p = 1.26e-05 to 3.03e-04). In quantile regression analysis, rs41272110 as a predictor shows no impact on Lp(a) (ß = -0.06 [-0.79; 0.68], p = 0.879), but in a joint model including both SNPs as predictors, rs41272110 is associated with markedly higher Lp(a) (ß = +9.40 mg/dL [6.45; 12.34], p = 4.07e-10). Similarly, rs41272110 shows no effect on CAD in UKB (HR = 1.01 [0.97; 1.04], p = 0.731), while rs41272110 carriers not carrying 4925G>A show an increased CAD risk (HR = 1.10 [1.04; 1.16], p = 6.9e-04). This group corresponds to 4% of the population. Adjustment for apolipoprotein(a) isoforms further modified the effect estimates markedly. CONCLUSIONS: This work emphasizes the complexity of the genetic regulation of Lp(a) and the importance to account for genetic subgroups in Lp(a) association studies and when interpreting genetic cardiovascular risk profiles.

An in-depth analysis of the mitochondrial phylogenetic landscape of Cambodia.

Cambodia harbours a variety of human aboriginal populations that have scarcely been studied in terms of genetic diversity of entire mitochondrial genomes. Here we present the matrilineal gene pool of 299 Cambodian refugees from three different ethnic groups (Cham, Khmer, and Khmer Loeu) deriving from 16 Cambodian districts. After establishing a DNA-saving high-throughput strategy for mitochondrial whole-genome Sanger sequencing, a HaploGrep based workflow was used for quality control, haplogroup classification and phylogenetic reconstruction. The application of diverse phylogenetic algorithms revealed an exciting picture of the genetic diversity of Cambodia, especially in relation to populations from Southeast Asia and from the whole world. A total of 224 unique haplotypes were identified, which were mostly classified under haplogroups B5a1, F1a1, or categorized as newly defined basal haplogroups or basal sub-branches of R, N and M clades. The presence of autochthonous maternal lineages could be confirmed as reported in previous studies. The exceptional homogeneity observed between and within the three investigated Cambodian ethnic groups indicates genetic isolation of the whole population. Between ethnicities, genetic barriers were not detected. The mtDNA data presented here increases the phylogenetic resolution in Cambodia significantly, thereby highlighting the need for an update of the current human mtDNA phylogeny.

Frequent LPA KIV-2 Variants Lower Lipoprotein(a) Concentrations and Protect Against Coronary Artery Disease.

BACKGROUND: Lipoprotein(a) (Lp(a)) concentrations are a major independent risk factor for coronary artery disease (CAD) and are mainly determined by variation in LPA. Up to 70% of the LPA coding sequence is located in the hypervariable kringle IV type 2 (KIV-2) region. It is hardly accessible by conventional technologies, but may contain functional variants. OBJECTIVES: This study sought to investigate the new, very frequent splicing variant KIV-2 4733G>A on Lp(a) and CAD. METHODS: We genotyped 4733G>A in the GCKD (German Chronic Kidney Disease) study (n = 4,673) by allele-specific polymerase chain reaction, performed minigene assays, identified proxy single nucleotide polymorphisms and used them to characterize its effect on CAD by survival analysis in UK Biobank (n = 440,234). Frequencies in ethnic groups were assessed in the 1000 Genomes Project. RESULTS: The 4733G>A variant (38.2% carrier frequency) was found in most isoform sizes. It reduces allelic expression without abolishing protein production, lowers Lp(a) by 13.6 mg/dL (95% CI: 12.5-14.7; P < 0.0001) and is the strongest variance-explaining factor after the smaller isoform. Splicing of minigenes was modified. Compound heterozygosity (4.6% of the population) for 4733G>A and 4925G>A, another KIV-2 splicing mutation, reduces Lp(a) by 31.8 mg/dL and most importantly narrows the interquartile range by 9-fold (from 42.1 to 4.6 mg/dL) when compared to the wild type. In UK Biobank 4733G>A alone and compound heterozygosity with 4925G>A reduced HR for CAD by 9% (95% CI: 7%-11%) and 12% (95% CI: 7%-16%) (both P < 0.001). Frequencies in ethnicities differ notably. CONCLUSIONS: Functional variants in the previously inaccessible LPA KIV-2 region cooperate in determining Lp(a) variance and CAD risk. Even a moderate but lifelong genetic Lp(a) reduction translates to a noticeable CAD risk reduction.

When the genome bluffs: a tandem duplication event during generation of a novel Agmo knockout mouse model fools routine genotyping.

BACKGROUND: Genome editing in mice using either classical approaches like homologous recombination or CRISPR/Cas9 has been reported to harbor off target effects (insertion/deletion, frame shifts or gene segment duplications) that lead to mutations not only in close proximity to the target site but also outside. Only the genomes of few engineered mouse strains have been sequenced. Since the role of the ether-lipid cleaving enzyme alkylglycerol monooxygenase (AGMO) in physiology and pathophysiology remains enigmatic, we created a knockout mouse model for AGMO using EUCOMM stem cells but unforeseen genotyping issues that did not agree with Mendelian distribution and enzyme activity data prompted an in-depth genomic validation of the mouse model. RESULTS: We report a gene segment tandem duplication event that occurred during the generation of an Agmo knockout-first allele by homologous recombination. Only low homology was seen between the breakpoints. While a single copy of the recombinant 18 kb cassette was integrated correctly around exon 2 of the Agmo gene, whole genome nanopore sequencing revealed a 94 kb duplication in the Agmo locus that contains Agmo wild-type exons 1-3. The duplication fooled genotyping by routine PCR, but could be resolved using qPCR-based genotyping, targeted locus amplification sequencing and nanopore sequencing. Despite this event, this Agmo knockout mouse model lacks AGMO enzyme activity and can therefore be used to study its physiological role. CONCLUSIONS: A duplication event occurred at the exact locus of the homologous recombination and was not detected by conventional quality control filters such as FISH or long-range PCR over the recombination sites. Nanopore sequencing provides a cost convenient method to detect such underrated off-target effects, suggesting its use for additional quality assessment of gene editing in mice and also other model organisms.

Investigation of a nonsense mutation located in the complex KIV-2 copy number variation region of apolipoprotein(a) in 10,910 individuals.

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, R2 = 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.

Profiling of Mitochondrial DNA Heteroplasmy in a Prospective Oral Squamous Cell Carcinoma Study.

While a shift in energy metabolism is essential to cancers, the knowledge about the involvement of the mitochondrial genome in tumorigenesis and progression in oral squamous cell carcinoma (OSCC) is still very limited. In this study, we evaluated 37 OSCC tumors and the corresponding benign mucosa tissue pairs by deep sequencing of the complete mitochondrial DNA (mtDNA). After extensive quality control, we identified 287 variants, 137 in tumor and 150 in benign samples exceeding the 1% threshold. Variant heteroplasmy levels were significantly increased in cancer compared to benign tissues (p = 0.0002). Furthermore, pairwise high heteroplasmy frequency difference variants (∆HF% > 20) with potential functional impact were increased in the cancer tissues (p = 0.024). Fourteen mutations were identified in the protein-coding region, out of which thirteen were detected in cancer and only one in benign tissue. After eight years of follow-up, the risk of mortality was higher for patients who harbored at least one ∆HF% > 20 variant in mtDNA protein-coding regions relative to those with no mutations (HR = 4.6, (95%CI = 1.3-17); p = 0.019 in primary tumor carriers). Haplogroup affiliation showed an impact on survival time, which however needs confirmation in a larger study. In conclusion, we observed a significantly higher accumulation of somatic mutations in the cancer tissues associated with a worse prognosis.