Lp(a) und KIV-2

Die KIV-2-Wiederholungen verlaengern oder verkuerzen das Apolipoprotein(a). Kurze Isoformen werden effizienter sezerniert und gehen tendenziell mit hoeheren Plasmaspiegeln einher, doch der Zusammenhang ist nicht streng: zusaetzliche funktionelle Varianten innerhalb und ausserhalb von KIV-2 koennen Isoformgroesse und Lp(a)-Spiegel voneinander entkoppeln. Da der hochrepetitive Bereich fuer die kurzgelesene Sequenzierung schwer zugaenglich ist, braucht es spezielle Verfahren statt einfacher Alignment-Zaehlung. Die Kopienzahl ist daher ein Anhaltspunkt fuer die Isoformgroesse, kein direkter Messwert des Plasmaspiegels.

Genome nutzt KILDA (KIv2 Length Determined from a kmer Analysis), eine alignmentfreie Nextflow-Pipeline, die die KIV-2-Kopienzahl direkt aus den Reads schaetzt. KILDA zaehlt mit jellyfish kurze Sequenzbausteine fester Laenge (k-mere, Standard k gleich 31) und vergleicht die mittlere Haeufigkeit KIV-2-spezifischer k-mere mit der mittleren Haeufigkeit von k-meren aus einer single-copy Normierungsregion im LPA-Gen. Das Verhaeltnis ergibt die geschaetzte Kopienzahl, wobei die Normierung tiefenbedingte Unterschiede zwischen Proben ausgleicht. In der von Genome verwendeten Konfiguration werden zusaetzlich drei Lp(a)-relevante rsIDs als k-mer-Signale gezaehlt. Das Ergebnis erscheint als technische Evidenz im Bericht zur Medizinischen Genomik, getrennt von Interpretation und Grenzen.

Der KIV-2-Block ist ein hochidentischer Tandem-Wiederholungsabschnitt (VNTR) von etwa 5,5 Kilobasen pro Kopie, weshalb sich kurzgelesene Reads vielen Positionen gleich gut zuordnen lassen und Mapper sie oft als mehrdeutig verwerfen oder falsch platzieren. Eine alignmentbasierte Zaehlung verliert dadurch genau dort an Genauigkeit, wo die Kopienzahl bestimmt werden soll. KILDA umgeht das, indem es nicht Positionen, sondern Vorkommen kurzer Sequenzbausteine zaehlt, die fuer KIV-2 spezifisch sind, und diese gegen eine single-copy Region normiert. Die Autoren berichten auf dem 1000-Genomes-Datensatz (n gleich 2459) eine Uebereinstimmung von R-Quadrat gleich 0,923 gegenueber dem DRAGEN-LPA-Caller und R-Quadrat gleich 0,962 gegenueber optisch kartierten Bionano-Allelen. Die Schaetzung bleibt eine relative Groesse und ersetzt weder eine direkte Lp(a)-Messung noch eine langgelesene Aufloesung des Repeats.